Mathematical medicine and biology : a journal of the IMA最新文献

Stochastically perturbed models, where the white noise type stochastic perturbations are proportional to the current system state, the most realistically describe real-life biosystems. However, such models essentially have no equilibrium states apart from one at the origin. This feature makes analysis of such models extremely difficult. Probably, the best result that can be found for such models is finding of accurate estimations of a region in the model phase space that serves as an attractor for model trajectories. In this paper, we consider a classical stochastically perturbed Lotka-Volterra model of competing or symbiotic populations, where the white noise type perturbations are proportional to the current system state. Using the direct Lyapunov method in a combination with a recently developed technique, we establish global asymptotic properties of this model. In order to do this, we, firstly, construct a Lyapunov function that is applicable to the both competing (and globally stable) and symbiotic deterministic Lotka-Volterra models. Then, applying this Lyapunov function to the stochastically perturbed model, we show that solutions with positive initial conditions converge to a certain compact region in the model phase space and oscillate around this region thereafter. The direct Lyapunov method allows to find estimates for this region. We also show that if the magnitude of the noise exceeds a certain critical level, then some or all species extinct via process of the stochastic stabilization ('stabilization by noise'). The approach applied in this paper allows to obtain necessary conditions for the extinction. Sufficient conditions for the extinction (that for this model occurs via the process that is known as the 'stochastic stabilization', or the 'stabilization by noise') are found applying the Khasminskii-type Lyapunov functions.

Drug resistance is a significant obstacle to effective cancer treatment. To gain insights into how drug resistance develops, we adopted a concept called fitness landscape and employed a phenotype-structured population model by fitting to a set of experimental data on a drug used for ovarian cancer, olaparib. Our modeling approach allowed us to understand how a drug affects the fitness landscape and track the evolution of a population of cancer cells structured with a spectrum of drug resistance. We also incorporated pharmacokinetic (PK) modeling to identify the optimal dosages of the drug that could lead to long-term tumor reduction. We derived a formula that indicates that maximizing variation in plasma drug concentration over a dosing interval could be important in reducing drug resistance. Our findings suggest that it may be possible to achieve better treatment outcomes with a drug dose lower than the levels recommended by the drug label. Acknowledging the current limitations of our work, we believe that our approach, which combines modeling of both PK and drug resistance evolution, could contribute to a new direction for better designing drug treatment regimens to improve cancer treatment.

Altay Prefecture, a typical arid region in northwestern China, has experienced the climate transition from warming-drying to warming-wetting since 1980s and has attracted widespread attention. Nonetheless, it is still unclear how climate change has influenced the distribution of vegetation in this region. In this paper, a reaction-diffusion model of the climate-vegetation system is proposed to study the impact of climate change (precipitation, temperature and carbon dioxide concentration) on vegetation patterns in Altay Prefecture. Our results indicate that the tendency of vegetation growth in Altay Prefecture improved gradually from 1985 to 2010. Under the current climate conditions, the increase of precipitation results in the change of vegetation pattern structures, and eventually vegetation coverage tends to be uniform. Moreover, we found that there exists an optimal temperature where the spot vegetation pattern structure remains stable. Furthermore, the increase in carbon dioxide concentration induces vegetation pattern transition. Based on four climate change scenarios of the Coupled Model Intercomparison Project Phase 6 (CMIP6), we used the power law range (PLR) to predict the optimal scenario for the sustainable development of the vegetation ecosystem in Altay Prefecture.

The emergence of multiple strains of SARS-COV-2 has made it complicated to predict and control the COVID-19 pandemic. Although some vaccines have been effective in reducing the severity of the disease, these vaccines are designed for a specific strain of the virus and are usually less effective for other strains. In addition, the waning of vaccine-induced immunity, reinfection of recovered people, and incomplete vaccination are challenging to the vaccination program. In this study, we developed a detailed model to describe the multi-strain transmission dynamics of COVID-19 under vaccination. We implemented our model to examine the impact of inter-strain transmission competition under vaccination on the critical outbreak indicators: hospitalized cases, undiagnosed cases, basic reproduction numbers, and the overtake-time by a new strain to the existing strain. In particular, our results on the dependence of the overtake-time on vaccination rates, progression-to-infectious rate, and relative transmission rates provide helpful information for managing a pandemic with circulating two strains. Furthermore, our results suggest that a reduction in the relative transmission rates and a decrease in vaccination dropout rates or an increase in vaccination rates help keep the reproduction number of both strains below unity and keep the number of hospitalized cases and undiagnosed cases at their lowest levels. Moreover, our analysis shows that the second and booster-dose vaccinations are useful for further reducing the reproduction number.

Structure-function relationships occur throughout the sciences. Motivated by optimization of such systems, we develop a framework for estimating the input modes from the singular value decomposition from the action of the forward operator alone. These can then be used to determine the input (structure) changes, which induce the largest output (function) changes. The accuracy of the estimate is determined by reference to the method of snapshots. The proposed method is demonstrated on several example problems, and finally used to approximate the optimal airway treatment set for a problem in respiratory physiology.

We employ the multiphase, moving boundary model of Byrne et al. (2003, Appl. Math. Lett., 16, 567-573) that describes the evolution of a motile, viscous tumour cell phase and an inviscid extracellular liquid phase. This model comprises two partial differential equations that govern the cell volume fraction and the cell velocity, together with a moving boundary condition for the tumour edge, and here we characterize and analyse its travelling-wave and pattern-forming behaviour. Numerical simulations of the model indicate that patterned solutions can be obtained, which correspond to multiple regions of high cell density separated by regions of low cell density. In other parameter regimes, solutions of the model can develop into a forward- or backward-moving travelling wave, corresponding to tumour growth or extinction, respectively. A travelling-wave analysis allows us to find the corresponding wave speed, as well as criteria for the growth or extinction of the tumour. Furthermore, a stability analysis of these travelling-wave solutions provides us with criteria for the occurrence of patterned solutions. Finally, we discuss how the initial cell distribution, as well as parameters related to cellular motion and cell-liquid drag, control the qualitative features of patterned solutions.