  Standard epidemiological models for COVID-19 employ variants of compartment
(SIR) models at local scales, implicitly assuming spatially uniform local
mixing. Here, we examine the effect of employing more geographically detailed
diffusion models based on known spatial features of interpersonal networks,
most particularly the presence of a long-tailed but monotone decline in the
probability of interaction with distance, on disease diffusion. Based on
simulations of unrestricted COVID-19 diffusion in 19 U.S cities, we conclude
that heterogeneity in population distribution can have large impacts on local
pandemic timing and severity, even when aggregate behavior at larger scales
mirrors a classic SIR-like pattern. Impacts observed include severe local
outbreaks with long lag time relative to the aggregate infection curve, and the
presence of numerous areas whose disease trajectories correlate poorly with
those of neighboring areas. A simple catchment model for hospital demand
illustrates potential implications for health care utilization, with
substantial disparities in the timing and extremity of impacts even without
distancing interventions. Likewise, analysis of social exposure to others who
are morbid or deceased shows considerable variation in how the epidemic can
appear to individuals on the ground, potentially affecting risk assessment and
compliance with mitigation measures. These results demonstrate the potential
for spatial network structure to generate highly non-uniform diffusion behavior
even at the scale of cities, and suggest the importance of incorporating such
structure when designing models to inform healthcare planning, predict
community outcomes, or identify potential disparities.


  This paper is concerned with a family of Reaction-Diffusion systems that we
introduced in [15], and that generalizes the SIR type models from epidemiology.
Such systems are now also used to describe collective behaviors.In this paper,
we propose a modeling approach for these apparently diverse phenomena through
the example of the dynamics of social unrest. The model involves two
quantities: the level of social unrest, or more generally activity, u, and a
field of social tension v, which play asymmetric roles. We think of u as the
actually observed or explicit quantity while v is an ambiant, sometimes
implicit, field of susceptibility that modulates the dynamics of u. In this
article, we explore this class of model and prove several theoretical results
based on the framework developed in [15], of which the present work is a
companion paper. We particularly emphasize here two subclasses of systems:
tension inhibiting and tension enhancing. These are characterized by
respectively a negative or a positivefeedback of the unrest on social tension.
We establish several properties for these classes and also study some
extensions. In particular, we describe the behavior of the system following an
initial surge of activity. We show that the model can give rise to many diverse
qualitative dynamics. We also provide a variety of numerical simulations to
illustrate our results and to reveal further properties and open questions.


  Though pandemic preparedness has been a focus of public health planning for
centuries, during which our understanding of infectious disease dynamics has
grown, our methodologies for managing outbreaks have remained relatively
unchanged. We propose leveraging this history to identify opportunities for
actual progress. We contrast current plans with historical outbreak control
measures and isolate how the complexities of a modern era yield additional
challenges in how best to anticipate and mitigate outbreaks. We analyze a
diversity of publicly available modern preparedness plans against the context
of a historically-based fictional outbreak control strategy described in
Defoe's A Journal of the Plague Year (published 1720). We identify themes in
preparedness planning that remain unchanged from historical settings even
though they continue to be actively evaluated in planning efforts. More
importantly, we isolate critical modern challenges in preparedness planning
that remain predominantly unsolved. These modern, unsolved issues offer best
avenues for meaningful improvement. Shifting our planning efforts to focus on
identified novel issues may greatly strengthen our local- to global- capacity
to deal with infectious threats.


  The paradigm for compartment models in epidemiology assumes exponentially
distributed incubation and removal times, which is not realistic in actual
populations. Commonly used variations with multiple exponentially distributed
variables are more flexible, yet do not allow for arbitrary distributions. We
present a new formulation, focussing on the SEIR concept that allows to include
general distributions of incubation and removal times. We compare the solution
to two types of agent-based model simulations, a spatially homogeneous one
where infection occurs by proximity, and a model on a scale-free network with
varying clustering properties, where the infection between any two agents
occurs via their link if it exists. We find good agreement in both cases.
Furthermore a family of asymptotic solutions of the equations is found in terms
of a logistic curve, which after a non-universal time shift, fits extremely
well all the microdynamical simulations. The formulation allows for a simple
numerical approach; software in Julia and Python is provided.


  The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute
respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected near 5 million
people and led to over 0.3 million deaths. Currently, there is no specific
anti-SARS-CoV-2 medication. New drug discovery typically takes more than ten
years. Drug repositioning becomes one of the most feasible approaches for
combating COVID-19. This work curates the largest available experimental
dataset for SARS-CoV-2 or SARS-CoV main protease inhibitors. Based on this
dataset, we develop validated machine learning models with relatively low root
mean square error to screen 1553 FDA-approved drugs as well as other 7012
investigational or off-market drugs in DrugBank. We found that many existing
drugs might be potentially potent to SARS-CoV-2. The druggability of many
potent SARS-CoV-2 main protease inhibitors is analyzed. This work offers a
foundation for further experimental studies of COVID-19 drug repositioning.


  In its more severe forms, COVID-19 progresses towards an excessive immune
response, leading to the systemic overexpression of proinflammatory cytokines
like IL6, mostly from the infected lungs. This cytokine storm can cause
multiple organ damage and death. Consequently, there is a pressing need to
identify therapies to treat and prevent severe symptoms during COVID-19. Based
on previous clinical evidence, we hypothesized that inhibiting T cell
co-stimulation by blocking CD80/86 could be an effective therapeutic strategy
against progression to severe proinflammatory states. To support this
hypothesis, we performed an analysis integrating blood transcriptional data we
generated from rheumatoid arthritis patients treated with abatacept -- a
CD80/86 co-stimulation inhibitor -- with the pathological features associated
with COVID-19, particularly in its more severe forms. We have found that many
of the biological processes that have been consistently associated with
COVID-19 pathology are reversed by CD80/86 co-stimulation inhibition, including
the downregulation of IL6 production. Also, analysis of previous
transcriptional data from blood of SARS-CoV-infected patients showed that the
response to abatacept has a very high level of antagonism to that elicited by
COVID-19. Finally, analyzing a recent single cell RNA-seq dataset from
bronchoalveolar lavage fluid cells from COVID-19 patients, we found a
significant correlation along the main elements of the C80/86 axis: CD86+/80+
antigen presenting cells, activated CD4+ T cells and IL6 production. Our
in-silico study provides additional support to the hypothesis that blocking of
the CD80/CD86 signaling axis may be protective of the excessive proinflammatory
state associated with COVID-19 in the lungs.


  The COVID-19 has caused more than three million infections and over two
hundred thousand deaths by April 20201. Limiting socioeconomic activities (SA)
is among the most adopted governmental mitigating efforts to combat the
transmission of the virus, though the degree varies dramatically among
different regimes2. This study aims to quantify the contribution from the SA
and weather conditions to the transmission of COVID-19 at global scale. Ruling
out the unobservable factors including medical facilities and other control
policies (MOC) through region-by-time fixed effects3,4, we show that the
limiting SA has a leading contribution to lower the reproductive number by
18.3%, while weather conditions, including ultraviolet, relative humidity, and
wind explain a smaller amount of variation. Temperature might have a
non-monotonic impact on the transmission. We further show that in developed
countries5 and China, the SA effect is more pronounced whereas the weather
effect is significantly downplayed possibly because people tend to stay indoors
most of the time with a controlled climate. We finally estimate the reduced
reproductive number and the population spared from infections due to
restricting SA at 40,964, 180,336, 174,494, in China, United States, and Europe
respectively. From late January to mid-April, all regions, except for China,
Australia, and south Korea show a steep upward trend of spared infections due
to restricting SA. US and Europe, in particular, show far steeper upward trends
of spared infections in the analyzed timeframe, signaling a greater risk of
reopening the economy too soon.


  I tentatively re-analyze data from two well-publicized studies on COVID-19,
namely the Charit\'{e} "viral load in children" and the Bonn "seroprevalence in
Heinsberg/Gangelt" study, from information available in the preprints. The
studies have the following in common:
  - They received worldwide attention and arguably had policy impact.
  - The thrusts of their findings align with the respective lead authors'
(different) public stances on appropriate response to COVID-19.
  - Tentatively, my reading of the Gangelt study neutralizes its thrust, and my
reading of the Charit\'{e} study reverses it.
  The exercise may aid in placing these studies in the literature. With all
caveats that apply to n=2 quickfire analyses based off preprints, one also
wonders whether it illustrates inadvertent effects of "researcher degrees of
freedom."


  Using a probability of novel encounter derived from a physical model, we
augment the SIR compartmental model for disease spread. Scenarios with the same
initial trajectories and identical $R_0$ values can diverge greatly depending
on the speed at which our circles of acquaintances grow stale - leading to
order of magnitude differences in final case counts. A momentum effect arises
from variation in the mean time since infection, and this feeds back into new
infection rate and faster decline in the late stages of an outbreak. Rapid
extinction of an outbreak can occur in the early stages, but once this
opportunity is missed the effect is diminished and then, only herd immunity can
help.


  The dynamics of molecular collisions in a macroscopic body are encoded by the
parameter Thermodynamic entropy - a statistical measure of the number of
molecular configurations that correspond to a given macrostate. Directionality
in the flow of energy in macroscopic bodies is described by the Second Law of
Thermodynamics: In isolated systems, that is systems closed to the input of
energy and matter, thermodynamic entropy increases. The dynamics of the lower
level interactions in populations of replicating organisms is encoded by the
parameter Evolutionary entropy, a statistical measure which describes the
number and diversity of metabolic cycles in a population of replicating
organisms. Directionality in the transformation of energy in populations of
organisms is described by the Fundamental Theorem of Evolution: In systems open
to the input of energy and matter, Evolutionary entropy increases, when the
energy source is scarce and diverse, and decreases when the energy source is
abundant and singular. This article shows that when rho to 0, and N to
infinity, where rho is the production rate of the external energy source, and N
denote the number of replicating units, evolutionary entropy, an organized
state of energy; and thermodynamic entropy, a randomized state of energy,
coincide. Accordingly, the Fundamental Theorem of Evolution, is a
generalization of the Second Law of Thermodynamics.


  By the end of July 2020, the COVID-19 pandemic had infected more than
seventeen million people and had spread to almost all countries worldwide. In
response, many countries all over the world have used different methods to
reduce the infection rate, such as including case isolation, the closure of
schools and universities, banning public events, and mostly forcing social
distancing, including local and national lockdowns. We use a Monte-Carlo (MC)
based algorithm to predict the virus infection rate for different population
densities using the most recent epidemic data in our work. We test the spread
of the Coronavirus using three different lockdown models, and eight various
combinations of constraints, which allow us to examine the efficiency of each
model and constraint. In this paper, we have tested three different time-cyclic
patterns of no-restrictions/lockdown patterns. This model's main prediction is
that a cyclic schedule of no-restrictions/lockdown that contains at least ten
days of lockdown for each time cycle can help control the virus infection. In
particular, this model reduces the infection rate when accompanied by social
distancing and complete isolation of symptomatic patients.


  The critical role of energy consumption in biological systems including T
cell discrimination process has been investigated in various ways. The kinetic
proofreading(KPR) in T cell recognition involving different levels of energy
dissipation influences functional outcomes such as error rates and specificity.
In this work, we study quantitatively how the energy cost influences error
fractions, sensitivity, specificity, kinetic speed in terms of Mean First
Passage Time(MFPT) and adaption errors. These provide the background to
adequately understand T cell dynamics. It is found that energy plays a central
role in the system that aims to achieve minimum error fractions and maximum
sensitivity and specificity with the fastest speed under our kinetic scheme for
which numerical values of kinetic parameters are specially chosen, but such a
condition can be broken with varying data. Starting with the application of
steady state approximation(SSA) to the evaluation of the concentration of each
complex produced associated with KPR, which is used to quantify various
observables, we present both analytical and numerical results in detail.


  The global COVID-19 pandemic (SARS-CoV-2 virus) is the defining health crisis
of our century. Due to the absence of vaccines and drugs that can help to fight
it, the world solution to control the spread has been to consider public social
distance measures that avoids the saturation of the health system. In this
context, we investigate a Model Predictive Control (MPC) framework to determine
the time and duration of social distancing policies. We use Brazilian data in
the period from March to May of 2020. The available data regarding the number
of infected individuals and deaths suffers from sub-notification due to the
absence of mass tests and the relevant presence of the asymptomatic
individuals. We estimate variations of the SIR model using an
uncertainty-weighted Least-Squares criterion that considers both nominal and
inconsistent-data conditions. Moreover, we add to our versions of the SIR model
an additional dynamic state variable to mimic the response of the population to
the social distancing policies determined by the government that affects the
speed of COVID-19 transmission. Our control framework is within a mixed-logical
formalism, since the decision variable is forcefully binary (the existence or
the absence of social distance policy). A dwell-time constraint is included to
avoid harsh shifting between these two states. Finally, we present simulation
results to illustrate how such optimal control policy would operate. These
results point out that no social distancing should be relaxed before mid August
2020. If relaxations are necessary, they should not be performed before the
beginning this date and should be in small periods, no longer than 25 days.
This paradigm would proceed roughly until January/2021. The second peak of
infections, which has a forecast to the beginning of October, can be reduced if
the periods of no-isolation days are shortened.


  Recent studies of cancer somatic mutation profiles seek to identify mutations
for targeted therapy in personalized medicine. Analysis of profiles, however,
is not trivial, as each profile is heterogeneous and there are multiple
confounding factors that influence the cause-and-effect relationships between
cancer genes such as cancer (sub)type, biological processes, total number of
mutations, and non-linear mutation interactions. Moreover, cancer is
biologically redundant, i.e., distinct mutations can result in the alteration
of similar biological processes, so it is important to identify all possible
combinatorial sets of mutations for effective patient treatment. To model this
phenomena, we propose the correlated zero-inflated negative binomial process to
infer the inherent structure of somatic mutation profiles through latent
representations. This stochastic process takes into account different, yet
correlated, co-occurring mutations using profile-specific negative binomial
dispersion parameters that are mixed with a correlated beta-Bernoulli process
and a probability parameter to model profile heterogeneity. These model
parameters are inferred by iterative optimization via amortized and stochastic
variational inference using the Pan Cancer dataset from The Cancer Genomic
Archive (TCGA). By examining the the latent space, we identify biologically
relevant correlations between somatic mutations.


  DNA strand displacement, where a single-stranded nucleic acid invades a DNA
duplex, is pervasive in genomic processes and DNA engineering applications. The
kinetics of strand displacement have been studied in bulk; however, the
kinetics of the underlying strand exchange were obfuscated by a slow
bimolecular association step. Here, we use a novel single-molecule Fluorescence
Resonance Energy Transfer (smFRET) approach termed the "fission" assay to
obtain the full distribution of first passage times of unimolecular strand
displacement. At a frame time of 4.4 ms, the first passage time distribution
for a 14-nt displacement domain exhibited a nearly monotonic decay with little
delay. Among the eight different sequences we tested, the mean displacement
time was on average 35 ms and varied by up to a factor of 13. The measured
displacement kinetics also varied between complementary invaders and between
RNA and DNA invaders of the same base sequence except for T$\rightarrow$U
substitution. However, displacement times were largely insensitive to the
monovalent salt concentration in the range of 0.25 M to 1 M. Using a
one-dimensional random walk model, we infer that the single-step displacement
time is in the range of $\sim 30 \mu s$ to $\sim 300 \mu s$ depending on the
base identity. The framework presented here is broadly applicable to the
kinetic analysis of multistep processes investigated at the single-molecule
level.


  Timely, creditable, and fine-granular case information is vital for local
communities and individual citizens to make rational and data-driven responses
to the COVID-19 pandemic. This paper presents CovidNet, a COVID-19 tracking
project associated with a large scale epidemic dataset, which was initiated by
1Point3Acres. To the best of our knowledge, the project is the only platform
providing real-time global case information of more than 4,124 sub-divisions
from over 27 countries worldwide with multi-language supports. The platform
also offers interactive visualization tools to analyze the full historical case
curves in each region. Initially launched as a voluntary project to bridge the
data transparency gap in North America in January 2020, this project by far has
become one of the major independent sources worldwide and has been consumed by
many other tracking platforms. The accuracy and freshness of the dataset is a
result of the painstaking efforts from our voluntary teamwork, crowd-sourcing
channels, and automated data pipelines. As of May 18, 2020, the project website
has been visited more than 200 million times and the CovidNet dataset has
empowered over 522 institutions and organizations worldwide in policy-making
and academic researches. All datasets are openly accessible for non-commercial
purposes at https://coronavirus.1point3acres.com via a formal request through
our APIs.


  We consider how the NHS COVID-19 application will initially calculate a risk
score for an individual based on their recent contact with people who report
that they have coronavirus symptoms.


  When the novel coronavirus disease SARS-CoV2 (COVID-19) was officially
declared a pandemic by the WHO in March 2020, the scientific community had
already braced up in the effort of making sense of the fast-growing wealth of
data gathered by national authorities all over the world. However, despite the
diversity of novel theoretical approaches and the comprehensiveness of many
widely established models, the official figures that recount the course of the
outbreak still sketch a largely elusive and intimidating picture. Here we show
unambiguously that the dynamics of the COVID-19 outbreak belongs to the simple
universality class of the SIR model and extensions thereof. Our analysis
naturally leads us to establish that there exists a fundamental limitation to
any theoretical approach, namely the unpredictable non-stationarity of the
testing frames behind the reported figures. However, we show how such bias can
be quantified self-consistently and employed to mine useful and accurate
information from the data. In particular, we describe how the time evolution of
the reporting rates controls the occurrence of the apparent epidemic peak,
which typically follows the true one in countries that were not vigorous enough
in their testing at the onset of the outbreak. The importance of testing early
and resolutely appears as a natural corollary of our analysis, as countries
that tested massively at the start clearly had their true peak earlier and less
deaths overall.


  Severe acute respiratory syndrome coronavirus two (SARS-CoV-2), the virus
responsible for the coronavirus disease 2019 (COVID-19) pandemic, represents an
unprecedented global health challenge. Consequently, a large amount of research
into the disease pathogenesis and potential treatments has been carried out in
a short time frame. However, developing novel drugs is a costly and lengthy
process, and is unlikely to deliver a timely treatment for the pandemic. Drug
repurposing, by contrast, provides an attractive alternative, as existing drugs
have already undergone many of the regulatory requirements. In this work we
used a combination of network algorithms and human curation to search
integrated knowledge graphs, identifying drug repurposing opportunities for
COVID-19. We demonstrate the value of this approach, reporting on eight
potential repurposing opportunities identified, and discuss how this approach
could be incorporated into future studies.


  The presence of a large number of infected individuals with few or no
symptoms is an important epidemiological difficulty and the main mathematical
feature of COVID-19. The A-SIR model, i.e. a SIR (Susceptible-Infected-Removed)
model with a compartment for infected individuals with no symptoms or few
symptoms was proposed by Giuseppe Gaeta, arXiv:2003.08720 [q-bio.PE] (2020). In
this paper we investigate a slightly generalized version of the same model and
propose a scheme for fitting the parameters of the model to real data using the
time series only of the deceased individuals. The scheme is applied to the
concrete cases of Lombardy, Italy and S\~ao Paulo state, Brazil, showing
different aspects of the epidemics. For each case we show that we may have good
fits to the data up to the present, but with very large differences in the
future behavior. The reasons behind such disparate outcomes are the uncertainty
on the value of a key parameter, the probability that an infected individual is
fully symptomatic, and on the intensity of the social distancing measures
adopted. This conclusion enforces the necessity of trying to determine the real
number of infected individuals in a population, symptomatic or asymptomatic.


  Latin America is experiencing severe impacts of the SARS-CoV-2 pandemic, but
poverty and weak public health institutions hamper gathering the kind of
refined data needed to inform classical SEIR models of epidemics. We present an
alternative approach that draws on advances in statistical ecology and
conservation biology to enhance the value of sparse data in projecting and
ameliorating epidemics. Our approach, leading to what we call a Stochastic
Epidemic Gompertz model, with few parameters can flexibly incorporate
heterogeneity in transmission within populations and across time. We
demonstrate that poverty has a large impact on the course of the pandemic,
across fourteen Latin American countries, and show how our approach provides
flexible, time-varying projections of disease risk that can be used to refine
public health strategies.


  Covid-19, a serious respiratory complications caused by SARS-CoV-2 has become
one of the global threat to human healthcare system. The present study
evaluated the possibility of plant originated approved 117 therapeutics against
the main protease protein (MPP), RNA-dependent RNA polymerase (RdRp) and spike
protein (S) of SARS-CoV-2 including drug surface analysis by using molecular
docking through drug repurposing approaches. The molecular interaction study
revealed that Rifampin (-16.3 kcal/mol) were topmost inhibitor of MPP where
Azobechalcone were found most potent plant therapeutics for blocking the RdRp
(-15.9 kcal /mol) and S (-14.4 kcal/mol) protein of SARS-CoV-2. After the
comparative analysis of all docking results, Azobechalcone, Rifampin,
Isolophirachalcone, Tetrandrine and Fangchinoline were exhibited as the most
potential inhibitory plant compounds for targeting the key proteins of
SARS-CoV-2. However, amino acid positions; H41, C145, and M165 of MPP played
crucial roles for the drug surface interaction where F368, L371, L372, A375,
W509, L514, Y515 were pivotal for RdRP. In addition, the drug interaction
surface of S proteins also showed similar patterns with all of its maximum
inhibitors. ADME analysis also strengthened the possibility of screened plant
therapeutics as the potent drug candidates against SARS-C with the highest drug
friendliness.


  Living systems produce copies of information-carrying molecules such as DNA
by assembling monomer units into finite-length oligomer (short polymer) copies.
We explore the role of initiation and termination of the copy process in the
thermodynamics of copying. By splitting the free-energy change of copy
formation into informational and chemical terms, we show that copy accuracy
plays no direct role in the overall thermodynamics. Instead, it is
thermodynamically costly to produce outputs that are more similar to the
oligomers in the environment than sequences obtained by randomly sampling
monomers. Copy accuracy can be thermodynamically neutral, or even favoured,
depending on the surroundings. Oligomer copying mechanisms can thus function as
information engines that interconvert chemical and information-based free
energy. Hard thermodynamic constraints on accuracy derived for infinite-length
polymers instead manifest as kinetic barriers experienced while the copy is
template-attached. These barriers are easily surmounted by shorter oligomers.


  Mass public quarantining, colloquially known as a lock-down, is a
non-pharmaceutical intervention to check spread of disease. This paper presents
ESOP (Epidemiologically and Socio-economically Optimal Policies), a novel
application of active machine learning techniques using Bayesian optimization,
that interacts with an epidemiological model to arrive at lock-down schedules
that optimally balance public health benefits and socio-economic downsides of
reduced economic activity during lock-down periods. The utility of ESOP is
demonstrated using case studies with VIPER
(Virus-Individual-Policy-EnviRonment), a stochastic agent-based simulator that
this paper also proposes. However, ESOP is flexible enough to interact with
arbitrary epidemiological simulators in a black-box manner, and produce
schedules that involve multiple phases of lock-downs.


  There are many hard-to-reconcile numbers circulating concerning Covid-19.
Using reports from random testing, the fatality ratio per infection is
evaluated and used to extract further information on the actual fraction of
infections and the success of their identification for different countries.


  We assess the impact of COVID-19 response measures implemented in Germany and
Switzerland on cumulative COVID-19-related hospitalization and death rates. Our
analysis exploits the fact that the epidemic was more advanced in some regions
than in others when certain lockdown measures came into force, based on
measuring health outcomes relative to the region-specific start of the epidemic
and comparing outcomes across regions with earlier and later start dates. When
estimating the effect of the relative timing of measures, we control for
regional characteristics and initial epidemic trends by linear regression
(Germany and Switzerland), doubly robust estimation (Germany), or synthetic
controls (Switzerland). We find for both countries that a relatively later
exposure to the measures entails higher cumulative hospitalization and death
rates on region-specific days after the outbreak of the epidemic, suggesting
that an earlier imposition of measures is more effective than a later one. For
Germany, we also evaluate curfews (as introduced in a subset of states) based
on cross-regional variation. We do not find any effects of curfews on top of
the federally imposed contact restriction that banned groups of more than 2
individuals. Finally, an analysis of mobility patterns in Switzerland shows an
immediate behavioral effect of the lockdown in terms of reduced mobility.


  This paper considers an explicit numerical scheme for solving the
mathematical model of the propagation of Covid-19 epidemic with undetected
infectious cases. We analyze the stability and convergence rate of the new
approach in $L^{\infty}$-norm. The proposed method is less time consuming.
Furthermore, the method is stable, at least second-order convergent and can
serve as a robust tool for the integration of general systems of ordinary
differential equations. A wide set of numerical evidences which consider the
case of Cameroon are presented and discussed.


  Most models of epidemic spread, including many designed specifically for
COVID-19, implicitly assume mass-action contact patterns and undirected contact
networks, meaning that the individuals most likely to spread the disease are
also the most at risk to receive it from others. Here, we review results from
the theory of random directed graphs which show that many important quantities,
including the reproduction number and the epidemic size, depend sensitively on
the joint distribution of in- and out-degrees ("risk" and "spread"), including
their heterogeneity and the correlation between them. By considering joint
distributions of various kinds, we elucidate why some types of heterogeneity
cause a deviation from the standard Kermack-McKendrick analysis of SIR models,
i.e., so-called mass-action models where contacts are homogeneous and random,
and some do not. We also show that some structured SIR models informed by
realistic complex contact patterns among types of individuals (age or activity)
are simply mixtures of Poisson processes and tend not to deviate significantly
from the simplest mass-action model. Finally, we point out some possible policy
implications of this directed structure, both for contact tracing strategy and
for interventions designed to prevent superspreading events. In particular,
directed graphs have a forward and backward version of the classic "friendship
paradox" -- forward edges tend to lead to individuals with high risk, while
backward edges lead to individuals with high spread -- such that a combination
of both forward and backward contact tracing is necessary to find
superspreading events and prevent future cascades of infection.


  Many countries have passed their first COVID-19 epidemic peak. Traditional
epidemiological models describe this as a result of non-pharmaceutical
interventions that pushed the growth rate below the recovery rate. In this new
phase of the pandemic many countries show an almost linear growth of confirmed
cases for extended time-periods. This new containment regime is hard to explain
by traditional models where infection numbers either grow explosively until
herd immunity is reached, or the epidemic is completely suppressed (zero new
cases). Here we offer an explanation of this puzzling observation based on the
structure of contact networks. We show that for any given transmission rate
there exists a critical number of social contacts, $D_c$, below which linear
growth and low infection prevalence must occur. Above $D_c$ traditional
epidemiological dynamics takes place, as e.g. in SIR-type models. When
calibrating our corresponding model to empirical estimates of the transmission
rate and the number of days being contagious, we find $D_c\sim 7.2$. Assuming
realistic contact networks with a degree of about 5, and assuming that lockdown
measures would reduce that to household-size (about 2.5), we reproduce actual
infection curves with a remarkable precision, without fitting or fine-tuning of
parameters. In particular we compare the US and Austria, as examples for one
country that initially did not impose measures and one that responded with a
severe lockdown early on. Our findings question the applicability of standard
compartmental models to describe the COVID-19 containment phase. The
probability to observe linear growth in these is practically zero.


  HIV-1 affects tens of millions of people worldwide. Current treatments often
involve a cocktail of antiretroviral drugs, which are effective in reducing the
virus and extending life spans. However, there is currently no FDA-approved
HIV-1 transcription inhibitor. Furthermore, there have only been a few attempts
to model the transcription process in HIV-1. In this work, we extend a novel
three-state model of HIV-1 transcription introduced in DeMarino et al. (2020)
that has been developed and validated against experimental data. After fitting
this model to in vitro data, significant differences in the transcription
process of HIV-1 in T-cells and macrophages have been observed. In particular,
the activation of the HIV-1 promoter in T-cells appears to take place rapidly
as the Tat protein approaches a critical threshold. In contrast, the same
process occurs smoother in macrophages.
  In this work, we carry out systematic mathematical analyses of the model to
complement experimental data fitting and sensitivity analysis performed
earlier. We derive explicit solutions of the model to obtain exact
transcription process decay rates for the original model and then study the
effect of nonlinearity on the system behavior, including the existence and the
local and global stability of the positive equilibrium. We were able to show
the stability of the positive steady state in limiting cases, with the global
stability in the general case remaining an open question.
  By modeling the effect of transcription-inhibiting drug therapy, we provide a
nontrivial condition for it to be effective in reducing viral load. Moreover,
our numerical simulations and analysis point out that the effect of the
transcription-inhibitor can be enhanced by synchronizing with standard
treatments, such as combination antiretroviral therapy, to allow the reduction
of total dosages and toxicity.


  We shall apply a generalized SEIR model to study the outbreak of COVID-19 in
Brazil. In particular, we would like to explain the projections of the increase
in the level of infection over a long period of time, overlapping large local
outbreaks in the most populous states in the country. A time-dependent dynamic
SEIR model inspired in a model previously used during the outbreak in China was
used to analyses the time trajectories of infected, recovered, and deaths. The
model has parameters that vary with time and are fitted considering a nonlinear
least-squares method. The simulations starting from April 8, 2020, concluded
that the time for a peak in Brazil will be in July 21, 2020 with total
cumulative infected cases around 982K people; in addition, an estimated total
death case will reach to 192K in the end. Besides that, Brazil will reach a
peak in terms of daily new infected cases and death cases around the middle of
July with 50K cases of infected and almost 6.0K daily deaths.


  DNA strand displacement (DSD) reactions have been used to construct chemical
reaction networks in which species act catalytically at the level of the
overall stoichiometry of reactions. These effective catalytic reactions are
typically realised through one or more of the following: many-stranded gate
complexes to coordinate the catalysis, indirect interaction between the
catalyst and its substrate, and the recovery of a distinct ``catalyst'' strand
from the one that triggered the reaction. These facts make emulation of the
out-of-equilibrium catalytic circuitry of living cells more difficult. Here, we
propose a new framework for constructing catalytic DSD networks: Active
Circuits of Duplex Catalysts (ACDC). ACDC components are all double-stranded
complexes, with reactions occurring through 4-way strand exchange. Catalysts
directly bind to their substrates, and and the ``identity'' strand of the
catalyst recovered at the end of a reaction is the same molecule as the one
that initiated it. We analyse the capability of the framework to implement
catalytic circuits analogous to phosphorylation networks in living cells. We
also propose two methods of systematically introducing mismatches within DNA
strands to avoid leak reactions and introduce driving through net base pair
formation. We then combine these results into a compiler to automate the
process of designing DNA strands that realise any catalytic network allowed by
our framework.


  The Susceptible-Exposed-Infected-Recovered (SEIR) epidemiological model is
one of the standard models of disease spreading. Here we analyse an extended
SEIR model that accounts for asymptomatic carriers, believed to play an
important role in COVID-19 transmission. For this model we derive a number of
analytic results for important quantities such as the peak number of
infections, the time taken to reach the peak and the size of the final affected
population. We also propose an accurate way of specifying initial conditions
for the numerics (from insufficient data) using the fact that the early time
exponential growth is well-described by the dominant eigenvector of the
linearized equations. Secondly we explore the effect of different intervention
strategies such as social distancing (SD) and testing-quarantining (TQ). The
two intervention strategies (SD and TQ) try to reduce the disease reproductive
number, $R_0$, to a target value $R^{\rm target}_0 < 1$, but in distinct ways,
which we implement in our model equations. We find that for the same $R^{\rm
target}_0 < 1$, TQ is more efficient in controlling the pandemic than SD.
However, for TQ to be effective, it has to be based on contact tracing and our
study quantifies the required ratio of tests-per-day to the number of new
cases-per-day. Our analysis shows that the largest eigenvalue of the linearised
dynamics provides a simple understanding of the disease progression, both pre-
and post- intervention, and explains observed data for many countries. We apply
our results to the COVID data for India to obtain heuristic projections for the
course of the pandemic, and note that the predictions strongly depend on the
assumed fraction of asymptomatic carriers.


  The success of DNA nanotechnology has been driven by the discovery of novel
structural motifs with a wide range of shapes and uses. We present a
comprehensive study of the T-motif, a 3-armed, planar, right-angled junction
that has been used in the self-assembly of DNA polyhedra and periodic
structures. The motif is formed through the interaction of a bulge loop in one
duplex and a sticky end of another. The polarity of the sticky end has
significant consequences for the thermodynamic and geometrical properties of
the T-motif: different polarities create junctions spanning different grooves
of the duplex. We compare experimental binding strengths with predictions of
oxDNA, a coarse-grained model of DNA, for various loop sizes. We find that,
although both sticky-end polarities can create stable junctions, junctions
resulting from 5$'$ sticky ends are stable over a wider range of bulge loop
sizes. We highlight the importance of possible coaxial stacking interactions
within the motif and investigate how each coaxial stacking interaction
stabilises the structure and favours a particular geometry.


  Key enzymatic processes in biology use the nonequilibrium error correction
mechanism called kinetic proofreading to enhance their specificity. Kinetic
proofreading typically requires several dedicated structural features in the
enzyme, such as a nucleotide hydrolysis site and multiple enzyme-substrate
conformations that delay product formation. Such requirements limit the
applicability and the adaptability of traditional proofreading schemes. Here,
we explore an alternative conceptual mechanism of error correction that
achieves delays between substrate binding and subsequent product formation by
having these events occur at distinct physical locations. The time taken by the
enzyme-substrate complex to diffuse from one location to another is leveraged
to discard wrong substrates. This mechanism does not require dedicated
structural elements on the enzyme, making it easier to overlook in experiments
but also making proofreading tunable on the fly. We discuss how tuning the
length scales of enzyme or substrate concentration gradients changes the
fidelity, speed and energy dissipation, and quantify the performance
limitations imposed by realistic diffusion and reaction rates in the cell. Our
work broadens the applicability of kinetic proofreading, and sets the stage for
the study of spatial gradients as a possible route to specificity.


  We consider the problem of distance estimation under the TKF91 model of
sequence evolution by insertions, deletions and substitutions on a phylogeny.
In an asymptotic regime where the expected sequence lengths tend to infinity,
we show that no consistent distance estimation is possible from sequence
lengths alone. More formally, we establish that the distributions of pairs of
sequence lengths at different distances cannot be distinguished with
probability going to one.


  Objective: Total 186 biologically important phenylpropanoids and polyketides
compounds from different Indian medicinal plants and dietary sources were
screened to filter potential compounds that bind at the active site of the
therapeutic target proteins of SARS-CoV-2. Method: The molecular docking
studies were carried out by using the Autodock Vina. The in silico ADMET and
drug-likeness properties of the compounds were predicted from SwissADME server.
Result: The molecular docking study of the 186 compounds with the therapeutic
target proteins (Mpro, PLpro, RdRp, SGp and ACE2) of SARS-CoV-2 resulted 40
compounds that bind at the active site with dock score above -8.0 kcal/mol.
Conclusion: Based on the in silico ADMET study and drug-likeness prediction of
40 compounds, we proposed petunidin, baicalein, cyanidin,
7-hydroxy-3',4'-methylenedioxyflavan, quercetin and ellagic acid among the 186
biologically important phenylpropanoids and polyketides as potential lead
compounds, which can further be investigated pharmacologically and clinically
to formulate therapeutic approaches for the COVID-19.


  We study Susceptible-Exposed-Asymptomatic-Infectious-Recovered (SEAIR)
epidemic spreading model of COVID-19. It captures two important characteristics
of the infectiousness of COVID-19: delayed start and its appearance before
onset of symptoms, or even with total absence of them. The model is
theoretically analyzed in continuous-time compartmental version and
discrete-time version on random regular graphs and complex networks. We show
analytically that there are relationships between the epidemic thresholds and
the equations for the susceptible populations at the endemic equilibrium in all
three versions, which hold when the epidemic is weak. We provide theoretical
arguments that eigenvector centrality of a node approximately determines its
risk to become infected.


  We propose a mathematical model to investigate the effects of
information-dependent vaccination behavior on meningitis transmission. The
information is represented by means of information index as early proposed in
(d'Onofrio et al., Theor. pop. biol., 2007). We perform a qualitative analysis
based on stability theory, focusing to the global stability of the disease free
equilibrium (DFE) and the related transcritical bifurcation taking place at the
threshold for the DFE. Finally, we assess the role of epidemiological and
information parameters in the model dynamics through numerical simulations. Our
simulations suggests that the impact of the human behavior critically depend on
the average information delay. For example, it can induce recurrent epidemics,
provided that transfer rate from the carrier to the infectious state is over a
threshold. Otherwise, the endemic equilibrium is (at least) locally stable.


  We introduce $\textit{scadnano}$ (https://scadnano.org) (short for
"scriptable cadnano"), a computational tool for designing synthetic DNA
structures. Its design is based heavily on cadnano, the most widely-used
software for designing DNA origami, with three main differences:
  1. scadnano runs entirely in the browser, with $\textit{no software
installation}$ required.
  2. scadnano designs, while they can be edited manually, can also be created
and edited by a $\textit{well-documented Python scripting library}$, to help
automate tedious tasks.
  3. The scadnano file format is $\textit{easily human-readable}$. This goal is
closely aligned with the scripting library, intended to be helpful when
debugging scripts or interfacing with other software. The format is also
somewhat more expressive than that of cadnano, able to describe a broader range
of DNA structures than just DNA origami.


  The actin cortex is a thin layer of actin, myosin, and actin binding proteins
that underlies the membrane of most animal cells. It is highly dynamic and can
undergo remodelling on time-scales of tens of seconds thanks to protein
turnover and myosin-mediated contractions. The cortex enables cells to resist
external mechanical stresses, controls cell shape, and allows cells to exert
forces on their neighbours. Thus, its mechanical properties are key to its
physiological function. Here, we give an overview of how cortex composition,
structure, and dynamics control cortex mechanics and cell shape. We use mitosis
as an example to illustrate how global and local regulation of cortex mechanics
give rise to a complex series of cell shape changes.


  The sterile insect technique consists in massive release of sterilized males
in the aim to reduce the size of mosquitoes population or even eradicate it. In
this work, we investigate the feasability of using the sterile insect technique
as a barrier against reinvasion. More precisely, we provide some numerical
simulations and mathematical results showing that performing the sterile insect
technique on a band large enough may stop reinvasion.


  We propose a mesoscale model structure for the coronavirus nucleocapsid,
assembled from the high resolution structures of the basic building blocks of
the N-protein, CryoEM imaging and mathematical constraints for an overall
quasi-spherical particle. The structure is a truncated octahedron that
accommodates two layers: an outer shell composed of triangular and quadrangular
lattices of the N-terminal domain and an inner shell of equivalent lattices of
coiled parallel helices of the C-terminal domain. The model is consistent with
the dimensions expected for packaging large viral genomes and provides a
rationale to interpret the apparent pleomorphic nature of coronaviruses.


  In this paper we develop a predictive model for the spread of COVID-19
infection at a provincial (i.e. EU NUTS-3) level in Italy by using official
data from the Italian Ministry of Health integrated with data extracted from
daily official press conferences of regional authorities and from local
newspaper websites. This integration is mainly concerned with COVID-19 cause
specific death data which are not available at NUTS-3 level from open official
data data channels. An adjusted time-dependent SIRD model is used to predict
the behavior of the epidemic, specifically the number of susceptible, infected,
deceased and recovered people. Predictive model performance is evaluated using
comparison with real data.


  Timely estimation of the current value for COVID-19 reproduction factor $R$
has become a key aim of efforts to inform management strategies. $R$ is an
important metric used by policy-makers in setting mitigation levels and is also
important for accurate modelling of epidemic progression. This brief paper
introduces a method for estimating $R$ from biased case testing data. Using
testing data, rather than hospitalisation or death data, provides a much
earlier metric along the symptomatic progression scale. This can be hugely
important when fighting the exponential nature of an epidemic. We develop a
practical estimator and apply it to Scottish case testing data to infer a
current (20 May 2020) $R$ value of $0.74$ with $95\%$ confidence interval
$[0.48 - 0.86]$.


  The impact of the outbreak of COVID-19 on health has been widely concerned.
Disease risk assessment, prediction, and early warning have become a
significant research field. Previous research suggests that there is a
relationship between air quality and the disease. This paper investigated the
impact of the atmospheric environment on the basic reproduction number (R$_0$)
in Australia, South Korea, and Italy by using atmospheric environment data,
confirmed case data, and the distributed lag non-linear model (DLNM) model
based on Quasi-Poisson regression. The results show that the air temperature
and humidity have lag and persistence on short-term R$_0$, and seasonal factors
have an apparent decorating effect on R$_0$. PM$_{10}$ is the primary pollutant
that affects the excess morbidity rate. Moreover, O$_3$, PM$_{2.5}$, and SO$_2$
as perturbation factors have an apparent cumulative effect. These results
present beneficial knowledge for correlation between environment and COVID-19,
which guiding prospective analyses of disease data.


  A model of optimal control of ventilation has recently been developed for
humans. This model highlights the importance of the localization of the
transition between a convective and a diffusive transport of respiratory gas.
This localization determines how ventilation should be controlled in order to
minimize its energetic cost at any metabolic regime. We generalized this model
to any mammal, based on the core morphometric characteristics shared by all
mammalian lungs and on their allometric scaling from the literature. Since the
main energetic costs of ventilation are related to convective transport, we
prove that, for all mammals, the localization of the shift from a convective
transport to a diffusive transport plays a critical role on keeping this cost
low while fulfilling the lung function. Our model predicts for the first time
the localization of this transition in order to minimize the energetic cost of
ventilation, depending on mammal mass and metabolic regime. From this optimal
localization, we are able to predict allometric scaling laws for both tidal
volumes and breathing rates, at any metabolic rate. We ran our model for the
three common metabolic rates -- basal, field and maximal -- and showed that our
predictions reproduce accurately experimental data available in the literature.
Our analysis supports the hypothesis that mammals allometric scaling laws of
tidal volumes and breathing rates at a given metabolic rate are driven by a few
core geometrical characteristics shared by mammalian lungs and by the physical
processes of respiratory gas transport.


  A recent work (DOI 10.1101/2020.05.06.20093310) indicated that temporarily
splitting larger populations into smaller groups can efficiently mitigate the
spread of SARS-CoV-2 virus. The fact that, soon afterwards, on May 15, 2020,
the two million people Slovenia was the first European country proclaiming the
end of COVID-19 epidemic within national borders may be relevant from this
perspective. Motivated by this evolution, in this paper we investigate the time
dynamics of coronavirus cases in Slovenia with emphasis on how efficient
various containment measures act to diminish the number of COVID-19 infections.
Noteworthily, the present analysis does not rely on any speculative theoretical
assumption; it is solely based on raw epidemiological data. Out of the results
presented here, the most important one is perhaps the finding that, while
imposing drastic curfews and travel restrictions reduce the infection rate
kappa by a factor of four with respect to the unrestricted state, they only
improve the \k{appa}-value by ~15 % as compared to the much bearable state of
social and economical life wherein (justifiable) wearing face masks and social
distancing rules are enforced/followed. Significantly for behavioral and social
science, our analysis of the time dependence \k{appa} = \k{appa}(t) may reveal
an interesting self-protection instinct of the population, which became
manifest even before the official lockdown enforcement.


  COVID-19 pandemic has reshaped our world in a timescale much shorter than
what we can understand. Particularities of SARS-CoV-2, such as its persistence
in surfaces and the lack of a curative treatment or vaccine against COVID-19,
have pushed authorities to apply restrictive policies to control its spreading.
As data drove most of the decisions made in this global contingency, their
quality is a critical variable for decision-making actors, and therefore should
be carefully curated. In this work, we analyze the sources of error in
typically reported epidemiological variables and usual tests used for
diagnosis, and their impact on our understanding of COVID-19 spreading
dynamics. We address the existence of different delays in the report of new
cases, induced by the incubation time of the virus and testing-diagnosis time
gaps, and other error sources related to the sensitivity/specificity of the
tests used to diagnose COVID-19. Using a statistically-based algorithm, we
perform a temporal reclassification of cases to avoid delay-induced errors,
building up new epidemiologic curves centered in the day where the contagion
effectively occurred. We also statistically enhance the robustness behind the
discharge/recovery clinical criteria in the absence of a direct test, which is
typically the case of non-first world countries, where the limited testing
capabilities are fully dedicated to the evaluation of new cases. Finally, we
applied our methodology to assess the evolution of the pandemic in Chile
through the Effective Reproduction Number $R_t$, identifying different moments
in which data was misleading governmental actions. In doing so, we aim to raise
public awareness of the need for proper data reporting and processing protocols
for epidemiological modelling and predictions.


  We demonstrate the ability of statistical data assimilation to identify the
measurements required for accurate state and parameter estimation in an
epidemiological model for the novel coronavirus disease COVID-19. Our context
is an effort to inform policy regarding social behavior, to mitigate strain on
hospital capacity. The model unknowns are taken to be: the time-varying
transmission rate, the fraction of exposed cases that require hospitalization,
and the time-varying detection probabilities of new asymptomatic and
symptomatic cases. In simulations, we obtain accurate estimates of undetected
(that is, unmeasured) infectious populations, by measuring the detected cases
together with the recovered and dead - and without assumed knowledge of the
detection rates. Given a noiseless measurement of the recovered population,
excellent estimates of all quantities are obtained using a temporal baseline of
101 days, with the exception of the time-varying transmission rate at times
prior to the implementation of social distancing. With low noise added to the
recovered population, accurate state estimates require a lengthening of the
temporal baseline of measurements. Estimates of all parameters are sensitive to
the contamination, highlighting the need for accurate and uniform methods of
reporting. The aim of this paper is to exemplify the power of SDA to determine
what properties of measurements will yield estimates of unknown parameters to a
desired precision, in a model with the complexity required to capture important
features of the COVID-19 pandemic.


  Background: Without proven effect treatments and vaccines, Social Distancing
is the key protection factor against COVID-19. Social distancing alone should
have been enough to protect again the virus, yet things have gone very
differently, with a big mismatch between theory and practice. What are the
reasons? A big problem is that there is no actual social distancing data, and
the corresponding people behavior in a pandemic is unknown. We collect the
world-first dataset on social distancing during the COVID-19 outbreak, so to
see for the first time how people really implement social distancing, identify
dangers of the current situation, and find solutions against this and future
pandemics.
  Methods: Using a sensor-based social distancing belt we collected social
distance data from people in Italy for over two months during the most critical
COVID-19 outbreak. Additionally, we investigated if and how wearing various
Personal Protection Equipment, like masks, influences social distancing.
  Results: Without masks, people adopt a counter-intuitively dangerous
strategy, a paradox that could explain the relative lack of effectiveness of
social distancing. Using masks radically changes the situation, breaking the
paradoxical behavior and leading to a safe social distance behavior. In
shortage of masks, DIY (Do It Yourself) masks can also be used: even without
filtering protection, they provide social distancing protection. Goggles should
be recommended for general use, as they give an extra powerful safety boost.
Generic Public Health policies and media campaigns do not work well on social
distancing: explicit focus on the behavioral problems of necessary mobility are
needed.


  The novel Corona Virus COVID-19 arrived on Australian shores around 25
January 2020. This paper presents a novel method of dynamically modeling and
forecasting the COVID-19 pandemic in Australia with a high degree of accuracy
and in a timely manner using limited data; a valuable resource that can be used
to guide government decision-making on societal restrictions on a daily and/or
weekly basis. The "partially-observable stochastic process" used in this study
predicts not only the future actual values with extremely low error, but also
the percentage of unobserved COVID-19 cases in the population. The model can
further assist policy makers to assess the effectiveness of several possible
alternative scenarios in their decision-making processes.


  We analyze an optimal control version of a simple SIRS epidemiology model.
The policy maker can adopt policies to diminish the contact rate between
infected and susceptible individuals, at a specific economic cost. The arrival
of a vaccine will immediately end the epidemic. Total or partial immunity is
modeled, while the contact rate can exhibit a (user specified) seasonality. The
problem is solved in a spreadsheet environment. A reasonable parameter
selection leads to optimal policies which are similar to those followed by
different countries. A mild response relying on eventually reaching a high
immunity level is optimal if ample health facilities are available. On the
other hand limited health care facilities lead to strict lock downs while
moderate ones allow a flattening of the curve approach.


  The Covid-19 outbreak of 2020 has required many governments to develop
mathematical-statistical models of the outbreak for policy and planning
purposes. This work provides a tutorial on building a compartmental model using
Susceptibles, Exposed, Infected, Recovered and Deaths status through time. A
Bayesian Framework is utilized to perform both parameter estimation and
predictions. This model uses interventions to quantify the impact of various
government attempts to slow the spread of the virus. Predictions are also made
to determine when the peak Active Infections will occur.


  Quantum computers can in principle solve certain problems exponentially more
quickly than their classical counterparts. We have not yet reached the advent
of useful quantum computation, but when we do, it will affect nearly all
scientific disciplines. In this review, we examine how current quantum
algorithms could revolutionize computational biology and bioinformatics. There
are potential benefits across the entire field, from the ability to process
vast amounts of information and run machine learning algorithms far more
efficiently, to algorithms for quantum simulation that are poised to improve
computational calculations in drug discovery, to quantum algorithms for
optimization that may advance fields from protein structure prediction to
network analysis. However, these exciting prospects are susceptible to "hype",
and it is also important to recognize the caveats and challenges in this new
technology. Our aim is to introduce the promise and limitations of emerging
quantum computing technologies in the areas of computational molecular biology
and bioinformatics.


  For practical reasons, many forecasts of case, hospitalization and death
counts in the context of the current COVID-19 pandemic are issued in the form
of central predictive intervals at various levels. This is also the case for
the forecasts collected in the COVID-19 Forecast Hub
(https://covid19forecasthub.org/). Forecast evaluation metrics like the
logarithmic score, which has been applied in several infectious disease
forecasting challenges, are then not available as they require full predictive
distributions. This article provides an overview of how established methods for
the evaluation of quantile and interval forecasts can be applied to epidemic
forecasts in this format. Specifically, we discuss the computation and
interpretation of the weighted interval score, which is a proper score that
approximates the continuous ranked probability score. It can be interpreted as
a generalization of the absolute error to probabilistic forecasts and allows
for a decomposition into a measure of sharpness and penalties for over- and
underprediction.


  The COVID-19 pandemic has caused major disturbance to human life. An
important reason behind the widespread social anxiety is the huge uncertainty
about the pandemic. A fundamental uncertainty is how many or what percentage of
people have been infected. There are published and frequently updated data on
various statistics of the pandemic, at local, country or global level. However,
due to various reasons, many cases were not included in those reported numbers.
We propose a structured approach for the estimation of the number of unreported
cases, where we distinguish cases that arrive late in the reported numbers and
those who had mild or no symptoms and thus were not captured by any medical
system at all. We use post-report data for the estimation of the former and
population matching to the latter. We estimate that the reported number of
infected cases in the US should be corrected by multiplying a factor of 220.54%
as of Apr 20, 2020, while the infection ratio out of the US population is
estimated to be 0.53%, implying a case mortality rate at 2.85% which is close
to the 3.4% suggested by the WHO in Mar 2020. Towards the end of the summer of
2020, the overall infection ratio of the US rises to 2.49% while the case
mortality decreases to 2.09%, and the ratio of asymptomatic cases out of all
infected cases reduces from the pre-summer 35-40% to around 20-25%.


  Kinetic proofreading is an error correction mechanism present in the
processes of the central dogma and beyond, and typically requires the free
energy of nucleotide hydrolysis for its operation. Though the molecular players
of many biological proofreading schemes are known, our understanding of how
energy consumption is managed to promote fidelity remains incomplete. In our
work, we introduce an alternative conceptual scheme called 'the piston model of
proofreading' where enzyme activation through hydrolysis is replaced with
allosteric activation achieved through mechanical work performed by a piston on
regulatory ligands. Inspired by Feynman's ratchet and pawl mechanism, we
consider a mechanical engine designed to drive the piston actions powered by a
lowering weight, whose function is analogous to that of ATP synthase in cells.
Thanks to its mechanical design, the piston model allows us to tune the 'knobs'
of the driving engine and probe the graded changes and trade-offs between
speed, fidelity and energy dissipation. It provides an intuitive explanation of
the conditions necessary for optimal proofreading and reveals the unexpected
capability of allosteric molecules to beat the Hopfield limit of fidelity by
leveraging the diversity of states available to them. The framework that we
built for the piston model can also serve as a basis for additional studies of
driven biochemical systems.


  Neuropeptide SF (NPSF) is a member of RFamide neuropeptides that play diverse
roles in central nervous system. Little is know about the effects of NPSF on
brain functions. Antidepressant-like effect of NPSF was studied in modified
mice FST. NPSF showed the antidepressant-like effects by decreasing the
immobility time and increasing the climbing and swimming time. Furthermore, the
involvement of the adrenergic, serotonergic, cholinergic or dopaminergic
receptors in the antidepressant-like effect of NPSF was studied in modified
mice FST. Mice were pretreated with a non-selective {\alpha}-adrenergic
receptor antagonist phenoxybenzamine, a \b{eta}-adrenergic receptor antagonist,
propranolol, a non-selective 5-HT2 serotonergic receptor antagonist,
cyproheptadine, nonselective muscarinic acetylcholine receptor antagonist,
atropine, or D2, D3, D4 dopamine receptor antagonist, haloperidol. The present
results confirmed that the antidepressant-like effect of NPSF is mediated, at
least in part, by an interaction of the {\alpha}-adrenergic, 5-HT2
serotonergic, muscarinic acetylcholine receptors and D2, D3, D4 dopamine
receptor in a modified mouse FST.


  In the article a virus transmission model is constructed on a simplified
social network. The social network consists of more than 2 million nodes, each
representing an inhabitant of Slovenia. The nodes are organised and
interconnected according to the real household and elderly-care center
distribution, while their connections outside these clusters are semi-randomly
distributed and fully-linked. The virus spread model is coupled to the disease
progression model. The ensemble approach with the perturbed transmission and
disease parameters is used to quantify the ensemble spread, a proxy for the
forecast uncertainty. The presented ongoing forecasts of COVID-19 epidemic in
Slovenia are compared with the collected Slovenian data. Results show that
infection is currently twice more likely to transmit within households/elderly
care centers than outside them. We use an ensemble of simulations (N = 1000) to
inversely obtain posterior distributions of model parameters and to estimate
the COVID-19 forecast uncertainty. We found that in the uncontrolled epidemic,
the intrinsic uncertainty mostly originates from the uncertainty of the virus
biology, i.e. its reproductive number. In the controlled epidemic with low
ratio of infected population, the randomness of the social network becomes the
major source of forecast uncertainty, particularly for the short-range
forecasts. Social-network-based models are thus essential for improving
epidemics forecasting.


  Prevention of protein aggregation and thus stabilization of proteins has
large biological and biotechnological implications. Here, we show that
inhibition of amyloid-like aggregates is possible in stoichiometric conjugates
of polymer surfactant and bovine serum albumin (BSA) chosen as a model protein.
We investigate using a combination of Thioflavin-T fluorescence spectroscopy,
dynamic light scattering and FTIR spectroscopy the aggregation behavior in
polymer surfactant modified and unmodified (native) BSA solutions. The
BSA-polymer surfactant conjugates are stable up to 5 days under aggregation
conditions, while native BSA forms amyloid fibrillar structures. Further,
DLS-based micro-rheology studies performed with heat-treated 100 to 200 {\mu}M
native BSA aggregates provided understanding of the equilibrium elastic and
viscous moduli over a very large frequency range, reaching MHz, which are
inaccessible using bulk rheology. Our results indicate that after 6 days of
aggregation conditions, elastic moduli showed values between 1.2 to 3.6 Pa
corresponding to an entanglement length ({\xi}) of 105 nm. Interestingly,
heating 200 {\mu}M native BSA solution at 65 degree C for 2 days in a plastic
Eppendorf resulted in self-standing films. These films exhibited strong
ThT-fluorescence intensity and a predominant \b{eta}-sheet secondary structure
from the FTIR studies, suggesting that self-standing microstructure resulted
from hierarchical self-assembly of amyloid fibrils.


  In this study three soil ecosystems, that differ in the type of management,
have been compared in the attempt to understand if and how anthropogenic action
affects them. The structure of the corresponding food webs was analyzed and
their robustness was calculated through the use of a dynamic model. With regard
to the topology, it has been found that the structure of all three networks is
small world. Furthermore, all three networks have a disassortative nature as
expected for foodwebs. The values of the clustering coefficient, of the
connectance and of the complexity, together with the calculation of the
robustness suggest that the ecosystem related to a fallowed pastures with low
pressure management is more robust than two ecosystems related to organic farms
subject to middle intensity management. If the evidence suggested by this study
was confirmed by further studies, the robustness shown by the networks could be
useful for evaluating, from an ecological point of view, the sustainability of
the agricultural practices to which the ecosystem is subject.


  Expanding the scope of graph-based, deep-learning models to noncovalent
protein-ligand interactions has earned increasing attention in structure-based
drug design. Modeling the protein-ligand interactions with graph neural
networks (GNNs) has experienced difficulties in the conversion of
protein-ligand complex structures into the graph representation and left
questions regarding whether the trained models properly learn the appropriate
noncovalent interactions. Here, we proposed a GNN architecture, denoted as
InteractionNet, which learns two separated molecular graphs, being covalent and
noncovalent, through distinct convolution layers. We also analyzed the
InteractionNet model with an explainability technique, i.e., layer-wise
relevance propagation, for examination of the chemical relevance of the model's
predictions. Separation of the covalent and noncovalent convolutional steps
made it possible to evaluate the contribution of each step independently and
analyze the graph-building strategy for noncovalent interactions. We applied
InteractionNet to the prediction of protein-ligand binding affinity and showed
that our model successfully predicted the noncovalent interactions in both
performance and relevance in chemical interpretation.


  We present a model analysis of the tumor and normal cell growth under the
influence of a carcinogenic agent, an immunomdulator (IM) and variable influx
of immune cells including relevant interactions. The tumor growth is
facilitated by carcinogens such as inorganic arsenic while the IM considered
here is black tea (Camellia sinesnsis). The model with variable influx of
immune cells is observed to have considerable advantage over the constant
influx model, and while the tumor cell population is greatly mitigated, normal
cell population remains above healthy levels. The evolution of normal and tumor
cells are computed from the proposed model and their local stabilities are
investigated analytically. Numerical simulations are performed to study the
long term dynamics and an estimation of the effects of various factors is made.
This helps in developing a balanced strategy for tumor mitigation without the
use of chemotherapeutic drugs that usually have strong side-effects.


  We consider a spatial (line) model for invasion of a population by a single
mutant with a stochastically selectively neutral fitness landscape, independent
from the fitness landscape for non-mutants. This model is similar to those
considered in Farhang-Sardroodi et al. [PLOS Comput. Biol., 13(11), 2017; J. R.
Soc. Interface, 16(157), 2019]. We show that the probability of mutant fixation
in a population of size $N$, starting from a single mutant, is greater than
$1/N$, which would be the case if there were no variation in fitness
whatsoever. In the small variation regime, we recover precise asymptotics for
the success probability of the mutant. This demonstrates that the introduction
of randomness provides an advantage to minority mutations in this model, and
shows that the advantage increases with the system size. We further demonstrate
that the mutants have an advantage in this setting only because they are better
at exploiting unusually favorable environments when they arise, and not because
they are any better at exploiting pockets of favorability in an environment
that is selectively neutral overall.


  A dynamical epidemic model optimized using genetic algorithm and cross
validation method to overcome the overfitting problem is proposed. The cross
validation procedure is applied so that available data are split into a
training subset used to fit the algorithm's parameters, and a smaller subset
used for validation. This process is tested on the countries of Italy, Spain,
Germany and South Korea before being applied to Algeria. Interestingly, our
study reveals an inverse relationship between the size of the training sample
and the number of generations required in the genetic algorithm. Moreover, the
enhanced compartmental model presented in this work is proven to be a reliable
tool to estimate key epidemic parameters and non-measurable asymptomatic
infected portion of the susceptible population in order to establish realistic
nowcast and forecast of epidemic's evolution. The model is employed to study
the COVID-19 outbreak dynamics in Algeria between February 25th and May 24th,
2020. The basic reproduction number and effective reproduction number on May
24th, after three months of the outbreak, are estimated to be 3.78 (95% CI
3.033-4.53) and 0.651 (95% CI 0.539-0.761) respectively. Disease incidence, CFR
and IFR are also calculated. Numerical programs developed for the purpose of
this study are made publicly accessible for reproduction and further use.


  In this paper a Bayesian SEIR model is studied to estimate the proportion of
the population infected with SARS-CoV-2, the virus responsible for COVID-19. To
capture heterogeneity in the population and the effect of interventions to
reduce the rate of epidemic spread, the model uses a time-varying contact rate,
whose logarithm has a Gaussian process prior. A Poisson point process is used
to model the occurrence of deaths due to COVID-19 and the model is calibrated
using data of daily death counts in combination with a snapshot of the the
proportion of individuals with an active infection, performed in Stockholm in
late March. The methodology is applied to regions in Sweden. The results show
that the estimated proportion of the population who has been infected is around
13.5% in Stockholm, by 2020-05-15, and ranges between 2.5% - 15.6% in the other
investigated regions. In Stockholm where the peak of daily death counts is
likely behind us, parameter uncertainty does not heavily influence the expected
daily number of deaths, nor the expected cumulative number of deaths. It does,
however, impact the estimated cumulative number of infected individuals. In the
other regions, where random sampling of the number of active infections is not
available, parameter sharing is used to improve estimates, but the parameter
uncertainty remains substantial.


  Physical unclonable functions (PUFs) are small circuits that are widely used
as hardware security primitives for authentication. These circuits can generate
unique signatures because of the inherent randomness in manufacturing and
process variations. This paper introduces molecular PUFs based on multiplexer
(MUX) PUFs using dual-rail representation. It may be noted that molecular PUFs
have not been presented before. Each molecular multiplexer is synthesized using
16 molecular reactions. The intrinsic variations of the rate constants of the
molecular reactions are assumed to provide inherent randomness necessary for
uniqueness of PUFs. Based on Gaussian distribution of the rate constants of the
reactions, this paper simulates intra-chip and inter-chip variations of linear
molecular MUX PUFs containing 8, 16, 32 and 64 stages. These variations are,
respectively, used to compute reliability and uniqueness. It is shown that, for
the rate constants used in this paper, although 8-state molecular MUX PUFs are
not useful as PUFs, PUFs containing 16 or higher stages are useful as molecular
PUFs. Like electronic PUFs, increasing the number of stages increases
uniqueness and reliability of the PUFs


  Currently, there is no effective antiviral drugs nor vaccine for coronavirus
disease 2019 (COVID-19) caused by acute respiratory syndrome coronavirus 2
(SARS-CoV-2). Due to its high conservativeness and low similarity with human
genes, SARS-CoV-2 main protease (M$^{\text{pro}}$) is one of the most favorable
drug targets. However, the current understanding of the molecular mechanism of
M$^{\text{pro}}$ inhibition is limited by the lack of reliable binding affinity
ranking and prediction of existing structures of M$^{\text{pro}}$-inhibitor
complexes. This work integrates mathematics and deep learning (MathDL) to
provide a reliable ranking of the binding affinities of 92 SARS-CoV-2
M$^{\text{pro}}$ inhibitor structures. We reveal that Gly143 residue in
M$^{\text{pro}}$ is the most attractive site to form hydrogen bonds, followed
by Cys145, Glu166, and His163. We also identify 45 targeted covalent bonding
inhibitors. Validation on the PDBbind v2016 core set benchmark shows the MathDL
has achieved the top performance with Pearson's correlation coefficient ($R_p$)
being 0.858. Most importantly, MathDL is validated on a carefully curated
SARS-CoV-2 inhibitor dataset with the averaged $R_p$ as high as 0.751, which
endows the reliability of the present binding affinity prediction. The present
binding affinity ranking, interaction analysis, and fragment decomposition
offer a foundation for future drug discovery efforts.


  SARS-CoV-2 causing COVID-19 disease has moved rapidly around the globe,
infecting millions and killing hundreds of thousands. The basic reproduction
number, which has been widely used and misused to characterize the
transmissibility of the virus, hides the fact that transmission is stochastic,
is dominated by a small number of individuals, and is driven by super-spreading
events (SSEs). The distinct transmission features, such as high stochasticity
under low prevalence, and the central role played by SSEs on transmission
dynamics, should not be overlooked. Many explosive SSEs have occurred in indoor
settings stoking the pandemic and shaping its spread, such as long-term care
facilities, prisons, meat-packing plants, fish factories, cruise ships, family
gatherings, parties and night clubs. These SSEs demonstrate the urgent need to
understand routes of transmission, while posing an opportunity that outbreak
can be effectively contained with targeted interventions to eliminate SSEs.
Here, we describe the potential types of SSEs, how they influence transmission,
and give recommendations for control of SARS-CoV-2.


  Using a beta-hairpin protein as a representative example of two-state
folders, we studied how the exploration of native-like states affects the
folding kinetics. It has been found that the first-passage time (FPT)
distributions are essentially single-exponential not only for the times to
overcome the free energy barrier that separates unfolded and native-like states
but also for the times to find the native state among the native-like ones. If
the protein explores native-like states for a time much longer than the time to
overcome the free energy barrier, which was found to be characteristic of high
temperatures, the resulting FPT distribution to reach the native state remains
close to exponential but the mean FPT (MFPT) is determined not by the height of
the free energy barrier but by the time to explore native-like states. The mean
time to overcome the free energy barrier is found to be in reasonable agreement
with the Kramers rate formula and generally far shorter than the MFPT to reach
the native state. The time to find the native state among native-like ones
increases with temperature, which explains the known U-shape dependence of the
MFPTs on temperature.


  Following the highly restrictive measures adopted by many countries for
combating the current pandemic, the number of individuals infected by
SARS-CoV-2 and the associated number of deaths is steadily decreasing. This
fact, together with the impossibility of maintaining the lockdown indefinitely,
raises the crucial question of whether it is possible to design an exit
strategy based on quantitative analysis. Guided by rigorous mathematical
results, we show that this is indeed possible: we present a robust numerical
algorithm which can compute the cumulative number of deaths that will occur as
a result of increasing the number of contacts by a given multiple, using as
input only the most reliable of all data available during the lockdown, namely
the cumulative number of deaths.


  All genes interact with other genes, and their additive effects and epistatic
interactions affect an organism's phenotype and fitness. Recent theoretical and
empirical work has advanced our understanding of the role of multi-locus
interactions in speciation. However, relating different models to one another
and to empirical observations is challenging. This review focuses on
multi-locus interactions that lead to reproductive isolation (RI) through
reduced hybrid fitness. We first review theoretical approaches and show how
recent work incorporating a mechanistic understanding of multi-locus
interactions recapitulates earlier models, but also makes novel predictions
concerning the build-up of RI. These include high variance in the build-up rate
of RI among taxa, the emergence of strong incompatibilities producing localised
barriers to introgression, and an effect of population size on the build-up of
RI. We then review recent experimental approaches to detect multi-locus
interactions underlying RI using genomic data. We argue that future studies
would benefit from overlapping methods like Ancestry Disequilibrium scans,
genome scans of differentiation and analyses of hybrid gene expression.
Finally, we highlight a need for further overlap between theoretical and
empirical work, and approaches that predict what kind of patterns multi-locus
interactions resulting in incompatibilities will leave in genome-wide
polymorphism data.


  The current outbreak of the coronavirus disease 2019 (COVID-19) is an
unprecedented example of how fast an infectious disease can spread around the
globe (especially in urban areas) and the enormous impact it causes on public
health and socio-economic activities. Despite the recent surge of
investigations about different aspects of the COVID-19 pandemic, we still know
little about the effects of city size on the propagation of this disease in
urban areas. Here we investigate how the number of cases and deaths by COVID-19
scale with the population of Brazilian cities. Our results indicate small towns
are proportionally more affected by COVID-19 during the initial spread of the
disease, such that the cumulative numbers of cases and deaths per capita
initially decrease with population size. However, during the long-term course
of the pandemic, this urban advantage vanishes and large cities start to
exhibit higher incidence of cases and deaths, such that every 1% rise in
population is associated with a 0.14% increase in the number of fatalities per
capita after about four months since the first two daily deaths. We argue that
these patterns may be related to the existence of proportionally more health
infrastructure in the largest cities and a lower proportion of older adults in
large urban areas. We also find the initial growth rate of cases and deaths to
be higher in large cities; however, these growth rates tend to decrease in
large cities and to increase in small ones over time.


  For decades, it was unknown how electron bifurcating systems in Nature
prevented energy-wasting short-circuiting reactions that have large driving
forces, so synthetic electron bifurcating molecular machines could not be
designed and built. The underpinning free energy landscapes for electron
bifurcation were also enigmatic. We predict that a simple and universal free
energy landscape enables electron bifurcation, and we show that it enables
high-efficiency bifurcation with limited short-circuiting (the EB-scheme). The
landscape relies on steep free energy slopes in the two redox branches to
insulate against short-circuiting without relying on nuanced changes in the
microscopic rate constants for the short-circuiting reactions. The EB-scheme
thus provides a blueprint for future campaigns to establish synthetic electron
bifurcating machines.


  Cells regulate themselves via dizzyingly complex biochemical processes called
signaling pathways. These are usually depicted as a network, where nodes
represent proteins and edges indicate their influence on each other. In order
to understand diseases and therapies at the cellular level, it is crucial to
have an accurate understanding of the signaling pathways at work. Since
signaling pathways can be modified by disease, the ability to infer signaling
pathways from condition- or patient-specific data is highly valuable. A variety
of techniques exist for inferring signaling pathways. We build on past works
that formulate signaling pathway inference as a Dynamic Bayesian Network
structure estimation problem on phosphoproteomic time course data. We take a
Bayesian approach, using Markov Chain Monte Carlo to estimate a posterior
distribution over possible Dynamic Bayesian Network structures. Our primary
contributions are (i) a novel proposal distribution that efficiently samples
sparse graphs and (ii) the relaxation of common restrictive modeling
assumptions. We implement our method, named Sparse Signaling Pathway Sampling,
in Julia using the Gen probabilistic programming language. Probabilistic
programming is a powerful methodology for building statistical models. The
resulting code is modular, extensible, and legible. The Gen language, in
particular, allows us to customize our inference procedure for biological
graphs and ensure efficient sampling. We evaluate our algorithm on simulated
data and the HPN-DREAM pathway reconstruction challenge, comparing our
performance against a variety of baseline methods. Our results demonstrate the
vast potential for probabilistic programming, and Gen specifically, for
biological network inference. Find the full codebase at
https://github.com/gitter-lab/ssps


  The COVID-19 pandemic led several countries to resort to social distancing,
the only known way to slow down the spread of the virus and keep the health
system under control. Here we use an individual based model (IBM) to study how
the duration, start date and intensity of quarantine affect the height and
position of the peak of the infection curve. We show that stochastic effects,
inherent to the model dynamics, lead to variable outcomes for the same set of
parameters, making it crucial to compute the probability of each result. To
simplify the analysis we divide the outcomes in only two categories, that we
call {best and worst scenarios. Although long and intense quarantine is the
best way to end the epidemic, it is very hard to implement in practice. Here we
show that relatively short and intense quarantine periods can also be very
effective in flattening the infection curve and even killing the virus, but the
likelihood of such outcomes are low. Long quarantines of relatively low
intensity, on the other hand, can delay the infection peak and reduce its size
considerably with more than 50% probability, being a more effective policy than
complete lockdown for short periods.


  Power-law behaviors are common in many disciplines, especially in network
science. Real-world networks, like disease spreading among people, are more
likely to be interconnected communities, and show richer power-law behaviors
than isolated networks. In this paper, we look at the system of two communities
which are connected by bridge links between a fraction $r$ of bridge nodes, and
study the effect of bridge nodes to the final state of the
Susceptible-Infected-Recovered model, by mapping it to link percolation. By
keeping a fixed average connectivity, but allowing different transmissibilities
along internal and bridge links, we theoretically derive different power-law
asymptotic behaviors of the total fraction of the recovered $R$ in the final
state as $r$ goes to zero, for different combinations of internal and bridge
link transmissibilities. We also find crossover points where $R$ follows
different power-law behaviors with $r$ on both sides when the internal
transmissibility is below but close to its critical value, for different bridge
link transmissibilities. All of these power-law behaviors can be explained
through different mechanisms of how finite clusters in each community are
connected into the giant component of the whole system, and enable us to pick
effective epidemic strategies and to better predict their impacts.


  We portray the evolution of the Covid-19 epidemic during the crisis of
March-April 2020 in the Paris area, by analyzing the medical emergency calls
received by the EMS of the four central departments of this area (Centre 15 of
SAMU 75, 92, 93 and 94). Our study reveals strong dissimilarities between these
departments. We show that the logarithm of each epidemic observable can be
approximated by a piecewise linear function of time. This allows us to
distinguish the different phases of the epidemic, and to identify the delay
between sanitary measures and their influence on the load of EMS. This also
leads to an algorithm, allowing one to detect epidemic resurgences. We rely on
a transport PDE epidemiological model, and we use methods from Perron-Frobenius
theory and tropical geometry.


  Current drug discovery is expensive and time-consuming. It remains a
challenging task to create a wide variety of novel compounds with desirable
pharmacological properties and cheaply available to low-income people. In this
work, we develop a generative network complex (GNC) to generate new drug-like
molecules based on the multi-property optimization via the gradient descent in
the latent space of an autoencoder. In our GNC, both multiple chemical
properties and similarity scores are optimized to generate and predict
drug-like molecules with desired chemical properties. To further validate the
reliability of the predictions, these molecules are reevaluated and screened by
independent 2D fingerprint-based predictors to come up with a few hundreds of
new drug candidates. As a demonstration, we apply our GNC to generate a large
number of new BACE1 inhibitors, as well as thousands of novel alternative drug
candidates for eight existing market drugs, including Ceritinib, Ribociclib,
Acalabrutinib, Idelalisib, Dabrafenib, Macimorelin, Enzalutamide, and
Panobinostat.


  Clustering of ligand-binding receptors of different types on thickened isles
of the cell membrane, namely lipid rafts, is an experimentally observed
phenomenon. Although its influence on cell response is deeply investigated, the
role of the coupling between mechanical processes and multiphysics involving
the active receptors and the surrounding lipid membrane during ligand-binding
has not yet been understood. Specifically, the focus of this work is on
G-protein-coupled receptors (GPCRs), which regulate specific cell processes
through chemical signalling pathways involving a synergistic balance between
the cyclic Adenosine Monophosphate (cAMP) produced by active GPCRs in the
intracellular environment and its efflux, mediated by Multidrug Resistance
Proteins (MRPs) transporters. This paper develops a multiphysics approach based
on the interplay among energetics, multiscale geometrical changes and mass
balance of species, i.e. active GPCRs and MRPs, including diffusion and
kinetics of binding and unbinding. Because the obtained energy depends upon
both the kinematics and the changes of species densities, balance of mass and
of linear momentum are coupled and govern the space-time evolution of the cell
membrane. The mechanobiology involving remodelling and change of lipid ordering
of the cell membrane allows to predict dynamics of transporters and active
receptors -in agreement with experimentally observed cAMP levels- and how the
latter trigger rafts formation and cluster on such sites. Within the current
scientific debate on Severe Acute Respiratory Syndrome CoronaVirus 2
(SARS-CoV-2) and on the basis of the ascertained fact that lipid rafts often
serve as an entry port for viruses, it is felt that approaches accounting for
strong coupling among mechanobiological aspects could even turn helpful in
better understanding membrane-mediated phenomena such as COVID-19 virus-cell
interaction.


  Aquatic biospheres reliant on oxygenic photosynthesis are expected to play an
important role on Earth-like planets endowed with large-scale oceans insofar as
carbon fixation (i.e., biosynthesis of organic compounds) is concerned. We
investigate the properties of aquatic biospheres comprising Earth-like biota
for habitable rocky planets orbiting Sun-like stars and late-type M-dwarfs such
as TRAPPIST-1. In particular, we estimate how these characteristics evolve with
the available flux of photosynthetically active radiation (PAR) and the ambient
ocean temperature ($T_W$), the latter of which constitutes a key environmental
variable. We show that many salient properties, such as the depth of the
photosynthesis zone and the net primary productivity (i.e., the effective rate
of carbon fixation), are sensitive to PAR flux and $T_W$ and decline
substantially when the former is decreased or the latter is increased. We
conclude by exploring the implications of our analysis for exoplanets around
Sun-like stars and M-dwarfs.


  We study the influence of the rate of the attainment of herd immunity (HI),
in the absence of an approved vaccine, on the vulnerable population. We
essentially ask the question: how hard the evolution towards the desired herd
immunity could be on the life of the vulnerables? We employ mathematical
modelling (chemical network theory) and cellular automata based computer
simulations to study the human cost of an epidemic spread and an effective
strategy to introduce HI. Implementation of different strategies to counter the
spread of the disease requires a certain degree of quantitative understanding
of the time dependence of the outcome. In this paper, our main objective is to
gather understanding of the dependence of outcome on the rate of progress of
HI. We generalize the celebrated SIR model (Susceptible-Infected-Removed) by
compartmentalizing the susceptible population into two categories- (i)
vulnerables and (ii) resilients, and study dynamical evolution of the disease
progression. We achieve such a classification by employing different rates of
recovery of vulnerables vis-a-vis resilients. We obtain the relative fatality
of these two sub-categories as a function of the percentages of the vulnerable
and resilient population, and the complex dependence on the rate of attainment
of herd immunity. Our results quantify the adverse effects on the recovery
rates of vulnerables in the course of attaining the herd immunity. We find the
important result that a slower attainment of the HI is relatively less fatal.
However, a slower progress towards HI could be complicated by many intervening
factors.