Mathematical modeling of viral infection and the immune response controlled by the circadian clock

Time of day affects how well the immune system responds to viral or bacterial infections. While it is well known that the immune system is regulated by the circadian clock, the dynamic origin of time-of-day-dependent immunity remains unclear. In this paper, we studied the circadian control of immune response upon infection of influenza A virus through mathematical modeling. Dynamic simulation analyses revealed that the time-of-day-dependent immunity was rooted in the relative phase between the circadian clock and the pulse of viral infection. The relative phase, which depends on the time the infection occurs, plays a crucial role in the immune response. It can drive the immune system to one of two distinct bistable states, a high inflammatory state with a higher mortality rate or a safe state characterized by low inflammation. The mechanism we found here also explained why the same species infected by different viruses has different time-of-day-dependent immunities. Further, the time-of-day-dependent immunity was found to be abolished when the immune system was regulated by an impaired circadian clock with decreased oscillation amplitude or without oscillations.

Temperature compensation through kinetic regulation in biochemical oscillators

Nearly all circadian clocks maintain a period that is insensitive to temperature changes, a phenomenon known as temperature compensation (TC). Yet, it is unclear whether there is any common feature among different systems that exhibit TC. From a general timescale invariance, we show that TC relies on the existence of certain period-lengthening reactions wherein the period of the system increases strongly with the rates in these reactions. By studying several generic oscillator models, we show that this counterintuitive dependence is nonetheless a common feature of oscillators in the nonlinear (far-from-onset) regime where the oscillation can be separated into fast and slow phases. The increase of the period with the period-lengthening reaction rates occurs when the amplitude of the slow phase in the oscillation increases with these rates while the progression speed in the slow phase is controlled by other rates of the system. The positive dependence of the period on the period-lengthening rates balances its inverse dependence on other kinetic rates in the system, which gives rise to robust TC in a wide range of parameters. We demonstrate the existence of such period-lengthening reactions and their relevance for TC in all four model systems we considered. Theoretical results for a model of the Kai system are supported by experimental data. A study of the energy dissipation also shows that better TC performance requires higher energy consumption. Our study unveils a general mechanism by which a biochemical oscillator achieves TC by operating in parameter regimes far from the onset where period-lengthening reactions exist.

Determination of Spin-Parity Quantum Numbers of X(2370) as 0^{-+} from J/ψ→γK_{S}^{0}K_{S}^{0}η^{'}

Based on (10087±44)×10^{6}  J/ψ events collected with the BESIII detector, a partial wave analysis of the decay J/ψ→γK_{S}^{0}K_{S}^{0}η^{'} is performed. The mass and width of the X(2370) are measured to be 2395±11(stat)_{-94}^{+26}(syst)  MeV/c^{2} and 188_{-17}^{+18}(stat)_{-33}^{+124}(syst)  MeV, respectively. The corresponding product branching fraction is B[J/ψ→γX(2370)]×B[X(2370)→f_{0}(980)η^{'}]×B[f_{0}(980)→K_{S}^{0}K_{S}^{0}]=(1.31±0.22(stat)_{-0.84}^{+2.85}(syst))×10^{-5}. The statistical significance of the X(2370) is greater than 11.7σ and the spin parity is determined to be 0^{-+} for the first time. The measured mass and spin parity of the X(2370) are consistent with the predictions of the lightest pseudoscalar glueball.

Observation of Structures in the Processes e^{+}e^{-}→ωχ_{c1} and ωχ_{c2}

We present measurements of the Born cross sections for the processes e^{+}e^{-}→ωχ_{c1} and ωχ_{c2} at center-of-mass energies sqrt[s] from 4.308 to 4.951 GeV. The measurements are performed with data samples corresponding to an integrated luminosity of 11.0  fb^{-1} collected with the BESIII detector operating at the Beijing Electron Positron Collider storage ring. Assuming the e^{+}e^{-}→ωχ_{c2} signals come from a single resonance, the mass and width are determined to be M=(4413.6±9.0±0.8)  MeV/c^{2} and Γ=(110.5±15.0±2.9)  MeV, respectively, which is consistent with the parameters of the well-established resonance ψ(4415). In addition, we also use one single resonance to describe the e^{+}e^{-}→ωχ_{c1} line shape and determine the mass and width to be M=(4544.2±18.7±1.7)  MeV/c^{2} and Γ=(116.1±33.5±1.7)  MeV, respectively. The structure of this line shape, observed for the first time, requires further understanding.

Coupled-Channel Analysis of the χ_{c1}(3872) Line Shape with BESIII Data

We perform a study of the χ_{c1}(3872) line shape using the data samples of e^{+}e^{-}→γχ_{c1}(3872), χ_{c1}(3872)→D^{0}D[over ¯]^{0}π^{0}, and π^{+}π^{-}J/ψ collected with the BESIII detector. The effects of the coupled channels and the off-shell D^{*0} are included in the parametrization of the line shape. The line shape mass parameter is obtained to be M_{X}=(3871.63±0.13_{-0.05}^{+0.06})  MeV. Two poles are found on the first and second Riemann sheets corresponding to the D^{*0}D[over ¯]^{0} branch cut. The pole location on the first sheet is much closer to the D^{*0}D[over ¯]^{0} threshold than the other, and is determined to be 7.04±0.15_{-0.08}^{+0.07}  MeV above the D^{0}D[over ¯]^{0}π^{0} threshold with an imaginary part -0.19±0.08_{-0.19}^{+0.14}  MeV.

Observation of the Anomalous Shape of X(1840) in J/ψ→γ3(π^{+}π^{-}) Indicating a Second Resonance Near pp[over ¯] Threshold

Using a sample of (10087±44)×10^{6}  J/ψ events, which is about 45 times larger than that was previously analyzed, a further investigation on the J/ψ→γ3(π^{+}π^{-}) decay is performed. A significant distortion at 1.84  GeV/c^{2} in the line shape of the 3(π^{+}π^{-}) invariant mass spectrum is observed for the first time, which could be resolved by two overlapping resonant structures, X(1840) and X(1880). The new state X(1880) is observed with a statistical significance larger than 10σ. The mass and width of X(1880) are determined to be 1882.1±1.7±0.7  MeV/c^{2} and 30.7±5.5±2.4  MeV, respectively, which indicates the existence of a pp[over ¯] bound state.

Study of the f_{0}(980) and f_{0}(500) Scalar Mesons through the Decay D_{s}^{+}→π^{+}π^{-}e^{+}ν_{e}

Using e^{+}e^{-} collision data corresponding to an integrated luminosity of 7.33  fb^{-1} recorded by the BESIII detector at center-of-mass energies between 4.128 and 4.226 GeV, we present an analysis of the decay D_{s}^{+}→π^{+}π^{-}e^{+}ν_{e}, where the D_{s}^{+} is produced via the process e^{+}e^{-}→D_{s}^{*±}D_{s}^{∓}. We observe the f_{0}(980) in the π^{+}π^{-} system and the branching fraction of the decay D_{s}^{+}→f_{0}(980)e^{+}ν_{e} with f_{0}(980)→π^{+}π^{-} measured to be (1.72±0.13_{stat}±0.10_{syst})×10^{-3}, where the uncertainties are statistical and systematic, respectively. The dynamics of the D_{s}^{+}→f_{0}(980)e^{+}ν_{e} decay are studied with the simple pole parametrization of the hadronic form factor and the Flatté formula describing the f_{0}(980) in the differential decay rate, and the product of the form factor f_{+}^{f_{0}}(0) and the c→s Cabibbo-Kobayashi-Maskawa matrix element |V_{cs}| is determined for the first time to be f_{+}^{f_{0}}(0)|V_{cs}|=0.504±0.017_{stat}±0.035_{syst}. Furthermore, the decay D_{s}^{+}→f_{0}(500)e^{+}ν_{e} is searched for the first time but no signal is found. The upper limit on the branching fraction of D_{s}^{+}→f_{0}(500)e^{+}ν_{e}, f_{0}(500)→π^{+}π^{-} decay is set to be 3.3×10^{-4} at 90% confidence level.

Observation of D^{+}→K_{S}^{0}a_{0}(980)^{+} in the Amplitude Analysis of D^{+}→K_{S}^{0}π^{+}η

We perform for the first time an amplitude analysis of the decay D^{+}→K_{S}^{0}π^{+}η and report the observation of the decay D^{+}→K_{S}^{0}a_{0}(980)^{+} using 2.93  fb^{-1} of e^{+}e^{-} collision data taken at a center-of-mass energy of 3.773 GeV with the BESIII detector. As the only W-annihilation-free decay among D to a_{0}(980) pseudoscalar, D^{+}→K_{S}^{0}a_{0}(980)^{+} is the ideal decay in extracting the contributions of the W-emission amplitudes involving a_{0}(980) and to study the final-state interactions. The absolute branching fraction of D^{+}→K_{S}^{0}π^{+}η is measured to be (1.27±0.04_{stat}±0.03_{syst})%. The branching fractions of intermediate processes D^{+}→K_{S}^{0}a_{0}(980)^{+} with a_{0}(980)^{+}→π^{+}η and D^{+}→π^{+}K[over ¯]_{0}^{*}(1430)^{0} with K[over ¯]_{0}^{*}(1430)^{0}→K_{S}^{0}η are measured to be (1.33±0.05_{stat}±0.04_{syst})% and (0.14±0.03_{stat}±0.01_{syst})%, respectively.

Observation of Significant Flavor-SU(3) Breaking in the Kaon Wave Function at 12<Q^{2}<25 GeV^{2} and Discovery of the Charmless Decay ψ(3770)→K_{S}^{0}K_{L}^{0}

We present cross sections for the reaction e^{+}e^{-}→K_{S}^{0}K_{L}^{0} at center-of-mass energies ranging from 3.51 to 4.95 GeV using data samples collected in the BESIII experiment, corresponding to a total integrated luminosity of 26.5  fb^{-1}. The ratio of neutral-to-charged kaon form factors at large momentum transfers (12 Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

The dynamic-process characterization and prediction of synthetic gene circuits by dynamic delay model

Differential equation models are widely used to describe genetic regulations, predict multicomponent regulatory circuits, and provide quantitative insights. However, it is still challenging to quantitatively link the dynamic behaviors with measured parameters in synthetic circuits. Here, we propose a dynamic delay model (DDM) which includes two simple parts: the dynamic determining part and the doses-related steady-state-determining part. The dynamic determining part is usually supposed as the delay time but without a clear formula. For the first time, we give the detail formula of the dynamic determining function and provide a method for measuring all parameters of synthetic elements (include 8 activators and 5 repressors) by microfluidic system. Three synthetic circuits were built to show that the DDM can notably improve the prediction accuracy and can be used in various synthetic biology applications.

Investigation of the ΔI=1/2 Rule and Test of CP Symmetry through the Measurement of Decay Asymmetry Parameters in Ξ^{-} Decays

Using (10087±44)×10^{6}  J/ψ events collected with the BESIII detector, numerous Ξ^{-} and Λ decay asymmetry parameters are simultaneously determined from the process J/ψ→Ξ^{-}Ξ[over ¯]^{+}→Λ(pπ^{-})π^{-}Λ[over ¯](n[over ¯]π^{0})π^{+} and its charge-conjugate channel. The precisions of α_{Λ0} for Λ→nπ^{0} and α[over ¯]_{Λ0} for Λ[over ¯]→n[over ¯]π^{0} compared to world averages are improved by factors of 4 and 1.7, respectively. The ratio of decay asymmetry parameters of Λ→nπ^{0} to that of Λ→pπ^{-}, ⟨α_{Λ0}⟩/⟨α_{Λ-}⟩, is determined to be 0.873±0.012_{-0.010}^{+0.011}, where the first and the second uncertainties are statistical and systematic, respectively. The ratio is smaller than unity more than 5σ, which signifies the existence of the ΔI=3/2 transition in Λ for the first time. Besides, we test for CP symmetry in Ξ^{-}→Λπ^{-} and in Λ→nπ^{0} with the best precision to date.

Observation of D_{s}^{+}→η^{'}μ^{+}ν_{μ}, Precision Test of Lepton Flavor Universality with D_{s}^{+}→η^{(')}l^{+}ν_{l}, and First Measurements of D_{s}^{+}→η^{(')}μ^{+}ν_{μ} Decay Dynamics

By analyzing 7.33  fb^{-1} of e^{+}e^{-} annihilation data collected at center-of-mass energies between 4.128 and 4.226 GeV with the BESIII detector, we report the observation of the semileptonic decay D_{s}^{+}→η^{'}μ^{+}ν_{μ}, with a statistical significance larger than 10σ, and the measurements of the D_{s}^{+}→ημ^{+}ν_{μ} and D_{s}^{+}→η^{'}μ^{+}ν_{μ} decay dynamics for the first time. The branching fractions of D_{s}^{+}→ημ^{+}ν_{μ} and D_{s}^{+}→η^{'}μ^{+}ν_{μ} are determined to be (2.235±0.051_{stat}±0.052_{syst})% and (0.801±0.055_{stat}±0.028_{syst})%, respectively, with precision improved by factors of 6.0 and 6.6 compared to the previous best measurements. Combined with the results for the decays D_{s}^{+}→ηe^{+}ν_{e} and D_{s}^{+}→η^{'}e^{+}ν_{e}, the ratios of the decay widths are examined both inclusively and in several ℓ^{+}ν_{ℓ} four-momentum transfer ranges. No evidence for lepton flavor universality violation is found within the current statistics. The products of the hadronic form factors f_{+,0}^{η^{(')}}(0) and the c→s Cabibbo-Kobayashi-Maskawa matrix element |V_{cs}| are determined. The results based on the two-parameter series expansion are f_{+,0}^{η}(0)|V_{cs}|=0.452±0.010_{stat}±0.007_{syst} and f_{+,0}^{η^{'}}(0)|V_{cs}|=0.504±0.037_{stat}±0.012_{syst}, which help to constrain present models on f_{+,0}^{η^{(')}}(0). The forward-backward asymmetries are determined to be ⟨A_{FB}^{η}⟩=-0.059±0.031_{stat}±0.005_{syst} and ⟨A_{FB}^{η^{'}}⟩=-0.064±0.079_{stat}±0.006_{syst} for the first time, which are consistent with the theoretical calculation.

Determination of the Σ^{+} Timelike Electromagnetic Form Factors

Based on data samples collected with the BESIII detector at the BEPCII collider, the process e^{+}e^{-}→Σ^{+}Σ[over ¯]^{-} is studied at center-of-mass energies sqrt[s]=2.3960, 2.6454, and 2.9000 GeV. Using a fully differential angular description of the final state particles, both the relative magnitude and phase information of the Σ^{+} electromagnetic form factors in the timelike region are extracted. The relative phase between the electric and magnetic form factors is determined to be sinΔΦ=-0.67±0.29(stat)±0.18(syst) at sqrt[s]=2.3960  GeV, ΔΦ=55°±19°(stat)±14°(syst) at sqrt[s]=2.6454  GeV, and 78°±22°(stat)±9°(syst) at sqrt[s]=2.9000  GeV. For the first time, the phase of the hyperon electromagnetic form factors is explored in a wide range of four-momentum transfer. The evolution of the phase along with four-momentum transfer is an important input for understanding its asymptotic behavior and the dynamics of baryons.

Landscape of T cell epitopes displays hot mutations of SARS-CoV-2 variant spikes evading cellular immunity

The continuous evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been accompanied by the emergence of viral mutations that pose a great challenge to existing vaccine strategies. It is not fully understood with regard to the role of mutations on the SARS-CoV-2 spike protein from emerging viral variants in T cell immunity. In the current study, recombinant eukaryotic plasmids were constructed as DNA vaccines to express the spike protein from multiple SARS-CoV-2 strains. These DNA vaccines were used to immunize BALB/c mice, and cross-T cell responses to the spike protein from these viral strains were quantitated using interferon-γ (IFN-γ) Elispot. Peptides covering the full-length spike protein from different viral strains were used to detect epitope-specific IFN-γ+ CD4+ and CD8+ T cell responses by fluorescence-activated cell sorting. SARS-CoV-2 Delta and Omicron BA.1 strains were found to have broad T cell cross-reactivity, followed by the Beta strain. The landscapes of T cell epitopes on the spike protein demonstrated that at least 30 mutations emerging from Alpha to Omicron BA.5 can mediate the escape of T cell immunity. Omicron and its sublineages have 19 out of these 30 mutations, most of which are new, and a few are inherited from ancient circulating variants of concerns. The cross-T cell immunity between SARS-CoV-2 prototype strain and Omicron strains can be attributed to the T cell epitopes located in the N-terminal domain (181-246 aa [amino acids], 271-318 aa) and C-terminal domain (1171-1273 aa) of the spike protein. These findings provide in vivo evidence for optimizing vaccine manufacturing and immunization strategies for current or future viral variants.

First Measurement of the Decay Asymmetry in the Pure W-Boson-Exchange Decay Λ_{c}^{+}→Ξ^{0}K^{+}

Based on 4.4  fb^{-1} of e^{+}e^{-} annihilation data collected at the center-of-mass energies between 4.60 and 4.70 GeV with the BESIII detector at the BEPCII collider, the pure W-boson-exchange decay Λ_{c}^{+}→Ξ^{0}K^{+} is studied with a full angular analysis. The corresponding decay asymmetry is measured for the first time to be α_{Ξ^{0}K^{+}}=0.01±0.16(stat)±0.03(syst). This result reflects the noninterference effect between the S- and P-wave amplitudes. The phase shift between S- and P-wave amplitudes has two solutions, which are δ_{p}-δ_{s}=-1.55±0.25(stat)±0.05(syst)  rad or 1.59±0.25(stat)±0.05(syst)  rad.

Corrigendum to "Improving cooperativity of transcription activators by oligomerization domains in mammalian cells" [Synth Syst Biotechnol 8 (1) (2023) 114-120]

[This corrects the article DOI: 10.1016/j.synbio.2022.12.003.].

Test of CP Symmetry in Hyperon to Neutron Decays

The quantum entangled J/ψ→Σ^{+}Σ[over ¯]^{-} pairs from (1.0087±0.0044)×10^{10}  J/ψ events taken by the BESIII detector are used to study the nonleptonic two-body weak decays Σ^{+}→nπ^{+} and Σ[over ¯]^{-}→n[over ¯]π^{-}. The CP-odd weak decay parameters of the decays Σ^{+}→nπ^{+} (α_{+}) and Σ[over ¯]^{-}→n[over ¯]π^{-} (α[over ¯]_{-}) are determined to be 0.0481±0.0031_{stat}±0.0019_{syst} and -0.0565±0.0047_{stat}±0.0022_{syst}, respectively. The decay parameter α[over ¯]_{-} is measured for the first time, and the accuracy of α_{+} is improved by a factor of 4 compared to the previous results. The simultaneously determined decay parameters allow the first precision CP symmetry test for any hyperon decay with a neutron in the final state with the measurement of A_{CP}=(α_{+}+α[over ¯]_{-})/(α_{+}-α[over ¯]_{-})=-0.080±0.052_{stat}±0.028_{syst}. Assuming CP conservation, the average decay parameter is determined as ⟨α_{+}⟩=(α_{+}-α[over ¯]_{-})/2=-0.0506±0.0026_{stat}±0.0019_{syst}, while the ratios α_{+}/α_{0} and α[over ¯]_{-}/α[over ¯]_{0} are -0.0490±0.0032_{stat}±0.0021_{syst} and -0.0571±0.0053_{stat}±0.0032_{syst}, where α_{0} and α[over ¯]_{0} are the decay parameters of the decays Σ^{+}→pπ^{0} and Σ[over ¯]^{-}→p[over ¯]π^{0}, respectively.

Measurement of Energy-Dependent Pair-Production Cross Section and Electromagnetic Form Factors of a Charmed Baryon

We study the process e^{+}e^{-}→Λ_{c}^{+}Λ[over ¯]_{c}^{-} at twelve center-of-mass energies from 4.6119 to 4.9509 GeV using data samples collected by the BESIII detector at the BEPCII collider. The Born cross sections and effective form factors (|G_{eff}|) are determined with unprecedented precision after combining the single and double-tag methods based on the decay process Λ_{c}^{+}→pK^{-}π^{+}. Flat cross sections around 4.63 GeV are obtained and no indication of the resonant structure Y(4630), as reported by Belle, is found. In addition, no oscillatory behavior is discerned in the |G_{eff}| energy dependence of Λ_{c}^{+}, in contrast to what is seen for the proton and neutron cases. Analyzing the cross section together with the polar-angle distribution of the Λ_{c}^{+} baryon at each energy point, the moduli of electric and magnetic form factors (|G_{E}| and |G_{M}|) are extracted and separated. For the first time, the energy dependence of the form factor ratio |G_{E}/G_{M}| is observed, which can be well described by an oscillatory function.

Determination of Spin-Parity Quantum Numbers for the Narrow Structure near the pΛ[over ¯] Threshold in e^{+}e^{-}→pK^{-}Λ[over ¯]+c.c

A narrow structure in the pΛ[over ¯] system near the mass threshold, named as X(2085), is observed in the process e^{+}e^{-}→pK^{-}Λ[over ¯] with a statistical significance greater than 20σ. Its spin and parity are determined for the first time to be J^{P}=1^{+} in an amplitude analysis, with a statistical significance greater than 5σ over other quantum numbers (0^{-},1^{-} and 2^{+}). The pole positions of X(2085) are measured to be M_{pole}=(2084_{-2}^{+4}±9)  MeV and Γ_{pole}=(58_{-3}^{+4}±25)  MeV, where the first uncertainties are statistical and the second ones are systematic. The analysis is based on the study of the process e^{+}e^{-}→pK^{-}Λ[over ¯] and uses the data samples collected with the BESIII detector at the center-of-mass energies sqrt[s]=4.008, 4.178, 4.226, 4.258, 4.416, and 4.682 GeV with a total integrated luminosity of 8.35  fb^{-1}.

Precise Measurement of the e^{+}e^{-}→D_{s}^{*+}D_{s}^{*-} Cross Sections at Center-of-Mass Energies from Threshold to 4.95 GeV

The process e^{+}e^{-}→D_{s}^{*+}D_{s}^{*-} is studied with a semi-inclusive method using data samples at center-of-mass energies from threshold to 4.95 GeV collected with the BESIII detector operating at the Beijing Electron Positron Collider. The Born cross sections of the process are measured for the first time with high precision in this energy region. Two resonance structures are observed in the energy-dependent cross sections around 4.2 and 4.4 GeV. By fitting the cross sections with a coherent sum of three Breit-Wigner amplitudes and one phase-space amplitude, the two significant structures are assigned masses of (4186.8±8.7±30) and (4414.6±3.4±6.1)  MeV/c^{2}, widths of (55±15±53) and (122.5±7.5±8.1)  MeV, where the first errors are statistical and the second ones are systematic. The inclusion of a third Breit-Wigner amplitude is necessary to describe a structure around 4.79 GeV.

First Experimental Study of the Purely Leptonic Decay D_{s}^{*+}→e^{+}ν_{e}

Using 7.33  fb^{-1} of e^{+}e^{-} collision data taken with the BESIII detector at the BEPCII collider, we report the first experimental study of the purely leptonic decay D_{s}^{*+}→e^{+}ν_{e}. Our data contain a signal of this decay with a statistical significance of 2.9σ. The branching fraction of D_{s}^{*+}→e^{+}ν_{e} is measured to be (2.1_{-0.9_{stat}}^{+1.2}±0.2_{syst})×10^{-5}, corresponding to an upper limit of 4.0×10^{-5} at the 90% confidence level. Taking the total width of the D_{s}^{*+} [(0.070±0.028)  keV] predicted with the radiative D_{s}^{*+} decay from the lattice QCD calculation as input, the decay constant of the D_{s}^{*+} is determined to be f_{D_{s}^{*+}}=(214_{-46_{stat}}^{+61}±44_{syst})  MeV, corresponding to an upper limit of 354 MeV at the 90% confidence level.

Search for Λ[over ¯]-Λ Baryon-Number-Violating Oscillations in the Decay J/ψ→pK^{-}Λ[over ¯]+c.c

We report on the first search for Λ[over ¯]-Λ oscillations in the decay J/ψ→pK^{-}Λ[over ¯]+c.c. by analyzing 1.31×10^{9}  J/ψ events accumulated with the BESIII detector at the BEPCII collider. The J/ψ events are produced using e^{+}e^{-} collisions at a center of mass energy sqrt[s]=3.097  GeV. No evidence for hyperon oscillations is observed. The upper limit for the oscillation rate of Λ[over ¯] to Λ hyperons is determined to be P(Λ)=[B(J/ψ→pK^{-}Λ+c.c.)/B(J/ψ→pK^{-}Λ[over ¯]+c.c.)]<4.4×10^{-6} corresponding to an oscillation parameter δm_{ΛΛ[over ¯]} of less than 3.8×10^{-18}  GeV at the 90% confidence level.

Determining subunit-subunit interaction from statistics of cryo-EM images: observation of nearest-neighbor coupling in a circadian clock protein complex

Biological processes are typically actuated by dynamic multi-subunit molecular complexes. However, interactions between subunits, which govern the functions of these complexes, are hard to measure directly. Here, we develop a general approach combining cryo-EM imaging technology and statistical modeling and apply it to study the hexameric clock protein KaiC in Cyanobacteria. By clustering millions of KaiC monomer images, we identify two major conformational states of KaiC monomers. We then classify the conformational states of (>160,000) KaiC hexamers by the thirteen distinct spatial arrangements of these two subunit states in the hexamer ring. We find that distributions of the thirteen hexamer conformational patterns for two KaiC phosphorylation mutants can be fitted quantitatively by an Ising model, which reveals a significant cooperativity between neighboring subunits with phosphorylation shifting the probability of subunit conformation. Our results show that a KaiC hexamer can respond in a switch-like manner to changes in its phosphorylation level.

Dynamical modelling of viral infection and cooperative immune protection in COVID-19 patients

Once challenged by the SARS-CoV-2 virus, the human host immune system triggers a dynamic process against infection. We constructed a mathematical model to describe host innate and adaptive immune response to viral challenge. Based on the dynamic properties of viral load and immune response, we classified the resulting dynamics into four modes, reflecting increasing severity of COVID-19 disease. We found the numerical product of immune system's ability to clear the virus and to kill the infected cells, namely immune efficacy, to be predictive of disease severity. We also investigated vaccine-induced protection against SARS-CoV-2 infection. Results suggested that immune efficacy based on memory T cells and neutralizing antibody titers could be used to predict population vaccine protection rates. Finally, we analyzed infection dynamics of SARS-CoV-2 variants within the construct of our mathematical model. Overall, our results provide a systematic framework for understanding the dynamics of host response upon challenge by SARS-CoV-2 infection, and this framework can be used to predict vaccine protection and perform clinical diagnosis.

Molecular mechanism for activation of the 26S proteasome by ZFAND5

Various hormones, kinases, and stressors (fasting, heat shock) stimulate 26S proteasome activity. To understand how its capacity to degrade ubiquitylated proteins can increase, we studied mouse ZFAND5, which promotes protein degradation during muscle atrophy. Cryo-electron microscopy showed that ZFAND5 induces large conformational changes in the 19S regulatory particle. ZFAND5's AN1 Zn-finger domain interacts with the Rpt5 ATPase and its C terminus with Rpt1 ATPase and Rpn1, a ubiquitin-binding subunit. Upon proteasome binding, ZFAND5 widens the entrance of the substrate translocation channel, yet it associates only transiently with the proteasome. Dissociation of ZFAND5 then stimulates opening of the 20S proteasome gate. Using single-molecule microscopy, we showed that ZFAND5 binds ubiquitylated substrates, prolongs their association with proteasomes, and increases the likelihood that bound substrates undergo degradation, even though ZFAND5 dissociates before substrate deubiquitylation. These changes in proteasome conformation and reaction cycle can explain the accelerated degradation and suggest how other proteasome activators may stimulate proteolysis.

Free energy dissipation enhances spatial accuracy and robustness of self-positioned Turing pattern in small biochemical systems

Accurate and robust spatial orders are ubiquitous in living systems. In 1952, Turing proposed a general mechanism for pattern formation exemplified by a reaction-diffusion model with two chemical species in a large system. However, in small biological systems such as a cell, the existence of multiple Turing patterns and strong noise can lower the spatial order. Recently, a modified reaction-diffusion model with an additional chemical species is shown to stabilize the Turing pattern. Here, we study non-equilibrium thermodynamics of this three-species reaction-diffusion model to understand the relationship between energy cost and the performance of self-positioning. By using computational and analytical approaches, we show that beyond the onset of pattern formation the positioning error decreases as energy dissipation increases. In a finite system, we find that a specific Turing pattern exists only within a finite range of total molecule number. Energy dissipation broadens this range, which enhances the robustness of Turing pattern against molecule number fluctuations in living cells. The generality of these results is verified in a realistic model of the Muk system underlying DNA segregation in Escherichia coli, and testable predictions are made for the dependence of the accuracy and robustness of the spatial pattern on the ATP/ADP ratio.

First Study of Reaction Ξ^{0}n→Ξ^{-}p Using Ξ^{0}-Nucleus Scattering at an Electron-Positron Collider

Using (1.0087±0.0044)×10^{10}  J/ψ events collected with the BESIII detector at the BEPCII storage ring, the process Ξ^{0}n→Ξ^{-}p is studied, where the Ξ^{0} baryon is produced in the process J/ψ→Ξ^{0}Ξ[over ¯]^{0} and the neutron is a component of the ^{9}Be, ^{12}C, and ^{197}Au nuclei in the beam pipe. A clear signal is observed with a statistical significance of 7.1σ. The cross section of the reaction Ξ^{0}+^{9}Be→Ξ^{-}+p+^{8}Be is determined to be σ(Ξ^{0}+^{9}Be→Ξ^{-}+p+^{8}Be)=(22.1±5.3_{stat}±4.5_{sys})  mb at the Ξ^{0} momentum of 0.818  GeV/c, where the first uncertainty is statistical and the second is systematic. No significant H-dibaryon signal is observed in the Ξ^{-}p final state. This is the first study of hyperon-nucleon interactions in electron-positron collisions and opens up a new direction for such research.

Measurements of Normalized Differential Cross Sections of Inclusive π^{0} and K_{S}^{0} Production in e^{+}e^{-} Annihilation at Energies from 2.2324 to 3.6710 GeV

Based on electron positron collision data collected with the BESIII detector operating at the BEPCII storage rings, the differential cross sections of inclusive π^{0} and K_{S}^{0} production as a function of hadron momentum, normalized by the total cross section of the e^{+}e^{-}→hadrons process, are measured at six center-of-mass energies from 2.2324 to 3.6710 GeV. Our results, which cover a relative hadron energy range from 0.1 to 0.9, significantly deviate from several theoretical calculations based on existing fragmentation functions.

A randomized trial of eribulin monotherapy versus eribulin plus anlotinib in patients with locally recurrent or metastatic breast cancer

BACKGROUND

Eribulin mesylate is a novel, nontaxane, microtubule dynamics inhibitor. In this study, we assessed the efficacy and safety of eribulin versus eribulin plus the oral small-molecule tyrosine kinase inhibitor anlotinib in patients with locally recurrent or metastatic breast cancer.

PATIENTS AND METHODS

In this single-center, open-label, phase II clinical study (NCT05206656) conducted in a Chinese hospital, patients with human epidermal growth factor receptor 2 (HER2)-negative, locally recurrent or metastatic breast cancer previously treated with anthracycline- or taxane-based chemotherapy were randomized (1 : 1) to receive eribulin alone or in combination with anlotinib. The primary efficacy endpoint was investigator-assessed progression-free survival (PFS).

RESULTS

From June 2020 to April 2022, a total of 80 patients were randomly assigned to either eribulin monotherapy or eribulin plus anlotinib combination therapy, with 40 patients in each group. The data cut-off was 10 August 2022. The median PFS was 3.5 months [95% confidence interval (CI) 2.8-5.5 months] for eribulin and 5.1 months (95% CI 4.5-6.9 months) for eribulin plus anlotinib (hazard ratio = 0.56, 95% CI 0.32-0.98; P = 0.04). The objective response rates were 32.5% versus 52.5% (P = 0.07), respectively, and disease control rates were 67.5% versus 92.5% (P = 0.01), respectively. Patients <50 years of age, with an Eastern Cooperative Oncology Group performance status score of 0, visceral metastasis, number of treatment lines of four or more, hormone receptor negative (triple-negative), and HER2 low expression appeared to benefit more from combined treatment. The most common adverse events in both groups were leukopenia (n = 28, 70.0%, patients in the eribulin monotherapy group versus n = 35, 87.5%, patients in the combination therapy group), aspartate aminotransferase elevations (n = 28, 70.0%, versus n = 35, 87.5%), neutropenia (n = 25, 62.5%, versus n = 31, 77.5%), and alanine aminotransferase elevations (n = 25, 62.5%, versus n = 30, 75.0%).

CONCLUSION

Eribulin plus anlotinib can be considered an alternative treatment option for HER2-negative locally advanced or metastatic breast cancer.

Precision Measurement of the Decay Σ^{+}→pγ in the Process J/ψ→Σ^{+}Σ[over ¯]^{-}

Using (10 087±44)×10^{6} J/ψ events collected with the BESIII detector, the radiative hyperon decay Σ^{+}→pγ is studied at an electron-positron collider experiment for the first time. The absolute branching fraction is measured to be (0.996±0.021_{stat}±0.018_{syst})×10^{-3}, which is lower than its world average value by 4.2 standard deviations. Its decay asymmetry parameter is determined to be -0.652±0.056_{stat}±0.020_{syst}. The branching fraction and decay asymmetry parameter are the most precise to date, and the accuracies are improved by 78% and 34%, respectively.

Coupled oscillator cooperativity as a control mechanism in chronobiology

Control of dynamical processes is vital for maintaining correct cell regulation and cell-fate decisions. Numerous regulatory networks show oscillatory behavior; however, our knowledge of how one oscillator behaves when stimulated by two or more external oscillatory signals is still missing. We explore this problem by constructing a synthetic oscillatory system in yeast and stimulate it with two external oscillatory signals. Letting model verification and prediction operate in a tight interplay with experimental observations, we find that stimulation with two external signals expands the plateau of entrainment and reduces the fluctuations of oscillations. Furthermore, by adjusting the phase differences of external signals, one can control the amplitude of oscillations, which is understood through the signal delay of the unperturbed oscillatory network. With this we reveal a direct amplitude dependency of downstream gene transcription. Taken together, these results suggest a new path to control oscillatory systems by coupled oscillator cooperativity.

Molecular mechanisms for activation of the 26S proteasome

Various hormones, kinases, and stressors (fasting, heat shock) stimulate 26S proteasome activity. To understand how its capacity to degrade ubiquitylated protein can increase, we studied ZFAND5, which promotes protein degradation during muscle atrophy. Cryo-electron microscopy showed that ZFAND5 induces large conformational changes in the 19S regulatory particle. ZFAND5's AN1 Zn finger interacts with the Rpt5 ATPase and its C-terminus with Rpt1 ATPase and Rpn1, a ubiquitin-binding subunit. Surprisingly, these C-terminal interactions are sufficient to activate proteolysis. With ZFAND5 bound, entry into the proteasome's protein translocation channel is wider, and ZFAND5 dissociation causes opening of the 20S gate for substrate entry. Using single-molecular microscopy, we showed that ZFAND5 binds ubiquitylated substrates, prolongs their association with proteasomes, and increases the likelihood that bound substrates undergo degradation, even though ZFAND5 dissociates before substrate deubiquitylation. These changes in proteasome conformation and reaction cycle can explain the accelerated degradation and suggest how other proteasome activators may stimulate proteolysis.

Plant-on-chip: Core morphogenesis processes in the tiny plant Wolffia australiana

A plant can be thought of as a colony comprising numerous growth buds, each developing to its own rhythm. Such lack of synchrony impedes efforts to describe core principles of plant morphogenesis, dissect the underlying mechanisms, and identify regulators. Here, we use the minimalist known angiosperm to overcome this challenge and provide a model system for plant morphogenesis. We present a detailed morphological description of the monocot Wolffia australiana, as well as high-quality genome information. Further, we developed the plant-on-chip culture system and demonstrate the application of advanced technologies such as single-nucleus RNA-sequencing, protein structure prediction, and gene editing. We provide proof-of-concept examples that illustrate how W. australiana can decipher the core regulatory mechanisms of plant morphogenesis.

A 3D Scalable Chamber-Array Chip for Digital LAMP

As an absolute quantification method at the single-molecule level, digital PCR (dPCR) offers the highest accuracy. In this work, we developed a 3D scalable chamber-array chip that multiplied the number of partitions by stacking chamber-array layers and realized digital loop-mediated isothermal amplification to quantify DNA molecules. It greatly increases the number of partitions to improve the performance of dPCR without increasing the chip size, the operation workflow complicity, and operation time. For the three-chamber-array-layer chip which contains 200,000 reactors of a 0.125 nL volume, it has been proved that the reagent filling and partition were finished within 3 min, and the whole detection could be finished within 1 h. The method demonstrated that it could be scalable to a six-chamber-array layer, which contains 400,000 reactors without increasing the size of the chip and the complication of filling/partition workflow but only takes an additional hour for scanning. Due to its potential for high throughput, low cost, and simple operation, our device may significantly expand the clinical application range of dPCR.

Structure of the Spring Viraemia of Carp Virus Ribonucleoprotein Complex Reveals Its Assembly Mechanism and Application in Antiviral Drug Screening

Spring viremia of carp virus (SVCV) is a highly pathogenic Vesiculovirus infecting the common carp, yet neither a vaccine nor effective therapies are available to treat spring viremia of carp (SVC). Like all negative-sense viruses, SVCV contains an RNA genome that is encapsidated by the nucleoprotein (N) in the form of a ribonucleoprotein (RNP) complex, which serves as the template for viral replication and transcription. Here, the three-dimensional (3D) structure of SVCV RNP was resolved through cryo-electron microscopy (cryo-EM) at a resolution of 3.7 Å. RNP assembly was stabilized by N and C loops; RNA was wrapped in the groove between the N and C lobes with 9 nt nucleotide per protomer. Combined with mutational analysis, our results elucidated the mechanism of RNP formation. The RNA binding groove of SVCV N was used as a target for drug virtual screening, and it was found suramin had a good antiviral effect. This study provided insights into RNP assembly, and anti-SVCV drug screening was performed on the basis of this structure, providing a theoretical basis and efficient drug screening method for the prevention and treatment of SVC. IMPORTANCE Aquaculture accounts for about 70% of global aquatic products, and viral diseases severely harm the development of aquaculture industry. Spring viremia of carp virus (SVCV) is the pathogen causing highly contagious spring viremia of carp (SVC) disease in cyprinids, especially common carp (Cyprinus carpio), yet neither a vaccine nor effective therapies are available to treat this disease. In this study, we have elucidated the mechanism of SVCV ribonucleoprotein complex (RNP) formation by resolving the 3D structure of SVCV RNP and screened antiviral drugs based on the structure. It is found that suramin could competitively bind to the RNA binding groove and has good antiviral effects both in vivo and in vitro. Our study provides a template for rational drug discovery efforts to treat and prevent SVCV infections.

Measurements of the Electric and Magnetic Form Factors of the Neutron for Timelike Momentum Transfer

We present the first measurements of the electric and magnetic form factors of the neutron in the timelike (positive q^{2}) region as function of four-momentum transfer. We explored the differential cross sections of the reaction e^{+}e^{-}→n[over ¯]n with data collected with the BESIII detector at the BEPCII accelerator, corresponding to an integrated luminosity of 354.6  pb^{-1} in total at twelve center-of-mass energies between sqrt[s]=2.0-2.95  GeV. A relative uncertainty of 18% and 12% for the electric and magnetic form factors, respectively, is achieved at sqrt[s]=2.3935  GeV. Our results are comparable in accuracy to those from electron scattering in the comparable spacelike region of four-momentum transfer. The electromagnetic form factor ratio R_{em}≡|G_{E}|/|G_{M}| is within the uncertainties close to unity. We compare our result on |G_{E}| and |G_{M}| to recent model predictions, and the measurements in the spacelike region to test the analyticity of electromagnetic form factors.

Observation of a New X(3872) Production Process e^{+}e^{-}→ωX(3872)

Using 4.7  fb^{-1} of e^{+}e^{-} collision data at center-of-mass energies from 4.661 to 4.951 GeV collected by the BESIII detector at the BEPCII collider, we observe the X(3872) production process e^{+}e^{-}→ωX(3872) for the first time. The significance is 7.8σ, including both the statistical and systematic uncertainties. The e^{+}e^{-}→ωX(3872) Born cross section and the corresponding upper limit at 90% confidence level at each energy point are reported. The line shape of the cross section indicates that the ωX(3872) signals may be from the decays of some nontrivial structures.

Mobile and Self-Sustained Data Storage in an Extremophile Genomic DNA

DNA has been pursued as a novel biomaterial for digital data storage. While large-scale data storage and random access have been achieved in DNA oligonucleotide pools, repeated data accessing requires constant data replenishment, and these implementations are confined in professional facilities. Here, a mobile data storage system in the genome of the extremophile Halomonas bluephagenesis, which enables dual-mode storage, dynamic data maintenance, rapid readout, and robust recovery. The system relies on two key components: A versatile genetic toolbox for the integration of 10-100 kb scale synthetic DNA into H. bluephagenesis genome and an efficient error correction coding scheme targeting noisy nanopore sequencing reads. The storage and repeated retrieval of 5 KB data under non-laboratory conditions are demonstrated. The work highlights the potential of DNA data storage in domestic and field scenarios, and expands its application domain from archival data to frequently accessed data.

Observation of Three Charmoniumlike States with J^{PC}=1^{--} in e^{+}e^{-}→D^{*0}D^{*-}π^{+}

The Born cross sections of the process e^{+}e^{-}→D^{*0}D^{*-}π^{+} at center-of-mass energies from 4.189 to 4.951 GeV are measured for the first time. The data samples used correspond to an integrated luminosity of 17.9  fb^{-1} and were collected by the BESIII detector operating at the BEPCII storage ring. Three enhancements around 4.20, 4.47, and 4.67 GeV are visible. The resonances have masses of 4209.6±4.7±5.9  MeV/c^{2}, 4469.1±26.2±3.6  MeV/c^{2}, and 4675.3±29.5±3.5  MeV/c^{2} and widths of 81.6±17.8±9.0  MeV, 246.3±36.7±9.4  MeV, and 218.3±72.9±9.3  MeV, respectively, where the first uncertainties are statistical and the second systematic. The first and third resonances are consistent with the ψ(4230) and ψ(4660) states, respectively, while the second one is compatible with the ψ(4500) observed in the e^{+}e^{-}→K^{+}K^{-}J/ψ process. These three charmoniumlike ψ states are observed in the e^{+}e^{-}→D^{*0}D^{*-}π^{+} process for the first time.

Precise programming of multigene expression stoichiometry in mammalian cells by a modular and programmable transcriptional system

Context-dependency of mammalian transcriptional elements has hindered the quantitative investigation of multigene expression stoichiometry and its biological functions. Here, we describe a host- and local DNA context-independent transcription system to gradually fine-tune single and multiple gene expression with predictable stoichiometries. The mammalian transcription system is composed of a library of modular and programmable promoters from bacteriophage and its cognate RNA polymerase (RNAP) fused to a capping enzyme. The relative expression of single genes is quantitatively determined by the relative binding affinity of the RNAP to the promoters, while multigene expression stoichiometry is predicted by a simple biochemical model with resource competition. We use these programmable and modular promoters to predictably tune the expression of three components of an influenza A virus-like particle (VLP). Optimized stoichiometry leads to a 2-fold yield of intact VLP complexes. The host-independent orthogonal transcription system provides a platform for dose-dependent control of multiple protein expression which may be applied for advanced vaccine engineering, cell-fate programming and other therapeutic applications.

In vitro activity of tetracycline analogs against multidrug-resistant and extensive drug resistance clinical isolates of Mycobacterium tuberculosis

BACKGROUND

Multidrug-resistant tuberculosis (MDR-TB) has become a big threaten to global health. The current strategy for treatment of MDR-TB and extensive drug resistant tuberculosis (XDR-TB) is with low efficacy and high side effect. While new drug is fundamental for cure MDR-TB, repurposing the Food and Drug Administration (FDA)-approved drugs represents an alternative soluation with less cost.

METHODS

The activity of 8 tetracycline-class antibiotics against mycobacterium tuberculosis (M.tb) were determined by Minimum Inhibitory Concentration (MIC) in vitro. A transposon M.smeg libraries was generated by using the Harm phage and then used to isolate the conditional growth mutants in doxycycline containing plate. Eleven mutants were isolated and genomic DNAs were extracted using the cetyltrimethyl ammonium bromide (CTAB) method and analyzed by whole genome sequencing.

RESULTS

We found that three of eight drugs efficiently inhibited mycobacteria growth under the peak plasma concentration in the human body. Further tests showed these three tetracycline analogs (demeclocycline, doxycycline and methacycline) had antimicrobial activity against seven clinical isolates, including MDR and XDR strains. Among them, Doxycycline had the lowest MICs in all mycobacteria strains tested in this study. By using a transposon library, we identify the insertion of transposon in two genes, porin and MshA, associatewith the resistant to doxycycline.

CONCLUSIONS

Our findings show that tetracycline analogs such as doxycycline, has bactericidal activity against not only drug sensitive M.tb, but also clinical MDR and XDR strains, provided proof of concept to repurpose doxycycline to fight MDR-TB and XDR-TB. Further investigations are warranted to clarify the underlying mechanism and optimize the strategy in combination with other anti-TB drugs.

Improving cooperativity of transcription activators by oligomerization domains in mammalian cells

Cooperative activation is critical for the applications of synthetic biology in mammalian cells. In this study, we have developed cooperative transcription factor by fusing oligomerization domain in mammalian cells. Firstly, we demonstrated that two oligomerized domains (CI434 and CI) successfully improved transcription factor cooperativity in bacterial cells but failed to increase cooperativity in mammalian cells, possibly because the additional mammalian activation domain disrupted their oligomerization capability. Therefore, we chose a different type of oligomerized domain (CarHC), whose ability to oligomerize is not dependent on its C-terminal domains, to fuse with a transcription factor (RpaR) and activation domain (VTR3), forming a potential cooperative transcription activator RpaR-CarH-VTR3 for mammalian regulatory systems. Compared with RpaR-VTR3, the cooperativity of RpaR-CarH-VTR3 was significantly improved with higher Hill coefficient and a narrower input range in the inducible switch system in mammalian cells. Moreover, a mathematical model based on statistical mechanics model was developed and the simulation results supported the hypothesis that the tetramer of the CarH domain in mammalian cells was the reason for the cooperative capacity of RpaR-CarH-VTR3.

Evidence for the Cusp Effect in η' Decays into ηπ^{0}π^{0}

Using a sample of 4.3×10^{5} η^{'}→ηπ^{0}π^{0} events selected from the ten billion J/ψ event dataset collected with the BESIII detector, we study the decay η^{'}→ηπ^{0}π^{0} within the framework of nonrelativistic effective field theory. Evidence for a structure at π^{+}π^{-} mass threshold is observed in the invariant mass spectrum of π^{0}π^{0} with a statistical significance of around 3.5σ, which is consistent with the cusp effect as predicted by the nonrelativistic effective field theory. After introducing the amplitude for describing the cusp effect, the ππ scattering length combination a_{0}-a_{2} is determined to be 0.226±0.060_{stat}±0.013_{syst}, which is in good agreement with theoretical calculation of 0.2644±0.0051.

Ni, Co, and Yb Cation Co-doping and Defect Engineering of FeOOH Nanorods as an Electrocatalyst for the Oxygen Evolution Reaction

Electrocatalytic water splitting is a feasible technology that can produce hydrogen from renewable sources. The oxygen evolution reaction (OER), which has a slower kinetics and higher overpotential than the hydrogen evolution reaction, is the bottleneck that limits the overall water splitting. It is essential to develop efficient OER catalysts to reduce the anode overpotential. Herein, Ni,Co,Yb-FeOOH nanorod arrays grown directly on a carbon cloth are synthesized by a simple one-step hydrothermal method. The doped Ni²⁺ and Co²⁺ can occupy Fe²⁺ and Fe³⁺ sites in FeOOH, increasing the concentration of oxygen vacancies (V_O), and the doped Yb³⁺ with a larger ionic radius can occupy the interstitial sites, which leads to more edge dislocations. V_O and edge dislocations greatly enrich the active sites in FeOOH/CC. In addition, density functional theory calculations confirm that doping of Ni²⁺, Co²⁺, and Yb³⁺ modulates the electronic structure of the main active Fe sites, bringing its d-band center closer to the Fermi level and reducing the Gibbs free energy change of the rate-determining step of the OER. When the current density reaches 10 mA cm^-2, the overpotential of Ni,Co,Yb-FeOOH/CC is only 230.9 mV, and the Tafel slope is 22.7 mV dec^-1. In particular, a mechanism of multi-cation doping synergistic interaction with the oxygen vacancy and edge dislocation to enhance the OER catalytic activity of the material is proposed.

The Topological Characteristics of Biological Ratio-Sensing Networks

Ratio sensing is a fundamental biological function observed in signal transduction and decision making. In the synthetic biology context, ratio sensing presents one of the elementary functions for cellular multi-signal computation. To investigate the mechanism of the ratio-sensing behavior, we explored the topological characteristics of biological ratio-sensing networks. With exhaustive enumeration of three-node enzymatic and transcriptional regulatory networks, we found that robust ratio sensing was highly dependent on network structure rather than network complexity. Specifically, a set of seven minimal core topological structures and four motifs were deduced to be capable of robust ratio sensing. Further investigations on the evolutionary space of robust ratio-sensing networks revealed highly clustered domains surrounding the core motifs which suggested their evolutionary plausibility. Our study revealed the network topological design principles of ratio-sensing behavior and provided a design scheme for constructing regulatory circuits with ratio-sensing behavior in synthetic biology.

Two-Layered Microfluidic Devices for High-Throughput Dynamic Analysis of Synthetic Gene Circuits in E. coli

Escherichia coli is a common chassis for synthetic gene circuit studies. In addition to the dose-response of synthetic gene circuits, the analysis of dynamic responses is also an important part of the future design of more complicated synthetic systems. Recently, microfluidic-based methods have been widely used for the analysis of gene expression dynamics. Here, we established a two-layered microfluidic platform for the systematic characterization of synthetic gene circuits (eight strains in eight different culture environments could be observed simultaneously with a 5 min time resolution). With this platform, both dose responses and dynamic responses with a high temporal resolution could be easily derived for further analysis. A controlled environment ensures the stability of the bacterial growth rate, excluding changes in gene expression dynamics caused by changes of the growth dilution rate. The precise environmental switch and automatic micrograph shooting ensured that there was nearly no time lag between the inducer addition and the data recording. We studied four four-node incoherent-feedforward-loop (IFFL) networks with different operators using this device. The experimental results showed that as the effect of inhibition increased, two of the IFFL networks generated pulselike dynamic gene expressions in the range of the inducer concentrations, which was different from the dynamics of the two other circuits with only a simple pattern of rising to the platform. Through fitting the dose-response curves and the dynamic response curves, corresponding parameters were derived and introduced to a simple model that could qualitatively explain the generation of pulse dynamics.

Study of the Semileptonic Decay Λ_{c}^{+}→Λe^{+}ν_{e}

The study of the Cabibbo-favored semileptonic decay Λ_{c}^{+}→Λe^{+}ν_{e} is reported using a 4.5  fb^{-1} data sample of e^{+}e^{-} annihilations collected at center-of-mass energies ranging from 4.600 GeV to 4.699 GeV with the BESIII detector at the BEPCII collider. The branching fraction of the decay is measured to be B(Λ_{c}^{+}→Λe^{+}ν_{e})=(3.56±0.11_{stat}±0.07_{syst})%, which is the most precise measurement to date. Furthermore, we perform an investigation of the internal dynamics in Λ_{c}^{+}→Λe^{+}ν_{e}. We provide the first direct comparisons of the differential decay rate and form factors with those predicted from lattice quantum chromodynamics (LQCD) calculations. Combining the measured branching fraction with a q^{2}-integrated rate predicted by LQCD, we determine |V_{cs}|=0.936±0.017_{B}±0.024_{LQCD}±0.007_{τ_{Λ_{c}}}.

Measurement of the Absolute Branching Fraction and Decay Asymmetry of Λ→nγ

The radiative hyperon decay Λ→nγ is studied using (10087±44)×10^{6}  J/ψ events collected with the BESIII detector operating at BEPCII. The absolute branching fraction of the decay Λ→nγ is determined to be (0.832±0.038_{stat}±0.054_{syst})×10^{-3}, which is a factor of 2.1 lower and 5.6 standard deviations different than the previous measurement. By analyzing the joint angular distribution of the decay products, the first determination of the decay asymmetry α_{γ} is reported with a value of -0.16±0.10_{stat}±0.05_{syst}.

Observation of an Isoscalar Resonance with Exotic J^{PC}=1^{-+} Quantum Numbers in J/ψ→γηη^{'}

Using a sample of (10.09±0.04)×10^{9}  J/ψ events collected with the BESIII detector operating at the BEPCII storage ring, a partial wave analysis of the decay J/ψ→γηη^{'} is performed. The first observation of an isoscalar state with exotic quantum numbers J^{PC}=1^{-+}, denoted as η_{1}(1855), is reported in the process J/ψ→γη_{1}(1855) with η_{1}(1855)→ηη^{'}. Its mass and width are measured to be (1855±9_{-1}^{+6})  MeV/c^{2} and (188±18_{-8}^{+3})  MeV, respectively, where the first uncertainties are statistical and the second are systematic, and its statistical significance is estimated to be larger than 19σ.