Development of a network based model to simulate the between-farm transmission of the porcine reproductive and respiratory syndrome virus.

Contact structure within a population can significantly affect the outcomes of infectious disease spread models. The objective of this study was to develop a network based simulation model for the between-farm spread of porcine reproductive and respiratory syndrome virus to assess the impact of contact structure on between-farm transmission of PRRS virus. For these farm level models, a hypothetical population of 500 swine farms following a multistage production system was used. The contact rates between farms were based on a study analyzing movement of pigs in Canada, while disease spread parameters were extracted from published literature. Eighteen distinct scenarios were designed and simulated by varying the mode of transmission (direct versus direct and indirect contact), type of index herd (farrowing, nursery and finishing), and the presumed network structures among swine farms (random, scale-free and small-world). PRRS virus was seeded in a randomly selected farm and 500 iterations of each scenario were simulated for 52 weeks. The median epidemic size by the end of the simulated period and percentage die-out for each scenario, were the key outcomes captured. Scenarios with scale-free network models resulted in the largest epidemic sizes, while scenarios with random and small-world network models resulted in smaller and similar epidemic sizes. Similarly, stochastic die-out percentage was least for scenarios with scale-free networks followed by random and small-world networks. Findings of the study indicated that incorporating network structures among the swine farms had a considerable impact on the spread of PRRS virus, highlighting the importance of understanding and incorporating realistic contact structures when developing infectious disease spread models for similar populations.

Discontinuous molecular dynamics (DMD) study of heteropolymer collapse in an explicit solvent.

In this study, we employ the discontinuous molecular dynamics simulation method to investigate the collapse properties of a single heteropolymer chain in an explicit solvent. Solvent density ρ, fraction of hydrophobic monomers n H (defined as the ratio of the number of hydrophobic monomers to the total number of monomers) and a hydrophobicity parameter λ (which controls the energy mismatch between the monomers and solvent particles) were systematically varied to examine their role in polymer collapse. The average static structure factor of the polymer was used to find the so-called θ-point characterizing the state of an ideal chain. Phase diagrams of ρ versus λ for the coil-globule transition were mapped out for different values of n H. Increasing the fraction of hydrophobic monomers n H, solvent density ρ, and hydrophobicity parameter λ were all shown to aid in stabilizing the globule phase. In an effort to explore scaling behaviour of the coil-globule phase diagram as a function of n H, and to investigate whether the phase boundaries for different n H collapsed on to one universal curve, we rescaled λ by n H (δ) λ; we determined δ = 1.72, in contrast to mean-field predictions of δ = 2.0.

Theoretical study of solvent effects on the coil-globule transition.

The coil-globule transition of a polymer in a solvent has been studied using Monte Carlo simulations of a single chain subject to intramolecular interactions as well as a solvent-mediated effective potential. This solvation potential was calculated using several different theoretical approaches for two simple polymer/solvent models, each employing hard-sphere chains and hard-sphere solvent particles as well as attractive square-well potentials between some interaction sites. For each model, collapse is driven by variation in a parameter which changes the energy mismatch between monomers and solvent particles. The solvation potentials were calculated using two fundamentally different methodologies, each designed to predict the conformational behavior of polymers in solution: (1) the polymer reference interaction site model (PRISM) theory and (2) a many-body solvation potential (MBSP) based on scaled particle theory introduced by Grayce [J. Chem. Phys. 106, 5171 (1997)]. For the PRISM calculations, two well-studied solvation monomer-monomer pair potentials were employed, each distinguished by the closure relation used in its derivation: (i) a hypernetted-chain (HNC)-type potential and (ii) a Percus-Yevick (PY)-type potential. The theoretical predictions were each compared to results obtained from explicit-solvent discontinuous molecular dynamics simulations on the same polymer/solvent model systems [J. Chem. Phys. 125, 194904 (2006)]. In each case, the variation in the coil-globule transition properties with solvent density is mostly qualitatively correct, though the quantitative agreement between the theory and prediction is typically poor. The HNC-type potential yields results that are more qualitatively consistent with simulation. The conformational behavior of the polymer upon collapse predicted by the MBSP approach is quantitatively correct for low and moderate solvent densities but is increasingly less accurate for higher densities. At high solvent densities, the PRISM-HNC and MBSP approaches tend to overestimate, while the PRISM-PY approach underestimates the tendency of the solvent to drive polymer collapse.

Discontinuous molecular dynamics for rigid bodies: applications.

Event-driven molecular dynamics simulations are carried out on two rigid-body systems which differ in the symmetry of their molecular mass distributions. First, simulations of methane in which the molecules interact via discontinuous potentials are compared with simulations in which the molecules interact through standard continuous Lennard-Jones potentials. It is shown that under similar conditions of temperature and pressure, the rigid discontinuous molecular dynamics method reproduces the essential dynamical and structural features found in continuous-potential simulations at both gas and liquid densities. Moreover, the discontinuous molecular dynamics approach is demonstrated to be between 3 and 100 times more efficient than the standard molecular dynamics method depending on the specific conditions of the simulation. The rigid discontinuous molecular dynamics method is also applied to a discontinuous-potential model of a liquid composed of rigid benzene molecules, and equilibrium and dynamical properties are shown to be in qualitative agreement with more detailed continuous-potential models of benzene. The few qualitative differences in the angular dynamics of the two models are related to the relatively crude treatment of variations in the discontinuous repulsive interactions as one benzene molecule rotates by another.

Discontinuous molecular dynamics for semiflexible and rigid bodies.

A general framework for performing event-driven simulations of systems with semiflexible or rigid bodies interacting under impulsive forces is outlined. The method consists of specifying a means of computing the free evolution of constrained motion, evaluating the times at which interactions occur, and determining the consequences of interactions on subsequent motion. Algorithms for computing the times of interaction events and carrying out efficient event-driven simulations are discussed. The semiflexible case and the rigid case differ qualitatively in that the free motion of a rigid body can be computed analytically and need not be integrated numerically.

Discontinuous molecular dynamics simulation study of polymer collapse.

Discontinuous molecular dynamics simulations were used to study the coil-globule transition of a polymer in an explicit solvent. Two different versions of the model were employed, which are differentiated by the nature of monomer-solvent, solvent-solvent, and nonbonded monomer-monomer interactions. For each case, a model parameter lambda determines the degree of hydrophobicity of the monomers by controlling the degree of energy mismatch between the monomers and solvent particles. We consider a lambda-driven coil-globule transition at constant temperature. The simulations are used to calculate average static structure factors, which are then used to determine the scaling exponents of the system in order to determine the theta-point values lambda(theta) separating the coil from the globule state. For each model we construct coil-globule phase diagrams in terms of lambda and the particle density rho. Additionally, we explore for each model the effects of varying the range of the attractive interactions on the phase boundary separating the coil and globule phases. The results are analyzed in terms of a simple Flory-type theory of the collapse transition.