## Extreme-value analysis in nano-biological systems: Applications and Implications

Kumiko Hayashi 1 , Nobumichi Takamatsu  and Shunki Takaramoto 1 1

1 Institute for Solid State Physics, The University of Tokyo, Kashiwano-ha 5-1-5, Kashiwa, Chiba 277-8581,

Japan

*Correspondence to Kumiko Hayashi hayashi@issp.u-tokyo.ac.jp

## ABSTRACT

Extreme value analysis (EVA) is a statistical method that studies the properties of extreme values of datasets, crucial for fields like engineering, meteorology, finance, insurance, and environmental science. EVA models extreme  events  using  distributions  such  as  FrŽchet,  Weibull,  or  Gumbel,  aiding  in  risk  prediction  and management. This review explores EVAÕs application to nanoscale biological systems. Traditionally, biological research focuses on average values from repeated experiments. However, EVA offers insights into molecular mechanisms by examining extreme data points. We introduce EVAÕs concepts with simulations and review its use in studying motor protein movements within cells, highlighting the importance of in vivo analysis due to the  complex  intracellular  environment.  We  suggest  EVA  as  a  tool  for  extracting  motor  proteinsÕ  physical properties in  vivo and  discuss  its  potential  in  other  biological  systems.  While  there  have  been  only  a  few applications of EVA to biological systems, it holds promise for uncovering hidden properties in extreme data, promoting its broader application in life sciences.

## Introduction

Extreme value analysis (EVA) is a branch of statistics that focuses on extreme values. It is the study of the statistical properties of particularly large or small values within a data set. EVA is widely applied in various fields where extreme phenomena hold significant importance (Coles 2001; Gilleland and Katz 2016), such as disaster prevention (de Haan and Ferreira 2006; Tippett et al. 2016), finance (Kratz 2017), safety estimation (Songchitruksa  and  Tarko  2006),  sports  (Einmahl  and  Magnus  2008;  Gembris  et  al.  2002),  human  lifespan (Dong et al. 2016; Rootzen and Zolud 2017), and the recent pandemic (Wong and Collins 2020). Recently, its applications in the biological data analysis has also become active (Basnayake et al. 2019; Tsuduki 2024). The central concepts of EVA involve identifying the largest or smallest values in a dataset, determining whether these extremes follow distributions like the FrŽchet distribution, Weibull distribution, or Gumbel distribution, and analyzing data points that exceed specific thresholds to estimate their distribution. EVA is utilized in a broad range of fields. For instance, in meteorology, it is used to predict and assess extreme weather events such as typhoons, floods, and droughts. In finance, it aids in risk management for extreme price fluctuations, such as market crashes or surges. In engineering, it is applied to analyze extreme stresses or loads for evaluating the durability and safety design of structures. EVA serves as a powerful tool for predicting the frequency and impact of extreme events, aiding in risk management and safety measures.

In this review paper, we discuss the application of EV A to the analysis of experimental data in nanoscale biological systems. By using EVA, it is possible to explore novel molecular mechanisms in life sciences through the  behaviour  of  extreme  values  in  the  data.  Rare  events  in  nanoscale  biological  systems,  such  as  protein misfolding  (Harada  et  al.  2015)  and  synaptic  delay  (Tsuduki  2024),  have  already  been  studied,  and  their prediction is considered important due to their relevance to diseases. In the following sections, we first learn the concepts of EVA through simple simulations using random numbers. Next, we introduce the application of EVA to the movement of motor proteins within cells, investigated in our recent research paper (Naoi et al. 2024). The physical properties of motor proteins, such as force and velocity, have been investigated by in-vitro singlemolecule experiments, in which the functions of motor proteins consisting of minimal complexes were analyzed in glass chambers (Brenner et al. 2020; Elshenawy et al. 2019; Gennerich et al. 2007; Hirakawa et al. 2000; Mallik et al. 2004; Schnitzer et al. 2000; Toba et al. 2006). However, because motor proteins function fully in the intracellular environment and are equipped with accessory proteins, the investigation of motor proteins in vivo is  as  important  as in-vitro single-molecule  experiments.  We  believe  that  EV A  can  be  useful  as  a  new information science tool for extracting the physical properties of motor proteins in the complex intracellular environment.

Finally,  we  discuss  the  potential  applications  of  EV A  in  nanoscale  biological  systems,  with  a  specific example of its application to the analysis of droplet sizes in liquid-liquid phase separation (LLPS) (Takaramoto and Inoue 2024). This application of EVA was inspired by Takaramoto and InoueÕs poster at the IUPAB 2024 Congress.  While  research  utilizing  EV A  in  nanoscale  biological  systems  is  still  limited,  it  is  believed  that undiscovered properties may lie hidden within the extreme values of experimental data.

## Simulation using random numbers

In  this  section,  the  analytical  methods  of  EV A  are  reviewed  through  a  simple  simulation  using  uniform random numbers ranging from 0 to 1 and Gaussian random numbers. In the case of random numbers ranging from 0 to 1, it is obvious that the finite maximum value is 1, while Gaussian random numbers do not have a finite maximum value. In other words, we examine two random numbers with different extreme properties.

One block of EVA is considered as a data set ( i.e. , a block) of elements 𝑀 , from which the largest value is selected. In meteorology, 𝑀 is often considered as one year, which is 365 days. The selection of the block maximum ( 𝑋#$% &amp; )  from the 𝑖 -th block (the 𝑖 -th data set) is repeated 𝑛 times to collect the block maximum

data set )𝑋#$%* &amp; ( 𝑖 = 1, ⋯ , 𝑛 ). Using the block maximum data set )𝑋#$%* &amp; , the return level plot is investigated (the plot was calculated by using the ismev and evd packages in R ( R Core Team 2018 )) (Fig. 1a and 1b). The two axes of the return level plot (Fig. 1a and 1b, bottom) represent the return period 𝑟 0 and return level 𝑧0 . For  a  given  probability 𝑝 , 𝑟 0 = -{log	(1 - 𝑝)	} &lt;= ,  and 𝑧0 is  defined  by  the  generalized  extreme  value distribution as 1 - 𝑝 = 𝐺(𝑧0 ) , where

$$G ( z _ { p } ) = \exp \left [ - \left \{ 1 + \xi \left ( \frac { z _ { p } - \mu } { \sigma } \right ) \right \} ^ { - 1 / \xi } \right ].$$

Note that 𝑧0 represents 𝑋 Q #$% &amp; ,  where )𝑋 Q &amp; #$%* is  the  rearranged  data  of )𝑋#$%* &amp; ,  such  that 𝑋 Q #$% = ≤ 𝑋 Q S #$% ≤ ⋯≤𝑋 Q T #$% . Roughly, 𝑟 0 is the sample number. (Note that we obtained parameters of the generalized extreme value distribution 𝜉 , 𝜇 , and 𝜎 in Eq. (1) by using the ismev and evd packages in R ( R Core Team 2018 )) . Here, the return level 𝑧0 refers to the magnitude of extreme values expected to occur once within the sample size we are considering, such as 100 years, 1000 people and so on. For example, when analyzing river water level data for environmental science case, its return level plot can provide information such as the maximum water level of a flood expected once every 100 years. The distribution characteristics of the extreme value data, such as the values of location ( 𝜇 ), scale ( 𝜎 ), and shape parameters ( 𝜉 ) of the extreme value distribution (the Fr é chet distribution, Weibull distribution, or Gumbel distribution) can be estimated from the behavior of the return level plot (Fig. 1c).

The shape parameters of the extreme value distribution 𝜉 is important because it characterizes the behavior of extreme values, determining whether a data set has a finite maximum value or not. We found that 𝜉	 &lt; 0 for the case of uniform random numbers, which have the finite maximum value (Fig. 1a), and 𝜉	~0 for the case of Gaussian random numbers, which do not have the finite maximum value (Fig. 1b). The existence of a finite maximum value is related to the value of the shape parameter 𝜉 . Therefore, we particularly focus on the value of 𝜉 in EVA.

Depending on the value of 𝜉 , it is classified into the Weibull distribution ( 𝜉 &lt; 0 ), Gumbel distribution ( 𝜉 = 0 ), or Fr é chet distribution ( 𝜉 &gt; 0 ) (Fig. 1c). In the case of a Weibull distribution, the finite extreme value (the maximum value) 𝑉 [% is proved to exist and is estimated using the following equation:

$$V _ { \text{ex} } = \mu - \sigma / \xi$$

It  has been known that the extreme values (the maximum values) of a uniform random number distribution converge to a Weibull distribution (Fig. 1a), and those of a Gaussian distribution to Gumbel distribution (Fig. 1b). In our simulations, we  obtained the set of parameters (𝜇, 𝜎, 𝜉) = (0.911 ± 0.003SE, 0.0793 ± 0.0034SE, -0.890 ± 0.012SE) for uniform random numbers ranging from 0 to 1. For the set of parameters (𝜇, 𝜎, 𝜉) , the quantity 𝜇 - 𝜎/𝜉 (Eq. (2)) was calculated to be 1.00 as the maximum number.

The  properties  of  the  maximum  value  of  the  simulation  examples  shown  in  Fig.  1  are  known  from  the beginning.  However, the  important  point is  that  the  extreme  value  distribution  of  any  data,  whose  extreme properties are not known, can be described by Eq. (1), based on the mathematical theorem. In unknown complex systems like nanoscale biological systems, it is a significant advantage to know the functional form (Eq. (1)) for their extreme value distributions. In the next section, we consider the application of this strength in EV A to intracellular motor proteins.

## Applications of EVA to cargo transport by motor proteins.

Motor protein is a general term for proteins that move and function using energy obtained from adenosine triphosphate (ATP) hydrolysis. Here, we focus on kinesin and dynein among such motor proteins, which are responsible for cargo transport and form the basis of intracellular material logistics (Fig. 2a). Particularly, we investigate anterograde transport by kinesin and retrograde transport by dynein in the axons of neurons. In the

logistics of neurons, both anterograde transport, which carries synaptic materials to the terminals, and retrograde transport, which collects waste materials, are important. Because logistics in long axons is particularly important for  neurons,  and  disruptions  in  this  process  are  often  associated  with  neurological  diseases  (GuedesDias, Holzbaur  2019;  Keefe  et  al.  2023),  the  application  of  EV A  to  axonal  transport  is  considered  to  be significant to help elucidate disease mechanisms in the future by increasing in vivo physical quantities we can estimate.

Cargo motion can be observed using fluorescence microscopy (Fig. 2a) (Hayashi et al. 2018; Naoi et al. 2024). Then  the  transport  velocity  is  measured  by  analyzing  the  recorded  movies  obtained  from  the  fluorescence microscopy. Note that fluorescence observation is a standard technique in biophysics, and this review will not discuss the methods of fluorescence observation. In in vivo cases, the measured velocity values varied affected by the size of a cargo and the number of motor proteins involved in the cargo transport (Fig. 3a). Especially since cargo sizes can differ by a factor of 10, the average velocity value reflects the cargo size rather than the performance of the motor proteins due to high viscosity in cells, unlike the situations in in vitro single-molecule experiments.  The  differences  in  the  mechanical  properties  of  the  two  different  motor  proteins,  kinesin  and dynein, are also less likely to be reflected in the velocity values, because the influence of cargo size on the velocity values is larger than the differences in the two motors. Then, the difference and true performance of the motor proteins can be considered to appear in the limit of small cargo sizes. Under this condition, it is believed that the difference and true performance of motor proteins can be accurately evaluated using EV A. In this review, in the following, we introduce the application of EV A to axonal transport in two different species: C. elegans worms and mice.

## C. elegans worms

The measurement method of transport velocity using fluorescence imaging in living C. elegans worms, was explained in the original paper (Naoi et al. 2024), noting that fluorescent proteins were labeled on the cargo for the fluorescence microscopy (Fig. 2a). Approximately 2,000 velocity data were collected from about 200 worms for anterograde and retrograde transport. One thing to be careful about is the fact that a measurement difficulty in live worms caused the small size of 𝑀 𝑀 = 10 ( for this experiment), which is typically set to of the order of 1000 to obtain a correct extreme value distribution, highlighting a key problem in the application of EV A to nanoscale biological systems, which is discussed in the Discussion section.

Using the block maximum dataset )𝑣 #$% * &amp; ( 𝑖 = 1, ⋯ , 𝑛 where 𝑛 =228 for anterograde transport and 𝑛 =217 for retrograde transport), the return-level plot was calculated (it was calculated concretely by using the ismev and  evd  packages  in  R ( R  Core  Team  2018 )) (Fig.  2b). 𝜉	 &lt; 0 for  anterograde  transport  and 𝜉	~0 for retrograde  transport.  For  anterograde  transport,  the  return  level  plot  shows  a  convergent  behavior  as 𝑟 0 becomes larger, a property specific to a Weibull distribution ( 𝜉	 &lt; 0 ), and the extreme value 𝑉 [% was estimated to be 4.0 ± 0.4 ∝ m/s using Eq. (2).

Because the return-level plot of the retrograde transport (Fig. 2b) shows 𝜉~	0 , and 𝑉 [% cannot be estimated from Eq. (2) for retrograde transport by dynein, noting that the finite maximum exists only for the case 𝜉	 &lt; 0 . In the following simulation section, we explore why dynein does not exhibit a behavior of Weibull distribution.

## Mice

To investigate a mammalian case, EVA was also applied to examine the velocities of synaptic cargo transport by motor proteins in mouse hippocampal neurons (Fig. 2c, left), as originally reported in our previous study (Hayashi  et  al  2021).  The  return-level  plot  with 𝜉 &lt; 0 was  observed  for  anterograde  transport.  Then  the maximum velocity for anterograde transport in mice hippocampal neurons was calculated to be 5.6 ± 1.4 ∝ m/s using Eq. (2). The maximum velocity was higher than that of the worms. Comparing the maximum velocities among  various  species  is  a  future  important  issue.  Return-level  plot  with 𝜉 &gt; 0 was  also  observed  for

retrograde transport (Fig. 2c, right). It is known that dynein exhibits different mechanical properties depending on the species, but it is expected to have similar mechanical properties in both worms and mice, judging from the similar behaviors of the return level plots.

## N umerical simulation

In our previous paper on the application of EVA to motor proteins (Naoi et al. 2024), the opposite signs of 𝜉 shown in the return level plots (Fig. 2b, c) was attributed to the different force (load)-velocity relationship between kinesin and dynein. The difference in the convexity of force (load)-velocity relationships (Fig. 3b), which have been clarified based on in-vitro single-molecule studies using optical tweezers, reflects the different mechanisms  underlying  the  walking  behavior  of  motor  proteins  on  microtubules,  although  the  biological meanings of convexity are not clearly understood yet. Due to high viscosity in cells, a motor protein experiences a load proportional to the cargo size transported by it. Therefore, we predicted that the behavior of velocity values in vivo would be influenced by the force (load)-velocity relationship. The observed extreme velocity values, which are close to those under the zero-load condition, are particularly related to the properties of the force (load)-velocity relationship under low load condition areas (green and yellow area depicted in Fig. 3b, right). The property that 𝜉~0 for the retrograde velocity data was attributed to the fact that the force (load)velocity relationship was concave. The steep velocity decrease for the case of the concave force (load)-velocity relationship in the low-load condition (green area in Fig. 3b) caused a major variation in the larger velocity values, and this behavior tends to generate a major variation in velocity. Then large gaps generated between 𝑣 fg#,#$% &amp; and 𝑣 fg#,#$% &amp;h= for a large 𝑖 in this case. This gap caused the non-convergence of 𝑣#$% &amp; and 𝜉~0 as a result. In the first place, the occurrence of small cargos, which are thought to generate high velocities, is a rare event. As a result, in a concave force (load)-velocity relationship, large velocity values are less likely to converge due to this rarity.

We can actually reproduce the behavior observed for anterograde and retrograde transport using a theoretical model of the force (load)-velocity relationship derived based on the mechanisms underlying ATP hydrolysis by motor  proteins  (Sasaki  et  al.  2018).  The  force  (load)-velocity  relationship  for  the  model  is  represented  as follows:

𝑣

$$\nu _ { \text{three} } ( F ) = ( k _ { 0 1 } - k _ { 0 2 } ) l \\ \frac { 1 } { k _ { 0 1 } } = \frac { 1 } { \kappa _ { 1 } } + \frac { 1 } { \lambda _ { 1 } } e ^ { d _ { 1 } F / k _ { B } T } \\ \frac { 1 } { k _ { 0 2 } } = \frac { 1 } { \lambda _ { 2 } } e ^ { - d _ { 2 } F / k _ { B } T } \quad \ ( 3 )$$

See reference (Sasaki et al. 2018) for the definitions and values of parameters. The differences in the model parameters  (Eq.  (3))  for  kinesin  and  dynein  were  resulted  in  the  different  convexities  of  the  force-velocity relationship. Fig. 3c (top) represents the )𝑣 fg# &amp; * (4000 data) obtained from the simulations using the model (Eq. (3)). )𝑣 fg#,#$% &amp; * were chosen from )𝑣 fg# &amp; * ( 𝑀 = 10 ). Using these )𝑣 fg#,#$% &amp; * ( 𝑛 = 400 ), we calculated the return-level plots for kinesin and dynein (Fig. 3c, bottom). We found that the tendency that 𝑧0 ( = 𝑣 fg#,#$% &amp; arranged in ascending order) did not converge for a large 𝑟 0 (Roughly, 𝑟 0 is the sample number (1 ≤ 𝑟 0 ≤ 400)) in the case of dynein, i.e .,  the  gaps created between 𝑣 fg#,#$% &amp; and 𝑣 fg#,#$% &amp;h= for a large 𝑖 . This is because a large velocity value is likely to be generated in the case of a concave force (load)-velocity relationship, owing to its steep slope in the load-sensitive regime.

## Discussion

Typically, a block size (the number of elements in a block, 𝑀 ) of around 1000 is used to obtain a correct extreme value distribution (Coles 2001). In comparison, our study uses an extremely small 𝑀 , around 10. This is due to the limited cargo transport observable in a single C. elegans worm, resulting in an 𝑀 value of around 10. It is difficult to collect a large number of samples within living organisms or cells. This is a fundamental issue  in  biological  systems.  In  our  research,  we  checked  that  while  the  parameters 𝜇 , 𝜎 ,  and 𝜉 showed dependency  on 𝑀 due  to  its  small  size,  the  qualitative  results  of 𝜉 &lt; 0  for  anterograde  and 𝜉 &gt; 0  for retrograde transport were independent of 𝑀 (Naoi et al. 2024). Due to the small sample size, it may be helpful to use other statistical methods, such as a bootstrap method, to cope with outliers of data when we apply EV A (Naoi et al. 2024).

There  is  one  thing  that  should  be  mentioned  regarding  the  conclusion  obtained  from  the  numerical simulation  depicted  in  Fig.  3c.  We  attributed  the  property 𝜉~0 in  the  extreme  value  data  for  retrograde transport (specifically the rare occurrence of large velocity values and the divergence observed in the return level plot at large 𝑟 0 values)  to  a  concave  force  (load)-velocity  relationship.  However, several mechanisms could  lead  to  the  occasional  large  velocity  values  in  retrograde  transport,  such  as  active  non-equilibrium fluctuations  originating  from  the  cellular  cytoskeleton  (Chaubet  et  al.  2020;  Mizuno  et  al.  2007)  and  the participation of multiple molecular motors (Fig. 3a) (Rai et al. 2013). In the future, we wish to further investigate the origin of large velocity values observed for retrograde transport.

To expand the application of EVA in nanoscale biological systems, it is crucial to consider applying EV A to systems entirely different from motor proteins. Inspired by the poster presentation by Takaramoto and Inoue at the IUPAB2024 Congress (Takaramoto and Inoue 2024), we became interested in the analysis of droplet sizes formed by liquid-liquid  phase  separation  (LLPS).  Within  cells, there  are  membrane-less  organelles  formed through LLPS, such as stress granules and nucleoli. These organelles are formed through local concentration and  separation  caused  by  interactions  between  specific  proteins  and RNA,  playing  critical  roles  in  cellular functions. While membrane-less organelles formed by LLPS are reversibly assembled and disassembled under normal conditions, in pathological states, this process can break down, leading to the formation of irreversible aggregates.  In  neurodegenerative  diseases  such  as  ALS  and  Parkinson's  disease,  abnormal  irreversible aggregation of proteins has been reported to cause neuronal cell death. Takaramoto and Inoue investigated the mechanism of liquid condensate formation by the Parkinson's disease-related protein α -synuclein ( α Syn). A previous study (Hoffmann et al. 2021) has shown that when α Syn coexists with the neuronal protein synapsin, α Syn  is  incorporated  into  the  droplets  formed  by  synapsin  and  can  undergo  phase  separation  at  lower concentrations than when α Syn forms droplets alone. Takaramoto and Inoue studied the detailed molecular mechanism of droplet formation by α Syn in the presence of synapsin. Takaramoto and Inoue confirmed droplet formation in solution of α Syn and synapsin ( in vitro experiment) (Fig. 1 of Reference (Takaramoto and Inoue 2024)), and obtained the droplet size distribution. From our preliminary application of EV A to the droplet size distribution, we found that this distribution follows the FrŽchet type, which does not have a finite maximum value (unpublished data). On the other hand, within neurons, a maximum value for droplet size is expected to exist due to the fixed cell size. It is important to compare both the droplet formation in the in vitro experiment and cellular environments. In addition to experiments, molecular simulations have emerged as powerful tools to elucidate  the  physicochemical  principles  of  biomolecular  condensates,  complementing  experimental approaches. At Joseph's keynote session at IUPAB2024, we learned the potential of LLPS-specific simulations (Mpipi)  (Joseph  et  al.  2021).  Pursuing  whether  the  maximum  value  of  LLPS  droplet  size  is  related  to  the mechanism of disease onset and describing life sciences through extreme values rather than averages is a crucial task for the future.

## Declaration

No funding was received for conducting this study.

## Conflict of interest

The authors declare no competing interests.

## Data availability

Data sharing is not applicable to this article as no new data were created.

## Reference

Basnayake K, Mazaud D, Bemelmans A, Rouach N, Korkotian E, Holcman D (2019) Fast calcium transients in dendritic spines driven by extreme statistics. PLoS Biol , 17(6): e2006202. https://doi.org/10.1371/journal.pbio.2006202

Brenner S, Berger F, Rao L, Nicholas MP, Gennerich A (2020) Force production of human cytoplasmic dynein is limited by its processivity. Sci Adv , 6 (15)  eaaz4295. : https://www.science.org/doi/10.1126/sciadv.aaz4295

Chaubet L, Chaudhary AR, Heris HK, Ehrlicher AJ, Hendricks AG (2020) Dynamic actin cross-linking governs the cytoplasm's transition to fluid-like behavior. Mol Biol Cell , 31 (16)  1744-1752. : https://doi.org/10.1091/mbc.E19-09-0504

Coles S (2001) An Introduction to Statistical Modeling of Extreme Values. Springer, London .

de Haan L, Ferreira A (2006) Extreme Value Theory. Springer .

Dong X, Milholland B, Vijg J (2016) Evidence for a limit to human lifespan. Nature , 538 (7624)  257-259. : https://www.nature.com/articles/nature19793

Einmahl JH, Magnus JR (2008) Records in Athletics Through Extreme-Value Theory. Journal of the American Statistical Association , 103: pp.1382-1391. https://doi.org/10.1198/016214508000000698

Elshenawy MM, Canty JT, Oster L, Ferro LS, Zhou Z, Blanchard SC, Yildiz A (2019) Cargo adaptors regulate stepping and force generation of mammalian dynein-dynactin. Nat Chem Biol , 15 (11)  1093-1101. : https://www.nature.com/articles/s41589-019-0352-0

Gembris D, Taylor JG, Suter D (2002) Trends and random fluctuations in athletics. Nature , 417 (6888)  506. : https://www.nature.com/articles/417506a

Gennerich A, Carter AP, Reck-Peterson SL, Vale RD (2007) Force-induced bidirectional stepping of cytoplasmic dynein. Cell , 131 (5)  952-965. https://doi.org/10.1016/j.cell.2007.10.016 :

Gilleland E, Katz RW (2016) extRemes 2.0: An Extreme Value Analysis Package. J Software , 72: 1 39. https://doi.org/10.18637/jss.v072.i08

Guedes-Dias P, Holzbaur ELF (2019) Axonal transport: Driving synaptic function. Science , 366 :199. https://www.science.org/doi/10.1126/science.aaw9997

Harada  R,  Takano  Y,  Baba  T,  Shigeta  Y  (2015)  Simple,  yet  powerful  methodologies  for  comformational sampling of proteins. Phys Chem Chem Phys, 17:6155-6173. https://pubs.rsc.org/en/content/articlehtml/2015/cp/c4cp05262e

Hayashi  K,  Miyamoto  MG,  Niwa  S  (2021)  Effects  of  dynein  inhibitor  on  the  number  of  motor  proteins transporting synaptic cargos. Biophys J , 120 :1605. https://doi.org/10.1016/j.bpj.2021.02.018

Hirakawa E, Higuchi H, Toyoshima YY (2000) Processive movement of single 22S dynein molecules occurs only at low ATP concentrations. Proc Natl Acad Sci U S A , 97 (6)  2533-2537. : https://www.pnas.org/doi/full/10.1073/pnas.050585297

| Hoffmann C, Sansevrino R, Morabito G, Logan C, Martin Vabulas R, Ulusoy A, Ganzella M, Milovanovic D (2021) J Mol Biol , 433 :166961. https://doi.org/10.1016/j.jmb.2021.166961                                                                                                         |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Joseph JA, Reinhardt A, Aguirre A, Chew PY, Russell KO, Espinosa JR, Garaizar A, Collepardo-Guevara R (2021) Physics-driven coarse-grained model for biomolecular phase separation with near-quantitative accurac y. Nat Comput Sci , 1:732. https://doi.org/10.1038/s43588-021-00155-3 |
| Keefe AJ, Gabrych DR, Zhu Y, Vocadlo DJ, Silverman MA (2023) Axonal Transport of Lysosomes Is Unaffected in Glucocerebrosidase-Inhibited iPSC-Derived Forebrain Neurons. eNeuro https://doi.org/10.1523/ENEURO.0079-23.2023                                                             |
| Kratz M (2017) Introduction to Extreme Value Theory. Applications Risk Analysis & Manegament. Matrix annals pp.591-636.                                                                                                                                                                 |
| Mizuno D, Tardin C, Schmidt CF, Mackintosh FC (2007) Nonequilibrium mechanics of active cytoskeletal networks. Science 2007, 315 (5810) : 370-373. https://www.science.org/doi/10.1126/science.1134404                                                                                  |
| Mallik R, Carter BC, Lex SA, King SJ, Gross SP (2004) Cytoplasmic dynein functions as a gear in response to load. Nature , 427 (6975) : 649-652. https://www.nature.com/articles/nature02293                                                                                            |
| Naoi T, Kagawa Y, Nagino K, Niwa S, Hayashi K (2024) Extreme-Value Analysis of Intracellular Cargo Transport by Motor Proteins. Commun Phys , 7 :50. https://doi.org/10.1038/s42005-024-01538-4                                                                                         |
| R Core Team (2018) R: Alanguage and environment for statistical computinh. R Foundation for Statistical Computing, Vienna, Austria. URLhttps://www.R-project.org/.                                                                                                                      |
| Rai AK, Rai A, Ramaiya AJ, Jha R, Mallik R (2013) Molecular adaptations allow dynein to generate large collective forces inside cells. Cell , 152 (1-2) : 172-182. https://doi.org/10.1016/j.cell.2012.11.044                                                                           |
| Rootzen H, Zolud D. Human life is unlimited. but short (2017) Extremes , 20: pp. 713-728. https://link.springer.com/article/10.1007/s10687-017-0305-5                                                                                                                                   |
| Sasaki K, Kaya M, Higuchi H(2018) AUnified Walking Model for Dimeric Motor Proteins. Biophys J , 115 (10) : 1981-1992. https://doi.org/10.1016/j.bpj.2018.09.032                                                                                                                        |
| Schnitzer MJ, Visscher K, Block SM (2000) Force production by single kinesin motors. Nat Cell Biol , 2 (10) : 718-723. https://www.nature.com/articles/ncb1000_718                                                                                                                      |
| Songchitruksa P, Tarko AP (2006) The extreme value theory approach to safety estimation. Accid Anal Prev , 38 (4) : 811-822. https://doi.org/10.1016/j.aap.2006.02.003                                                                                                                  |
| Takaramoto S, Inoue K(2024) Interaction Mechanism of α -Synuclein with Synapsin in the Liquid Condensates. the collection of abstracts of the IUPAB2024 congress, 25P-067, Biophys and Physicobiol . https://doi.org/10.2142/biophysico.bppb-v21.e2001                                  |
| Tippett MK, Lepore C, Cohen JE (2016) More tornadoes in the most extreme U.S. tornado outbreaks. Science , 354 (6318) : 1419-1423. https://www.science.org/doi/10.1126/science.aah7393                                                                                                  |
| Toba S, Watanabe TM, Yamaguchi-Okimoto L, Toyoshima YY, Higuchi H(2006) Overlapping hand-over-hand mechanism of single molecular motility of cytoplasmic dynein. Proc Natl Acad Sci USA , 103 (15) : 5741-5745. https://www.pnas.org/doi/full/10.1073/pnas.0508511103                   |
| Tsuduki S (2024) Extreme value statistics of nerve transmission delay. PLoS One . https://doi.org/10.1371/journal.pone.0306605                                                                                                                                                          |
| Wong F, Collins JJ (2020) Evidence that coronavirus superspreading is fat-tailed. Proc Natl Acad Sci U S A ,                                                                                                                                                                            |

Wong F, Collins JJ (2020) Evidence that coronavirus superspreading is fat-tailed. Proc Natl Acad Sci U S A , 117 (47)  29416-29418. https://www.pnas.org/doi/full/10.1073/pnas.2018490117 :

Fig. 1 Simulation using random numbers. The cases for uniform random numbers ranging from 0 to 1 ( 𝑀 = 10 , 𝑛 = 500 )  ( a )  and  Gaussian  random  numbers  ( 𝑀 = 1000 , 𝑛 = 500 )  ( b ).  Each  panel  represents  the  counts, cumulative  probability  and  return  level  plot  for  the  block  maximum.  The  shape  parameter 𝜉	 &lt; 0 for  the uniform random numbers, and 𝜉	~0 for  the  Gaussian  random numbers. Red lines represent 𝐺(𝑋#$%) (Eq. (1)). ( c ) Examples of extreme value distributions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01007v2-figure-593f2adf-5a67-43f1-963e-7b4e85b1c295.jpg)

fluorescence protein

Vmax

Fig.  2  Schematics  of  synaptic  cargo  transport  by  motor  proteins  in  the  axon  of  a  neuron. a A  cargo  is anterogradely transported by KIF1A (kinesin) and retrogradely by cytoplasmic dynein. b , c Return level plot and histogram of the block maximum velocity data of synaptic cargo transport by motor proteins. Results for anterograde (red) and retrograde (blue) transport, in the cases of the motor neurons of C. elegans worms (Naoi et al. 2024) ( b ), and the hippocampal neurons of mice (Hayashi et al. 2021) ( c ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01007v2-figure-f5f25bcd-57c5-4008-b086-3afae87c15f3.jpg)

Fig. 3 Physical parameters of intracellular cargo transport by motor proteins. a The variety in the values of cargo transport resulted from the number of motor proteins and cargo sizes. b Load dependence of velocity in convex and concave cases. An existence of a cargo can be a load in the axons of neurons, because of high viscosity in the cells. c Simulation results for the model (Eq. (3)).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01007v2-figure-c90688c2-2c5c-4932-b774-3109d58cff24.jpg)
Kumiko Hayashi Nobumichi Takamatsu Shunki Takaramoto
Extreme-value analysis in nano-biological systems: Applications and Implications
Extreme value analysis (EVA) is a statistical method that studies the
properties of extreme values of datasets, crucial for fields like engineering,
meteorology, finance, insurance, and environmental science. EVA models extreme
events using distributions such as Fr\'echet, Weibull, or Gumbel, aiding in
risk prediction and management. This review explores EVA's application to
nanoscale biological systems. Traditionally, biological research focuses on
average values from repeated experiments. However, EVA offers insights into
molecular mechanisms by examining extreme data points. We introduce EVA's
concepts with simulations and review its use in studying motor protein
movements within cells, highlighting the importance of in vivo analysis due to
the complex intracellular environment. We suggest EVA as a tool for extracting
motor proteins' physical properties in vivo and discuss its potential in other
biological systems. While there have been only a few applications of EVA to
biological systems, it holds promise for uncovering hidden properties in
extreme data, promoting its broader application in life sciences.