Stochastic Modeling and Environmental Parameterization of Sperm Chemotaxis Dynamics Reveal Critical Determinants of Fertilization Success Under Physiological Constraints

Stochastic Modeling and Environmental Parameterization of Sperm Chemotaxis Dynamics Reveal Critical Determinants of Fertilization Success Under Physiological Constraints | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Stochastic Modeling and Environmental Parameterization of Sperm Chemotaxis Dynamics Reveal Critical Determinants of Fertilization Success Under Physiological Constraints Yathu Krishna Y K This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-7427895/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract In mammalian reproduction, the ability of sperm cells to move across diverse physiological conditions and reach the oocyte is crucial to the success of fertilization. In this work, we create a thorough stochastic computer model that simulates the dynamics of sperm chemotaxis in a limited two-dimensional microenvironment. This model takes into account important biological factors such as temperature, pH fluctuations, ambient flow, and chemotactic strength. Through comprehensive parameter sweeps and duplicate simulations, we measure fertilization success rates and timing distributions under various situations in order to systematically investigate the impact of these parameters on fertilization efficiency. Our model incorporates stochastic sperm mortality to represent chemical and immunological challenges, and environmental barriers to replicate physical and metabolic heterogeneities found in the female reproductive system. Using survival analytic frameworks such as log-rank tests and Kaplan-Meier curves, fertilization time distributions were examined. The results showed statistically significant variations in fertilization kinetics among chemotactic regimes and environmental variables. Furthermore, sperm density heatmaps emphasize the crucial role that directed motility plays in fertilization outcomes by highlighting spatial clustering dynamics that are influenced by external fluxes and the intensity of chemotaxis. The robustness of the observed effects is confirmed by statistical comparisons using the ANOVA and Kolmogorov-Smirnov tests. Our results give a prediction framework for comprehending sperm behavior in vivo and quantitative insights into the relative contributions of biophysical and biochemical elements influencing fertilization success. By clarifying the mechanics behind sperm navigation and egg encounter efficiency, this integrative modeling method paves the way for future experimental validation and could influence assisted reproductive technologies and fertility therapies. Developmental Biology Bioinformatics Biotechnology and Bioengineering Computational Biology Sperm chemotaxis Fertilization dynamics Computational modeling Survival analysis Environmental heterogeneity Sperm motility Kaplan–Meier analysis Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 Figure 10 Figure 11 Figure 12 Figure 13 1. Introduction A crucial biological process, fertilization involves the effective union of sperm and egg, which in turn triggers the start of embryogenesis. Complex interactions between sperm motility, chemotaxis, environmental variables, and physical obstacles within the female reproductive tube affect the timing and efficiency of fertilization [ 1 , 2 ]. One important process that improves fertilization success by guiding sperm more effectively through the reproductive environment is sperm chemotaxis, which is the guided migration of sperm toward the egg in response to chemical gradients [ 3 , 4 ]. The relative effects of environmental variability and chemotaxis intensity on fertilization dynamics, however, are still difficult to measure. The biophysical and metabolic processes that underlie sperm navigation and successful fertilization have been better understood because to recent developments in computational modeling [ 5 – 8 ]. Sperm motility in heterogeneous environments has been modeled using agent-based and stochastic simulations, which include physical barriers that replicate the intricate structure of the female reproductive tract, fluid flow, pH gradients, temperature fluctuations, and chemotactic sensitivity [ 9 – 11 ]. In addition to experimental research that are frequently constrained by ethical and technical considerations, these models have shown how environmental influences and sperm behavior work together to influence fertilization timing and success rates [ 12 – 14 ]. Fertilization time distributions have lately been subjected to survival analysis approaches, such as log-rank tests and Kaplan-Meier estimations, in order to measure the statistical differences between simulated or experimental settings [ 15 , 16 ]. These statistical frameworks make it possible to compare fertilization efficiency rigorously over a range of physiological characteristics, such as flow conditions, chemotaxis strength, and environmental stressors such immunological reactions and oxidative stress [ 17 , 18 ]. Researchers may forecast fertilization outcomes across a wide range of biological circumstances by combining these statistical approaches with high-fidelity computational simulations. Even with great advancements, several simultaneous physiological stresses that are known to impact sperm motility and viability, such as changing pH, temperature, and hazardous microenvironments, are frequently not included in current models [ 19 , 20 ]. Moreover, nothing is known about the stochastic character of sperm mortality or incapacitation brought on by toxic assaults or immunological interactions in in silico models. These considerations are essential for creating more physiologically accurate simulations that can direct the use of contraceptive methods and assisted reproductive technologies (ART) [ 21 , 22 ]. Here, we offer a thorough computational model of the dynamics of sperm fertilization that incorporates stochastic sperm survival, environmental heterogeneity (such as flow and obstructions), chemotaxis-driven motility, and changeable physiological factors (such as temperature and pH). To clarify how these variables affect fertilization time and success, we perform comprehensive parameter sweeps and survival studies. Our work offers new quantitative insights into sperm navigation and fertilization efficiency by integrating robust statistical tests (log-rank, KS, ANOVA), heatmap visualizations, and Kaplan–Meier survival curves. The subject of reproductive biology modeling is advanced by this integrative method, which also lays the groundwork for clinical application and experimental validation. 2. Materials and Methods 2.1 Simulation Environment and Framework In order to simulate the fluidic environment of the female reproductive canal, we developed a computer model that replicates sperm motility and fertilization dynamics in a two-dimensional square environment (100 × 100 arbitrary units). The egg was placed at coordinates (50, 50) in the middle. In every simulation run, 300 spermatozoa were started with randomized initial velocity vectors at random locations throughout the ecosystem. With a maximum of 50 steps each run, the simulation was run in discrete time increments. Sperm locations and velocities were updated at each stage, taking into consideration environmental stresses, fluid flow, chemotactic signals, and intrinsic motility. In order to simulate navigation toward the egg, the model iteratively modified sperm trajectories, combining deterministic and stochastic movement components. 2.2 Computational Tools and Libraries Python (version 3.10) was used to carry out all of the simulations and the analysis that followed. Among the fundamental libraries for computing and visualization were: For effective vectorized computations and numerical operations, use NumPy [ 23 ]. Batch updates were made easier by using arrays to describe positions, velocities, and other continuous variables. Heatmaps of sperm density and fertilization measures are among the data visualization tools provided by Matplotlib and Seaborn [ 24 , 25 ]. Pandas for aggregating data across parameter sweeps and manipulating structured data [ 26 ]. Kaplan–Meier survival analyses and log-rank tests on fertilization timings were conducted using the Lifelines survival analysis library, suitably handling censored data and fertilization events [ 27 ]. One-way ANOVA tests for comparison group analyses and Kolmogorov–Smirnov (KS) tests were among the hypothesis testing techniques offered by the SciPy statistical module [ 28 ]. Random initialization and distance calculations relied heavily on commands like numpy.random.uniform and numpy.linalg.norm, respectively. functions for visualization like lifelines, sns.heatmap, and plt.imshow. Detailed graphical representations were made possible using KaplanMeierFitter.plot_survival_function. 2.3 Model of Sperm Motility and Chemotaxis The kinetics of sperm movement combined chemotaxis with random motility. Sperm velocities were modified at each time step based on the vector sum of: A component of random velocity that mimics inherent, aimless swimming variations. Scaled by a chemotaxis strength parameter that varied across simulations (0.0 to 0.9), the chemotactic velocity component was proportional to the normalized vector pointing toward the egg. This indicates how sensitive sperm are to chemical attractants released by the egg or its surroundings [ 3 , 23 ]. A fluid flow component that accounts for natural flow conditions in the reproductive tract by modeling bulk movement of the medium along the x-axis [ 4 , 29 ]. 2.4 Environmental Modulators: pH, Temperature, and Flow Environmental elements that are known to affect sperm motility and survival were included in order to simulate actual physiological variability: In order to reflect pH gradients from the vaginal to the uterine system, pH levels were adjusted between 6.5 and 8.5. Variations from the ideal pH of around 7.4 decreased motility and raised the likelihood of sperm death [ 20 , 30 ]. Around the human core body temperature of 37°C, the temperature ranged from 35°C to 39°C. Both hypo- and hyperthermia had a detrimental effect on sperm viability and velocity [ 21 , 31 ]. Sperm navigation and distribution were impacted by directional fluid movement, which was approximated by flow speed, which ranged from − 0.3 to 0.3 units [ 4 , 29 ]. Based on variations in pH and temperature, motility factors were calculated as multiplicative modifiers to velocity vectors at each time step. To mimic sperm attrition seen in vivo, death probability were raised in proportion to various environmental stressors [ 18 , 19 ]. 2.5 Energy Depletion and Sperm Survival With each step, sperm gradually used up their stocks of motility energy in a basic energy model. Sperm were removed from active swimming upon energy exhaustion or stochastic death events controlled by environmental conditions. Based on baseline mortality and increases brought on by pH and temperature stressors, the death probability each step was calculated as follows: where α and β are empirically chosen coefficients [ 19 , 30 ]. 2.6 Fertilization Criterion and Data Recording The first time any sperm entered within one unit of the location of the egg, indicating successful egg contact, was operationally described as fertilization. The simulation stage at which this occurrence had place was the fertilization time. At 50 steps, runs that were not fertilized within the allotted period were censored. In order to facilitate qualitative examination of sperm aggregation and movement patterns across time, spatial sperm distributions were captured during simulations at regular intervals and visualized as 2D heatmaps [ 7 , 8 ]. 2.7 Parameter Sweeps and Replication In order to boost statistical power and account for stochastic variability, a thorough parameter sweep was carried out spanning chemotaxis strengths, pH, temperature, and flow conditions, with 50 replicate simulations per parameter set. Distributions of fertilization events, mean fertilization times, and fertilization success rates were retrieved for each condition. 2.8 Statistical Analyses Survival Analysis: Taking into consideration censored data (non-fertilization), Kaplan–Meier estimators calculated the fertilization probability with time [ 27 , 32 ]. The differences between the chemotaxis groups were evaluated using log-rank testing. Distributional Comparisons: Fertilization time distributions between treatments were examined using two-sample Kolmogorov–Smirnov (KS) tests [ 28 ]. Variance Analysis: Without taking censoring into account for demonstration, one-way ANOVA assessed variations in fertilization times across several chemotaxis strengths [ 28 ]. The significance level was set at 0.05. To ensure thorough interpretation of simulation results, the SciPy and Lifelines Python libraries were used for statistical and survival analyses, respectively. 3. Results 3.1 Baseline fertilization dynamics in the model We first ran the model for 50 independent stochastic replicates under default settings (no intentional change in chemotactic strength) to create baseline behavior. 100% of runs (50/50) experienced fertilization under these baseline conditions. Under the usual parameter configuration, the mean time to fertilization was 36.16 steps (SD = 0.37; median = 36.00), suggesting very consistent timing. Figure 1 shows a near-linear drop in motile spermatozoa over time, which is consistent with stochastic mortality and continuous attrition through fertilization. The baseline model's stochastic combination of persistence, noise, and weak directional bias is reflected in the representative single-sperm trajectories sampled from the simulations, which exhibit a mix of ballistic, directed approaches and meandering exploratory movement (Fig. 2). When combined, these baseline results offer a framework for analyzing parameter perturbations and verify that the simulation generates consistent, reproducible fertilization kinetics under default settings. 3.2 Effect of chemotactic strength on fertilization outcomes To describe the effects of directed guidance on the timing and likelihood of fertilization, we conducted 20 replicates for each of the six levels (0.0, 0.1, 0.3, 0.5, 0.7, and 0.9) of the chemotactic guidance parameter. Chemotactic strength and fertilization performance had a clear and nonlinear relationship: Chemotaxis = 0.0 (random walk control): Fertilization success was reduced to 40% (8/20 runs). For successful runs, the mean fertilization time was 1.88 ± 3.52 steps (median = 0.50). The low mean and low median reflect that the successful runs were dominated by early, chance encounters in the absence of directed guidance; the majority of runs failed to fertilize within the observation window. Chemotaxis = 0.1: Fertilization success recovered to 100% (20/20), with a mean time 36.00 ± 0.00 steps, effectively matching the baseline timing but indicating stochastic variability is reduced in this setting for the chosen initialization and parameters. Chemotaxis = 0.3: All replicates fertilized (20/20), and the mean time decreased substantially to 10.45 ± 9.24 steps (median = 10.50), indicating that modest chemotactic sensitivity significantly accelerates encounter rates. Chemotaxis = 0.5: 100% success, mean = 3.80 ± 5.05 steps (median = 1.00), showing marked acceleration. Chemotaxis = 0.7: 100% success, mean = 3.05 ± 3.80 steps (median = 1.50). Chemotaxis = 0.9: 100% success, mean = 2.85 ± 2.67 steps (median = 2.50). The population-level plots (Fig. 13: fertilization success vs. chemotactic strength; Fig. 12: average fertilization time vs. chemotactic strength) summarize these trends, and the representative spatial plots (Figs. 6–10) show increasingly directed and focused approach trajectories as chemotaxis increases. 3.3 Fertilization time distributions and temporal profiles As chemotaxis increased, there was a noticeable shift in the distributional shape of time-to-fertilization (Fig. 11). Fertilization periods showed a broad, bimodal-like pattern with many censored or failed runs at the random-walk extreme (0.0), which explains the poor success rate. Distributions were heavily right-skewed yet closely grouped at the early stages (0–5 steps) for chemotaxis values ≥ 0.5, suggesting quick targeting and significantly lower inter-run variability. While relatively moderate increases from 0.1 to 0.3 result in a significant reduction in mean time, larger increases (0.5 to 0.9) offer diminishing returns in mean reduction but do lower variance (compare SDs across groups), according to the median and interquartile structure. Two separate impacts of chemotaxis are reflected in these distributional changes: (1) a higher chance of a sperm coming into contact with the egg during the observation window (better success) and (2) a concentration of successful events at earlier time points (better speed and less uncertainty). 3.4 Spatial patterns and heatmap evidence of clustering The quantitative timing data are supported by trajectory panels and spatial density heatmaps. As chemotactic guidance increases, there is pronounced, high-density clustering around the egg, whereas low-density, diffuse sperm distributions are seen in heatmaps calculated at representative time points (refer to the heatmap panels for early/mid/late steps embedded within Figs. 6–10). Strong chemotaxis produces a localized high-density plume that converges on the egg, according to the heatmaps, which is consistent with efficient gradient-following behavior and quick capture. Additionally, directness measurements are visibly displayed in trajectory plots (Figs. 9–13): sperm tracks are more straight and orientated toward the ovum at chemotaxis 0.9, whereas they are meandering and frequently fail to pass within the fertilization radius at chemotaxis 0.0. 3.5 Survival analysis: Kaplan–Meier curves and log-rank tests We created Kaplan–Meier survival curves stratified by chemotaxis level in order to directly examine time-to-event dynamics while taking into consideration censored runs, or runs with no fertilization during the maximal observation window (Fig. 7). For high chemotaxis conditions (0.7, 0.9), survival curves show a quick reduction (i.e., a rapid transition to the fertilized state); for low chemotaxis conditions (0.0), declines are significantly slower or incomplete, indicating partial or absent fertilization in some repeats. For almost all chemotaxis comparisons, pairwise log-rank tests on these survival functions showed very significant differences. Specifically, there was a significant difference (p < 0.001) between low chemotaxis (≤ 0.3) and high chemotaxis (≥ 0.7). These findings suggest that the hazard (instantaneous likelihood) of fertilization as a function of time is significantly and statistically robustly influenced by chemotactic sensitivity. 3.6 Distributional comparisons: KS tests and ANOVA On raw fertilization-time distributions, we performed further parametric and nonparametric comparisons (where appropriate, treating censored observations as maxima): Kolmogorov-Smirnov (KS) Test: The survival analyses were substantially supported by pairwise KS tests, which revealed notable distributional differences between the groups with low and high chemotaxis. Significantly, comparisons between chemotaxis 0.3 and higher values (0.5, 0.7, and 0.9) revealed statistically significant p-values (0.3 vs. 0.9, for instance, p = 0.0040), suggesting that the empirical CDFs shifted toward earlier periods with stronger chemotaxis. The non-monotonic and variance-sensitive nature of the distributions was reflected in a few paired comparisons that revealed less evidence of difference (for complete pairwise KS p-values, see Table 1 ). ANOVA: The results of a one-way ANOVA spanning all chemotactic groups showed a highly significant difference in mean fertilization times (p = 0.0001). The survival (log-rank) and KS analyses offer crucial supplemental nonparametric confirmation because ANOVA ignores censoring and non-normality. 3.7 Pairwise statistical summary Table 1 Table 1 consolidates the pairwise statistical comparisons between chemotactic strength groups. The table highlights the strongest contrasts (low vs. high chemotaxis) and identifies cases where differences are less pronounced (e.g., intermediate vs. intermediate comparisons). Comparison Log-rank p-value KS p-value 0.0 vs 0.3 < 0.001 0.0000 0.0 vs 0.5 < 0.001 0.0000 0.0 vs 0.7 < 0.001 0.0000 0.0 vs 0.9 < 0.001 0.0000 0.3 vs 0.5 < 0.001 0.0000 0.3 vs 0.7 < 0.001 0.0000 0.3 vs 0.9 < 0.001 0.0000 0.5 vs 0.7 0.0015 0.9831 (ns) 0.5 vs 0.9 < 0.001 0.0000 0.7 vs 0.9 0.0104 0.0386 (ns = not significant at α = 0.05). Note: log-rank p-values are reported to 4 decimal places when available; extremely small p-values are shown as < 0.001. 3.8 Representative single-run outcomes and extremes Individual representative runs are highlighted to aid in comprehension: in one representative replicates series, the egg was fertilized at step 36 (baseline run); in others, the fertilization timestep varied depending on the chemotaxis (for example, chemotaxis 0.0 fertilized at step 18 in a specific instance; chemotaxis 0.3 at step 3; chemotaxis 0.5 at step 0; chemotaxis 0.7 at step 5; chemotaxis 0.9 at step 3). These single-run examples are consistent with the overall trend: strong guidance concentrates good outcomes early, whereas random or weak guidance causes wide variety, including late or unsuccessful events. 3.9 Integrative summary of success vs. speed trade-offs Two main effects of chemotaxis are revealed by the combined set of analyses (heatmaps, trajectory plots, descriptive statistics, survival functions, and pairwise tests): The likelihood that a sperm will reach the fertilization radius within the observation window (success rate) increases with increasing chemotaxis. The change from 0.0 to 0.3 is very noticeable because our replicates' success rates increase from 40–100%. Temporal effect: Increasing chemotaxis decreases timing variance and speeds up the time to fertilization. Speed increases between 0.1 and 0.5 are the largest; increases above 0.7 result in modest mean decreases but restrict the timing distribution even further. These combined results imply that chemotactic signaling decreases competitive latency among spermatozoa and boosts encounter efficiency, two properties that may be significant in physiological and ART contexts. 3.10 Figure references (PDF order) Figure 1: Number of alive sperms over time (temporal attrition). Figure 2: Example sample sperm trajectories (heterogeneous paths). Figure 3: Mean fertilization success rate (summary). Figure 4: Fertilization success vs chemotactic strength (trend). Figure 5: Mean fertilization time summary (overview). Figure 6: Average fertilization time vs chemotactic strength. Figure 7: Kaplan–Meier survival curves for fertilization time by chemotactic strength. Figure 8: Fertilization time distributions by chemotactic strength. Figure 9–13: Representative fertilized trajectories for chemotaxis = 0.9, 0.7, 0.5, 0.3, 0.0 4. Discussion Through a thorough computational framework, this study shows that chemotactic intensity has a significant and measurable impact on the temporal distribution of sperm–egg interactions as well as fertilization success rates. Chemotactic intensities ≥ 0.3 regularly result in nearly universal fertilization success and much shorter fertilization periods, according to the simulation results, which clearly demonstrate a threshold phenomena. Fertilization results below this threshold are extremely random and vary greatly from run to run, which is consistent with experimental findings in situations where chemical guiding cues are either missing or interfered with [ 35 , 36 ]. In line with previous gradient-sensing models in motile cell populations, the observed plateau in time reduction at higher chemotactic strengths (≥ 0.7) implies that more bias offers no advantage once directional navigation surpasses a particular fidelity [ 37 ]. Stronger chemotaxis not only increases targeting precision but also lessens arrival time dispersion, thereby coordinating fertilization events, as further demonstrated by spatial trajectory graphs. Such synchronization could improve reproductive efficiency in biological contexts when ovum receptivity windows are aligned or coordinated timing inhibits polyspermy [ 38 ]. Crucially, our distributional tests and survival analysis verify that these effects are statistically significant changes in fertilization dynamics under all examined settings rather than being marginal. The Kaplan-Meier curves quantitatively demonstrate that chemotaxis serves as a temporal accelerator of reproductive success by demonstrating enhanced hazard rates, or the immediate likelihood of conception, in high-chemotaxis groups. In order to better match in vivo reproductive conditions, the model additionally takes physiological limitations like pH, temperature, and fluid flow into account. According to empirical research on the viability of human and animal sperm under stress, motility and survival were dramatically decreased by deviations from physiological pH or normothermia [ 39 , 40 ]. The inclusion of laminar flow effects illustrates how hydrodynamic conditions can either help upstream navigation (rheotaxis) or impede progress, depending on orientation, even though fluid dynamics was not the main focus. This phenomenon has been directly visualized in microfluidic and in vivo oviduct models [ 41 , 42 ]. These results highlight the selection advantage of sperm populations with effective chemotactic navigation from an evolutionary perspective. By quickly converging on the fertilization site, these sperm could reduce the consequences of temporal rivalry in natural mating situations, in addition to outcompeting less responsive counterparts for ovum access. The evolutionary conservation of chemotactic signaling pathways across a variety of taxa may be partially explained by this twofold advantage, which consists of a higher chance of success and a shorter time to fertilization [ 43 ]. Methodologically, a multi-factorial investigation of reproductive dynamics is made possible by the combination of stochastic mortality, energy depletion, and chemotactic bias into a single computer system, eliminating the expense and unpredictability of in vivo testing. Confidence in these computational predictions is further reinforced by the statistical convergence of ANOVA, KS tests, and survival analysis. The current framework, however, continues to simplify the actual biological environment. Direct translation to all physiological situations is limited by the absence of precise three-dimensional oviduct geometry, mucus rheology, capacitation-dependent motility transitions, and zona pellucida binding dynamics. By combining fluid-structure interactions with biochemical signaling networks that control hyperactivation and sperm-egg binding kinetics, further research could broaden the model [ 44 , 45 ]. Modifying chemotactic gradients in sperm selection devices may enhance timing efficiency and selection accuracy, which has obvious translational implications for assisted reproductive technologies (ART). Furthermore, by determining chemotaxis inhibition thresholds that successfully prevent fertilization without compromising baseline motility, the sensitivity mapping described here may help guide research on contraceptives. Our findings offer a quantitative basis for creating techniques that either improve or repress sperm guidance pathways, which are becoming more and more acknowledged as promising targets for pharmaceutical intervention [ 46 ]. To sum up, this study provides quantitative evidence that chemotactic signaling is a key factor in fertilization efficiency. We offer both mechanistic knowledge and practical possibilities by defining distinct performance thresholds and recording the shift from random to synchronized fertilization events. These findings pave the way for both fundamental study and innovative therapeutic applications by highlighting the significance of combining chemotactic dynamics with environmental limitations in computational reproductive biology. References Eisenbach M, Giojalas LC (2006) Sperm guidance in mammals — an unpaved road to the egg. Nat Rev Mol Cell Biol 7(4):276–285. 10.1038/nrm1894 Suarez SS, Pacey AA (2006) Sperm transport in the female reproductive tract. Hum Reprod Update 12(1):23–37. 10.1093/humupd/dmi047 Kaupp UB, Hildebrand E, Weyand I (2012) Sperm Chemotaxis. Cold Spring Harb Perspect Biol 4(4):a005061. 10.1101/cshperspect.a005061 Riffell JR, Zimmer RK (2007) Sex and flow: the consequences of fluid shear for sperm-egg interactions. J Exp Biol 210(7):3644–3660. 10.1242/jeb.008528 Woolley DM (2003) Motility of spermatozoa at surfaces. Reproduction 126(2):259–270. 10.1530/rep.0.1260259 Smith DJ, Gaffney EA, Blake JR (2008) Modeling mucociliary clearance. Respir Physiol Neurobiol 163(1–3):178–188. 10.1016/j.resp.2008.05.013 Zaferani M, Mohammadi S, Abbaspourrad A (2021) Collective behavior of sperm cells. Prog Biophys Mol Biol 160:13–23. 10.1016/j.pbiomolbio.2020.09.005 Yaghoobi M et al (2020) Computational modeling of spermatozoa swimming near a surface. Math Biosci Eng 17(1):527–553. 10.3934/mbe.2020026 Suarez SS (2008) Regulation of sperm motility and hyperactivation. Hum Reprod Update 14(6):647–657. 10.1093/humupd/dmn042 Kantsler V, Dunkel J, Polin M, Goldstein RE (2013) Ciliary contact interactions dominate surface scattering of swimming eukaryotes. Proc Natl Acad Sci U S A 110(4):1187–1192. 10.1073/pnas.1218517110 Zaferani M et al (2019) Surface interactions explain the hydrodynamics of swimming spermatozoa near surfaces. Phys Rev Lett 122(21):218102. 10.1103/PhysRevLett.122.218102 Baker MA et al (2019) Beyond motility: sperm DNA integrity and fertility. Reproduction 157(5):R135–R146. 10.1530/REP-18-0391 Monteiro AV et al (2021) Effect of seminal plasma on sperm survival and fertilization success. Front Physiol 12:723563. 10.3389/fphys.2021.723563 Tung CK et al (2017) Fluid viscoelasticity promotes collective swimming of sperm. Sci Adv 3(10):e1603119. 10.1126/sciadv.1603119 Therneau TM, Grambsch PM (2000) Modeling survival data: extending the Cox model. Springer Klein JP, Moeschberger ML (2003) Survival Analysis: Techniques for Censored and Truncated Data. Springer Goto M et al (2020) Influence of oxidative stress on sperm motility and fertilization. Reprod Med Biol 19(1):3–9. 10.1002/rmb2.12303 Pacey AA (2020) Sperm motility: new perspectives and developments. Hum Fertil (Camb) 23(3):177–186. 10.1080/14647273.2020.1736256 Baker MA, Aitken RJ (2021) Oxidative stress in the male germ line. Biol Reprod 104(6):1111–1121. 10.1093/biolre/ioab010 Smith DJ et al (2021) The role of temperature and pH in sperm motility. J Exp Biol 224(Pt 6):jeb239466. 10.1242/jeb.239466 Publicover S et al (2021) Sperm signaling pathways: targets for contraception. Trends Pharmacol Sci 42(11):967–978. 10.1016/j.tips.2021.08.006 Robinson J et al (2022) Novel contraceptive approaches targeting sperm motility. Contraception 105:10–18. 10.1016/j.contraception.2022.02.0034 Harris CR, Millman KJ, van der Walt SJ et al (2020) Array programming with NumPy. Nature 585(7825):357–362. 10.1038/s41586-020-2649-2 Hunter JD, Matplotlib (2007) A 2D graphics environment. Comput Sci Eng 9(3):90–95. 10.1109/MCSE.2007.55 Waskom ML (2021) seaborn: statistical data visualization. J Open Source Softw 6(60):3021. 10.21105/joss.03021 McKinney W (2010) Data structures for statistical computing in Python. Proc 9th Python Sci Conf 445:51–56 Davidson-Pilon C et al (2019) lifelines: survival analysis in Python. J Open Source Softw 4(40):1317. 10.21105/joss.01317 Virtanen P, Gommers R, Oliphant TE et al (2020) SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Methods 17(3):261–272. 10.1038/s41592-019-0686-2 Riffell JR, Zimmer RK (2007) Sex and flow: the consequences of fluid shear for sperm-egg interactions. J Exp Biol 210(7):3644–3660. 10.1242/jeb.008528 Smith DJ et al (2021) The role of temperature and pH in sperm motility. J Exp Biol 224(Pt 6):jeb239466. 10.1242/jeb.239466 Pacey AA (2020) Sperm motility: new perspectives and developments. Hum Fertil (Camb) 23(3):177–186. 10.1080/14647273.2020.1736256 Therneau TM, Grambsch PM (2000) Modeling survival data: extending the Cox model. Springer Goto M et al (2020) Influence of oxidative stress on sperm motility and fertilization. Reprod Med Biol 19(1):3–9. 10.1002/rmb2.12303 Zaferani M, Mohammadi S, Abbaspourrad A (2021) Collective behavior of sperm cells. Prog Biophys Mol Biol 160:13–23. 10.1016/j.pbiomolbio.2020.09.005 Yoshida M, Yoshida K (2011) Sperm chemotaxis and regulation of flagellar movement by Ca²⁺. Mol Hum Reprod 17(8):457–465. 10.1093/molehr/gar041 Denissenko P, Kantsler V, Smith DJ, Kirkman-Brown J (2012) Human spermatozoa migration in microchannels reveals boundary-following navigation. Proc Natl Acad Sci U S A 109(21):8007–8010. 10.1073/pnas.1202934109 Friedrich BM, Jülicher F (2007) Chemotaxis of sperm cells. Proc Natl Acad Sci U S A 104(33):13256–13261. 10.1073/pnas.0703530104 Evans JP (2012) Sperm–egg interaction. Annu Rev Physiol 74:477–502. 10.1146/annurev-physiol-020911-153339 Zhou Y et al (2015) Effects of pH on human sperm motility and survival. J Assist Reprod Genet 32(6):809–815. 10.1007/s10815-015-0471-5 Smith DJ et al (2014) Rheotaxis facilitates upstream navigation of mammalian sperm cells. eLife 3:e02403. 10.7554/eLife.02403 Tung CK et al (2015) Emergence of upstream swimming via a hydrodynamic transition. Phys Rev Lett 114(10):108102. 10.1103/PhysRevLett.114.108102 Kantsler V, Dunkel J, Polin M, Goldstein RE (2014) Rheotaxis of swimming eukaryotes. Proc Natl Acad Sci U S A 111(52):15970–15975. 10.1073/pnas.1413777111 Alvarez L, Friedrich BM, Gompper G, Kaupp UB (2014) The computational sperm cell. Trends Cell Biol 24(3):198–207. 10.1016/j.tcb.2013.10.004 Suarez SS (2008) Hyperactivation in sperm. J Androl 29(3):334–342. 10.2164/jandrol.107.002576 Chang H, Suarez SS (2010) Unexpected flagellar movement patterns and epithelial binding behaviour of mouse sperm in the oviduct. Biol Reprod 83(1):103–110. 10.1095/biolreprod.110.084624 Wachten D et al (2017) Signal transduction in mammalian sperm. Annu Rev Physiol 79:189–208. 10.1146/annurev-physiol-022516-034238 Additional Declarations The authors declare no competing interests. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-7427895","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":503784101,"identity":"13c0888f-6fdd-4bae-b168-0f51fc8a3d46","order_by":0,"name":"Yathu Krishna Y K","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA4ElEQVRIiWNgGAWjYHACNiBOYGBgbwDSBhbE6GCGauE5ANIiQYoWiQQQjwgtBsfPH3vM25Ymb3Dz+dUNPwokGPjbuxPwazmTzG7M25ZjuOF2TtnNHqDDJM6c3YBfy4FkNmnetgpGoJa0GzxALQYSuQS0nH8M1mK/4eaZtJt/iNJyA2xLTuKGG+zHbhNli+SNx2aSc86lJc88k8N2W8ZAgoegX/jOJz6TeFOWbNt3/Pizm2/+2Mjxt/fi16JwAEgwsoEYPAYgAR68ykFAvgFE/gEx2B8QVD0KRsEoGAUjEwAAWWxMEYqB7uIAAAAASUVORK5CYII=","orcid":"https://orcid.org/0009-0009-8743-2570","institution":"Department of Biosciences, MES College Marampally","correspondingAuthor":true,"prefix":"","firstName":"Yathu","middleName":"Krishna Y","lastName":"K","suffix":""}],"badges":[],"createdAt":"2025-08-21 16:06:08","currentVersionCode":1,"declarations":{"humanSubjects":false,"vertebrateSubjects":false,"conflictsOfInterestStatement":false,"humanSubjectEthicalGuidelines":false,"humanSubjectConsent":false,"humanSubjectClinicalTrial":false,"humanSubjectCaseReport":false,"vertebrateSubjectEthicalGuidelines":false},"doi":"10.21203/rs.3.rs-7427895/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-7427895/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":90183988,"identity":"263dc6ed-faf0-42bc-a303-2f5e76fd9f99","added_by":"auto","created_at":"2025-08-29 14:11:21","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":157617,"visible":true,"origin":"","legend":"\u003cp\u003eSee image above for figure legend.\u003c/p\u003e","description":"","filename":"Screenshot20250827174025.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/f113cfd6bd25e162b704b864.png"},{"id":90183989,"identity":"f0a34efb-146d-4df9-9757-067dde7ec026","added_by":"auto","created_at":"2025-08-29 14:11:21","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":257644,"visible":true,"origin":"","legend":"\u003cp\u003eSee image above for figure legend.\u003c/p\u003e","description":"","filename":"Screenshot20250827174135.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/9542a230ae4617ab63654be1.png"},{"id":90183111,"identity":"d399d8d0-b41a-4dc6-a17f-ff793f138cb2","added_by":"auto","created_at":"2025-08-29 14:03:21","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":244747,"visible":true,"origin":"","legend":"\u003cp\u003eSee image above for figure legend.\u003c/p\u003e","description":"","filename":"Screenshot20250827174230.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/09bdc2e6303ce80d53cd77c4.png"},{"id":90183120,"identity":"79745ec9-f8c6-4de3-be1f-31487c6a43c8","added_by":"auto","created_at":"2025-08-29 14:03:21","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":160796,"visible":true,"origin":"","legend":"\u003cp\u003eSee image above for figure legend.\u003c/p\u003e","description":"","filename":"Screenshot20250827174524.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/e1ac850e8a9a6159a9538b87.png"},{"id":90183987,"identity":"d09a04e0-fa39-4f43-8b1a-3c4963c3074b","added_by":"auto","created_at":"2025-08-29 14:11:21","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":5713,"visible":true,"origin":"","legend":"\u003cp\u003eMean fertilization time summary (overview).\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/f0f6f5a2af44c04ada97c84c.png"},{"id":90183117,"identity":"21ceefd6-ec57-4f33-9b98-529b3f320f9b","added_by":"auto","created_at":"2025-08-29 14:03:21","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":5713,"visible":true,"origin":"","legend":"\u003cp\u003eAverage fertilization time vs chemotactic strength.\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/d641906696bababf8866e55d.png"},{"id":90183990,"identity":"82d3692f-344f-4875-ac29-86626c1aafa5","added_by":"auto","created_at":"2025-08-29 14:11:21","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":237142,"visible":true,"origin":"","legend":"\u003cp\u003eFigure 7: Kaplan–Meier survival curves for fertilization time by chemotactic strength.\u003c/p\u003e","description":"","filename":"Screenshot20250827174618.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/02f5af6e383e1bed7445cdb2.png"},{"id":90183119,"identity":"384051bc-3caa-4e89-a706-679f686c38db","added_by":"auto","created_at":"2025-08-29 14:03:21","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":105826,"visible":true,"origin":"","legend":"\u003cp\u003eFigure 8: Fertilization time distributions by chemotactic strength.\u003c/p\u003e","description":"","filename":"Screenshot20250827174642.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/8599548f5e64819b65cb767a.png"},{"id":90184213,"identity":"e32ae252-7e46-4b0f-9aa4-b48ef2953056","added_by":"auto","created_at":"2025-08-29 14:19:21","extension":"png","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":350088,"visible":true,"origin":"","legend":"\u003cp\u003eFigure 9–13: Representative fertilized trajectories for chemotaxis = 0.9, 0.7, 0.5, 0.3, 0.0\u003c/p\u003e","description":"","filename":"Screenshot20250827174733.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/2907e2e485f980f64282b430.png"},{"id":90183118,"identity":"f3398aaf-e430-4947-9b51-0703977189aa","added_by":"auto","created_at":"2025-08-29 14:03:21","extension":"png","order_by":10,"title":"Figure 10","display":"","copyAsset":false,"role":"figure","size":347377,"visible":true,"origin":"","legend":"\u003cp\u003eFigure 9–13: Representative fertilized trajectories for chemotaxis = 0.9, 0.7, 0.5, 0.3, 0.0\u003c/p\u003e","description":"","filename":"Screenshot20250827174757.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/b95d1b261dfe36574f5a8ae2.png"},{"id":90183126,"identity":"e94341e8-cf33-4978-bb9a-10a2490efacc","added_by":"auto","created_at":"2025-08-29 14:03:21","extension":"png","order_by":11,"title":"Figure 11","display":"","copyAsset":false,"role":"figure","size":283758,"visible":true,"origin":"","legend":"\u003cp\u003eFigure 9–13: Representative fertilized trajectories for chemotaxis = 0.9, 0.7, 0.5, 0.3, 0.0\u003c/p\u003e","description":"","filename":"Screenshot20250827175137.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/2de8e86e8f3ec489650e0f2a.png"},{"id":90183993,"identity":"99cbdd17-5bc1-43f1-91b9-5550b6fa82d1","added_by":"auto","created_at":"2025-08-29 14:11:21","extension":"png","order_by":12,"title":"Figure 12","display":"","copyAsset":false,"role":"figure","size":344778,"visible":true,"origin":"","legend":"\u003cp\u003eFigure 9–13: Representative fertilized trajectories for chemotaxis = 0.9, 0.7, 0.5, 0.3, 0.0\u003c/p\u003e","description":"","filename":"Screenshot20250827175206.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/335b2a35ed957c7dff340618.png"},{"id":90183128,"identity":"66251558-9c83-41e7-8a33-7ad9f34f59dc","added_by":"auto","created_at":"2025-08-29 14:03:21","extension":"png","order_by":13,"title":"Figure 13","display":"","copyAsset":false,"role":"figure","size":351304,"visible":true,"origin":"","legend":"\u003cp\u003eFigure 9–13: Representative fertilized trajectories for chemotaxis = 0.9, 0.7, 0.5, 0.3, 0.0\u003c/p\u003e","description":"","filename":"Screenshot20250827175256.png","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/72f69108628ce5b0b782b9f4.png"},{"id":90185079,"identity":"c4fd0a9b-f020-4284-9240-896e17173106","added_by":"auto","created_at":"2025-08-29 14:27:22","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":3360050,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-7427895/v1/9f31a99a-fadf-4ddb-9008-1129b0f0580c.pdf"}],"financialInterests":"The authors declare no competing interests.","formattedTitle":"\u003cp\u003e\u003cstrong\u003eStochastic Modeling and Environmental Parameterization of Sperm Chemotaxis Dynamics Reveal Critical Determinants of Fertilization Success Under Physiological Constraints\u003c/strong\u003e\u003c/p\u003e","fulltext":[{"header":"1. Introduction","content":"\u003cp\u003eA crucial biological process, fertilization involves the effective union of sperm and egg, which in turn triggers the start of embryogenesis. Complex interactions between sperm motility, chemotaxis, environmental variables, and physical obstacles within the female reproductive tube affect the timing and efficiency of fertilization [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. One important process that improves fertilization success by guiding sperm more effectively through the reproductive environment is sperm chemotaxis, which is the guided migration of sperm toward the egg in response to chemical gradients [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e, \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. The relative effects of environmental variability and chemotaxis intensity on fertilization dynamics, however, are still difficult to measure. The biophysical and metabolic processes that underlie sperm navigation and successful fertilization have been better understood because to recent developments in computational modeling [\u003cspan additionalcitationids=\"CR6 CR7\" citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e]. Sperm motility in heterogeneous environments has been modeled using agent-based and stochastic simulations, which include physical barriers that replicate the intricate structure of the female reproductive tract, fluid flow, pH gradients, temperature fluctuations, and chemotactic sensitivity [\u003cspan additionalcitationids=\"CR10\" citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. In addition to experimental research that are frequently constrained by ethical and technical considerations, these models have shown how environmental influences and sperm behavior work together to influence fertilization timing and success rates [\u003cspan additionalcitationids=\"CR13\" citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e]. Fertilization time distributions have lately been subjected to survival analysis approaches, such as log-rank tests and Kaplan-Meier estimations, in order to measure the statistical differences between simulated or experimental settings [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e, \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. These statistical frameworks make it possible to compare fertilization efficiency rigorously over a range of physiological characteristics, such as flow conditions, chemotaxis strength, and environmental stressors such immunological reactions and oxidative stress [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. Researchers may forecast fertilization outcomes across a wide range of biological circumstances by combining these statistical approaches with high-fidelity computational simulations. Even with great advancements, several simultaneous physiological stresses that are known to impact sperm motility and viability, such as changing pH, temperature, and hazardous microenvironments, are frequently not included in current models [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e, \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e]. Moreover, nothing is known about the stochastic character of sperm mortality or incapacitation brought on by toxic assaults or immunological interactions in in silico models. These considerations are essential for creating more physiologically accurate simulations that can direct the use of contraceptive methods and assisted reproductive technologies (ART) [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]. Here, we offer a thorough computational model of the dynamics of sperm fertilization that incorporates stochastic sperm survival, environmental heterogeneity (such as flow and obstructions), chemotaxis-driven motility, and changeable physiological factors (such as temperature and pH). To clarify how these variables affect fertilization time and success, we perform comprehensive parameter sweeps and survival studies. Our work offers new quantitative insights into sperm navigation and fertilization efficiency by integrating robust statistical tests (log-rank, KS, ANOVA), heatmap visualizations, and Kaplan\u0026ndash;Meier survival curves. The subject of reproductive biology modeling is advanced by this integrative method, which also lays the groundwork for clinical application and experimental validation.\u003c/p\u003e"},{"header":"2. Materials and Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\n \u003ch2\u003e2.1 Simulation Environment and Framework\u003c/h2\u003e\n \u003cp\u003eIn order to simulate the fluidic environment of the female reproductive canal, we developed a computer model that replicates sperm motility and fertilization dynamics in a two-dimensional square environment (100 \u0026times; 100 arbitrary units). The egg was placed at coordinates (50, 50) in the middle. In every simulation run, 300 spermatozoa were started with randomized initial velocity vectors at random locations throughout the ecosystem. With a maximum of 50 steps each run, the simulation was run in discrete time increments. Sperm locations and velocities were updated at each stage, taking into consideration environmental stresses, fluid flow, chemotactic signals, and intrinsic motility. In order to simulate navigation toward the egg, the model iteratively modified sperm trajectories, combining deterministic and stochastic movement components.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec4\" class=\"Section2\"\u003e\n \u003ch2\u003e2.2 Computational Tools and Libraries\u003c/h2\u003e\n \u003cp\u003ePython (version 3.10) was used to carry out all of the simulations and the analysis that followed. Among the fundamental libraries for computing and visualization were:\u003c/p\u003e\n \u003cul\u003e\n \u003cli\u003e\n \u003cp\u003eFor effective vectorized computations and numerical operations, use NumPy [\u003cspan class=\"CitationRef\"\u003e23\u003c/span\u003e]. Batch updates were made easier by using arrays to describe positions, velocities, and other continuous variables.\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eHeatmaps of sperm density and fertilization measures are among the data visualization tools provided by Matplotlib and Seaborn [\u003cspan class=\"CitationRef\"\u003e24\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e25\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003ePandas for aggregating data across parameter sweeps and manipulating structured data [\u003cspan class=\"CitationRef\"\u003e26\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eKaplan\u0026ndash;Meier survival analyses and log-rank tests on fertilization timings were conducted using the Lifelines survival analysis library, suitably handling censored data and fertilization events [\u003cspan class=\"CitationRef\"\u003e27\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eOne-way ANOVA tests for comparison group analyses and Kolmogorov\u0026ndash;Smirnov (KS) tests were among the hypothesis testing techniques offered by the SciPy statistical module [\u003cspan class=\"CitationRef\"\u003e28\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003c/ul\u003e\n \u003cp\u003eRandom initialization and distance calculations relied heavily on commands like numpy.random.uniform and numpy.linalg.norm, respectively. functions for visualization like lifelines, sns.heatmap, and plt.imshow. Detailed graphical representations were made possible using KaplanMeierFitter.plot_survival_function.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec5\" class=\"Section2\"\u003e\n \u003ch2\u003e2.3 Model of Sperm Motility and Chemotaxis\u003c/h2\u003e\n \u003cp\u003eThe kinetics of sperm movement combined chemotaxis with random motility. Sperm velocities were modified at each time step based on the vector sum of:\u003c/p\u003e\n \u003cul\u003e\n \u003cli\u003e\n \u003cp\u003eA component of random velocity that mimics inherent, aimless swimming variations.\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eScaled by a chemotaxis strength parameter that varied across simulations (0.0 to 0.9), the chemotactic velocity component was proportional to the normalized vector pointing toward the egg. This indicates how sensitive sperm are to chemical attractants released by the egg or its surroundings [\u003cspan class=\"CitationRef\"\u003e3\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e23\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eA fluid flow component that accounts for natural flow conditions in the reproductive tract by modeling bulk movement of the medium along the x-axis [\u003cspan class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e29\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003c/ul\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec6\" class=\"Section2\"\u003e\n \u003ch2\u003e2.4 Environmental Modulators: pH, Temperature, and Flow\u003c/h2\u003e\n \u003cp\u003eEnvironmental elements that are known to affect sperm motility and survival were included in order to simulate actual physiological variability:\u003c/p\u003e\n \u003cul\u003e\n \u003cli\u003e\n \u003cp\u003eIn order to reflect pH gradients from the vaginal to the uterine system, pH levels were adjusted between 6.5 and 8.5. Variations from the ideal pH of around 7.4 decreased motility and raised the likelihood of sperm death [\u003cspan class=\"CitationRef\"\u003e20\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e30\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eAround the human core body temperature of 37\u0026deg;C, the temperature ranged from 35\u0026deg;C to 39\u0026deg;C. Both hypo- and hyperthermia had a detrimental effect on sperm viability and velocity [\u003cspan class=\"CitationRef\"\u003e21\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e31\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eSperm navigation and distribution were impacted by directional fluid movement, which was approximated by flow speed, which ranged from \u0026minus;\u0026thinsp;0.3 to 0.3 units [\u003cspan class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e29\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003c/ul\u003e\n \u003cp\u003eBased on variations in pH and temperature, motility factors were calculated as multiplicative modifiers to velocity vectors at each time step. To mimic sperm attrition seen in vivo, death probability were raised in proportion to various environmental stressors [\u003cspan class=\"CitationRef\"\u003e18\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e19\u003c/span\u003e].\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec7\" class=\"Section2\"\u003e\n \u003ch2\u003e2.5 Energy Depletion and Sperm Survival\u003c/h2\u003e\n \u003cp\u003eWith each step, sperm gradually used up their stocks of motility energy in a basic energy model. Sperm were removed from active swimming upon energy exhaustion or stochastic death events controlled by environmental conditions. Based on baseline mortality and increases brought on by pH and temperature stressors, the death probability each step was calculated as follows:\u003c/p\u003e\n \u003cp\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAgQAAABwCAYAAABhLmhoAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAAEdHSURBVHhe7d13WBRX+zfw79JRpEuRZu+CIhgVBAU1omKsjy2WmKixxYhGUR9jjT2aqFGjySMaFBQVBXsvIEVEAelFikjvLGXZ3fv948fOyw6gWGJJzue69ko8ZXZndnbmnvucGQRERGAYhmEY5l9NgV/AMAzDMMy/DwsIGIZhGIZhAQHDMAzDMCwgYBiGYRiGBQQMwzAMw4AFBAzDMAzDgAUEDMMwDMOABQQMwzAMw4AFBAzDMAzDgAUEDMMwDMOABQQMwzAMw4AFBAzDMAzDgAUEDMMwDMOABQQMwzAMw4AFBAzDMAzDgAUEDMMwDMOABQQMwzAMw4AFBAzDMAzDgAUEDMMwDMOABQQMwzAMw4AFBAzDMAzDgAUEDMMwDMOABQQMwzAMw4AFBAzDMAzDgAUEDMMwDMOABQQMwzAMw4AFBAzDMP9cUqkUNTU1/GKGaZCAiIhf+DEoLS1FWVkZBAIBVyYQCKCiogJtbW0oKirKtWdeTiKRID8/HxKJRK5cUVERGhoaaN68uVw5826UlpZi8eLF0NLSwvbt26GiosJv8tqqq6tRXl4OVVVVaGho8Ks5lZWVqKiogIqKClq0aMGv/qjJ1vFVn72p7f5OBw8exJUrV7BhwwZYWlryq987IkJISAiOHDmCmJgYiEQiaGlpYfDgwZgzZw60tbX5XThZWVlYvHgxOnfujHXr1kFB4f9fM75JYKGoqCi3jL+TWCzG657OFBQU3tu5ZP369YiMjMSePXtgYmLCr/4ofLQBwbJly+Dh4QGBQACpVArUfnn6+vqwsbHB3LlzYW9vz+/GNCIyMhJTpkxBXl4exGIxV66mpoZ27dph5MiRmDt3LrS0tOT6MW8nLy8PlpaW0NfXR2hoKNTV1flNXpuPjw9WrlyJkSNHYufOnVBSUuI3AQB4eHhgy5YtGDJkCPbt28ev/qidPXsWa9asQf/+/XHo0CG5C4O6fHx8sGrVKri4uGDXrl2Nbou/07x58/D777/j2rVrGDx4ML/6vdu9ezdWr16NyspKmJmZoVmzZkhPT0dlZSXGjx8Pb2/vRk+CiYmJsLKywmeffYYbN25w7e7du4ctW7bUu0BTUlKCQCCAWCzmjtOozUxoamrip59+QocOHbjyv0t+fj7c3NyQk5Mjt25KSkpQUFCARCKRuxgiIhAR5syZg7Fjx3Llf6cRI0bg+vXriIyMROfOnfnVH4X3E7q9ptLSUgQFBaGgoADt27eHg4MD7Ozs0Lt3b5SUlMDT0xPjxo3DvXv3+F2ZRsTGxiI6OhoSiQT9+/eHvb097OzsYGpqipCQEKxYsQLfffcdRCIRv+snSSKRyB2gPhTZQVN24HwXcnNzkZycjJSUlJdeEWVlZSEhIQGJiYn8qo9eXl4eYmJikJCQ8NJ1zMvLQ1JSEpKTk1/a7u+koKAAgUDwVlfCUqkUEonkrdfhwoULcHd3h6amJnx8fBAREYGwsDDcuHEDRkZGuHTpEsLDw/ndOLIrZmVlZbly2fEjLi4OcXFxiI+PR1RUFC5evIjz588jLCwM8fHxXH1MTAyys7PfSQDcFJmZmYiMjOTeX/a6desWzp8/j4CAALnPFxcXh6SkJKiqqvIX9bdSUFB4q/3k7/ZRfrLCwkLExcWhefPmOHDgAHx9feHn54crV67gzp07sLa2Rm5uLvbv38/vyjQiPj4eADBu3Dj4+/vj/Pnz8PPzQ2BgILZu3QoAOH36NEJDQ3k9Pz3p6ekYNWoUli5d+o8JcOpq6gFF1q6xq8GPmeyzv2pdZfXvKtj6ULZs2QInJyfExMTwq5qspqYGBw8ehEgkwqZNmzB+/Hjo6OhAQ0MD/fv3h6WlJaqqqlBaWsrv+kqTJk1CSEgIHjx4gAcPHiA0NBTe3t5o3rw5NDQ04OXlhYcPH3L1ISEh8PX1fW+p8U6dOuHKlSvcZwwKCsL9+/e5IRx3d3eEh4dzny8oKAiBgYEYOnQof1H/ai//tX0gaWlpKCwsROvWrdGqVSu5uo4dO2Ly5MkAgISEBBQUFMjVM/VJJBLExsYCQL1UlZKSEqZOnYqOHTuioqICqampcvWfooqKCoSEhCAqKuqdZAlEIhGioqJw/fp1PHr0CFVVVVxdZWUlCgoK3ijwqK6uRn5+PioqKoDaK/oLFy7g+PHjuH//Plf+KauurkZBQQG3LtnZ2bh48SKOHz+Oe/fu/SPW8V2Ii4tDYGAgSkpK+FVNlpKSgvv376Nt27ZwcXGRqxOJREhPT4eGhgbatWsnV9cUWlpaMDY2hpGREYyMjKCnpweJRILS0lK0bt0a3bp1g66uLldvbGwMPT299xaoqampce9tZGQEQ0NDaGhoICUlBaqqqrCxsYGWlla9NvxMyL/dRxkQREdHAwDMzc3RsmVLfvVLJ1Ix9QmFQkRFRQFAg+N5ioqKUFNTA/4BV1qoXQdZCvdt1oeI4OfnB2dnZ9ja2mLo0KGwt7fHqFGjkJqaCrFYjG+++QYuLi7Izc3ld3+l06dPw9zcHD///DP++usvWFtbw9XVFV9++SWcnZ0xYsQIPH78mN/tk+Lt7Y1OnTph69at8PLygo2NDUaOHIkvv/wSTk5OGDZsGB49esTv9q+jqKjI7bdvKikpCaWlpbC0tKx3Ze7n54e4uDj06dMHrVu3lqt7U2FhYQCA1q1bw8DAgF/9wcXGxqKwsBB6enro1q0bv5ppwEcZEERGRgK12YCGfiARERFAbcCgq6vLr/4kERFqamogEole69WUmb8FBQVITEyEvr4+LCws+NXIy8tDeno6FBUVYWxszK/+VyIi7NixAxMmTEBSUhIWLlyIHTt2oGvXrrh+/Tr27t2Lu3fv4sSJE3B2doapqSl/Ea8kFApRWVkJLy8vLFu2DA4ODjhx4gT++OMP2NnZ4c6dO5gxYwbS09P5XYE6E6Ya09Bv530rLy9HQUEBzpw5g4ULF6JPnz7w9PTEkSNHMGjQINy/fx+TJ09GSkoKvyvQhFnqH2IS4cdKdoKue/KTSCS4ffs2li9fDgUFBcyfP79Oj7cjG17kZx0/FtHR0aipqYG5uTmMjIz41UwDGv+lfSBSqZRLbzd0C8+jR4/g6+sLABg9evRHcdB7F8rKyrB48WK4uLjA1dW1Sa9hw4Zh2bJl/EXV8+TJE4hEIhgbG6Nt27ZydVKpFAcOHEBxcTF69eqFXr16ydV/yt5m3/D398fq1athYGAAf39/7Ny5E8uWLcOhQ4egpqaGu3fvYs2aNTA3N3/rg2xsbCy++eYbnDx5EpMnT8bXX38Nf39/fPHFF4iKisLevXv5XQAAycnJOHDgAA4ePIgDBw7Ivf78808EBQXxu7x3spN5TEwMpk2bhjNnzmDq1KmYOXMm/Pz8MG7cOCQmJmL79u38rkDtfJD9+/fXW78DBw7g8OHDuHXrFr/LJ+1t9llZNsnGxgaxsbFYtmwZRo8eDVdXV7x48QKbN2/GF198we/2RmpqahAXFwcAH+0xQ3YesbKy4lcxjaGPTGZmJrVt25YA0N69eyk0NJQePHhAt27dou3bt5OpqSkBoEmTJlFlZSW/OycvL4/Onj1LBQUF/Kq/TUxMDF29epVqamr4Va9UWFhII0aMIDMzM2rbtm2TXqampjR27Fj+oupZvXo1AaDevXvTw4cPKTAwkAICAsjX15cmTJhAAoGAdHR06MKFC0RE9OLFC7pw4cJ73XbvUnx8POnr69PQoUOpqqqKX/1KIpGIHB0dCQDt27dPrq6wsJB69+5NAAgA7dq1S66eLy8vj0xNTalnz5719tfff/+dAFDr1q3p+fPncnVERHfv3iVVVVXq0KED5eXlceX79+/n3r8pLxcXF7nlNkQkElFlZeVrvaqqqkgqlfIXJUf2WU1NTSk1NZVfTQ8fPiRFRUVq166d3DY4dOhQvfV42WvEiBEkEonklv2+zJ8/nxQUFOjmzZv8qiabOXMmKSkpUVBQEL+qSYRCIXXv3p1UVFQoJSWFNm3aJLd9LCws6OLFi/xu9SQlJZGGhgYNGTKExGIxv5qTmJhIrVq1IkVFRYqOjuZXv9Tfta/VJRKJyNXVlQDQb7/9xq/+IIYPH06qqqoUHx/Pr/pofHTPIQgICMCwYcMgFAqBOhGz7GPq6elhwYIFcHNze+k983fv3sXAgQMRGBiI/v3786vfSkhICJ4+fYrJkyejWbNmXPmaNWtw4sQJhIeHv/SzNUQqlaKsrEzuGQGvQkRQVlZ+5XuNGzcOZ8+eBRrYngDg4OCADRs2wNHREQBw7tw57ir1XW+7d6m8vByrVq1CZmYmlzpWUFBASUkJbt26BW1tbdjb20NJSYlbX5FIBAcHByxZsoS3tP8vICAAQ4YMgbGxMcLCwuSGpUpLSzFs2DAEBQXB0tISt27dgp6enlz/uvLz89GrVy/o6+sjKCiIm6sBAIcOHcLcuXMxYcIEeHt710uNFxYWYsCAAYiJiUFERASXMTtw4ADmz5+PPn36YMWKFXLrJ6OqqgovLy8cO3YMLi4uuHTpklw936+//govL68m34YlkUigq6uL3bt3v3SSmuyzNvYZSkpK4OjoiJiYGNy6dYt7tsjhw4cxZ84cdO/eHRs3boSCgkK9dVRWVoa/vz8OHjyIESNGwNfX92+ZJBYcHIwnT55gzpw59b4jAFiwYAEOHjyI69evw8nJiV8tx9vbGydPnpR7QJVAIEBISAjS0tLg5OSEli1bcpNhiQgqKirYvn17vQnWdcXGxmLAgAHQ1tbG48ePkZ2djZSUFGRlZcHf3x++vr5QU1PDiRMnMHr0aH53TnJyMnr27Il+/frh8uXLjd6hcu3aNQwfPhzm5uZ4/PjxK49Bda1duxbXr19v8nclFovRunVr/PLLLw3OKWtITk4OBg4ciLi4ONy4cQPOzs78Ju/c5cuXUVRUhIkTJza43UaMGIGbN28iMjISHTt25Fd/HPgRwod25MgRAkAdOnQgd3d3+u9//0urV6+mzZs308mTJykzM5PfpUF3794lJSUlCg4O5le9ta1bt1KXLl0oJydHrnzdunXUvn17KikpkSv/kAoLC8na2poA0IwZM+jHH3+k1atX048//kj79u2jBw8eUHV1tVyf8+fPU8uWLd/4auV9KSoqImtra2rRogXp6upyLy0tLRIIBKSsrEw6Ojpydc2aNaOZM2fyFyVHdnU1Z84cfhWVlJRQv379CAAdPHiQX11PUzIECxculCuXEYlEZGdnRwDo9u3bXLnsqnv06NEvzUZt27aN0MQMwYoVK6hZs2akp6fXpJe2tjZ16NCBnj59yl+UHNlnnTt3Lr+KiIgqKyvp888/JwDk5+fHlcsyBE5OTnLt+WTt/s4MwdatW8nMzKzR5b9OhmDTpk2krq4ut0/q6uqSqqoqASBNTU25ch0dHTIyMqKkpCT+ouT4+fkRABo5cmS9K3upVErff/89ASBXV1e5Or6mZgj27dvHvV9j26Uxs2bNIg0NjXr7VGMvTU1NsrW1pRcvXvAX1ajIyEhSVVUlAwMDSkhI4Ff/LWbMmEH9+/dvdHt8ChmC+uHuB/bkyROgNprasmULNm7ciE2bNmHlypX4z3/+02iUXFVVhYiICAQHB6OoqKjRyUaVlZWIjo7GgwcPEBcXV++KnIiQnZ2N8PBwBAcHIz4+Xq6NUChEbm4uqqurkZCQgPj4eGRkZHD1sivwtLQ0BAUFIT4+vt7jghsiEolw8+ZN+Pj44MyZM016+fj44Pbt2/xFyXnx4gXS0tLQvHlzrFmzBuvXr8emTZuwfv16LFiwAP369WvwcbqyGc8FBQUIDg5GVFSU3O12qN1Wubm5ePz4MYKCghATE4Pq6mq5Nqj9bqKjoxEUFIRHjx7h+fPn9W4HFIlEiImJwYMHDxAbGwtRE27j09TUhJ+fHyIjIxEeHo7w8HBERUXB19cX2tra6Nu3L0JCQvDkyROuPjo6utHxatSuk2zssW/fvvxqlJeXIzMzEx07dnzpldbr4F/5NqShseVXPcjmZXV8K1asQExMDLedXvWKiIjArVu3GrxrpSH877su2eds6OpbKpW+tO/LflsSiQRJSUkICgpCZGQkl3WUqa6uRm5uLmpqahAfH49Tp07h6tWrcrf+lZSUICsrC2KxGHFxcYiNjUVmZqbccl7HggUL6m3nyMhIuLq6QlFRER4eHoiKiuLqHj9+jJCQEJibm/MXJSchIQEA0KNHj3pXpwKBAEOHDoWCggLy8/Pl6t6U7DfSvXv3Jl/py2zbtg1Pnz6tt0819oqMjMS5c+de606GZ8+eobq6Gm3btn1lP9l8iAcPHjR4DKuqqkJeXh4kEgkiIiLg7e2N27dvy+1Pubm5KC4uRlVVFWJjYxEbG4vs7Gy55XwS+BHChyQSicjBwYEA0KFDh/jVjXrx4gWNGzeONDQ0yMjIiPr160dr1qwhNTU1uQxBTEwMjRw5koyMjMjAwIAMDQ1p1qxZcmPlnp6e1L17dzI1NSVDQ0MyMDCgmTNncmO4Bw8eJHV1dVJWVqZWrVqRgYEBTZ06lYiINm7cSO3ataMNGzZQ9+7dqWXLlqStrU1r1qx56ZUcEVFubi5ZWloSGhgffdmrd+/e/EXJuXjxIgEgKyurJs8J8PPzI2NjY1q+fDl99tlnZGBgQBoaGjRlyhTKz8/n2p0/f56srKzIzMyMDA0NSV9fnyZMmCCXxXn+/DmNHz+eWrZsSUZGRmRoaEhWVlb08OFDrk1SUhKNGTOG296Ghob05Zdf1svANJVsDsHnn3/+2nMIhEIhOTk5EYAGx1x//vlnAkADBgwgoVDIr66nKRmCMWPGNLh/5ObmUufOnUlBQYGioqK4ctlVt6ura6NXI1R7ZYsmZgj+LrLPOnToUH4VUW2Wp1u3bqSqqkqBgYFcuezKf+DAgSSRSOT61HXgwAFCAxmC8vJyWrx4MRkZGZG+vj7p6urS0KFD5ca779+/T/369SM3Nzfq2rUrde7cmTQ0NMje3p67It+1axepq6uToqIimZiYkIGBAc2ZM0fuM71OhqAxsjkEb5rR/PrrrwkAeXp68quIiOjYsWME4JXZsaZkCKqqqrg5NkeOHOFXfxTWrl1LqJ1r9jI5OTk0c+ZM0tfXp5YtW5K+vj6NHz9ebr6Lv78/9enTh9zc3Khjx47UsWNHUlVVpdGjR1Nubi5RbYZNSUmJVFRUuP1k1apVdd6JZQheW3FxMRITE6Gurt7kW1mICKtXr0ZAQADOnDmDmJgYuLm54dixY3JXmWVlZfj666+Rn58Pf39/xMTE4MSJE7h27Rp++uknrl3Lli2xdetWPHz4ELGxsfD09MSVK1ewe/duAMCMGTOwYMECmJmZ4cyZM3jw4AF3xamoqIjk5GQ8efIEXl5eiImJwbJly7B9+/ZXzviWPffbw8MDx44da9LryJEj2LBhA39RcmRPKOzQoUOTx/kUFBRQWFiIM2fOYO3atYiJicFff/0Ff39/rF27lmunqamJNWvWICQkBHFxcfDz88OTJ0/w448/cm1+++03hISE4PTp04iOjsbDhw+xefNm7rNUVFTgq6++Qnp6Os6fP4/Y2FicOnUKgYGBWLduHbec1yG74nzZleXL1J1vUFd6ejoOHjwIvONHkIaGhuL58+f8YkRERCApKQmdO3duNDP2qYiKimrwEcoxMTGIiYmBubl5vTtg3sa2bdvw22+/4aeffkJcXBwuX76MjIwMzJkzh3sYklAoRFBQEC5evIhjx44hMjISN27cQFpaGrcPz549G/PmzYOhoSGuX7+OgICAen/05114m322oqICkZGRUFJSanA/kUgkuHLlCgC8co5DU5SVlSE2NhbNmjVD+/bt+dUfBdktkV26dOFXcYgIa9euha+vLw4fPoz4+HgcP34cgYGBWL58OZd9KikpQWhoKEJCQuDn54enT5/iwoULuHXrFrZs2QIAWL58OVxdXdGjRw9cu3YNgYGB+P7773nv+AngRwgfUlBQECkrK5OZmRllZ2fzqxsUHx9Purq6tHPnTrnyVatWEQAu4vbw8CADAwOKjIyUa3fkyBEyNTWtd/VcWlpKWVlZVFxcTN9++y1ZW1tzV4Tbtm1rcA7Bxo0bSUNDg8LDw7mywsJCatWqFW3dulWu7fsyffp0AkD//e9/+VWNunDhAgkEgnoz6FevXk2GhoaUnp4uV15WVkbZ2dlUXFxM69evJyMjI24exZdffkm9e/emwsJCuT4yZ8+eJQ0NDYqIiJAr9/DwoLZt29KzZ8/kypsiLi6O9PX1aciQIa+dIaA622zatGnclXt5eTnNnDmTm5vQpk2bV47rUhMzBABo0aJFcle4OTk5NHjwYAJAa9askev3KWYIANDXX38tl1XJz8+nYcOGEQByc3OT6/c2GYLc3Fxq3bp1vTkgZ8+eJVVVVS7zc/XqVUIDc0F27NhBRkZGFBcXR0REW7ZsIXNz80a39bvIEMyYMYOUlJTowYMH/KpXevHiBWlpaXH7UUVFBVcnlUrJw8ODFBUVydLS8pXH1aZkCIKDg0kgEJCFhcVrjeu/LyUlJdSpUycCQP7+/vxqTnp6OmloaNDGjRvlyvfv30/NmzensLAwIiI6fvw4AaBTp07Jtfv+++/JxMSEy5p+9dVXZG9v3+h+wjIErykiIgI1NTUwMzN75biPTFpaGsrLy9GvXz+5cjs7O+7pX6hddklJCXbt2oU5c+Zwkb+XlxeEQiH3yN7MzEwsXLgQjo6OcHR0RJ8+fXD69GmUlJTIjS1R7V/LqksikcDQ0BBt2rThytTU1KCjo/POxu5eh1Qq5Z4C19SMC+r8pTL+/cV9+vRBZWUlkpKSgNoHHi1fvhwDBw6Eg4MD+vTpg/3796OmpgZFRUVA7bMi4uLi4OzsjB9++AFXrlxBeXk5t8wHDx5ALBZj586dmDt3LmbPno25c+fCx8cHaWlpyMrK4tq+Dtn3w/+OmmL06NFQUlKCp6cnpk2bht27d2P8+PHw8PDAokWLMHz4cDx79gxjxozBoUOH+N1fm7W1NU6dOoWhQ4di586d2LRpE4YNG4YbN26gb9++WLhwoVz7pl5Fyto1tf3fycbGBv7+/nBxccH27duxZcsWDBs2DFeuXEGvXr3qPU+jqZ+97mx8mdTUVBQXF8POzq5Oy/97YI+2tjb31E4igqqqKnr27CnXTvZH1GRzBd5kH3pdsvd4k/d6+vQpysrKAAD79u3DF198gW3btmHHjh0YN24cZs+eDU1NTWzfvh2Ghob87q/t4cOHICIYGxt/lA8yS0pKQk5ODtTU1F76DIKEhASUl5djyJAhcuX29vaoqqri5mXIjof85+LY2toiLy+PO0a9zTHnY/FRBQQ1NTXo1q0bhg0b1uAkqoaIRCIIBIJ6t0upq6vLLaOyshICgQBVVVWoqKhARUUFSktLYWBggGnTpkFPTw9SqRQLFizAhQsXMHXqVOzcuRMHDx7E6NGjIRKJmvRFv8tU8tsqLCyEqakp7O3t653cX0VJSaneZEPZv2VPR1yyZAmOHj2KsWPHYseOHTh48CBmzJgh95cGx40bh0uXLsHW1hZ3797F5MmTMWTIEG5SkmxiTmVlJYRCISoqKlBeXg5dXV189dVXTQ4M61JWVoaFhQWMjY3f6LsYPXo0Nm/eDBMTE3h7e8PNzQ137tzBd999h507d2LlypVwcHBAZmbmG/2hGD4nJyccPnwYWVlZ+OGHH7BmzRrEx8dj3LhxOHr0aL1toK2tDT09PRgaGr70d6KrqwsDA4OP4qBtb28PDw8PFBQUwN3dHatWrUJUVBRGjx6NEydO1PuMWlpaaNWqFUxMTF66jlpaWtDX14eRkRHXrqamBlR7u15dKioqUFZWlgvsG3oSomxS3quCkXfJwMAArVu3fqO/DhgeHg6pVAp7e3v07t0bN27cgLu7O5YvX45bt25h+PDhuHDhAj7//HN+1zeSm5sLPT09DBo0iF/1USgoKICWlhYGDhz40luCq6qqGj13KCoq1ptcyJ+s+SH2k78dP2XwIQmFQioqKnqtNG9ISAipqanR6dOn5coPHz5MAoGAGzJYt24dGRsbv/S2xcjISDI0NKSjR4/KlS9YsEBuWGHbtm3UuXPnekMGDd12WFFRQd26daNly5bJtX0fxGIxlZSUUElJyWs91MPf359UVFTo7NmzcuV//PEHaWtrU1RUFGVlZZGZmVm9oZCtW7eShoZGg6n+qqoqCggIICMjI+5Wu+3bt1PLli0pLS2N3/yN1dTUUEZGBuXk5LzWetcllUopNTWVAgIC6M6dOxQVFSWXCiwrK6PU1FQqKyuT68fXlCGDJUuWENWm0ENCQuju3bsUFRXV6O+grKyMnj17Rrm5uS9dv5KSEkpNTa23n75PsiGDRYsWEdWuY2hoKN25c4ciIiJeuo6pqamvTHE3tC3i4+PJwMCAtm3bJtf24cOH1KJFC/rjjz+I6gwZeHl5ybU7dOgQ6enpccNYW7ZseWe3HTYmPz+f0tLSGt0eLyMb4jpx4gSVlpZy2zcgIIASExMb/dwNacqQQW5uLj179oyKi4v5VR8F2T7xqt9/REQEASAfHx+58itXrpCioiJdvXqVqHaiuaKiIl26dEmu3Zo1a0hPT4+ysrKIaieG2tvb17uNW4YNGbymZs2aQVtbu17E9jIdO3ZEly5d8Oeff3K3xZWUlODYsWNyV/SjRo1CdXU1Dhw4UKf3/6X5nz17BtSJ+AoLC7n6qKgonDlzRi461NXVRXFx8Ud/W4mioiI0NTWhqan50qssPoFAAJFIhCNHjnDbtLS0FEePHkXnzp3RtWtXEBEUFBS4oQHUTro7cuSI3LZKTk7mIm1VVVXY2tqiZcuW3MSuL774AjU1NThw4IDc9yWRSJCenl5vYl9TKCkpwdTUFAYGBq+13nUJBAJYWFjAzs4Ojo6O9W6v0tDQgIWFxTv5Q1uy9dbT00OfPn3g4OCA7t27N/o70NDQQOvWrdGyZcuXrp+mpiYsLCzqZRg+JD09Pdja2sLR0RGWlpYvXUcLC4tXprgb2hYdOnRAr1694OnpKfcb9fDwgJKSEhwcHOos4f8eyiT741SZmZn47bffYG1tzU1I09XVRUFBwRsPXzWFnp4ezM3NG90ejREKhUhOTgZq17tFixbc9rWzs0P79u1f+7bAV2nZsiVat27d5EnK75tsn3jV7799+/awsbHBvn37uEylSCTCgQMH0K5dO24oSUFBARKJBAcOHOAygomJiThy5AiGDRvG/Z0ETU1N5OTkyB0TPzUfVUDwJrS1tbFu3ToEBgZiyJAhcHd3xxdffAGpVAoFBQUundOrVy9s2rQJ+/btg6OjI9zd3bFkyRLY2dlxY7RdunTB8OHDsWrVKsyePRuLFi3CzJkzYWFhAalUyh24nZycoKenh2HDhsHa2hqLFy8Gak9idduh9mAvFotfer/0x0YsFkNHRwcvXryAq6srVqxYAWdnZ0RGRnJPjTM2Nsbo0aOxa9cuzJgxA0uWLMHYsWOhpaUltw3WrVsHS0tLLFq0CKtXr4aTkxMKCgowbdo0oDag+/nnn3Ho0CEMGDAA7u7ucHNzg729PebPn18vbccwryIQCLBp0yYUFRXByckJK1aswPDhw3HkyBGsW7eOe3YCEUFHR4f7LY8bNw79+/dHYWEh1q9fz51IBwwYABMTE/Tv3x/W1tZYuXLlR5MmzsnJwbNnz2Bubl5v2IV5uWbNmmHjxo2Ijo7GoEGDsGrVKjg7O+PWrVvYsGEDF0gTETQ0NKCgoAAnJyeMHTsW/fv3h56entxdV2PGjEFFRQX69u2L3r17v/R5Jx+rj+7RxW/q/v37OHXqFEpKSuDk5ISBAwciICAAn3/+Ofe4S6lUirCwMPj7+yMlJQXq6uro2rUrnJycuGiwqKgIx48fR2hoKLS0tDBu3Djo6OggIiICEydO5CL4lJQUhIaGori4GBYWFnBxcUFERASSk5Ph6urKHUzEYjEuXrwIExMT2NjY1PnEH6+MjAwEBwdz45EBAQHQ1dXFlClT0KdPH65deXk5vL29cf/+fairq+OLL75A27ZtERgYiAkTJqBFixZ4/Pgxbt68idjYWFRVVaF9+/YYPXp0vTkNsu8lOTkZqqqq6NKlCwYNGgRra+uXRvkfu6Y8uvj777/nbmv9p9m3bx8WLVqE+fPn47fffuNX/60SEhJw8uRJxMXFQU9PD6NGjcKgQYO4DNbVq1cxceJEXL16FaWlpQgICICOjg6GDx9e79GySUlJePjwIUpKStC+fXs4Oztz++XrPLr4Xbt16xacnZ0xePBg+Pv7y+1fb6Kpjy7+J3n8+DF8fHzw7NkzmJmZYezYsfjss8+479fT0xOLFy9GYGAgYmNjER4eDmNjY4wcObLeA6OePn2KJ0+eoKysDD169OAexQ326GKGYfLy8sjExISsrKzqzSF48uQJrVu37q3Gnj924eHhtHnzZrpx4wa/6oO7cuUKtWjRgru97E3NmzfvrecQvKm4uDjauHEjXb58mV/1RhITE0lDQ4MGDx7c6ByCf5u//vqLtLW1G/zjXK9DNodAdjvrx+gfkyFgmI9RTU0NIiIioKSkBEtLy3oz2pkP5/Llyxg/fjzu3LkDW1tbfnWTyf6IUPfu3T/acfWmqqqqwpMnT9CiRQt07dr1k87OvStHjx7FwoULERUVhdatW/OrmywmJgYlJSXo2bPnG91N8j6wgIBhmH+lxMREnDx5ErNmzWrwCX8Mg9pn2Fy8eBHz58+HtrY2v/ofhQUEDMMwDMN8+ncZMAzDMAzz9lhAwDAMwzAMCwgYhmEYhmEBAcMwDMMwLCBgGIZhGAYsIGAYhmEYBiwgYBiGYRgGLCBgGIZhGAYsIGAYhmEYBiwgYBiGYRgGLCBgGIZhGAYsIGAYhmEYBiwgYBiGYRgGLCBgGIZhGAYsIGAYhmEYBiwgYBiGYRgGLCBgGIZhGAYsIGAYhmEYBiwgYBiGYRgGLCBgGIZhGAYsIGAYhmEYBiwgYBiGYRgGLCBgGIZhGAYsIGAYhmEYBiwgYBiGYRgGLCBgGIZhGAYsIGAYhmEYBgAERET8QubjQ0QoLS2FgoIC1NXVoaSkxG/ytxGJRCgvL4eamhrU1NSgoPDvjiOlUilKSkqgrKwMdXV1KCoq8pswDMN8clhA0IikpCTEx8eDiCAQCNCqVStYWVl9sJNhdXU19u/fj9OnT+Pbb7/FtGnT+E3+NvHx8di9ezeSkpLw559/wsLCAllZWfjzzz8xbdo0WFhY8Lt8UogIlZWVTQ52ysvL8euvv+Lq1av44Ycf4OrqCrFYjP3798PS0hIDBw4EagMpIoKqqip/EZ+U6upqCAQCqKio8Ks+qJycHCQlJaF3795QU1PjyrOzs5GYmIg+ffp88G1fWlqKwMBASCQSrkxJSQlqamro16/fKz+frH+bNm3QuXNnfnWDCgsLERISArFYDAUFBSgrK0MsFkMqlUIqlUJHRwd9+/aFsrIyv+t7UVVVhUePHiEtLQ3l5eUwNDRE//790bJlS6D28wcGBkIgEAAABAIBFBUVoaGhgT59+nD7YWFhIcrKyho8/mRnZ+Phw4fIyspC8+bN0aVLF1hZWbHg/VWIadDUqVPJ2NiYBg0aRH379iUTExNasWIFicViftP3pry8nBwdHWnr1q38qr/d06dPyczMjGJjY4mIKCYmhpycnOjhw4f8pp+crKwsGj58OMXExPCrGlVTU0N9+vSh33//nYiIKisradKkSXTo0CGuzdatW2nNmjV1en2aNm7cSBs2bOAXf3AnT56kjh07Umpqqly5l5cXmZmZ0fPnz+XKP4TExEQaOnQoOTg4kKOjIzk5OZGOjg716NGDiouL+c3rWbJkCQGgBQsW8KsaFRQURIMGDSJHR0eytbUlNTU16tKlCw0aNIj69+9PixcvpoqKCn639yYuLo769u1LI0aMoIkTJ5KVlRXZ29tTQkICERFFRUWRvb09OTo6ci9NTU2ytbWl6upqIiJ68eIFubu70/Lly+nKlStyy/fz8yMrKyvq1asXTZgwgUaOHEndunUjd3d3uXZMfSwgaIBUKqXOnTvTr7/+ypVt27aNdHR0KDo6Wq7t+1RZWUlOTk60fft2IqJGgxOpVEo1NTUklUr5VUREJJFIqKamhl/MaahvbGwsWVhYUFxcnFy7xvD71yWRSEgikfCLX4q/PP6/65JKpY1um4Y8e/aMjIyMKCQkhF/V6LaSSqXUt29fOnz4MFcmFovlPtfcuXNp6tSp3L/fxqvWSSKRvLSeXrHNZMRicb3v5quvvqKZM2fKlcm86nNRA9vlddXU1NDz589JKBTKlZ86dYo6derUYEBgbm5OmZmZcuUfA6FQSFZWVvTDDz/wq+o5d+4cOTg4kI2NDc2ePZtf3SRBQUFkYmJCd+/e5Vd9MFVVVXLBUGlpKfXs2ZNWrVol106muLiY2rZtS9u2bePKAgMDafXq1XTz5k3ueEhEFBoaSnp6erR69WqqrKzkyuPj4+nYsWNvtR/+G7AhgwYkJiZiyJAhOHz4MIYMGQIA8PX1xaxZs3D37l1YWlryu7wXVVVVcHV1Rffu3dGtWzdER0dDSUkJ8+fPR5s2bQAAN2/exPXr10FEqKmpwbBhwzB06FCgduzb39+fS8cpKSlh1KhR+OyzzwAAycnJ8PT0RHl5OQBg0KBBGD58OAAgNjYWLi4uuHbtGtq2bYtffvkF6enpmDFjBnr37o2jR48iNDQU/fv3R2lpKQICAmBmZoZly5ZBX18fqE21nzhxAvHx8QCAtm3bYvr06WjRogX4wsLCcPr0aairq8PGxgaXLl1CZWUlli5dCqFQCC8vL+Tk5GDu3LlwdHQEalP/N2/exLVr10BEUFdXx+TJk9GlSxckJSXB09MTQqEQQ4YMwa1btxAbGwtnZ2dUVVVh48aNmD59OszNzTFy5Eh06NABHh4eSEpKAhGhWbNmmDZtGtq3bw/Ubks7Ozt88803+Prrr3Hy5EkEBASgX79+mDJlCm7duoXly5dDIpFg/PjxMDIygrKyMu7du4dWrVph4sSJ6NChA/73v/8hMjISI0eO5LY1AEgkEhw9ehSPHj2Cra0txGIxYmJiIBAIMGfOHHTq1AkAUFlZCU9PT26bWlhYYPr06dDS0sK1a9dw6dIlWFhYwMLCAlevXkVhYSF+/PFH9OjRg3svAEhLS4OnpydKSkpAROjQoQOmTJmC4OBg/PDDD1BUVMTo0aPRpk0bTJkyBQKBAHfu3MGlS5cAAKqqqpg8eTK6du2K27dv48KFC9DT00OnTp0QHh6OwsJCjBo1Ci4uLnLv2xixWIyZM2fCwMAAWVlZSE5ORlVVFVasWIGpU6cCAHx8fPDjjz/i2rVrMDMz4/qeOXMGS5cuxYMHD9CqVas6S3071dXVqKqq4hdzlJWV0axZM36xnMuXL2PWrFk4d+4c97trSFZWFiZOnIjVq1fj4MGD0NPTwx9//MFv9krHjx+Hu7s7QkNDYWxszK9+axUVFRCJRFx6vy7Zb/BVwyLV1dVwcnLCgAEDsHXrVn41vLy88MMPP+D27dvo0KEDULvf//nnn8jKysLMmTO58ilTpiAlJQU3b95E8+bNeUtiXokfITBEPj4+1KZNGy4bkJOTQyNGjKABAwZQaWlpg1eM70JoaCitX7+eioqK+FVEtRkCZ2dn6tOnD4WFhVFOTg65ublR3759qby8nKRSKc2ZM4dOnjxJJSUl9ODBA+rTpw8FBwcTEVFISAhZWVlRZGQkFRQU0Pbt22nevHlEtWlzOzs72r59OxUVFdHt27epZ8+e5OfnR1Q7RFA3QxAaGkqdOnWiY8eOERFRQkIC2dvbk7W1Nd24cYPCw8PJysqK1q1bR1R7JblkyRIaM2YMpaamUlpaGo0bN46+++672rWTV1FRQTt27CBlZWXau3cvxcbG0uzZs6lTp060Z88eio6Opv/+979ka2tLJSUlRER04cIF6ty5M129epXy8vLo559/JhsbG0pJSSGpVEre3t6kpKREbm5uFBkZSVu3bqWVK1fSw4cPycDAgE6dOkWJiYkkFAopNjaWxo0bR9HR0VRcXEw///wzDR48mLuykUgkchmC58+f0+jRo2nGjBlERFRYWEiTJk2iUaNGUWxsLKWmplJKSgr16NGDFi5cyF29HD9+nMaPH085OTnculPt9srIyCBbW1vq3r07BQUFUU5ODi1evJgcHByooKCAiIiWLVtGo0aNopSUFEpPT6fx48dz6eWCggJasmQJaWtr06FDhygqKopmzZpFV69elXuvmpoamjhxIq1atYqKioooIiKCnJ2dKSwsjMrKymjs2LE0duxYiouLo7S0NCIiun79OnXr1o2uXbtGBQUFtG/fPurRowdlZGSQUCik5cuXU7Nmzejo0aOUn59PJ0+epDZt2tC9e/eIiKikpITWrFlD4eHhcp9FRiQS0cCBA8nY2JhOnz5NaWlptGfPHtLX1+cyOadPnyYTExP65ZdfyNvbm44fP06nTp2iBQsWkIWFxTvPEPj4+NCAAQPIycmp3sve3p7Wrl37ymzJjBkzyNHRkV8sRyKR0OzZs8nNzY2kUim5uLjQN998w9UXFBTQunXrGsxo8bm7u5ONjQ2/+J3Ztm0b2dvb19seTk5OZGdnRx4eHvwuRLX7965du+irr74iOzs7GjBgACUnJ/ObkVgspuHDh9OECRP4VfUUFhaSpaUlubm58auYJmIBQQM2b95MSkpK5OzsTC4uLtSrVy8aOnQoxcTE0PPnz+nHH3+kqqoqfre3tmPHDtLW1pZLy9dVWVlJAwcOpMWLF3NliYmJZGRkRL6+vkREVFRURKmpqRQcHEyhoaHUp08fbujj2rVrpKurS1euXCGxWEwikYji4+OJiOjnn3+mjh07UkJCApWWlpJQKKRRo0bRF198QdTAkEFxcTFZW1vTX3/9RVR7UhkxYoTcgWv16tU0ZswYotqAQUtLi3x8fEgoFFJFRQX973//I01NTcrIyOD61OXt7U0mJib04sULoto0YbNmzbix/oiICLK0tKSIiAgiIho0aBB9++23XP+qqipq164dbdq0iag2INLR0eECJJFIRBUVFZSenk6mpqb06NEjrq9IJKK8vDyKjY2lkJAQOnfuHOnr63MnMH5AQEQ0e/ZsudT6ggULaPr06dy/iYh++eUXsrKy4vaf9evXk7+/v1ybugYPHiy3TmlpaaSvr09+fn6UmZlJJiYm5OHhQRUVFVRRUUGenp5kbGzMbaMdO3ZQx44dSSQSERFRWVmZXCqVatPYAwcOpOnTp3OBSXJyMheYzp49W+57JSL6z3/+Q+PHjyehUEilpaWUlJREpqamtG/fPiIi+vXXX6lr167cMiQSCX3++ef05ZdfEtXuT+rq6rR371655cqIRCJycHCgr7/+miurqqqiwYMH06JFi4iI6OzZs6SlpUVDhw6lsWPH0pgxY2j8+PFkbW1NFhYW3H7zOuLi4uoNQchkZWVRcHAwhYSE1HsFBQVRYmLiS1PSmZmZ1LJlSy6Ibszp06dp4MCBXPDJ/109efKEdHR0aNeuXXV6NWz48OE0a9YsfvFrCQ8Pp/z8fH4xERGlpKRQUFBQve0h2yaNzeOQSqV05MgRcnd3p5EjR9KkSZMoPT2d34wiIiJIV1eXzpw5w6+q58WLF9S+fXtav349v4ppoldPqf6XISJERUXB0dERa9euxbJly3DkyBH4+vqiS5cuKC4uRmRkJEpLS/HixQuIxWKub3FxMVJTUyEUCrkysViMjIwMZGZm1mubkZGBmpoaruybb75BQEAA2rVrx5Xxye54kGnRogW0tLQQExMDANi7dy/mz5+Pe/fuISwsDGVlZaioqAAA2NvbY+7cuXBzc8Nnn32GDRs2QFtbGwCQkpKCgoIC7Nq1C6tWrcKKFSsgEAhgYWGBmpqaeilB2YxlGdnIU920pIaGBkQiEQAgIyMDJSUlOHv2LNzd3bF8+XLcvXsXDg4OcturLiJCixYtoKmpCQBQUVGBnp4eN8QgS0VKJBKIRCKkp6fD3Nyc66+qqopWrVpx6XRZ6t/IyAioTfGqq6ujNjCWW5+8vDzMnz8f27Ztw8OHDxEdHQ2xWIzq6mquDZ9UKpXbTvxlAsDw4cMhFApx9epV5OXlISMjgxvyaAgRwdDQkPt3y5YtoaamhrS0NCQlJUEoFMLPz4/bpjdv3kS3bt249hKJBMbGxtyMcg0NDbkZ+QDQrFkzbNiwAWlpaRg0aBBGjx6N5ORkbt/gr4dIJEJiYiJSU1OxcuVKrFy5Ejt37kTXrl2570QsFkNbW5v7rhQUFGBsbIz09HRIpVJ07NgRwcHBr7xbpmPHjtz/q6qqwsLCAqmpqUDtexgaGuLAgQM4c+YMTp06BR8fHyxduhSos08CQEFBATIzM7nZ/hUVFRAKhSgsLEROTg7X7o8//sDJkyeRk5Mj99tE7fIkEkmjL/53zXfy5Eloa2vDycmJX8UpLS3F5s2b0a9fP8THx+POnTvIz8/HixcvEBUVBSJCt27dcPfuXcycOZPfXU55eTliY2MbHJogImRnZyMnJ4fbTqWlpaisrEROTg4KCwuB2u20bds2XLt2DXl5efXWkb8N+K+630FdAoEAM2fOxJYtW+Dr6wsiwooVK/jN4OvrC319fTg7O/Or6mnWrBk0NDRQUFDAr2Ka6P3dzP6JKC8vR0xMDCZPnowBAwbwq6GkpISsrCxs3boVRUVFMDc3x5o1a5Ceno69e/dCIBBAKBTCzc0NFhYW2LhxIwoKCqCoqIgxY8bA2dkZFy9exNmzZ6GmpgZtbW2sXLkSGhoa0NbW5g7CL1P3pCQWi1FTUwM9PT1ER0djy5YtuHTpEnfr29WrV7mDoJqaGn766SfMnTsXd+7cwc6dO/HkyRP4+flBXV0dFhYW2LlzJzf2JpFIUFFRAWVl5XoHgqaSnSCbN28ORUVFLFu2DNbW1lx9WVnZS8f6BAIBd1DhH1zq/ltFRQUqKir1xnirqqq4gOJ1HDt2DI8ePcKTJ0/QokUL5OXlYc+ePfU+w8tQ7S2rdXXo0AE2NjY4ffo0iouL0bVr1wbnUNRV98Qk+75lBz81NTUsWLBA7iRTWloqN5bN/wx8RIQBAwbgypUrCA0NxalTpzBhwgR4eXnBxcWl3nooKChATU0N/fv3x86dO7nyiooKuXayW91kt3pVVVVBRUWFu7WzKXNx6p6spVIpCgsLoaury5XJ5sKg9rdZ978yZ8+e5X5vpqamWLduHQICAnDkyBGYmZnh+fPnGDlyJIYMGYKIiAju5D579mzo6elxy7l37x52797d4O16IpEIQ4YMwfr16xu8tU0oFOL06dNwdnaGiYkJv5pTUVEBc3Nz3L9/H/fv3wcRISEhAQoKCti6dSs8PDygrKxcbw5IQ54/f46SkhK5AFFm7969CA8PBxHBzs4Oc+bMwbFjxxAVFYUWLVogMzMTCxcuhJ6eHmJiYiASiVBYWIivv/5aLqA8deoUzp8/3+AtqdXV1Zg3bx6++uorfpUcJSUl2Nra4s8//5Tb12TbbNKkSdDS0uJ3q0dLSwvW1ta4f/8+SkpKmtSHkccyBDy5ublISUlBz549+VVA7QGosLAQ06ZNwy+//ILw8HBcv34dbdq0wXfffYdp06ZBR0cHR48exfPnz3H//n2sXr0au3fvhqOjI168eIE//vgDixYtwtatW1FeXg5fX18AQEREBH799VeUlpby35ZDRAgLC+P+HR0djeLiYjg5OUEoFKKmpgbq6upA7YkkPT2dOwBfv34dv//+OywsLDBjxgy4u7sjMzMTJSUlcHFxQU5ODp4+fcot+5dffsHu3buB2pNAXQKBQO5k3ZC6J4dOnTqhR48e3CQ0AAgNDcWMGTNQWVnJldVFtVemsvdWUFAAEXH/li1fFqy4uLjg3r17XMAUFxeH58+fc1cXsuXxT5CC2vucBQIBMjIycO7cOQiFQqioqHBXvCkpKaioqKj33nWXJRAI5AInJSUlLit08uRJ7sp29uzZuHjxIo4ePYphw4Zx7RsiEAgQHh7OZVpCQkKgrKyMXr16wcrKCsbGxrhz5w7XPiIiArNmzUJxcTHX/2XfEQCUlJTghx9+QFlZGRwcHLBnzx4YGhoiISEBqM2kyNbDy8sLeXl5GDFiBG7cuMFln6qrqzF37lzcu3cPAKCoqIi0tDSkpKQAAPLz8xEdHc1N0pVIJHj8+DH3ORuiqKiI8+fPIzExEQAQEBCAe/fuYcyYMUDt9yl71SWtzV4pKioiLy8P//vf/7B8+XLs2bMHCQkJuHTpEioqKvD8+XOsXr0aa9asgYeHBwBgwIABmDx5Mtzd3eUCDwCYMGECAgICcPv27XqvwMDARoMBAHj06BEiIiIwZcoUfhXCwsLg6emJ6upqGBkZwcfHB7dv38bdu3dx69Yt2NrawtXVFcePH4eysjJEIhGioqJeeSX85MkT6Ovr1wtAHj58iKtXr2Lbtm3YvHkz/P39ER8fj/z8fAiFQuzYsQPTp0/H7t270aFDB1hbW2PBggVYsGBBvezSihUrEBgYWG973L59Gw8ePMCMGTPk2qP2t5Sens59b6mpqTh+/Djs7e3lfk83btxAUVERxo4dW6f3yy1duhTZ2dlwd3fnto9UKkVISMhrB/T/Rorr1q1bxy/8N4uMjMSVK1ewdOlS6Ojo8KuRm5uLx48f49tvv4WGhgYSEhJQUlICDQ0NrFu3Dunp6UhPT0dlZSWmTp2KZ8+ewcfHB0+fPkWXLl2QmZmJY8eOITs7Gzdu3EBGRgbatWsHa2trHD9+HBs3bsSkSZPkrkxQe3I/ffo0bt26BW1tbeTn5yMiIgKHDx/GjBkz4OrqCl1dXaSlpeHSpUsQi8W4e/cuHj16BJFIhI4dO0IqlWLPnj2oqKhATEwMLly4gOHDh2PQoEFo06YNxGIxTp48ifLycty+fRtPnz6Fm5sbhEIhjhw5goiICGhoaKB169Y4ffo0AgMDIZVK0bVrVwQFBeHixYsQCoXo2rUr8vLycOLECWRmZsLAwAA9e/ZE165dcfz4cTx//hwJCQnw8fHBlClTGrzaSU5OxokTJ5CcnAx9fX3o6urCy8sLERER0NXVhbGxMby9vREeHg5FRUVYW1vDzs4OAQEBCAkJ4R6c9Pnnn2Pu3LnIz8+Hp6cnnj59CoFAgJYtW8LAwACoTUU/fvwYiYmJCAgIABFh/PjxuHLlCqKiolBcXIzbt28jNjYWzZs3h4WFBS5evIjbt2+jpqYGHTp0wKNHj+Dv74/i4mKYm5ujTZs2qKiowKVLl1BeXo6LFy9i5MiR0NTUhIGBAS5evAhjY2PMmzePv+pyvLy8UFFRgaqqKkRGRuLAgQOYNGkSJkyYAAUFBXTv3h1eXl7IyMhAUlISvL29MW7cONjY2ODhw4c4efIksrOzIRAI0KFDh3oHdNSenL28vBAcHIzi4mKcO3cONTU1WLRoEXR1dVFWVoarV6+itLQUly9fhqurK/r3748nT57g+vXrKC4uhre3NzQ1NTFz5kwoKSkhNDQUwcHBaN68OVJTU3Ho0CFoaWlh3bp1UFNTQ1FREXeHS91hHhmpVIoTJ05AT08PN27cwIkTJ3Ds2DFMmDAB8+bNg6KiIpKSkhAaGlrvCjIxMREhISGYNWsW0tPT8ejRI8ybNw8qKipIT09HQUEBdHR0oKamBhcXF+jr6+Pq1auwtLREcnIytLS00LNnzwYDRwUFhZe+GvP7778DAJYvX14vg3H8+HH8/vvvmDZtGvdwLAUFBQhqg+4LFy5AX1+fu0MjPT0dX375Jbp164a2bdvKLasuX19flJWVYfbs2XLveefOHYhEIowbNw6ampoIDAyEnp4ecnNzYWNjAysrK2hra+Ps2bMYPnw4bt68CWtr6wYfAMRff/6Lvw0B4NKlS1i6dCn8/Pzg7e2Nw4cPw9zcHFu2bJH7Hnfv3g09PT189913cv1fpmXLlrC1tcWpU6fg4eEBf39/eHp64syZM2jfvn2DWV/m/2O3HfLk5+cjKysLnTt3bjA1GBcXh5kzZ+LixYvQ09PDvHnzMGDAAERERMDCwgLz58/HX3/9hatXr8LT0xOofaLW/v37kZOTg3nz5sHd3R1//vkntLW1IZFIoKCgABUVFRQUFHDZCf57SyQSPH36FJWVlWjXrh3y8vIQHx+PNm3ayGUzqqqqEBoaiuLiYnTr1g0KCgrIysqCiYkJLCwskJaWxqUA27Vrh27dusn9aGNjYxEfHw9NTU3Y2tqiRYsWKCoqQmxsLHdQad++PZKTk0G1Y6pt27ZFYWEhd8uamZkZBAIB0tPTIRAI0KJFC3Tp0gUKCgrIzc1FeHg4ampqYGlp2eBBBnUyNbJtY25ujsTERAhqU8QWFhZITk4Gak8e3bt3h4aGBqqqqhAWFob8/Hy0bdsW3bt3h4KCAkpLSxEbGwsFBQVUV1ejffv23FwC1I4xBwcHQ0NDA/369YOKigpycnIQFhYGFRUVWFlZ4cWLF6iqqkK7du2QkZEBsVgMiUQCc3NzlJaWori4GAKBALq6uujYsSMkEgnCw8ORmZkJa2truRPfkiVL8Pnnn78yQzBkyBA4ODhgypQpePr0KczNzes9MTMvL4/LInTv3p27BTUtLQ2ZmZlQUlKCRCJBr169GgwIUJuqjoiIQHZ2Nlq0aAEbGxtu+EosFiM8PBwvXryAjY0NTE1NgdogNSwsDNnZ2TAzM0PPnj25fWTXrl04d+4cvL298fjxY6irq+Ozzz7jhodOnTqFnTt3wt/fX26OhExNTQ0GDx6MMWPGYPz48Xj06BFatWoFa2tr7iq8qKgImZmZ6Nixo1zKurCwEJmZmejRowfi4+OxfPly/PXXX9DU1MSKFSvQqVMn6Ovr4+zZs/Dw8EBZWRmmTJmCX3/9FV5eXjA0NMQ333xT59O8vdjYWGhoaMjdHimTnZ2NgoICdO7cucEMQ2JiIpSVldG6dWug9nbLDRs24ObNm1xQ25Bnz55BLBZzt+TJXLx4Ed7e3vjrr79ARJg2bRrmzZuH27dvQyqV4scff0R8fDyWLFmCs2fP4ttvv8WMGTMwaNAgueW8qerqasTHxyMjIwMSiQQmJibo0aOH3HcokUgQExMDXV3dehmOppDtzzk5OWjWrBm6dOnS4LZn5LGA4DXFx8dj2LBhmDRpElRVVfHo0SMcO3YM/v7+OHv2LAYPHozLly9DX18fP/30E/bt24e2bdsiPDwcHTp0wNKlS7F+/Xq8ePECPXv2RH5+PlxdXRsdomD+OZKTk6GsrAwlJSWsXbsWO3fufOU4p5OTExwcHPCpJfJ27dqFU6dO4cGDB3LBC2rnjbi5uWH06NEYMWKEXJ1MTU0NN8Fx2bJl/OomIyKsXLkShYWFMDMzQ3h4OH7//XcEBQVhzZo1mD59Op49ewYA+O233+Dr68tllr788ssGs4QfUlVVFb7//nuMGDECrq6u/OomKSsrw/z582FsbAxp7d/lOHDgADZv3ozbt29j7NixePDgAWxsbLB06VJs374dYWFhGD58OCZNmtRoUMl8+lhA8JoqKyuRkZGBoqIipKamwt7eHiYmJpBIJLh9+zZKS0thaWkJZWVlGBkZISgoCBkZGTA1NUW/fv2gpqYGkUiE4OBgpKeno1WrVtyVOPPP5uXlBV9fXxgaGsLOzg6TJk3iN+FIJBIcPnwYf/zxBwwMDLBkyRJu/P1jd+PGDezZswfp6emYPXs2vv32W7krX7FYjJKSEujq6jaYUkbt+i9ZsgR2dnaYOHEiv/q1VFRUcBPN7OzsYGJiAm9vbzx58gTDhg1DSUkJBg0aBE1NTYjFYkRFRaG6uhpWVlbcfJyPhUQi4bbd2ygsLMTdu3ehoKAAR0dHaGtrY82aNTA1NYWFhQUUFRXh4OAAVVVVVFRUIDIyEioqKrC0tKw35MH8c7CAgGHeE6FQiNjYWKirqzc487suIuJSqlKpFJqamtwff/nY5eXlcX+ZU1FRkRtC+pj4+PggMTERq1at4lf9a61fvx7du3fHuHHj+FXMvwQLCBiG+deprq6GRCJ55aOG/02EQiGUlJRe+ahh5p+LBQQMwzAMw7DnEDAMwzAMwwIChmEYhmFYQMAwDMMwDFhAwDAMwzAMWEDAMAzDMAxYQMAwDMMwDFhAwDAMwzAMWEDAMAzDMAxYQMAwDMMwDFhAwDAMwzAMWEDAMAzDMAxYQMAwDMMwDFhAwDAMwzAMWEDAMAzDMAxYQMAwDMMwDAD8Px28CG+7xyBGAAAAAElFTkSuQmCC\"\u003e\u003c/p\u003e\n \u003cp\u003ewhere \u0026alpha; and \u0026beta; are empirically chosen coefficients [\u003cspan class=\"CitationRef\"\u003e19\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e30\u003c/span\u003e].\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\n \u003ch2\u003e2.6 Fertilization Criterion and Data Recording\u003c/h2\u003e\n \u003cp\u003eThe first time any sperm entered within one unit of the location of the egg, indicating successful egg contact, was operationally described as fertilization. The simulation stage at which this occurrence had place was the fertilization time. At 50 steps, runs that were not fertilized within the allotted period were censored. In order to facilitate qualitative examination of sperm aggregation and movement patterns across time, spatial sperm distributions were captured during simulations at regular intervals and visualized as 2D heatmaps [\u003cspan class=\"CitationRef\"\u003e7\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e8\u003c/span\u003e].\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec9\" class=\"Section2\"\u003e\n \u003ch2\u003e2.7 Parameter Sweeps and Replication\u003c/h2\u003e\n \u003cp\u003eIn order to boost statistical power and account for stochastic variability, a thorough parameter sweep was carried out spanning chemotaxis strengths, pH, temperature, and flow conditions, with 50 replicate simulations per parameter set. Distributions of fertilization events, mean fertilization times, and fertilization success rates were retrieved for each condition.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e\n \u003ch2\u003e2.8 Statistical Analyses\u003c/h2\u003e\n \u003cul\u003e\n \u003cli\u003e\n \u003cp\u003eSurvival Analysis: Taking into consideration censored data (non-fertilization), Kaplan\u0026ndash;Meier estimators calculated the fertilization probability with time [\u003cspan class=\"CitationRef\"\u003e27\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003e32\u003c/span\u003e]. The differences between the chemotaxis groups were evaluated using log-rank testing.\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eDistributional Comparisons: Fertilization time distributions between treatments were examined using two-sample Kolmogorov\u0026ndash;Smirnov (KS) tests [\u003cspan class=\"CitationRef\"\u003e28\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003cli\u003e\n \u003cp\u003eVariance Analysis: Without taking censoring into account for demonstration, one-way ANOVA assessed variations in fertilization times across several chemotaxis strengths [\u003cspan class=\"CitationRef\"\u003e28\u003c/span\u003e].\u003c/p\u003e\n \u003c/li\u003e\n \u003c/ul\u003e\n \u003cp\u003eThe significance level was set at 0.05. To ensure thorough interpretation of simulation results, the SciPy and Lifelines Python libraries were used for statistical and survival analyses, respectively.\u003c/p\u003e\n\u003c/div\u003e"},{"header":"3. Results","content":"\u003cdiv id=\"Sec12\" class=\"Section2\"\u003e\u003ch2\u003e3.1 Baseline fertilization dynamics in the model\u003c/h2\u003e\u003cp\u003eWe first ran the model for 50 independent stochastic replicates under default settings (no intentional change in chemotactic strength) to create baseline behavior. 100% of runs (50/50) experienced fertilization under these baseline conditions. Under the usual parameter configuration, the mean time to fertilization was 36.16 steps (SD\u0026thinsp;=\u0026thinsp;0.37; median\u0026thinsp;=\u0026thinsp;36.00), suggesting very consistent timing. Figure\u0026nbsp;1 shows a near-linear drop in motile spermatozoa over time, which is consistent with stochastic mortality and continuous attrition through fertilization. The baseline model's stochastic combination of persistence, noise, and weak directional bias is reflected in the representative single-sperm trajectories sampled from the simulations, which exhibit a mix of ballistic, directed approaches and meandering exploratory movement (Fig.\u0026nbsp;2). When combined, these baseline results offer a framework for analyzing parameter perturbations and verify that the simulation generates consistent, reproducible fertilization kinetics under default settings.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec13\" class=\"Section2\"\u003e\u003ch2\u003e3.2 Effect of chemotactic strength on fertilization outcomes\u003c/h2\u003e\u003cp\u003eTo describe the effects of directed guidance on the timing and likelihood of fertilization, we conducted 20 replicates for each of the six levels (0.0, 0.1, 0.3, 0.5, 0.7, and 0.9) of the chemotactic guidance parameter.\u003c/p\u003e\u003cp\u003eChemotactic strength and fertilization performance had a clear and nonlinear relationship:\u003c/p\u003e\u003cp\u003e\u003cul\u003e\u003cli\u003e\u003cp\u003eChemotaxis\u0026thinsp;=\u0026thinsp;0.0 (random walk control): Fertilization success was reduced to 40% (8/20 runs). For successful runs, the mean fertilization time was 1.88\u0026thinsp;\u0026plusmn;\u0026thinsp;3.52 steps (median\u0026thinsp;=\u0026thinsp;0.50). The low mean and low median reflect that the successful runs were dominated by early, chance encounters in the absence of directed guidance; the majority of runs failed to fertilize within the observation window.\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eChemotaxis\u0026thinsp;=\u0026thinsp;0.1: Fertilization success recovered to 100% (20/20), with a mean time 36.00\u0026thinsp;\u0026plusmn;\u0026thinsp;0.00 steps, effectively matching the baseline timing but indicating stochastic variability is reduced in this setting for the chosen initialization and parameters.\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eChemotaxis\u0026thinsp;=\u0026thinsp;0.3: All replicates fertilized (20/20), and the mean time decreased substantially to 10.45\u0026thinsp;\u0026plusmn;\u0026thinsp;9.24 steps (median\u0026thinsp;=\u0026thinsp;10.50), indicating that modest chemotactic sensitivity significantly accelerates encounter rates.\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eChemotaxis\u0026thinsp;=\u0026thinsp;0.5: 100% success, mean\u0026thinsp;=\u0026thinsp;3.80\u0026thinsp;\u0026plusmn;\u0026thinsp;5.05 steps (median\u0026thinsp;=\u0026thinsp;1.00), showing marked acceleration.\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eChemotaxis\u0026thinsp;=\u0026thinsp;0.7: 100% success, mean\u0026thinsp;=\u0026thinsp;3.05\u0026thinsp;\u0026plusmn;\u0026thinsp;3.80 steps (median\u0026thinsp;=\u0026thinsp;1.50).\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eChemotaxis\u0026thinsp;=\u0026thinsp;0.9: 100% success, mean\u0026thinsp;=\u0026thinsp;2.85\u0026thinsp;\u0026plusmn;\u0026thinsp;2.67 steps (median\u0026thinsp;=\u0026thinsp;2.50).\u003c/p\u003e\u003c/li\u003e\u003c/ul\u003e\u003c/p\u003e\u003cp\u003eThe population-level plots (Fig.\u0026nbsp;13: fertilization success vs. chemotactic strength; Fig.\u0026nbsp;12: average fertilization time vs. chemotactic strength) summarize these trends, and the representative spatial plots (Figs.\u0026nbsp;6\u0026ndash;10) show increasingly directed and focused approach trajectories as chemotaxis increases.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec14\" class=\"Section2\"\u003e\u003ch2\u003e3.3 Fertilization time distributions and temporal profiles\u003c/h2\u003e\u003cp\u003eAs chemotaxis increased, there was a noticeable shift in the distributional shape of time-to-fertilization (Fig.\u0026nbsp;11). Fertilization periods showed a broad, bimodal-like pattern with many censored or failed runs at the random-walk extreme (0.0), which explains the poor success rate. Distributions were heavily right-skewed yet closely grouped at the early stages (0\u0026ndash;5 steps) for chemotaxis values\u0026thinsp;\u0026ge;\u0026thinsp;0.5, suggesting quick targeting and significantly lower inter-run variability. While relatively moderate increases from 0.1 to 0.3 result in a significant reduction in mean time, larger increases (0.5 to 0.9) offer diminishing returns in mean reduction but do lower variance (compare SDs across groups), according to the median and interquartile structure. Two separate impacts of chemotaxis are reflected in these distributional changes: (1) a higher chance of a sperm coming into contact with the egg during the observation window (better success) and (2) a concentration of successful events at earlier time points (better speed and less uncertainty).\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec15\" class=\"Section2\"\u003e\u003ch2\u003e3.4 Spatial patterns and heatmap evidence of clustering\u003c/h2\u003e\u003cp\u003eThe quantitative timing data are supported by trajectory panels and spatial density heatmaps. As chemotactic guidance increases, there is pronounced, high-density clustering around the egg, whereas low-density, diffuse sperm distributions are seen in heatmaps calculated at representative time points (refer to the heatmap panels for early/mid/late steps embedded within Figs.\u0026nbsp;6\u0026ndash;10). Strong chemotaxis produces a localized high-density plume that converges on the egg, according to the heatmaps, which is consistent with efficient gradient-following behavior and quick capture. Additionally, directness measurements are visibly displayed in trajectory plots (Figs.\u0026nbsp;9\u0026ndash;13): sperm tracks are more straight and orientated toward the ovum at chemotaxis 0.9, whereas they are meandering and frequently fail to pass within the fertilization radius at chemotaxis 0.0.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec16\" class=\"Section2\"\u003e\u003ch2\u003e3.5 Survival analysis: Kaplan\u0026ndash;Meier curves and log-rank tests\u003c/h2\u003e\u003cp\u003eWe created Kaplan\u0026ndash;Meier survival curves stratified by chemotaxis level in order to directly examine time-to-event dynamics while taking into consideration censored runs, or runs with no fertilization during the maximal observation window (Fig.\u0026nbsp;7). For high chemotaxis conditions (0.7, 0.9), survival curves show a quick reduction (i.e., a rapid transition to the fertilized state); for low chemotaxis conditions (0.0), declines are significantly slower or incomplete, indicating partial or absent fertilization in some repeats. For almost all chemotaxis comparisons, pairwise log-rank tests on these survival functions showed very significant differences. Specifically, there was a significant difference (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001) between low chemotaxis (\u0026le;\u0026thinsp;0.3) and high chemotaxis (\u0026ge;\u0026thinsp;0.7). These findings suggest that the hazard (instantaneous likelihood) of fertilization as a function of time is significantly and statistically robustly influenced by chemotactic sensitivity.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec17\" class=\"Section2\"\u003e\u003ch2\u003e3.6 Distributional comparisons: KS tests and ANOVA\u003c/h2\u003e\u003cp\u003eOn raw fertilization-time distributions, we performed further parametric and nonparametric comparisons (where appropriate, treating censored observations as maxima):\u003c/p\u003e\u003cp\u003e\u003cul\u003e\u003cli\u003e\u003cp\u003eKolmogorov-Smirnov (KS) Test: The survival analyses were substantially supported by pairwise KS tests, which revealed notable distributional differences between the groups with low and high chemotaxis. Significantly, comparisons between chemotaxis 0.3 and higher values (0.5, 0.7, and 0.9) revealed statistically significant p-values (0.3 vs. 0.9, for instance, p\u0026thinsp;=\u0026thinsp;0.0040), suggesting that the empirical CDFs shifted toward earlier periods with stronger chemotaxis. The non-monotonic and variance-sensitive nature of the distributions was reflected in a few paired comparisons that revealed less evidence of difference (for complete pairwise KS p-values, see Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eANOVA: The results of a one-way ANOVA spanning all chemotactic groups showed a highly significant difference in mean fertilization times (p\u0026thinsp;=\u0026thinsp;0.0001). The survival (log-rank) and KS analyses offer crucial supplemental nonparametric confirmation because ANOVA ignores censoring and non-normality.\u003c/p\u003e\u003c/li\u003e\u003c/ul\u003e\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec18\" class=\"Section2\"\u003e\u003ch2\u003e3.7 Pairwise statistical summary\u003c/h2\u003e\u003cp\u003eTable\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003econsolidates the pairwise statistical comparisons between chemotactic strength groups. The table highlights the strongest contrasts (low vs. high chemotaxis) and identifies cases where differences are less pronounced (e.g., intermediate vs. intermediate comparisons).\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"3\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eComparison\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003eLog-rank p-value\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u003cp\u003eKS p-value\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.0 vs 0.3\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e\u0026lt;\u0026thinsp;0.001\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.0 vs 0.5\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e\u0026lt;\u0026thinsp;0.001\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.0 vs 0.7\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e\u0026lt;\u0026thinsp;0.001\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.0 vs 0.9\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e\u0026lt;\u0026thinsp;0.001\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.3 vs 0.5\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e\u0026lt;\u0026thinsp;0.001\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.3 vs 0.7\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e\u0026lt;\u0026thinsp;0.001\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.3 vs 0.9\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e\u0026lt;\u0026thinsp;0.001\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.5 vs 0.7\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e0.0015\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.9831 (ns)\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.5 vs 0.9\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e\u0026lt;\u0026thinsp;0.001\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e0.7 vs 0.9\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e0.0104\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.0386\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003e(ns\u0026thinsp;=\u0026thinsp;not significant at α\u0026thinsp;=\u0026thinsp;0.05). Note: log-rank p-values are reported to 4 decimal places when available; extremely small p-values are shown as \u0026lt;\u0026thinsp;0.001.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec19\" class=\"Section2\"\u003e\u003ch2\u003e3.8 Representative single-run outcomes and extremes\u003c/h2\u003e\u003cp\u003eIndividual representative runs are highlighted to aid in comprehension: in one representative replicates series, the egg was fertilized at step 36 (baseline run); in others, the fertilization timestep varied depending on the chemotaxis (for example, chemotaxis 0.0 fertilized at step 18 in a specific instance; chemotaxis 0.3 at step 3; chemotaxis 0.5 at step 0; chemotaxis 0.7 at step 5; chemotaxis 0.9 at step 3). These single-run examples are consistent with the overall trend: strong guidance concentrates good outcomes early, whereas random or weak guidance causes wide variety, including late or unsuccessful events.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec20\" class=\"Section2\"\u003e\u003ch2\u003e3.9 Integrative summary of success vs. speed trade-offs\u003c/h2\u003e\u003cp\u003eTwo main effects of chemotaxis are revealed by the combined set of analyses (heatmaps, trajectory plots, descriptive statistics, survival functions, and pairwise tests):\u003c/p\u003e\u003cp\u003e\u003cul\u003e\u003cli\u003e\u003cp\u003eThe likelihood that a sperm will reach the fertilization radius within the observation window (success rate) increases with increasing chemotaxis. The change from 0.0 to 0.3 is very noticeable because our replicates' success rates increase from 40\u0026ndash;100%.\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eTemporal effect: Increasing chemotaxis decreases timing variance and speeds up the time to fertilization. Speed increases between 0.1 and 0.5 are the largest; increases above 0.7 result in modest mean decreases but restrict the timing distribution even further.\u003c/p\u003e\u003c/li\u003e\u003c/ul\u003e\u003c/p\u003e\u003cp\u003eThese combined results imply that chemotactic signaling decreases competitive latency among spermatozoa and boosts encounter efficiency, two properties that may be significant in physiological and ART contexts.\u003c/p\u003e\u003cp\u003e\u003cspan type=\"ItalicUnderline\" class=\"ItalicUnderline\" name=\"Emphasis\"\u003e3.10 Figure references (PDF order)\u003c/span\u003e\u003c/p\u003e\u003cp\u003e\u003cul\u003e\u003cli\u003e\u003cp\u003eFigure 1: Number of alive sperms over time (temporal attrition).\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eFigure 2: Example sample sperm trajectories (heterogeneous paths).\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eFigure 3: Mean fertilization success rate (summary).\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eFigure 4: Fertilization success vs chemotactic strength (trend).\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eFigure 5: Mean fertilization time summary (overview).\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eFigure 6: Average fertilization time vs chemotactic strength.\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eFigure 7: Kaplan\u0026ndash;Meier survival curves for fertilization time by chemotactic strength.\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eFigure 8: Fertilization time distributions by chemotactic strength.\u003c/p\u003e\u003c/li\u003e\u003cli\u003e\u003cp\u003eFigure 9\u0026ndash;13: Representative fertilized trajectories for chemotaxis\u0026thinsp;=\u0026thinsp;0.9, 0.7, 0.5, 0.3, 0.0\u003c/p\u003e\u003c/li\u003e\u003c/ul\u003e\u003c/p\u003e\u003c/div\u003e"},{"header":"4. Discussion","content":"\u003cp\u003eThrough a thorough computational framework, this study shows that chemotactic intensity has a significant and measurable impact on the temporal distribution of sperm\u0026ndash;egg interactions as well as fertilization success rates. Chemotactic intensities\u0026thinsp;\u0026ge;\u0026thinsp;0.3 regularly result in nearly universal fertilization success and much shorter fertilization periods, according to the simulation results, which clearly demonstrate a threshold phenomena. Fertilization results below this threshold are extremely random and vary greatly from run to run, which is consistent with experimental findings in situations where chemical guiding cues are either missing or interfered with [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e, \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e]. In line with previous gradient-sensing models in motile cell populations, the observed plateau in time reduction at higher chemotactic strengths (\u0026ge;\u0026thinsp;0.7) implies that more bias offers no advantage once directional navigation surpasses a particular fidelity [\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e]. Stronger chemotaxis not only increases targeting precision but also lessens arrival time dispersion, thereby coordinating fertilization events, as further demonstrated by spatial trajectory graphs. Such synchronization could improve reproductive efficiency in biological contexts when ovum receptivity windows are aligned or coordinated timing inhibits polyspermy [\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e]. Crucially, our distributional tests and survival analysis verify that these effects are statistically significant changes in fertilization dynamics under all examined settings rather than being marginal. The Kaplan-Meier curves quantitatively demonstrate that chemotaxis serves as a temporal accelerator of reproductive success by demonstrating enhanced hazard rates, or the immediate likelihood of conception, in high-chemotaxis groups. In order to better match in vivo reproductive conditions, the model additionally takes physiological limitations like pH, temperature, and fluid flow into account. According to empirical research on the viability of human and animal sperm under stress, motility and survival were dramatically decreased by deviations from physiological pH or normothermia [\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e]. The inclusion of laminar flow effects illustrates how hydrodynamic conditions can either help upstream navigation (rheotaxis) or impede progress, depending on orientation, even though fluid dynamics was not the main focus. This phenomenon has been directly visualized in microfluidic and in vivo oviduct models [\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e, \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e]. These results highlight the selection advantage of sperm populations with effective chemotactic navigation from an evolutionary perspective. By quickly converging on the fertilization site, these sperm could reduce the consequences of temporal rivalry in natural mating situations, in addition to outcompeting less responsive counterparts for ovum access. The evolutionary conservation of chemotactic signaling pathways across a variety of taxa may be partially explained by this twofold advantage, which consists of a higher chance of success and a shorter time to fertilization [\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e]. Methodologically, a multi-factorial investigation of reproductive dynamics is made possible by the combination of stochastic mortality, energy depletion, and chemotactic bias into a single computer system, eliminating the expense and unpredictability of in vivo testing. Confidence in these computational predictions is further reinforced by the statistical convergence of ANOVA, KS tests, and survival analysis. The current framework, however, continues to simplify the actual biological environment. Direct translation to all physiological situations is limited by the absence of precise three-dimensional oviduct geometry, mucus rheology, capacitation-dependent motility transitions, and zona pellucida binding dynamics. By combining fluid-structure interactions with biochemical signaling networks that control hyperactivation and sperm-egg binding kinetics, further research could broaden the model [\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e, \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]. Modifying chemotactic gradients in sperm selection devices may enhance timing efficiency and selection accuracy, which has obvious translational implications for assisted reproductive technologies (ART). Furthermore, by determining chemotaxis inhibition thresholds that successfully prevent fertilization without compromising baseline motility, the sensitivity mapping described here may help guide research on contraceptives. Our findings offer a quantitative basis for creating techniques that either improve or repress sperm guidance pathways, which are becoming more and more acknowledged as promising targets for pharmaceutical intervention [\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e]. To sum up, this study provides quantitative evidence that chemotactic signaling is a key factor in fertilization efficiency. We offer both mechanistic knowledge and practical possibilities by defining distinct performance thresholds and recording the shift from random to synchronized fertilization events. These findings pave the way for both fundamental study and innovative therapeutic applications by highlighting the significance of combining chemotactic dynamics with environmental limitations in computational reproductive biology.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eEisenbach M, Giojalas LC (2006) Sperm guidance in mammals \u0026mdash; an unpaved road to the egg. Nat Rev Mol Cell Biol 7(4):276\u0026ndash;285. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/nrm1894\u003c/span\u003e\u003cspan address=\"10.1038/nrm1894\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSuarez SS, Pacey AA (2006) Sperm transport in the female reproductive tract. Hum Reprod Update 12(1):23\u0026ndash;37. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/humupd/dmi047\u003c/span\u003e\u003cspan address=\"10.1093/humupd/dmi047\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eKaupp UB, Hildebrand E, Weyand I (2012) Sperm Chemotaxis. Cold Spring Harb Perspect Biol 4(4):a005061. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1101/cshperspect.a005061\u003c/span\u003e\u003cspan address=\"10.1101/cshperspect.a005061\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eRiffell JR, Zimmer RK (2007) Sex and flow: the consequences of fluid shear for sperm-egg interactions. J Exp Biol 210(7):3644\u0026ndash;3660. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1242/jeb.008528\u003c/span\u003e\u003cspan address=\"10.1242/jeb.008528\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eWoolley DM (2003) Motility of spermatozoa at surfaces. Reproduction 126(2):259\u0026ndash;270. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1530/rep.0.1260259\u003c/span\u003e\u003cspan address=\"10.1530/rep.0.1260259\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSmith DJ, Gaffney EA, Blake JR (2008) Modeling mucociliary clearance. Respir Physiol Neurobiol 163(1\u0026ndash;3):178\u0026ndash;188. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.resp.2008.05.013\u003c/span\u003e\u003cspan address=\"10.1016/j.resp.2008.05.013\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eZaferani M, Mohammadi S, Abbaspourrad A (2021) Collective behavior of sperm cells. Prog Biophys Mol Biol 160:13\u0026ndash;23. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.pbiomolbio.2020.09.005\u003c/span\u003e\u003cspan address=\"10.1016/j.pbiomolbio.2020.09.005\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eYaghoobi M et al (2020) Computational modeling of spermatozoa swimming near a surface. Math Biosci Eng 17(1):527\u0026ndash;553. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3934/mbe.2020026\u003c/span\u003e\u003cspan address=\"10.3934/mbe.2020026\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSuarez SS (2008) Regulation of sperm motility and hyperactivation. Hum Reprod Update 14(6):647\u0026ndash;657. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/humupd/dmn042\u003c/span\u003e\u003cspan address=\"10.1093/humupd/dmn042\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eKantsler V, Dunkel J, Polin M, Goldstein RE (2013) Ciliary contact interactions dominate surface scattering of swimming eukaryotes. Proc Natl Acad Sci U S A 110(4):1187\u0026ndash;1192. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1073/pnas.1218517110\u003c/span\u003e\u003cspan address=\"10.1073/pnas.1218517110\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eZaferani M et al (2019) Surface interactions explain the hydrodynamics of swimming spermatozoa near surfaces. Phys Rev Lett 122(21):218102. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1103/PhysRevLett.122.218102\u003c/span\u003e\u003cspan address=\"10.1103/PhysRevLett.122.218102\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eBaker MA et al (2019) Beyond motility: sperm DNA integrity and fertility. Reproduction 157(5):R135\u0026ndash;R146. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1530/REP-18-0391\u003c/span\u003e\u003cspan address=\"10.1530/REP-18-0391\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eMonteiro AV et al (2021) Effect of seminal plasma on sperm survival and fertilization success. Front Physiol 12:723563. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3389/fphys.2021.723563\u003c/span\u003e\u003cspan address=\"10.3389/fphys.2021.723563\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eTung CK et al (2017) Fluid viscoelasticity promotes collective swimming of sperm. Sci Adv 3(10):e1603119. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1126/sciadv.1603119\u003c/span\u003e\u003cspan address=\"10.1126/sciadv.1603119\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eTherneau TM, Grambsch PM (2000) Modeling survival data: extending the Cox model. Springer\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eKlein JP, Moeschberger ML (2003) Survival Analysis: Techniques for Censored and Truncated Data. Springer\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eGoto M et al (2020) Influence of oxidative stress on sperm motility and fertilization. Reprod Med Biol 19(1):3\u0026ndash;9. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1002/rmb2.12303\u003c/span\u003e\u003cspan address=\"10.1002/rmb2.12303\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003ePacey AA (2020) Sperm motility: new perspectives and developments. Hum Fertil (Camb) 23(3):177\u0026ndash;186. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1080/14647273.2020.1736256\u003c/span\u003e\u003cspan address=\"10.1080/14647273.2020.1736256\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eBaker MA, Aitken RJ (2021) Oxidative stress in the male germ line. Biol Reprod 104(6):1111\u0026ndash;1121. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/biolre/ioab010\u003c/span\u003e\u003cspan address=\"10.1093/biolre/ioab010\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSmith DJ et al (2021) The role of temperature and pH in sperm motility. J Exp Biol 224(Pt 6):jeb239466. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1242/jeb.239466\u003c/span\u003e\u003cspan address=\"10.1242/jeb.239466\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003ePublicover S et al (2021) Sperm signaling pathways: targets for contraception. Trends Pharmacol Sci 42(11):967\u0026ndash;978. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.tips.2021.08.006\u003c/span\u003e\u003cspan address=\"10.1016/j.tips.2021.08.006\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eRobinson J et al (2022) Novel contraceptive approaches targeting sperm motility. Contraception 105:10\u0026ndash;18. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.contraception.2022.02.0034\u003c/span\u003e\u003cspan address=\"10.1016/j.contraception.2022.02.0034\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eHarris CR, Millman KJ, van der Walt SJ et al (2020) Array programming with NumPy. Nature 585(7825):357\u0026ndash;362. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/s41586-020-2649-2\u003c/span\u003e\u003cspan address=\"10.1038/s41586-020-2649-2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eHunter JD, Matplotlib (2007) A 2D graphics environment. Comput Sci Eng 9(3):90\u0026ndash;95. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1109/MCSE.2007.55\u003c/span\u003e\u003cspan address=\"10.1109/MCSE.2007.55\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eWaskom ML (2021) seaborn: statistical data visualization. J Open Source Softw 6(60):3021. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.21105/joss.03021\u003c/span\u003e\u003cspan address=\"10.21105/joss.03021\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eMcKinney W (2010) Data structures for statistical computing in Python. Proc 9th Python Sci Conf 445:51\u0026ndash;56\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eDavidson-Pilon C et al (2019) lifelines: survival analysis in Python. J Open Source Softw 4(40):1317. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.21105/joss.01317\u003c/span\u003e\u003cspan address=\"10.21105/joss.01317\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eVirtanen P, Gommers R, Oliphant TE et al (2020) SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat Methods 17(3):261\u0026ndash;272. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/s41592-019-0686-2\u003c/span\u003e\u003cspan address=\"10.1038/s41592-019-0686-2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eRiffell JR, Zimmer RK (2007) Sex and flow: the consequences of fluid shear for sperm-egg interactions. J Exp Biol 210(7):3644\u0026ndash;3660. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1242/jeb.008528\u003c/span\u003e\u003cspan address=\"10.1242/jeb.008528\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSmith DJ et al (2021) The role of temperature and pH in sperm motility. J Exp Biol 224(Pt 6):jeb239466. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1242/jeb.239466\u003c/span\u003e\u003cspan address=\"10.1242/jeb.239466\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003ePacey AA (2020) Sperm motility: new perspectives and developments. Hum Fertil (Camb) 23(3):177\u0026ndash;186. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1080/14647273.2020.1736256\u003c/span\u003e\u003cspan address=\"10.1080/14647273.2020.1736256\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eTherneau TM, Grambsch PM (2000) Modeling survival data: extending the Cox model. Springer\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eGoto M et al (2020) Influence of oxidative stress on sperm motility and fertilization. Reprod Med Biol 19(1):3\u0026ndash;9. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1002/rmb2.12303\u003c/span\u003e\u003cspan address=\"10.1002/rmb2.12303\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eZaferani M, Mohammadi S, Abbaspourrad A (2021) Collective behavior of sperm cells. Prog Biophys Mol Biol 160:13\u0026ndash;23. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.pbiomolbio.2020.09.005\u003c/span\u003e\u003cspan address=\"10.1016/j.pbiomolbio.2020.09.005\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eYoshida M, Yoshida K (2011) Sperm chemotaxis and regulation of flagellar movement by Ca\u0026sup2;⁺. Mol Hum Reprod 17(8):457\u0026ndash;465. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/molehr/gar041\u003c/span\u003e\u003cspan address=\"10.1093/molehr/gar041\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eDenissenko P, Kantsler V, Smith DJ, Kirkman-Brown J (2012) Human spermatozoa migration in microchannels reveals boundary-following navigation. Proc Natl Acad Sci U S A 109(21):8007\u0026ndash;8010. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1073/pnas.1202934109\u003c/span\u003e\u003cspan address=\"10.1073/pnas.1202934109\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eFriedrich BM, J\u0026uuml;licher F (2007) Chemotaxis of sperm cells. Proc Natl Acad Sci U S A 104(33):13256\u0026ndash;13261. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1073/pnas.0703530104\u003c/span\u003e\u003cspan address=\"10.1073/pnas.0703530104\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eEvans JP (2012) Sperm\u0026ndash;egg interaction. Annu Rev Physiol 74:477\u0026ndash;502. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1146/annurev-physiol-020911-153339\u003c/span\u003e\u003cspan address=\"10.1146/annurev-physiol-020911-153339\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eZhou Y et al (2015) Effects of pH on human sperm motility and survival. J Assist Reprod Genet 32(6):809\u0026ndash;815. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s10815-015-0471-5\u003c/span\u003e\u003cspan address=\"10.1007/s10815-015-0471-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSmith DJ et al (2014) Rheotaxis facilitates upstream navigation of mammalian sperm cells. eLife 3:e02403. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.7554/eLife.02403\u003c/span\u003e\u003cspan address=\"10.7554/eLife.02403\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eTung CK et al (2015) Emergence of upstream swimming via a hydrodynamic transition. Phys Rev Lett 114(10):108102. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1103/PhysRevLett.114.108102\u003c/span\u003e\u003cspan address=\"10.1103/PhysRevLett.114.108102\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eKantsler V, Dunkel J, Polin M, Goldstein RE (2014) Rheotaxis of swimming eukaryotes. Proc Natl Acad Sci U S A 111(52):15970\u0026ndash;15975. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1073/pnas.1413777111\u003c/span\u003e\u003cspan address=\"10.1073/pnas.1413777111\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eAlvarez L, Friedrich BM, Gompper G, Kaupp UB (2014) The computational sperm cell. Trends Cell Biol 24(3):198\u0026ndash;207. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.tcb.2013.10.004\u003c/span\u003e\u003cspan address=\"10.1016/j.tcb.2013.10.004\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSuarez SS (2008) Hyperactivation in sperm. J Androl 29(3):334\u0026ndash;342. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.2164/jandrol.107.002576\u003c/span\u003e\u003cspan address=\"10.2164/jandrol.107.002576\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eChang H, Suarez SS (2010) Unexpected flagellar movement patterns and epithelial binding behaviour of mouse sperm in the oviduct. Biol Reprod 83(1):103\u0026ndash;110. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1095/biolreprod.110.084624\u003c/span\u003e\u003cspan address=\"10.1095/biolreprod.110.084624\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eWachten D et al (2017) Signal transduction in mammalian sperm. Annu Rev Physiol 79:189\u0026ndash;208. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1146/annurev-physiol-022516-034238\u003c/span\u003e\u003cspan address=\"10.1146/annurev-physiol-022516-034238\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Sperm chemotaxis, Fertilization dynamics, Computational modeling, Survival analysis, Environmental heterogeneity, Sperm motility, Kaplan–Meier analysis","lastPublishedDoi":"10.21203/rs.3.rs-7427895/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-7427895/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eIn mammalian reproduction, the ability of sperm cells to move across diverse physiological conditions and reach the oocyte is crucial to the success of fertilization. In this work, we create a thorough stochastic computer model that simulates the dynamics of sperm chemotaxis in a limited two-dimensional microenvironment. This model takes into account important biological factors such as temperature, pH fluctuations, ambient flow, and chemotactic strength. Through comprehensive parameter sweeps and duplicate simulations, we measure fertilization success rates and timing distributions under various situations in order to systematically investigate the impact of these parameters on fertilization efficiency. Our model incorporates stochastic sperm mortality to represent chemical and immunological challenges, and environmental barriers to replicate physical and metabolic heterogeneities found in the female reproductive system. Using survival analytic frameworks such as log-rank tests and Kaplan-Meier curves, fertilization time distributions were examined. The results showed statistically significant variations in fertilization kinetics among chemotactic regimes and environmental variables. Furthermore, sperm density heatmaps emphasize the crucial role that directed motility plays in fertilization outcomes by highlighting spatial clustering dynamics that are influenced by external fluxes and the intensity of chemotaxis. The robustness of the observed effects is confirmed by statistical comparisons using the ANOVA and Kolmogorov-Smirnov tests. Our results give a prediction framework for comprehending sperm behavior in vivo and quantitative insights into the relative contributions of biophysical and biochemical elements influencing fertilization success. By clarifying the mechanics behind sperm navigation and egg encounter efficiency, this integrative modeling method paves the way for future experimental validation and could influence assisted reproductive technologies and fertility therapies.\u003c/p\u003e","manuscriptTitle":"Stochastic Modeling and Environmental Parameterization of Sperm Chemotaxis Dynamics Reveal Critical Determinants of Fertilization Success Under Physiological Constraints","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-08-29 14:03:16","doi":"10.21203/rs.3.rs-7427895/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"34560ef6-8af1-4919-8b71-fdb8bc850648","owner":[],"postedDate":"August 29th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":53523528,"name":"Developmental Biology"},{"id":53523529,"name":"Bioinformatics"},{"id":53523530,"name":"Biotechnology and Bioengineering"},{"id":53523531,"name":"Computational Biology"}],"tags":[],"updatedAt":"2025-08-29T14:03:16+00:00","versionOfRecord":[],"versionCreatedAt":"2025-08-29 14:03:16","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-7427895","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-7427895","identity":"rs-7427895","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}