Research article Special Issues

Chiral selection in wrapping, crossover, and braiding of DNA mediated by asymmetric bend-writhe elasticity

  • Received: 28 August 2015 Accepted: 29 October 2015 Published: 06 November 2015
  • Wrapping, crossover, and braiding of DNA are the motifs of fundamental interest in genome packaging, gene regulation, and enzyme recognition. This study explores elastic mechanisms for the selection of chirality in wrapping, crossover, and braiding of DNA based on a coarse-grained model. The DNA model consists of two elastic chains that mutually intertwine in a right-handed manner forming a double-stranded helix with the distinction between major and minor grooves. Although individual potential energy functions of the DNA model have no asymmetry in terms of left and right twist, the model as a whole exhibits an asymmetric propensity to writhe in the left direction upon bending due to the right-handed helical geometry. Monte Carlo simulations of this model suggest that DNA has a propensity to prefer left-handed wrapping around a spherical core particle and also around a uniform rod due to the asymmetric elastic coupling between bending and writhing. This result indicates an elastic origin of the uniform left-handed wrapping of DNA in nucleosomes and also has implications on the wrapping of double-stranded DNA around rod-like molecules. Monte Carlo simulations of the DNA model also suggest that two juxtaposed DNA molecules can braid each other spontaneously under moderate attractive interactions with the preference for left-handed braiding due to the asymmetric coupling between bending and writhing. This result suggests the importance of asymmetric elasticity in the selection of chirality in braiding of a pair of DNA molecules.

    Citation: Tomohiro Yanao, Sosuke Sano, Kenichi Yoshikawa. Chiral selection in wrapping, crossover, and braiding of DNA mediated by asymmetric bend-writhe elasticity[J]. AIMS Biophysics, 2015, 2(4): 666-694. doi: 10.3934/biophy.2015.4.666

    Related Papers:

    [1] Chen Liang, Hai-Feng Huo, Hong Xiang . Modelling mosquito population suppression based on competition system with strong and weak Allee effect. Mathematical Biosciences and Engineering, 2024, 21(4): 5227-5249. doi: 10.3934/mbe.2024231
    [2] Yang Li, Jia Li . Stage-structured discrete-time models for interacting wild and sterile mosquitoes with beverton-holt survivability. Mathematical Biosciences and Engineering, 2019, 16(2): 572-602. doi: 10.3934/mbe.2019028
    [3] Rocio Caja Rivera, Shakir Bilal, Edwin Michael . The relation between host competence and vector-feeding preference in a multi-host model: Chagas and Cutaneous Leishmaniasis. Mathematical Biosciences and Engineering, 2020, 17(5): 5561-5583. doi: 10.3934/mbe.2020299
    [4] Maria Guadalupe Vazquez-Peña, Cruz Vargas-De-León, Jorge Velázquez-Castro . Global stability for a mosquito-borne disease model with continuous-time age structure in the susceptible and relapsed host classes. Mathematical Biosciences and Engineering, 2024, 21(11): 7582-7600. doi: 10.3934/mbe.2024333
    [5] Hongyong Zhao, Yangyang Shi, Xuebing Zhang . Dynamic analysis of a malaria reaction-diffusion model with periodic delays and vector bias. Mathematical Biosciences and Engineering, 2022, 19(3): 2538-2574. doi: 10.3934/mbe.2022117
    [6] Yan-Xia Dang, Zhi-Peng Qiu, Xue-Zhi Li, Maia Martcheva . Global dynamics of a vector-host epidemic model with age of infection. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1159-1186. doi: 10.3934/mbe.2017060
    [7] Mayra Núñez-López, Jocelyn A. Castro-Echeverría, Jorge X. Velasco-Hernández . Dynamic interaction between transmission, within-host dynamics and mosquito density. Mathematical Biosciences and Engineering, 2025, 22(6): 1364-1381. doi: 10.3934/mbe.2025051
    [8] Liming Cai, Shangbing Ai, Guihong Fan . Dynamics of delayed mosquitoes populations models with two different strategies of releasing sterile mosquitoes. Mathematical Biosciences and Engineering, 2018, 15(5): 1181-1202. doi: 10.3934/mbe.2018054
    [9] Karen R. Ríos-Soto, Baojun Song, Carlos Castillo-Chavez . Epidemic spread of influenza viruses: The impact of transient populations on disease dynamics. Mathematical Biosciences and Engineering, 2011, 8(1): 199-222. doi: 10.3934/mbe.2011.8.199
    [10] Bruce Pell, Amy E. Kendig, Elizabeth T. Borer, Yang Kuang . Modeling nutrient and disease dynamics in a plant-pathogen system. Mathematical Biosciences and Engineering, 2019, 16(1): 234-264. doi: 10.3934/mbe.2019013
  • Wrapping, crossover, and braiding of DNA are the motifs of fundamental interest in genome packaging, gene regulation, and enzyme recognition. This study explores elastic mechanisms for the selection of chirality in wrapping, crossover, and braiding of DNA based on a coarse-grained model. The DNA model consists of two elastic chains that mutually intertwine in a right-handed manner forming a double-stranded helix with the distinction between major and minor grooves. Although individual potential energy functions of the DNA model have no asymmetry in terms of left and right twist, the model as a whole exhibits an asymmetric propensity to writhe in the left direction upon bending due to the right-handed helical geometry. Monte Carlo simulations of this model suggest that DNA has a propensity to prefer left-handed wrapping around a spherical core particle and also around a uniform rod due to the asymmetric elastic coupling between bending and writhing. This result indicates an elastic origin of the uniform left-handed wrapping of DNA in nucleosomes and also has implications on the wrapping of double-stranded DNA around rod-like molecules. Monte Carlo simulations of the DNA model also suggest that two juxtaposed DNA molecules can braid each other spontaneously under moderate attractive interactions with the preference for left-handed braiding due to the asymmetric coupling between bending and writhing. This result suggests the importance of asymmetric elasticity in the selection of chirality in braiding of a pair of DNA molecules.


    Japanese encephalitis viral disease (JE) was first documented in 1871 in Japan [1]. Japanese encephalitis (JE) is the main cause of viral encephalitis in many countries of Asia and the western Pacific. Twenty four countries in the WHO South-East Asia and Western Pacific regions have Japanese encephalitis virus (JEV) transmission risk, which includes more than 3 billion people [1,2]. JEV is transmitted to humans through bites from infected mosquitoes of the Culex species (mainly Culex tritaeniorhynchus). The virus exists in a transmission cycle between mosquitoes, pigs, and/or water birds and is transmitted to humans by infected mosquito bite. JE is predominantly found in rural and periurban settings [3,4]. Each year there are nearly 68,000 clinical cases of JE globally, with approximately 13,600 to 20,400 deaths [1]. JE primarily affects children. Most adults in endemic countries have natural immunity after childhood infection, but individuals of any age may be affected. Most people infected with JE do not have symptoms or have only mild symptoms. However, a small percentage of infected people develop inflammation of the brain (encephalitis), with symptoms including sudden onset of headache, high fever, disorientation, coma, tremors, and convulsions. Approximately 1 in 250 infections results in severe clinical illness. Severe disease is characterized by rapid onset of high fever, headache, neck stiffness, disorientation, coma, seizures, spastic paralysis, and ultimately death [1,2]. The case fatality rate for the disease can be as high as 30% among those with disease symptoms. In addition, 20–30% of those who survive suffer permanent neuropsychiatric sequelae [3].

    Leptospirosis is a bacterial disease that affects humans and animals. Leptospirosis is considered to be the most widespread zoonotic disease in the world [5]. It is caused by bacteria of the genus Leptospira [6]. The bacteria are spread through the urine of infected animals, which can get into water or soil and can survive there for weeks to months. Humans can become infected through contact with urine (or other body fluids, except saliva) from infected animals; also, through contact with water, soil, or food contaminated with the urine of infected animals [7]. In humans, leptospirosis may occur in two phases, where about 10% of infected people move to the 2nd phase [5]. The first phase causes a wide range of symptoms, some of which may be mistaken for other diseases: fever, chills, headache, muscle aches, vomiting, or diarrhea. The second phase, if it occurs, is more severe; the person may have kidney damage, liver failure, meningitis, respiratory distress, and even death [6,8]. It is estimated that more than 1 million cases occur worldwide annually, including almost 60,000 deaths [8].

    In the last couple of decades, many ecological and field studies, along with a few mathematical ones, have been conducted to understand the dynamics of JEV transmission and to find some control measures to reduce JEV prevalence in humans. In 2001, a study used a probabilistic model of pathogen transmission to investigate various control measures for JEV transmission in humans [9]. Outcomes of the study show that a combination of control measures of similar effect (strategies to reduce vector population, or strategies to reduce vector-human interactions) is more effective compared to the combination of control measures of different effects (strategy to reduce vector population and strategy to reduce human-vector interaction). In 2014, Khan et al. studied the dynamics of JEV transmission in a pig population in northwest Bangladesh [10]. This study developed an SEIR model to understand transmission dynamics in pigs, and to estimate the potential impact of pig vaccination. Their results found that the prevalence of JE in pig populations can be reduced by up to 89% when 75% of susceptible pigs are vaccinated each year. The next year, in 2015, Lord et al. performed a study to rethink JEV transmission among hosts and vectors [11]. They suggested using a mathematical model parameterized with data to quantify the relative roles of potential species in JEV transmission.

    Very little research has been done, to the best of our knowledge, to understand the dynamics of leptospirosis in cattle and humans. In 2017, Chadsuthi et al. investigated the leptospirosis prevalence in livestock and humans in Thailand for 2010–2015 [12]. They tested humans, buffaloes, cattle, pigs, and analyzed collected data. Their analysis found livestock more susceptible to leptospirosis infection compared to humans. Later, in 2018, another study was done to understand the spread of leptospirosis in lambs in New Zealand. Here, researchers used a simple SI model to predict conditions under which the disease would persist in the lamb population [13]. Analysis of this study suggested that increasing the leptospira death rate in farms can reduce infection in livestock, and eventually in humans.

    There are many countries in Asia where both JE and leptospirosis are prevalent.Both diseases are endemic throughout southeast Asia from India to Japan and thus they are geographically correlated. Taking this reality into account, this study is developed to investigate the prevalence of both diseases in humans. Cattle contribute to leptospirosis infections in humans, acting as a source of leptospira. However, the presence of cattle in a domestic or peridomestic setting can be considered as an additional host for JE vectors, where humans act as the primary host. Hence, the presence of cattle increases host richness in the setting, which might eventually reduce the combined burden of these two diseases for humans. Earlier studies showed that an additional host (host richness) in a setting helps to reduce human infections of vector-borne diseases if certain conditions are maintained [14,15,16,17,18]. Johnson and Thieltges in 2010 identified that host diversity may help to reduce human infections depending on the relative abundance of additional host(s) relative to the focal host [14]. In 2013, Miller and Huppert proved that species diversity in host populations can amplify or can dilute disease prevalence depending on vectors' preference of host [15]. Recently in 2020, Zahid and Kribs showed that an additional host (dog, reservoir host) in a community can help in reducing the prevalence of visceral leishmaniasis in humans depending on the additional host's irritability to vector bites [17]. In another study, the same authors established that the presence of an additional host (chicken, incompetent host) in a domestic setting might help to reduce the prevalence of Chagas disease in humans if the distance between the two host populations (humans and chickens) is maintained at a certain range [18].

    The presence of cattle helps to reduce human-vector interactions by attracting mosquitoes towards them from humans [19]. Cattle have no contribution to the transmission of JE [20] whereas they act as a source of leptospirosis infections in humans. So, in terms of infections, the presence of cattle has two opposite influences on human health – it helps to reduce JE prevalence in humans and contributes to leptospirosis risk in humans. The goal of this study is to understand the dynamics of both diseases in each population, and eventually to understand if the presence of cattle in a setting, where JE and leptospirosis both are prevalent, is helpful to reduce the combined burden of JE and leptospirosis in humans. To answer this question we compare two different peridomestic settings: a setting involving cattle with humans, pigs, and mosquitoes, and the other setting involving humans, pigs, and mosquitoes without cattle. To understand disease dynamics, and to compare the proposed two settings, we develop SIR models for cattle, pig and human populations, and an SI model for the mosquito population. This study estimates the total number of $ DALY $s (Disability Adjusted Life Year) generated from JEV and (in the presence of cattle) leptospirosis infections for both settings and compare them to identify the better scenario. To the best of our knowledge, this is the first work ever to consider JE and leptospirosis infections together. More broadly, this extends the well-studied question in disease ecology of the impact of an additional host to a multi-pathogen context.

    The goal of this work is to understand the dynamics of Japanese encephalitis virus (JEV) and leptospirosis infections together and to estimate the impact of cattle presence on human infections of these diseases. Therefore, we select a location for this study where both diseases are prevalent. We choose Malkangiri, the southern district of Odisha (Orissa) State in India. Malkangiri has a history of JE outbreaks since 2009 where the cases are documented as Acute Encephalitis Syndrome (AES) [21]. This AES can include infections other than JE. Studies showed that AES cases are not always cases of JE infections, they can be cases of the severe form (the 2nd phase) of leptospirosis infections [22,23,24]. This fact, along with the high prevalence of leptospirosis in cattle in neighbor districts [25,26], illustrates the presence of human leptospirosis cases in our study location.

    We use an SIR model to understand the dynamics of JEV among pig, mosquito, and human populations; also, to understand the dynamics of leptospirosis between cattle and human populations. JEV mainly exists in a transmission cycle between mosquitoes and pigs, and humans get the infection of JE from bites of infected mosquitoes ($ I_m $) [1,2]. However, infected humans do not infect feeding mosquitoes due to the lack of sufficient viremia. So, humans are dead-end hosts for JEV [1,2,3,27]. Also, cattle have no role in the maintenance of JEV in nature [20]. In addition to these facts, we assume mosquito population density depends on breeding sites, not on cattle population density. Thus, humans and cattle have no contribution to JEV transmission. We consider total number of pig bites as $ \lambda_r = \frac{b_mN_m}{k_rN_r+k_cN_c+N_h}k_rN_r $ and total number of human bites as $ \lambda_h = \frac{b_mN_m}{k_rN_r+k_cN_c+N_h}N_h $, where $ N_m $, $ N_r $, $ N_c $, and $ N_h $ represent sizes of mosquito, pig, cattle, and human populations respectively, $ b_mN_m $ represents total mosquito bites. Here, $ k_r $ and $ k_c $ are defined based on mosquitoes' feeding preference among hosts (see Table 2). In our model development, we consider the vertical transmission of JEV in pigs, and mosquitoes. All these ideas are illustrated in Figures 13, and collectively system (1).

    $ dScdt=μcNcμcpcLc(βlcLc+μc)ScdLcdt=βlcLcSc+μcpcLc(γlc+μc)LcdRcdt=γlcLcμcRc
    $
    (1a)
    $ dSrdt=μrNrμrprIrλrβmrImNm1NrSrμrSrdIrdt=λrβmrImNm1NrSr+μrprIr(γr+μr)IrdRrdt=γrIrμrRrdSmdt=μmNmμmpmIm(λrβrmIrNr1Nm+μm)SmdImdt=λrβrmIrNr1NmSm+μmpmImμmIm
    $
    (1b)
    $ dSShdt=ΛhλhβmhImNm1NhSShβlhLcSShμhSShdEShdt=λhβmhImNm1NhSSh(βlhLc+γe)ESh(μh+αe)EShdRShdt=γeEShβlhLcRShμhRShdSL1hdt=βlhLcSShλhβmhImNm1NhSL1h(γl1+rl+μh)SL1hdEL1hdt=βlhLcESh+λhβmhImNm1NhSL1h(γl1+μh+αe+γe+rl)EL1hdRL1hdt=βlhLcRSh+γeEL1h(rl+γl1+μh)RL1hdSL2hdt=rlSL1hλhβmhImNm1NhSL2h(γl2+αl+μh)SL2hdEL2hdt=rlEL1h+λhβmhImNm1NhSL2h(γl2+γe+μh+αl+αe)EL2hdRL2hdt=rlRL1h+γeEL2h(γl2+μh+αl)RL2hdSRhdt=γl1SL1h+γl2SL2hλhβmhImNm1NhSRhμhSRhdERhdt=γl1EL1h+γl2EL2h+λhβmhImNm1NhSRh(γe+αe+μh)ERhdRRhdt=γl1RL1h+γl2RL2h+γeERhμhRRh
    $
    (1c)
    Figure 1.  Flow diagram for pig and mosquito population.
    Figure 2.  Flow diagram for cattle population.
    Figure 3.  Flow diagram for human population.

    where $ \lambda_r = \frac{b_mN_m}{k_rN_r+k_cN_c+N_h}k_rN_r $ $ \rm and $ $ \lambda_h = \frac{b_mN_m}{k_rN_r+k_cN_c+N_h}N_h. $

    Leptospirosis spreads through the urine of infected cattle or pigs and is transmitted through the contact of contaminated (by the urine of infected animal) surfaces or water [7]. Here, we ignore leptospirosis infections through the pig because pig's population size ($ N_r $) is very small compared to the cattle population size ($ N_c $), and pigs also recover and become immune much faster than cattle [28,29,30,31,32]. Therefore, in our model cattle are the only source of Leptospira, which causes the disease. Hence, cattle and humans are getting infections just due to the presence of cattle (Figures 2 and 3). Here, we consider that Leptospira can be transmitted in cattle through vertical transmission also.

    Susceptible humans get leptospirosis infections and move to $ SL1_h $ from $ SS_h $. Most leptospirosis infected humans recover from the first phase, while the remaining move to the second phase of the infection ($ SL2_h $). Finally, patients from the second phase move to $ SR_h $ upon their recovery. The $ SR_h $ compartment represents people who already recovered from leptospirosis infection; however, they are still susceptible to JEV. A portion of susceptible humans get JE infections by infected mosquito bites before leptospirosis infection and move to $ ES_h $. Some of these infected humans get leptospirosis infection before their recovery from JE and move to $ EL1_h $. The remaining population of $ ES_h $ move directly to the JE recovered class ($ RS_h $). Similar to $ SL1_h $, most people from $ EL1_h $ move to $ EL2_h $ while others move to $ ER_h $. Recovered people from JEV infections ($ RS_h $) get leptospirosis infections and move to $ RL1_h $. Then, similar to $ SL1_h $ and $ EL1_h $, most of them move to $ RR_h $ after recovery from leptospirosis, while others enter the second phase of infections ($ RL2_h $). Finally, people recover from $ RL2_h $ and move to $ RR_h $, which contains people who recovered from infections of both diseases. All model variables are summarized in Table 1.

    Table 1.  Model variables with definition.
    Variable Definition
    $ S_c $, $ S_r $, $ S_m $ Susceptible cattle, pigs, mosquitoes
    $ L_c, I_r, I_m $ Infected cattle, pigs, mosquitoes
    $ R_c $, $ R_r $ Recovered cattle, pigs
    $ SS_h $ Humans susceptible to both JE and leptospirosis
    $ ES_h $ Humans infected with JE, but susceptible to leptospirosis
    $ RS_h $ Humans recovered from JE and susceptible to leptospirosis
    $ SL1_h $ Humans susceptible to JE, but infected with 1st phase of leptospirosis
    $ EL1_h $ Humans infected with JE and 1st phase of leptospirosis
    $ RL1_h $ Humans recovered from JE and infected with 1st phase of leptospirosis
    $ SL2_h $ Humans susceptible to JE, but infected with 2nd phase of leptospirosis
    $ EL2_h $ Humans infected with JE and 2nd phase of leptospirosis
    $ RL2_h $ Humans recovered from JE and infected with 2nd phase of leptospirosis
    $ SR_h $ Humans susceptible to JE, but recovered from leptospirosis
    $ ER_h $ Humans infected with JE, but susceptible to leptospirosis
    $ RR_h $ Humans recovered from both JE and leptospirosis

     | Show Table
    DownLoad: CSV

    In our model development, we do not consider vaccination even though some JE endemic areas are under a vaccination program. This simplifying assumption is reasonable because many rural areas, where cattle are very common in houses, are not under vaccination program coverage. We also do not consider the seasonal effect on mosquitoes' abundance since the objective of this study is related to overall disease burden, not to the detailed short-term dynamics. Furthermore, we consider constant populations for cattle, pigs, and mosquitoes. Even though mosquito abundance varies seasonally, assuming a constant mosquito population is reasonable here because our focus is on time-average prevalence and incidence, which smoothes out seasonal effects but preserves the overall number of cases. For disease transmission, we choose standard incidence for JEV, because infected mosquitoes ($ I_m $) are free to bite any individual among our host populations. On the other hand, we use mass-action incidence for Leptospira transmission, because the new infections depend on the availability of infected cattle ($ L_c $).

    The World Bank recorded the life expectancy of people at birth in India in 2017 as 69.165 years [33]. Taking the reciprocal of this value we get the natural death rate for humans as $ \mu_h = 3.96\times10^{-5}\rm/day $. As per the 2014 census by the local health department, Malkangiri district had a population of 641,385 in 2014 with 109,483 households [34]. However, according to the World Bank data rural population in India in 2014 and 2018 was 876,035,725 and 892,321,651 respectively [35]. We use the ratio of these two populations to estimate the population of Malkangiri district in 2018, which gives us the population as 653,309 ($ N_h $) with 111,518 households. In our literature review, we did not find any value for the human recruitment rate. Hence, we use one of our equations (1st in the human subsystem (1c)) in a disease-free state which gives the equation $ \Lambda_h = \mu_hN_h $. Then we use this equation, and our already estimated values of $ \mu_h $ and $ N_h $, to get $ \Lambda_h = 25.87\rm \; human/day $.

    We did not find any documented data for the size of cattle and pig populations ($ N_c $ and $ N_r $ respectively) in Malkangiri district. However, we found one documented data source which mentioned that the numbers of humans, cattle and pigs in Korkunda block of Malkangiri district in 2014 were 143,867, 80,583 (75,772+4,811), and 10,007 respectively in 29,667 households [21]. We use these numbers to estimate cattle and pigs per household, and then use our estimated number of households (111,518) to calculate numbers of total cattle and total pigs in the district in 2018, which gives $ N_c = \frac{80,583}{29,667}\times111,518 = 302,911 $ and $ N_r = \frac{10,007}{29,667}\times111,518 = 37,616 $. Now we need to know the feeding pattern of JEV vectors (value of $ k_c $ and $ k_r $) to estimate $ \lambda_r $ and $ \lambda_h $. In the absence of this data in our study region, we use data from a study in a different part of India assuming them representative of mosquitoes' host preferences in the study area. That study estimated the feeding preference for JE vectors in the southern part of India as 46.4% on cattle, 4.8% on pigs, and 1.5% on humans [36]. Using these values we have $ k_c = \frac{46.4}{1.5} $ human/cattle and $ k_r = \frac{4.8}{1.5} $ human/pig.

    To estimate $ \lambda_r $ and $ \lambda_h $, we also need to have the size of JEV carrier mosquito populations ($ N_m $), along with their biting rate ($ b_m $). A field study in rural villages of Western Yunnan Province of China was done in 2013 to estimate the abundance of mosquitoes in an Asian rural setting [37]. Researchers collected mosquitoes from two households in each of four studied villages in each month for a 12 month-study. The total collection of mosquitoes was 85,307. It is equivalent to $ \frac{85,307}{(2\times4)\times12} $ = 888 mosquitoes/household. So, for our study we have a total of (111,518 $ \times $ 888) = 99,027,984 mosquitoes for a 12-months period which includes both high- and low-abundance seasons. Another study carried out in some villages of Malkangiri district showed that 69.1% of available mosquitoes are vectors for JEV [21]. Hence, the total number of JEV vectors we have is $ N_m $ = 99,027,984 $ \times $ 69.1% = 68,428,337. We found documented data for the biting rate of JEV vectors which give $ b_m $ = 0.25 bite/mosquito-day (mean of 0.2 and 0.3) [38]. Anyway, not all mosquito bites belong to these three hosts. Also, blood-meal analysis showed that some mosquitoes bite more than one host [36]. Hence, we estimated that about 70% of mosquito bites are distributed among cattle, pig and human population. Thus, $ \lambda_r = \frac{0.25 \times 68,428,337 \times70\%}{\frac{4.8}{1.5} \times37,616 + \frac{46.4}{1.5} \times302,911 + 653,309}\times(\frac{4.8}{1.5}\times37,616 $) = 142,102 bite (total)/day and $ \lambda_h = \frac{0.25 \times68,428,337 \times70\%}{\frac{4.8}{1.5}\times37,616 + \frac{46.4}{1.5} \times302,911+653,309} \times653,309 $ = 771,250 bite (total)/day.

    In 2000, the Indian Council of Medical Research carried out a study which found the transmission and infection rate for JEV vectors [39]. The study showed that vectors have 82% infection probability from pigs; however, the transmission probabilities of infections from mosquitoes to humans and mosquitoes to pigs are 31.6 and 63.5% respectively. Hence, we have $ \beta_{rm} $ = 0.82 mosquito/bite, $ \beta_{mr} $ = $ 0.635\rm\; pig/bite $, $ \rm\; $ and $ \beta_{mh} = 0.316\rm\; human/bite $. A separate review study on the expansion of JEV carried out in 2009 which mentioned that pigs maintain enough viremia to infect mosquitoes for up to 4 days [27]. Taking the reciprocal of this value we get $ \gamma_r = \frac{1}{4} $/day. For JE vectors' mortality rate, we found a study that estimated their life expectancy as 59.8 days [40], and this gives us $ \gamma_m = \frac{1}{59.8} $/day. Humans' mean recovery period during the 2012 JE outbreak in Odisha was 7 days [41], which gives us $ \gamma_e = \frac{1}{7} $/day. For the estimation of $ \mu_r $, we found that domestic pigs have an average lifespan of 6–10 years, but, that can be shorter due to certain problems [42]. Thus, instead of the median (8 years) of the range, we choose 7 years as the life expectancy for pig population which eventually estimated $ \mu_r = \frac{1}{7\times365}\rm/day = 3.91\times10^{-4} $/day.

    A review study from 2014 mentioned the mean durations for the acute phase (1st phase, in our model) and immune phase (2nd phase, in our model) as 10 and 14 days respectively [43]. Taking reciprocals of these values we get $ \gamma_{l1} = \frac{1}{10}\rm\; /day $ and $ \gamma_{l2} = \frac{1}{14}\rm\; /day $. Another study in coastal south India estimated that 51 out of 202 (25.24%) leptospirosis patients move to the 2nd phase [44]. Hence, we use this result and the ratio $ \frac{r_l}{\gamma_{l1}+r_l+\mu_h} = \frac{51}{202} $ to estimate the transfer rate from phase 1 to phase 2, and found $ r_l = 0.033 $/day. We found a couple of documented data for cattle recovery period of leptospirosis which are inconsistent in values [30,31,32]. Hence, based on those available values we choose that cattle can shed leptospira in their urine for 60 weeks, which gives $ \gamma_{lc} = \frac{1}{60\times7}\rm/day = 2.38\times10^{-3} $/day. In 2016, a different disease in cattle in Odisha state was studied which had the highest cattle age group of 6.5–7.5 years; however, they did not mention the highest or average age [45]. So, we take 8 years as the lifespan for cattle in Odisha, which gives us $ \mu_c = \frac{1}{8\times365}\rm/day = 3.42\times10^{-4} $/day.

    We did not find enough evidence for vertical transmission of leptospirosis in cattle ($ p_c $) and of JEV in pigs ($ p_r $), and therefore we assume $ p_c = 0\rm $ and $ p_r = 0 $. However, we found 4% effective vertical transmission for JEV vectors, which gave $ p_m = 0.04 $ [46]. In 2016, in Malkangiri district 37 died out of 175 JE patients [34]. Using this ratio and the relation $ \frac{\alpha_e}{\alpha_e+\gamma_e+\mu_h} = \frac{37}{175} $, we estimated $ \alpha_e = 0.0383/ $day. The case fatality rate for leptospirosis infection is 7.69% [47]. Also, we already know that about 25.24% leptospirosis cases move to the 2nd phase [44]. Hence, the case fatality rate among patients of 2nd phase is $ \frac{7.69}{25.24}\times 100 $, which implies $ \alpha_l = 3.13\times10^{-2} $/day.

    During our literature review, we did not find any value for infection rates of leptospirosis in cattle and in humans ($ \beta_{lc} $ and $ \beta_{lh} $ respectively). Hence, we follow the method used in [18,48] to estimate $ \beta_{lc} $ and $ \beta_{lh} $which gives $ \beta_{lc} = \frac{\gamma_{lc}+\mu_c}{\left(1-\left(\frac{\gamma_{lc}+\mu_c}{\mu_c}\right)y_c\right)N_c} $, $ \beta_{lh} = \frac{\mu_hN_hy_h}{\left(\left(1+\frac{r_l}{\gamma_{l2}+\mu_h+\alpha_l}\right)\frac{\Lambda_h}{\gamma_{l1}+\mu_h+r_l}-y_hN_h\right)N_cy_c} $ where $ y_c $ and $ y_h $ are prevalence of leptospirosis in cattle and in humans respectively. We did not find any documented data for $ y_c $ in Malkangiri. However, we found one study from 2013 which estimated leptospirosis prevalence as 20.7% in cattle in a neighbour district of Malkangiri [26]. We choose $ y_c = \frac{20.7\%}{2} = 11.35\% $ since the sample collections of this study were not random, rather were mostly from the villages with a history of abortions and other disorders. Also, we did not find documented data for leptospirosis prevalence in humans in Malkangiri. However, we found human and cattle prevalence of leptospirosis in another prevalent location of leptospirosis, northeast Thailand. We estimated leptospirosis prevalence in cattle for this area to be $ \frac{34+2+58}{130+183+238} = \frac{94}{551} = 17\% $ using the data from a study published in 2016 [49]. A different study estimated 12.5 annual cases of human leptospirosis per 100,000 people [50]. We used the ratio of estimated leptospirosis prevalence in cattle for our study area (11.35%) and of northeast Thailand (17%), and the annual cases of human leptospirosis in northeast Thailand (12.5 per 100,000) to calculate annual cases of human leptospirosis for our study area, and found $ \frac{12.5\times11.35}{17} = 8.35 $ cases per 100,000 people, which is equivalent to 0.00835%. Also, we have a weighted recovery period as $ \frac{(10+14)\times25.24\%+10\times74.76\%}{100\%} = 13.53 $ days, which gives us $ \frac{365}{13.53} = 27 $ as the number of generations in a year for infected people. Finally, we got $ y_h = \frac{0.00835\%}{27} = 3.09\times10^{-6} $. Using these prevalence values and our expressions for leptospirosis infection rates, we got $ \beta_{lc} = 9.66\times10^{-8} $/day-cattle and $ \beta_{lh} = 9.11\times10^{-12} $/day-cattle. All model parameters are summarized in Table 2.

    Table 2.  Summary of model parameters.
    Par. Definition Value Units Ref.
    $ \Lambda_h $ recruitment rate for humans 25.87 human/day This study
    $ p_c $ probability of vertical transmissions of leptospirosis in cattle 0 - This study
    $ p_r $ probability of vertical transmissions of JE in pigs 0 - This study
    $ p_m $ probability of vertical transmissions of JE in mosquitoes 0.04 - This study
    $ \beta_{lc} $ infection rate of leptospirosis to cattle $ 9.3\times10^{-8} $ 1/day-cattle This study
    $ \beta_{lh} $ infection rate of leptospirosis to humans $ 9.12\times10^{-12} $ 1/day-cattle This study
    $ \beta_{rm} $ infection probability of JE to mosquito from biting pigs 0.82 mosquito/bite [39]
    $ \beta_{mr} $ infection probability of JE to pigs from mosquitoes' bite 0.635 pig/bite [39]
    $ \beta_{mh} $ infection probability of JE to humans from mosquitoes' bite 0.316 human/bite [39]
    $ \gamma_{e} $ recovery rate from JE $ \frac{1}{7} $ 1/day [41]
    $ \gamma_{lc} $ recovery rate from leptospirosis for cattle $ 2.38\times10^{-3} $ 1/day This study
    $ \gamma_{l1} $ recovery rate from 1st phase of leptospirosis (humans) $ \frac{1}{10} $ 1/day [43]
    $ \gamma_{l2} $ recovery rate from 2nd phase of leptospirosis (humans) $ \frac{1}{14} $ 1/day [43]
    $ r_l $ transfer rate from phase 1 to phase 2 of leptospirosis 0.033 1/day This study
    $ \gamma_r $ recovery rate from JE for pigs $ \frac{1}{4} $ 1/day [27]
    $ \mu_c $ natural death rate for cattle $ 3.42\times10^{-4} $ 1/day This study
    $ \mu_r $ natural death rate for pigs $ 3.91\times10^{-4} $ 1/day This study
    $ \mu_m $ natural death rate for mosquitoes $ \frac{1}{59.8} $ 1/day [40]
    $ \mu_h $ natural death rate for humans $ 3.96\times10^{-5} $ 1/day [33]
    $ \alpha_e $ JE-induced death rate for humans 0.0383 1/day This study
    $ \alpha_l $ leptospirosis-induced death rate (during 2nd phase) for humans $ 3.13\times10^{-2} $ 1/day This study
    $ b_m $ mosquito's biting rate 0.25 bite/mosquito-day [38]
    $ N_c $ cattle population size of the study area 302,911 cattle This study
    $ N_m $ mosquito population size of the study area 68,428,337 mosquitoes This study
    $ N_r $ pig population size of the study area 37,616 pigs This study
    $ N_h $ human population size of the study area 653,309 humans This study
    $ k_c $ vector's feeding preference (humans over cattle) $ \frac{46.4}{1.5} $ human/cattle This study
    $ k_r $ vector's feeding preference (humans over pigs) $ \frac{4.8}{1.5} $ human/pig This study
    $ \lambda_r $ encounter rate between pigs and mosquitoes 142,102 bite (total)/day This study
    $ \lambda_h $ encounter rate between humans and mosquitoes 771,250 bite (total)/day This study

     | Show Table
    DownLoad: CSV

    We begin with the general model where we consider the presence of cattle which includes both JE and leptospirosis diseases. Later, we consider the scenario without cattle; this special scenario has JEV infection only. Finally, we estimate total disease burdens for these two scenarios and compare them to understand the impact of cattle on disease burden in JE prevalent areas. The goal of this study is actually to understand how the presence of domestic animals, here cattle, helps to reduce disease burden.

    First, we analyse the cattle system (1a), which is decoupled from the other two subsystems – the pig-mosquito system (1b), and the human system (1c). We set all the equations of this system equal to zero and solve them for state variables to find equilibria values. After performing some basic algebra we get the disease free equilibrium ($ DFE $) as $ S_c = N_c, L_c = 0, R_c = 0 $, and the endemic equilibrium ($ EE $) as

    $S^*_c = \frac{\gamma_{lc}+\mu_c(1-p_c)}{\beta_{lc}},\; L^*_c = \left(N_c-\frac{\gamma_{lc}+\mu_c(1-p_c)}{\beta_{lc}}\right)\frac{\mu_c}{\gamma_{lc}+\mu_c},\; R^*_c = \left(N_c-\frac{\gamma_{lc}+\mu_c(1-p_c)}{\beta_{lc}}\right)\frac{\gamma_{lc}}{\gamma_{lc}+\mu_c}.$

    Then we use the next generation method [51] to obtain the basic reproduction number ($ BRN $) which gives $ R_{0l} = \frac{\beta_{lc}N_c+\mu_cp_c}{\gamma_{lc}+\mu_c} $. We find that the $ EE $ exists under the condition $ R_{0l} > 1 $. Here the $ EE $ is globally asymptotic stable if and only if $ R_{0l} > 1 $, otherwise the $ DFE $ is globally stable (see Appendix 1 for details).

    Before we begin our discussion on the pig-mosquito subsystem, here we consider a simplifying assumption that the two diseases do not significantly affect the magnitude of the human population, so that we can replace $ N_h(t) $ by $ N_h(0) $ in the denominator of $ \lambda_r $ and $ \lambda_h $, in order to simplify the analysis. This assumption is supported by our literature review: as per the study [34], JEV infection causes $ \frac{37}{641,385}\times 100,000 = 5.77 $ annual deaths per 100,000 people. Also, leptospirosis infection causes $ 12.5\times7.69 = 0.96 $ annual death per 100,000 people [47,50]. These two estimations give us a total of 44 $ ( = \frac{653,309}{100,000}\times(5.77+0.96)) $ annual deaths due to both diseases, which is 0.0067% of the total human population. The assumption is also supported by later numerical analysis. Under this assumption, mosquito's biting rates to pigs, and humans ($ \frac{\lambda_r}{N_r} $ and $ \frac{\lambda_h}{N_h} $ respectively) become constants, and the vector-reservoir (pig-mosquito) subsystem is decoupled from the human subsystem. The pig-mosquito subsystem is already decoupled from the cattle subsystem. Now, following the same approaches, like the cattle system, for the pig-mosquito system, we get the $ DFE $ as $ S_r = N_r, I_r = 0, R_r = 0, S_m = N_m, I_m = 0, $ and the $ EE $ as

    $ Sr=λrβrm+Nmμm(1pm)μr+γrμrλrβrmNr(λrβmrNrμr+γrμr+(γr+μr(1pr)))(γr+μr(1pr)),Ir=λrβrmNrλrβmrNmμm(1pm)(γr+μr(1pr))λrβrmNr(λrβmrNrμr+γrμr+(γr+μr(1pr))),Rr=λrβrmNrλrβmrNmμm(1pm)(γr+μr(1pr))λrβrmNr(λrβmrNrμr+γrμr+(γr+μr(1pr)))γrμr,Sm=λrβmrμr+γrμr+Nr(γr+μr(1pr))λrβmrNm(λrβrmNm+μm(1pm)μr+γrμr)(1pm)μm,Im=λrβrmNmλrβmrNr(1pm)μm(γr+μr(1pr))λrβmrNm(λrβrmNm+μm(1pm)μr+γrμr).
    $
    (2)

    We get the $ BRN $ for the pig-mosquito system as

    $ R_{0e} = \frac{1}{2}\left(\left(\frac{\mu_r}{\mu_r+\gamma_r}p_r+p_m\right)+\sqrt{\left(\frac{\mu_r}{\mu_r+\gamma_r}p_r-p_m\right)^2+4\frac{\lambda_r\beta_{mr}}{\mu_m N_m}\frac{\lambda_r\beta_{rm}}{\left(\mu_r+\gamma_r\right)N_r}} \right). $

    Here the $ EE $ exists and is locally asymptotically stable (LAS) if and only if $ R_{0e} > 1 $, otherwise the $ DFE $ is LAS (see Appendix 2 for details). Next, we develop two strong Lyapunov functions to show that the $ DFE $ is globally stable when $ R_{0e} < 1 $ (see Appendix 3 for details) and the $ EE $ is globally stable if and only if $ R_{0e} > 1 $ (see Appendix 2 for details).

    The cattle subsystem and the pig-mosquito subsystem are decoupled from the human subsystem respectively by the model and the assumption that the magnitude of the human population is insignificantly affected by the effect of JE and leptospirosis. As these two subsystems go to equilibrium, the limiting system for the system (1) is given by the human subsystem with $ L_c\; {\rm{and}}\; I_m $ replaced by their equilibrium values (see system (3)). Here $ N_m $ is constant which makes the term $ \frac{I^*_m}{N_m} $ to be constant and $ \frac{\lambda_{h0}}{N_h} $ ($ = \frac{b_mN_m}{k_rN_r+k_cN_c+N_{h0}} $ = constant) is also constant. Thus, the resulting system (system (3)) is linear and has only one equilibrium (see Eq (17) in Appendix 4) where the equilibrium values are functions of $ L^*_c\; {\rm{and}}\; I^*_m $. Hence, the behavior of the limiting system follows the behavior of our cattle and pig-mosquito subsystems. Now, a theorem by Horst R. Thieme [52] leads us to the fact that the behavior of the entire system is asymptotic to the behavior of the limiting system, which is linear. Eventually, the entire system is governed by the behavior of the cattle and the pig–mosquito subsystems.

    Next, we extend our analysis to the entire system, which includes the cattle, the pig–mosquito, and the human subsystems. We obtain the invasion reproduction number ($ IRN $) for both diseases following the approach by Mitchell & Kribs [53]. Each of these $ IRN $s is found to be the same as the corresponding $ BRN $s. Thus, both of the $ IRN $s are independent of the presence of the other disease. This result is easily comprehensible since neither of these diseases' transmissions is affected by the other. Our preceding analyses found two different values for both $ L_c^* $ and $ I_m^* $; hence, their all possible combinations produces four different equilibria for the entire system. These four equilibria are – $ DFE $, endemic in leptospirosis and free of JE, endemic in JE and free of leptospirosis, and endemic in both diseases respectively. So, for the entire system, we have four different scenarios – (ⅰ) $ R_{0l} < 1\; {\rm {and}}\; R_{0e} < 1 $ (ⅱ) $ R_{0l} > 1\; {\rm {and}}\; R_{0e} < 1 $, (ⅲ) $ R_{0l} < 1\; {\rm {and}}\; R_{0e} > 1 $, and (ⅳ) $ R_{0l} > 1\; {\rm {and}}\; R_{0e} > 1 $.

    Now we consider the scenario without cattle which eliminates the possibility of leptospirosis infections in humans. Here the $ BRN $ for JE ($ \mathcal{R}_{0e} $) is different compared to the $ BRN $ ($ R_{0e} $) from the general model since the $ \lambda_r $ in $ R_{0e} $ is a function of cattle population size. The absence of cattle reduces the human system to a three-dimensional system, which gives the $ EE $ for the human system as

    $ SSh=Λhμh+λh0βmhNhImNm,ESh=λh0βmhNhImNm(μh+γe+αe)(Λhμh+λh0βmhNhImNm),RSh=λh0βmhNhImNmγeμh(μh+γe+αe)(Λhμh+λh0βmhNhImNm)
    $

    where $ \lambda_{h0} = \frac{b_mN_m}{k_rN_r+N_{h0}}N_h $ (here $ N_{h0} = N_h(0) $). This special case has $ DFE $, and the only endemic equilibrium ($ EE $) when $ \mathcal{R}_{0e} > 1 $. This $ EE $ is completely governed by $ I^*_m $ (see Eq (2) for $ I^*_m $).

    $ dSShdt=Λhβmhλh0NhImNmSShβlhLcSShμhSShdEShdt=βmhλh0NhImNmSSh(βlhLc+γe)ESh(μh+αe)EShdRShdt=γeEShβlhLcRShμhRShdSL1hdt=βlhLcSShβmhλh0NhImNmSL1h(γl1+rl+μh)SL1hdEL1hdt=βlhLcESh+βmhλh0NhImNmSL1h(γl1+μh+αe+γe+rl)EL1hdRL1hdt=βlhLcRSh+γeEL1h(rl+γl1+μh)RL1hdSL2hdt=rlSL1hβmhλh0NhImNmSL2h(γl2+αl+μh)SL2hdEL2hdt=rlEL1h+βmhλh0NhImNmSL2h(γl2+γe+μh+αl+αe)EL2hdRL2hdt=rlRL1h+γeEL2h(γl2+μh+αl)RL2hdSRhdt=γl1SL1h+γl2SL2hβmhλh0NhImNmSRhμhSRhdERhdt=γl1EL1h+γl2EL2h+βmhλh0NhImNmSRh(γe+αe+μh)ERhdRRhdt=γl1RL1h+γl2RL2h+γeERhμhRRh
    $
    (3)

    where $ \lambda_{h0} = \frac{b_mN_m}{k_rN_r+k_cN_c+N_{h0}}N_h $ (here $ N_{h0} = N_h(0) $).

    To better understand the impact of cattle, some numerical analyses are done based on our parameter estimations. Using our estimated parameter values, we find $ R_{0l} = 10.35\; {\rm{and}}\; R_{0e} = 1.008 $ respectively in the presence of cattle. However, our estimations in the absence of cattle found $ \mathcal{R}_{0e} = 12.97 $. Based on our estimated parameter values, in the presence of cattle we estimate 72 and 228 annual cases for leptospirosis and JE respectively, whereas the annual deaths are 6 and 48 respectively. In the absence of cattle, we estimate 9407 and 1988 number of annual cases and annual deaths due to JEV infections. Next, we perform numerical analysis to estimate total cases and the BRN of JE and leptospirosis as the average number of cattle per household varies and plot the related graphs (Figure 4). Figure 4(a) shows that the JEV infection decreases as the number of cattle increases; however, the human leptospirosis incidence increases with the cattle; Figure 4(b) illustrates that BRN for JE exponentially decreases, and for leptospirosis linearly increases as the number of average cattle increases. So, the presence of cattle has positive and negative impacts on human health.

    Figure 4.  The effect of cattle on infections and BRNs.

    To quantify the impact of cattle on human health, we calculate the total disease burden for our two settings – with cattle and without cattle. In this calculation we follow an approach similar to the approach of the study [54] where they followed WHO guidelines [55]. To calculate the total disease burden they estimated total $ DALY $ (Disability Adjusted Life Year) where one $ DALY $ means the loss of one year of healthy life. The $ DALY $ is composed of $ YLD $ (Years Lost due to Disability) resulting from infections and $ YLL $ (Years of Life Lost) caused by disease-induced premature deaths which are calculated using formulas $ YLD = I\times DW\times L_1 $ and $ YLL = D\times L_2 $, where $ I $ and $ D $ represent the total number of infections (cases) and total number of deaths respectively. Here $ L_1 = $ average duration (in years) of illness which is the reciprocal of the recovery rate for survival cases and of the death rate for non-survival cases, and $ L_2 = $ standard life expectancy at age of death (in years) = average life expectancy of the community – average age at premature death and the term $ DW $ is the disability weight for diseases which ranges from 0 (perfect health) to 1(death). It can be thought of as the proportional reduction in perfect health due to any adverse health condition caused by infections. So, the total burden of disease can be represented as the sum of $ YLD $ and $ YLL $, which gives

    $ DALY = I\times DW\times L_1+D\times L_2. $

    In this work, we study two different diseases, JE and leptospirosis. Hence, we will have different values of $ I, D, L_1, L_2 $ and of $ DW $ for JE, and leptospirosis. Here the quantities $ L_1 $s are the reciprocals of the $ \gamma_{l1}, \gamma_{l2} $ for the first and second phases of leptospirosis respectively, and the reciprocal of the $ \gamma_e $ for JE. To estimate $ L_2 $s, we need to know the average ages of infections which are 40.48 [44] and 5 years [41] for leptospirosis and JEV infections respectively. The average age of infections and the average age of deaths are the same for leptospirosis; however, they are different for JE (5 and 3 years respectively [41]). The values of $ DW $ differs between leptospirosis and JE; it also differs between phases of leptospirosis. Furthermore, the $ DW $ has different values for survival and non-survival. All these different values of $ DW $ are taken from the study [56] based on the severity (mild or severe) of illness and on the onward complications upon recovery. Therefore, using the above formula for $ DALY $ we calculate total $ DALY $s separately for JE and leptospirosis infections and add them to have total $ DALY $ cause from both diseases combined.

    Based on our parameter estimations, the above formula and values together calculate 4396 $ DALY $s and 157,642 $ DALY $s for the setting with cattle and without cattle respectively. These results clearly show that the presence of cattle has a huge positive impact on the disease burdens in JE-prevalent areas. Next, we calculate total $ DALY $s varying the average number of cattle per household to understand how the variation in the average number of cattle per household changes the annual disease burden. The results of these calculations are portrayed in Figure 5, which evidently illustrates that the presence of cattle is unconditionally helpful in terms of the annual disease burden. However, there is an optimum value for the average number of cattle per household to ensure the minimum disease burden which is 2.75. Maintaining cattle beyond this threshold value the annual disease burden increases slowly which is generated from the infections of leptospirosis. On the other hand, keeping fewer cattle per household than the optimum number causes a sharp rise in the annual disease burden due to a significant increase in JE infections.

    Figure 5.  Relation between avg cattle/household and annual disease burden (in $ DALY $s).

    To the best of our knowledge, this is the first work to understand the joint dynamics of JE and leptospirosis. Our qualitative analysis exhibits the classic threshold behavior of the system. The quantitative analysis shows that the presence of cattle reduces total disease burden even though it introduces leptospirosis. It happens because the number of JE cases is reduced by the cattle presence at a higher rate than the rate at which humans get leptospirosis infections due to the presence of cattle. Also, the case fatality rate (CFR) of JE infections is higher than the CFR of leptospirosis infections. The disability weight is higher for JE infections than for the leptospirosis infections, and the average age of JEV infections is very low compared to the average age of leptospirosis infections. Even though the presence of cattle is always helpful, there is an optimum average number of cattle per household which minimizes the annual disease burden.

    Our analysis found that the presence of cattle helps to reduce the annual disease burden in JE-prevalent areas because of the vector's (mosquitoes') feeding preference. Mosquitoes prefer to bite cattle at a much higher proportion compared to the proportion to bite humans or pigs. Thus, cattle reduce mosquito-human interactions and eventually reduce human cases of JEV infections. This result is consistent with the result of an ecological study by Miller and Huppert [15] where they identified that vectors' feeding preference decides whether or not host richness is helpful to dilute disease prevalence. Also, our results found that the JE prevalence in humans decreases as the cattle population increases. This finding is consistent with the finding of another ecological study on vector-borne diseases by Johnson and Thieltges [14] where they found that the strength of the dilution effects depends on the relative abundance of the dilution host to the focal host.

    As noted in the model development, we chose not to consider seasonality in this study because our focus was primarily with overall disease burden, not with detailed disease dynamics. Future work incorporating seasonal effects would provide a more accurate picture of the short-term disease dynamics. In particular, mosquito abundance is strongly influenced by temperature and precipitation, which are seasonal, so that mosquito population density is more properly periodic. Analysis of such a model could identify times of year when one disease is a more acute burden than the other.

    Future work could also benefit from techniques to reduce the dimensionality of large systems such as that studied here. In [57], Yu developed a technique, later applied in other studies, for reducing the dimension of a population density model for mosquitoes by considering the number of sterile mosquitoes to be a known function (typically piecewise constant in applications) rather than a state variable. This technique could be applied to other systems where a subpopulation is a known function.

    Our study has other potential limitations. In our model, we assumed that the total mosquito biting rate depends only on mosquito density, which should not have any significant impact on our results since the mosquito population is not limited by the available host populations, but rather by available breeding sites. In our analysis, we assumed that the biting rates to pigs and humans were independent of any fluctuations in the host (cattle, pig and human) populations, but the only such fluctuations (disease-induced deaths among humans) represent such a small proportion of the human population that this assumption should also have a negligible effect. We also assumed a homogeneous distribution among households and those cattle close enough to households to affect the infection risks for JEV and leptospirosis. Large cattle farms farther away from homes than pigs are, are not likely to affect those risks (except to employees). However, studies have shown that heterogeneity does affect disease risk [58,59]. In our case, heterogeneity is the heterogeneous distribution of household cattle per household.

    All authors declare that there is no conflicts of interest in this paper.

    Proposition 1: The endemic equilibrium of the cattle subsystem is locally stable if and only if $ R_{0l} > 1 $.

    Proof: We find the Jacobian matrix of the cattle system at the $ EE $ which gives the characteristic equation

    $ \left(\lambda^2+\mu_cR_{0l}\lambda+\mu_c(\gamma_{lc}+\mu_c)(R_{0l}-1)\right)(\lambda+\mu_c) = 0. $

    Here, $ \lambda = -\mu_c $ is negative. Next, we use Routh-Hurwitz Criteria to check if the real part of remaining eigenvalues are negative. As per Routh-Hurwitz Criteria, the coefficients of $ \lambda $ and the constant term need to be positive to ensure the negativity of real parts of eigenvalues. Here it is evident that the coefficient of $ \lambda $ is always positive and the constant term is positive when $ R_{0l} > 1 $.

    Hence, the $ EE $ of the cattle subsystem is locally stable if and only if $ R_{0l} > 1. $

    Proposition 2: The endemic equilibrium of the cattle subsystem is globally stable when $ R_{0l} > 1 $, otherwise the $ DFE $ is globally stable.

    Proof: The first two equations of the cattle subsystem are

    $ Sc=μcNcμcpcLc(βlcLc+μc)ScLc=βlcLcSc+μcpcLc(γlc+μc)Lc.
    $
    (4)

    It is evident from the system (4) that $ {S'}_c > 0 \; {\rm{if}}\; S_c = 0 \; {\rm{and}}\; {L'}_c = 0 \; {\rm{if}} \; L_c = 0 $, which means the system (4) is well posed. Now, the system (4) gives

    ${S'}_c+{L'}_c = {(S_c+L_c)'} = \mu_cN_c-\mu_cS_c-(\mu_c+\gamma_{lc})L_c$.

    From this relation it is evident that $ {(S_c+L_c)'}\le 0 $ if either $ S_c\ge N_c $, or $ L_c\ge N_c $ which shows all solutions are bounded.

    Now we define $ \beta(S_c, L_c) = \frac{1}{S_cL_c} $ to use Dulac's Criterion to check if the $ EE $ approaches a limit cycle. Here we assume $ F(S_c, L_c) = {S'}_c $ and $ G(S_c, L_c) = {L'}_c $ which give

    $ Sc(βF)+Lc(βG)=μcNcLcS2c+μcpcScS2c=μcS2c(pcNcLc)0
    $
    (5)

    So, no limit cycle exists.

    Hence, by Poincaré-Bendixson Theorem all solutions approach the $ EE $ when $ R_{0l} > 1 $ which means the $ EE $ is globally stable if and only if $ R_{0l} > 1 $, otherwise the $ DFE $ is globally stable.

    Proposition 3: The endemic equilibrium of the pig–mosquito subsystem is locally stable if and only if $ R_{0e} > 1 $.

    Proof: We can reduce our system to a system of three equations since pig and mosquito populations are constant. We replace $ S_r = N_r-(I_r+R_r) \; {\rm{and}}\; S_m = N_m-I_m $ in the last equation of the subsystem (1a) and 2nd equation of the subsystem (1b), and get the reduced pig-mosquito system as

    $ Ir=dIrdt=λrβmrImNm1Nr(Nr(Ir+Rr))(μr(1pr)+γr)IrRr=dRrdt=γrIrμrRrIm=dImdt=λrβrmIrNr1Nm(NmIm)μm(1pm)Im
    $
    (6)

    Next, we find the Jacobian matrix of the reduced system (6) at the $ EE $ which gives the characteristic equation

    $ \lambda^3+A_1\lambda^2+A_2\lambda+A_3 = 0. $

    Here,

    $ A1=(a1+b2+c3),A2=a1b2a2b1+a1c3a3c1+b2c3,A3=a3b2c1+a2b1c3a1b2c3
    $

    where

    $ a1=λrβmrNmNrIm(γr+μr(1pr)),a2=λrβmrNmNrIm,a3=λrβrmNmNr(NrIrRr),b1=γr,b2=μr,c1=λrβrmNmNr(NmIm),c3=λrβrmNmNrIrμm(1pm).
    $

    Next, as per Routh-Hurwitz Criteria, we need to verify $ A_1, A_2 > 0 \; {\rm{and}}\; A_1A_2 > A_3 $ to check if the real parts of all eigenvalues are negative. Here, $ A_1 $ is clearly positive.

    Now,

    $ A2=a1b2a2b1+a1c3a3c1+b2c3=a1b2a2b1+(λrβmrNmNrIm(γr+μr(1pr)))(λrβrmNmNrIrμm(1pm))λrβmrNmNr(NrIrRr)λrβrmNmNr(NmIm)+b2c3=a1b2a2b1+b2c3+λrβmrNmNrIm(λrβrmNmNrIr+μm(1pm))+(γr+μr(1pr))(λrβrmNmNrIr+μm(1pm))λrβmrNmNrλrβrmNmNrSrSm=a1b2a2b1+b2c3+λrβmrNmNrIm(λrβrmNmNrIr+μm(1pm))+λrβrmNmNrIr(γr+μr(1pr))+(γr+μr(1pr))μm(1pm)λrβmrNmNrλrβrmNmNrNmNrλrβmrNmNrλrβrm(γr+μr(1pr))μm(1pm)(usingSrandSmvaluesfromEq(2))=a1b2a2b1+b2c3+λrβmrNmNrIm(λrβrmNmNrIr+μm(1pm))+λrβrmNmNrIr(γr+μr(1pr))A2>0,sinceeachtermhereispositive.
    $

    Finally,

    $ A1A2A3=(a1+b2+c3)(a1b2a2b1+a1c3a3c1+b2c3)(a3b2c1+a2b1c3a1b2c3)=a1(a1b2a2b1+a1c3a3c1)a1b2c3b2(a1b2a2b1+a1c3+b2c3)+b2a3c1c3(a1b2+a1c3a3c1+b2c3)+c3a2b1(a3b2c1+a2b1c3a1b2c3)=a21b2+a1a2b1a21c3+a1a3c1a1b22+a2b1b2a1b2c3b22c3a1b2c3a1c23+a3c1c3b2c23=a21b2+a1a2b1a1b22+a2b1b2a1b2c3b22c3a1b2c3b2c23+(a1a3c1+a3c1c3a21c3a1c23)=B1+B2
    $
    (7)

    where

    $ B1=a21b2+a1a2b1a1b22+a2b1b2a1b2c3b22c3a1b2c3b2c23>0andB2=a1a3c1+a3c1c3a21c3a1c23=(a1+c3)(a3c1a1c3)=(a1+c3)(λrβmrNmNrIm(λrβrmNmNrIr+μm(1pm))+λrβrmNmNrIr(γr+μr(1pr))))B2>0,since(a1+c3)<0.
    $

    Thus, Eq (7) gives $ A_1A_2-A_3 > 0\; \implies\; A_1A_2 > A_3 $.

    Hence, the endemic equilibrium for the pig–mosquito subsystem is locally stable if and only if $ R_{0e} > 1 $.

    Proposition 4: The endemic equilibrium of the pig–mosquito subsystem is globally stable if and only if $ R_{0e} > 1 $.

    Proof: Here we use a alternative threshold quantity $ \mathcal{R} $ to understand how the value of $ R_{0e} $ affects the stability of the system which is defined as

    $ R=λrβrmλrβmrNrNm1μm(1pm)(γr+μr(1pr))
    $
    (8)

    where $ \mathcal{R} > 1\; {\rm{if \;and\; only\; if}}\; R_{0e} > 1. $

    Thus, $ I_r^*, R_r^*, \; {\rm{and}}\; I_m^* $ from Eq (2) can be written as

    $ Ir=(R21)Nmμm(1pm)(γr+μr(1pr))λrβrmNr(λrβmrNrμr+γrμr+(γr+μr(1pr))),Rr=(R21)Nmμm(1pm)(γr+μr(1pr))λrβrmNr(λrβmrNrμr+γrμr+(γr+μr(1pr)))γrμr,Im=(R21)Nrμm(1pm)(γr+μr(1pr))λrβmrNm(λrβrmNm+μm(1pm)μr+γrμr).
    $
    (9)

    Now, we define our Lyapunov function as

    $ V(Ir,Rr,Im)=IrlnIrIr+(IrIr)+RrlnRrRr+(RrRr)+ImlnImIm+(ImIm)
    $
    (10)

    which gives

    $ V(Ir,Rr,Im)=(1IrIr)Ir+(1RrRr)Rr+(1ImIm)Im
    $
    (11)

    From Eqs (10) and (11) it is evident that

    $ V(Ir,Rr,Im)=0andV(Ir,Rr,Im)=0for(Ir,Rr,Im)=(Ir,Rr,Im)ifandonlyifR>1.
    $

    Now, we will show $ V(I_r, R_r, I_m) > 0 $ for $ (I_r, R_r, I_m)\neq (I_r^*, R_r^*, I_m^*) $.

    Equation (10) gives

    $ VIr=1IrIr.Thus,VIr=0Ir=Ir.
    $
    $ Similarly,VRr=0Rr=RrandVIm=0Im=Im.
    $

    Hence, the only critical point is $ (I_r^*, R_r^*, I_m^*) $.

    Equation (10) also gives

    $ 2VI2r=IrI2r,2VR2r=RrR2r,and2VI2m=ImI2m.
    $

    Now,

    $ D1=2VI2r=IrI2rD1(Ir,Rr,Im)=1Ir,D2=|2VI2r2VIrRr2VRrIr2VR2r|D2(Ir,Rr,Im)=|1Ir001Rr|=1IrRr,andD3=|2VI2r2VIrRr2VIrIm2VRrIr2VR2r2VRrIm2VImIr2VImRr2VI2m|D3(Ir,Rr,Im)=|1Ir0001Rr0001Im|=1IrRrIm.
    $

    Since $ \mathcal{R} > 1 $ gives $ I_r^*, R_r^*, \; {\rm{and}}\; I_m^* > 0 $, we have $ D_1(I_r^*, R_r^*, I_m^*) > 0 $, $ D_2(I_r^*, R_r^*, I_m^*) > 0 $, and $ D_3(I_r^*, R_r^*, I_m^*) > 0 $.

    Thus, the function $ V(I_r, R_r, I_m) $ has the local minimum at the only critical point $ (I_r^*, R_r^*, I_m^*) $.

    So, $ V(I_r, R_r, I_m) $ has the global minimum at $ (I_r^*, R_r^*, I_m^*) $ and $ V(I_r^*, R_r^*, I_m^*) = 0 $.

    $ {\rm{Hence,}}\; V(I_r, R_r, I_m) > 0 \; {\rm{for}}\; (I_r, R_r, I_m)\ne (I_r^*, R_r^*, I_m^*). $

    To complete the proof, we only need to show $ {V'}(I_r, R_r, I_m) < 0 $ for $ (I_r, R_r, I_m)\neq (I_r^*, R_r^*, I_m^*) $.

    From Eq (11),

    $ {V'}(I_r, R_r, I_m) = \left(1-\frac{I_r^*}{I_r}\right){I_r'}+\left(1-\frac{R_r^*}{R_r}\right){R_r'}+\left(1-\frac{I_m^*}{I_m}\right){I_m'}, \; {\rm{here}}, I_r^*, R_r^*, I_m^* > 0 \; . $

    Thus,

    if $ I_r > I_r^*, \; {\rm{then}}\; (1-\frac{I_r^*}{I_r}) > 0 $; so, we need to show $ {I_r'} < 0 \; {\rm{to\; have}}\; \left(1-\frac{I_r^*}{I_r}\right){I_r'} < 0 $ and

    if $ I_r < I_r^*, \; {\rm{then}}\; (1-\frac{I_r^*}{I_r}) < 0 $; so, we need to show $ {I_r'} > 0 \; {\rm{to\; have}}\; \left(1-\frac{I_r^*}{I_r}\right){I_r'} < 0 $.

    Now, from Eq (6),

    $ Ir=λrβmrImNmNr(NrRr)(λrβmrImNmNr+(γr+μr(1pr)))Ir.
    $

    Here, $ \lambda_r\beta_{mr}\frac{I_m}{N_mN_r}+\left(\gamma_r+\mu_r(1-p_r)\right) > 0\; {\rm{and}}\; I_r^* > 0 \; ({\rm{for}}\; \mathcal{R} > 1).\; {\rm{So, }}\; {I_r'} < 0\; {\rm{if}}\; I_r > I_r^*\; {\rm{and}}\; {I_r'} $ $ > 0\; {\rm{if}}\; I_r < I_r^* $. Thus,

    $ (1IrIr)Ir<0forIrIr.
    $
    (12)

    Equation (6) also gives,

    $ Rr=γrIrμrRr.
    $

    Here, $ {R_r^*} > 0 \; ({\rm{for}}\; \mathcal{R} > 1). \; {\rm{So, }}\; {R_r'} < 0 \; {\rm{if}}\; R_r > R_r^*\; {\rm{and}}\; {R_r'} > 0 \; {\rm{if}}\; R_r < R_r^* $. Thus,

    $ (1RrRr)Rr<0forRrRr.
    $
    (13)

    Again, from Eq (6), $ {I'}_m = \lambda_r\beta_{rm}\frac{I_r}{N_r}-\left(\lambda_r\beta_{rm}\frac{I_r}{N_rN_m}+\mu_m(1-p_m)\right)I_m $.

    Here, $ \lambda_r\beta_{rm}\frac{I_r}{N_rN_m}+\mu_m(1-p_m) > 0\; {\rm{and}}\; {I_m^*} > 0 \; ({\rm{for}}\; \mathcal{R} > 1). \; {\rm{So}}\; , {I_m'} < 0 \; {\rm{if}}\; I_m > I_m^*\; {\rm{and}}\; {I_m'} > 0 $ $ \; {\rm{if}}\; I_m < I_m^* $. Thus,

    $ (1ImIm)Im<0forImIm.
    $
    (14)

    Thus, combining Eqs (12)–(14) we have,

    $ {V'}(I_r, R_r, I_m) < 0 \; {\rm{for}}\; (I_r, R_r, I_m)\ne (I_r^*, R_r^*, I_m^*) \; {\rm{when}}\; \mathcal{R} > 1. $

    Finally, if and only if $ \mathcal{R} > 1 $ we have,

    $ V(I_r, R_r, I_m) = 0 \; {\rm {and}}\; {V'}(I_r, R_r, I_m) = 0 \; {\rm{for}}\; (I_r, R_r, I_m) = (I_r^*,R_r^*,I_m^*). $
    $ {\rm{Also,}}\; V(I_r, R_r, I_m) > 0 \; {\rm{and}}\; {V'}(I_r, R_r, I_m) < 0 \; {\rm{for}}\; (I_r, R_r, I_m)\ne(I_r^*,R_r^*,I_m^*). $

    Hence, the endemic equilibrium for the pig-mosquito subsystem is globally stable if and only if $ R_{0e} > 1 $ (since $ \mathcal{R} > 1 \Leftrightarrow R_{0e} > 1 $).

    Proposition 5: The $ DFE $ of the pig–mosquito subsystem is globally stable when $ R_{0e} < 1 $.

    Proof: We define our Lyapunov function as

    $ V(Ir,Rr,Im)=μm(1pm)λrβmrNmIr+Im+(μm(1pm)(μr(1pr)+γr)λrβmrNmλrβrmNr)Rrγr
    $
    (15)
    $ =μm(1pm)λrβmrNmIr+Im+1γr(1R21)Rr
    $

    which gives

    $ V(Ir,Rr,Im)=μm(1pm)λrβmrNmIr+Im+(μm(1pm)(μr(1pr)+γr)λrβmrNmλrβrmNr)Rrγr=μm(1pm)Ir+RrNrμrγr(μm(1pm)(μr(1pr)+γr)λrβmrNmλrβrmNr)RrλrβrmNrImNm=μm(1pm)Ir+RrNrλrβrmNrImNmμrγr(1R21)Rr.
    $
    (16)

    Here we use our alternative threshold quantity $ \mathcal{R} $ (see Eq (8)) to understand if the function $ V $ satisfies all the conditions of a strong Lyapunov.

    Now, the alternative threshold value $ \mathcal{R} $, Eqs (15) and (16) easily verify that

    $ V(0,0,0) = 0 \; {\rm{and}}\; V(I_r, R_r, I_m) > 0 \; {\rm{for}}\; (I_r, R_r, I_m)\ne (0,0,0) \; {\rm{when}}\; \mathcal{R} < 1. $
    $ {\rm{Also,}} \; {V'}(0,0,0) = 0 \; {\rm{and}}\; {V'}(I_r, R_r, I_m) < 0 \; {\rm{for}} \; (I_r, R_r, I_m)\ne (0,0,0) \; {\rm{when}}\; \mathcal{R} < 1. $

    Hence, the $ DFE $ for the pig-mosquito system is globally stable when $ R_{0e} < 1 $ (since $ \mathcal{R} < 1 \Leftrightarrow R_{0e} < 1 $).

    $ SSh=ΛhβlhLc+μh+λhβmhNhImNmESh=λhβmhNhImNmβlhLc+μh+γe+αeSShRSh=λhβmhNhImNmγe(βlhLc+μh)(βlhLc+μh+γe+αe)SShSL1h=βlhLcλhβmhNhImNm+γl1+rl+μhSShEL1h=βlhLcλhβmhNhImNmγl1+γe+rl+μh+αe(1βlhLc+γe+μh+αe+1λhβmhNhImNm+γl+rl+μh)SShRL1h=βlhLcλhβmhNhImNmγeγl1+rl+μh(1(βlhLc+μh)(βlhLc+μh+γe+μe)+1γl1+rl+μh+γe+αe)SShSL2h=rlβlhLc(λhβmhNhImNm+γl2+αl+μh)(λhβmhNhImNm+γl1+rl+μh)SShEL2h=rlEL1h+λhβmhNhImNmSL2hγl2+αl+γe+αe+μhRL2h=rlRL1h+γeEL2hγl2+αl+μhSRh=βlhLc(γl1+rlγl2(λhβmhNhImNm+γl2+μh+αl))(λhβmhNhImNm+μh)(λhβmhNhImNm+γl1+rl+μh)SShERh=γl1EL1h+γl2EL2h+λhβmhNhImNmSRhγe+αe+μhRRh=γl1RL1h+γl2RL2h+γeERhμh
    $
    (17)
    [1] Watson JD, Baker TA, Bell SP, et al. (2008) Molecular Biology of the Gene, Sixth Edition, San Francisco: Pearson Education.
    [2] Alberts B, Johnson A, Lewis J, et al. (2008) Molecular Biology of the Cell, Fifth Edition, New York: Garland Science.
    [3] Nishizaka T, Yagi T, Tanaka Y, et al. (1993) Right-handed rotation of an actin filament in an in vitro motile system. Nature 361: 269-271. doi: 10.1038/361269a0
    [4] Tee YH, Shemesh T, Thiagarajan V, et al. (2015) Cellular chirality arising from the self-organization of the actin cytoskeleton. Nat Cell Biol 17: 445-457. doi: 10.1038/ncb3137
    [5] Bouchet-Marquis C, Zuber B, Glynn AM, et al. (2007) Visualization of cell microtubules in their native state. Biol Cell 99: 45-53. doi: 10.1042/BC20060081
    [6] Sumino Y, Nagai KH, Shikata Y, et al. (2012) Large-scale vortex lattice emerging from collectively moving microtubules. Nature 483: 448-452. doi: 10.1038/nature10874
    [7] Orgel JPRO, Irving TC, Miller A, et al. (2006) Microfibrillar structure of type I collagen in situ. Proc Natl Acad Sci U S A 103: 9001-9005. doi: 10.1073/pnas.0502718103
    [8] Grason GM (2015) Colloquium: Geometry and optimal packing of twisted columns and filaments. Rev Mod Phys 87: 401-419. doi: 10.1103/RevModPhys.87.401
    [9] Wan LQ, Ronaldson K, Guirguis M, et al. (2013) Micropatterning of cells reveals chiral morphogenesis. Stem Cell Res Therapy 4:24. doi: 10.1186/scrt172
    [10] Schiessel H (2003) The physics of chromatin. J Phys Condens Matter 15: R699-R774. doi: 10.1088/0953-8984/15/19/203
    [11] Arya G, Schlick T (2006) Role of histone tails in chromatin folding revealed by a mesoscopic oligonucleosome model. Proc Natl Acad Sci U S A 103: 16236-16241. doi: 10.1073/pnas.0604817103
    [12] Arya G, Zhang Q, Schlick T (2006) Flexible histone tails in a new mesoscopic oligonucleosome model. Biophys J 91: 133-150. doi: 10.1529/biophysj.106.083006
    [13] Müller O, Kepper N, Schöpflin R (2014) Changing chromatin fiber conformation by nucleosome repositioning. Biophys J 107: 2141-2150. doi: 10.1016/j.bpj.2014.09.026
    [14] Maeshima K, Imai R, Tamura S (2014) Chromatin as dynamic 10-nm fibers, Chromosoma 123: 225-237.
    [15] Ricci MA, Manzo C, García-Parajo MF, et al. (2015) Chromatin fibers are formed by heterogeneous groups of nucleosomes in vivo. Cell 160: 1145-1158. doi: 10.1016/j.cell.2015.01.054
    [16] Cherstvy AG, Everaers R (2006) Layering, bundling, and azimuthal orientations in dense phases of nucleosome core particles. J Phys Condens Matter 18: 11429-11442. doi: 10.1088/0953-8984/18/50/003
    [17] Luger K, Mäder AW, Richmond RK, et al. (1997) Crystal structure of the nucleosome core particle at 2.8 Å resolution. Nature 389: 251-260.
    [18] Richmond TJ, Davey CA (2003) The structure of DNA in the nucleosome core. Nature 423: 145-150. doi: 10.1038/nature01595
    [19] Schiessel H (2006) The nucleosome: A transparent, slippery, sticky and yet stable DNA-protein complex. Eur Phys J E 19: 251-262. doi: 10.1140/epje/i2005-10049-y
    [20] Tolstorukov MY, Colasanti AV, McCandlish DM, et al. (2007) A Novel Roll-and-Slide Mechanism of DNA Folding in Chromatin: Implications for Nucleosome Positioning. J Mol Biol 371: 725-738. doi: 10.1016/j.jmb.2007.05.048
    [21] Olson WK, Zhurkin VB (2011) Working the kinks out of nucleosomal DNA. Curr Opin Struct Biol 21: 348-357. doi: 10.1016/j.sbi.2011.03.006
    [22] Higuchi Y, Sakaue T, Yoshikawa K (2010) Torsional effect on the wrapping transition of a semiflexible polymer around a core as a model of nucleosome. Phys Rev E 82: 031909. doi: 10.1103/PhysRevE.82.031909
    [23] Boroudjerdi H, Naji A, Netz R (2014) Global analysis of the ground-state wrapping conformation of a charged polymer on an oppositely charged nano-sphere. Eur Phys J E 37: 21. doi: 10.1140/epje/i2014-14021-6
    [24] Korolev N, Berezhnoy NV, Dong EK, et al. (2012) A universal description for the experimental behavior of salt-(in)dependent oligocation-induced DNA condensation. Nucleic Acids Res 40: 2807-2821. doi: 10.1093/nar/gks214
    [25] Kenzaki H, Takada S (2015) Partial unwrapping and histone tail dynamics in nucleosome revealed by coarse-grained molecular simulations. PLoS Comput Biol 8: e1004443.
    [26] Walter JC, Baiesi M, Barkema G, et al. (2013) Unwinding relaxation dynamics of polymers. Phys Rev Lett 110: 068301. doi: 10.1103/PhysRevLett.110.068301
    [27] Walter JC, Baiesi M, Carlon E, et al. (2014) Unwinding dynamics of a helically wrapped polymer. Macromolecules 47: 4840-4846. doi: 10.1021/ma500635h
    [28] Zheng M, Jagota A, Semke ED, et al. (2003) DNA-assisted dispersion and separation of carbon nanotubes. Nature Mater 2: 338-342. doi: 10.1038/nmat877
    [29] Nakashima N, Okuzono S, Murakami H, et al. (2003) DNA dissolves single-walled carbon nanotubes in water. Chem Lett 32: 456-457. doi: 10.1246/cl.2003.456
    [30] Heller DA, Jeng ES, Yeung TK, et al. (2006) Optical detection of DNA conformational polymorphism on single-walled carbon nanotubes. Science 311: 508-511. doi: 10.1126/science.1120792
    [31] Manohar S, Tang T, Jagota A (2007) Structure of homopolymer DNA-CNT hybrids. J Phys Chem C 111: 17835-17845. doi: 10.1021/jp071316x
    [32] Alegret N, Santos E, Rodríguez-Fortea A, et al. (2012) Disruption of small double stranded DNA molecules on carbon nanotubes: A molecular dynamics study. Chem Phys Lett 525-526: 120-124. doi: 10.1016/j.cplett.2011.12.077
    [33] Marko JF (1997) Supercoiled and braided DNA under tension. Phys Rev E 55: 1758-1772. doi: 10.1103/PhysRevE.55.1758
    [34] Stone MD, Bryant Z, Crisona NJ, et al. (2003) Chirality sensing by Escherichia coli topoisomerase IV and the mechanism of type II topoisomerases. Proc Natl Acad Sci U S A 100: 8654-8659. doi: 10.1073/pnas.1133178100
    [35] Charvin G, Vologodskii A, Bensimon D, et al. (2005) Braiding DNA: Experiments, simulations, and models. Biophys J 88: 4124-4136. doi: 10.1529/biophysj.104.056945
    [36] Timsit Y, Várnai P (2010) Helical chirality: a link between local interactions and global topology in DNA. PLoS One 5: e9326-e9326.
    [37] Várnai P, Timsit Y (2010) Differential stability of DNA crossovers in solution mediated by divalent cations. Nucleic Acids Res 38: 4163-4172. doi: 10.1093/nar/gkq150
    [38] Timsit Y (2013) DNA self-assembly: from chirality to evolution. Int J Mol Sci 14: 8252-8270. doi: 10.3390/ijms14048252
    [39] Cherstvy AG (2008) DNA cholesteric phases: The role of DNA molecular chirality and DNA-DNA electrostatic interactions. J Phys Chem B 112: 12585-12595. doi: 10.1021/jp801220p
    [40] Kornyshev AA, Leikin S (2001) Sequence recognition in the pairing of DNA duplexes. Phys Rev Lett 86: 3666-3669. doi: 10.1103/PhysRevLett.86.3666
    [41] Kornyshev AA, Leikin S, Lee DJ, et al. (2007) Structure and interactions of biological helices. Rev Mod Phys 79: 943-996. doi: 10.1103/RevModPhys.79.943
    [42] Cortini R, Kornyshev AA, Lee DJ, et al. (2011) Electrostatic braiding and homologous pairing of DNA double helices. Biophys J 101: 875-884. doi: 10.1016/j.bpj.2011.06.058
    [43] Cortini R, Lee DJ, Kornyshev AA (2012) Chiral electrostatics breaks the mirror symmetry of DNA supercoiling. J Phys Condens Matter 24: 162203. doi: 10.1088/0953-8984/24/16/162203
    [44] Lee DJ, Cortini R, Korte AP, et al. (2013) Chiral effects in dual-DNA braiding. Soft Matt 41: 9833-9848.
    [45] Lee DJ (2014) Undulations in a weakly interacting mechanically generated molecular braid under tension. J Phys Condens Matter 26: 245101. doi: 10.1088/0953-8984/26/24/245101
    [46] Lee DJ (2015) Collapse and coexistence for a molecular braid with an attractive interaction component subject to mechanical forces. J Phys Condens Matter 27: 145101. doi: 10.1088/0953-8984/27/14/145101
    [47] Gore J, Bryant Z, Nöllmann M, et al. (2006) DNA overwinds when stretched. Nature 442: 836-839. doi: 10.1038/nature04974
    [48] Lionnet T, Joubaud S, Lavery R, et al. Wringing out DNA. (2006) Phys Rev Lett 96: 178102.
    [49] Strick TR, Allemand JF, Bensimon D, et al. (1998) Behavior of Supercoiled DNA. Biophys J 74: 2016-2028. doi: 10.1016/S0006-3495(98)77908-1
    [50] Allemand JF, Bensimon D, Lavery R, et al. (1998) Stretched and overwound DNA forms a Pauling-like structure with exposed bases. Proc Natl Acad Sci U S A 95: 14152-14157. doi: 10.1073/pnas.95.24.14152
    [51] Neukirch S (2004) Extracting DNA Twist Rigidity from Experimental Supercoiling Data. Phys Rev Lett 93: 198107. doi: 10.1103/PhysRevLett.93.198107
    [52] Besteman K, Hage S, Dekker NH, et al. (2007) Role of Tension and Twist in Single-Molecule DNA Condensation. Phys Rev Lett 98: 058103. doi: 10.1103/PhysRevLett.98.058103
    [53] Marko JF (2007) Torque and dynamics of linking number relaxation in stretched supercoiled DNA. Phys Rev E 76: 021926. doi: 10.1103/PhysRevE.76.021926
    [54] Marko JF, Siggia ED (1994) Bending and twisting elasticity of DNA, Macromolecules 27: 981-988.
    [55] Kulić IM, Andrienko D, Deserno M (2004) Twist-bend instability for toroidal DNA condensates. Europhys Lett 67: 418-424. doi: 10.1209/epl/i2004-10076-x
    [56] Yanao T, Yoshikawa K (2014) Chiral symmetry breaking of a double-stranded helical chain through bend-writhe coupling. Phys Rev E 89: 062713.
    [57] Kornyshev AA, Leikin S (2013) Helical structure determines different susceptibilities of dsDNA, dsRNA, and tsDNA to counterion-induced condensation. Biophys J 104: 2031-2041. doi: 10.1016/j.bpj.2013.03.033
    [58] Lee D (2014) Self-consistent treatment of electrostatics in molecular DNA braiding through external forces. Phys Rev E 89: 062711.
    [59] Savelyev A, Papoian GA (2009) Molecular renormalization group coarse-graining of polymer chains: Application to double-stranded DNA. Biophys J 96: 4044-4052. doi: 10.1016/j.bpj.2009.02.067
    [60] Knotts AT, Rathore N, Schwartz DC, et al. (2007) A coarse grain model for DNA. J Chem Phys 126: 084901. doi: 10.1063/1.2431804
    [61] Ouldridge TE, Louis AA, Doye JPK (2011) Structural, mechanical, and thermodynamic properties of a coarse-grained DNA model. J Chem Phys 134: 085101. doi: 10.1063/1.3552946
    [62] Snodin BEK, Randisi F, Mosayebi M, et al. (2015) Introducing improved structural properties and salt dependence into a coarse-grained model of DNA. J Chem Phys 142: 234901. doi: 10.1063/1.4921957
    [63] Patriciu A, Chirikjian GS, Pappu RV (2004) Analysis of the conformational dependence of mass-metric tensor determinants in serial polymers with constraints. J Chem Phys 121: 12708-12720. doi: 10.1063/1.1821492
    [64] Yanao T, Yoshikawa K (2008) Elastic origin of chiral selection in DNA wrapping. Phys Rev E 77: 021904. doi: 10.1103/PhysRevE.77.021904
    [65] Geggier S, Vologodskii A (2010) Sequence dependence of DNA bending rigidity. Proc Natl Acad Sci U S A 107: 15421-15426. doi: 10.1073/pnas.1004809107
    [66] Metropolis N, Rosenbluth AW, Rosenbluth MN, et al. (1953) Equation of State Calculations by Fast Computing Machines. J Chem Phys 21: 1087-1092. doi: 10.1063/1.1699114
    [67] Allen MP, Tildesley DJ (1987) Computer simulation of liquids, Oxford: Oxford Science Publications.
    [68] Humphrey W, Dalke A, Schulten K (1996) VMD: Visual Molecular Dynamics. J Mol Graphics 14: 33-38. doi: 10.1016/0263-7855(96)00018-5
  • This article has been cited by:

    1. Muhammad U. Riandi, Mara Ipa, Tri Wahono, Joni Hendri, Triwibowo A. Garjito, 2023, Chapter 83, 978-94-6463-111-1, 914, 10.2991/978-94-6463-112-8_83
    2. Troy A. Laidlow, Erin S. Johnston, Ruth N. Zadoks, Michael Walsh, Mafalda Viana, Kerrie E. Wiley, Balbir B. Singh, Francesco Baldini, Himani Dhanze, Cameron Webb, Victoria J. Brookes, Fedor Korennoy, Scoping Review of Japanese Encephalitis Virus Transmission Models, 2025, 2025, 1865-1674, 10.1155/tbed/9880670
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6695) PDF downloads(1147) Cited by(1)

Figures and Tables

Figures(13)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog