A systems perspective of host-pathogen interactions: predicting disease outcome in tuberculosis

The complex web of interactions between the host immune system and the pathogen determines the outcome of any infection. A computational model of this interaction network, which encodes complex interplay among host and bacterial components, forms a useful basis for improving the understanding of pathogenesis, in filling knowledge gaps and consequently to identify strategies to counter the disease. We have built an extensive model of the Mycobacterium tuberculosis host-pathogen interactome, consisting of 75 nodes corresponding to host and pathogen molecules, cells, cellular states or processes. Vaccination effects, clearance efficiencies due to drugs and growth rates have also been encoded in the model. The system is modelled as a Boolean network. Virtual deletion experiments, multiple parameter scans and analysis of the system's response to perturbations, indicate that disabling processes such as phagocytosis and phagolysosome fusion or cytokines such as TNF-alpha and IFN-gamma, greatly impaired bacterial clearance, while removing cytokines such as IL-10 alongside bacterial defence proteins such as SapM greatly favour clearance. Simulations indicate a high propensity of the pathogen to persist under different conditions.

A systems perspective of host-pathogen interactions: predicting disease outcome in tuberculosis

The complex web of interactions between the host immune system and the pathogen determines the outcome of any infection. A computational model of this interaction network, which encodes complex interplay among host and bacterial components, forms a useful basis for improving the understanding of pathogenesis, in filling knowledge gaps and consequently to identify strategies to counter the disease. We have built an extensive model of the Mycobacterium tuberculosis host-pathogen interactome, consisting of 75 nodes corresponding to host and pathogen molecules, cells, cellular states or processes. Vaccination effects, clearance efficiencies due to drugs and growth rates have also been encoded in the model. The system is modelled as a Boolean network. Virtual deletion experiments, multiple parameter scans and analysis of the system's response to perturbations, indicate that disabling processes such as phagocytosis and phagolysosome fusion or cytokines such as TNF-alpha and IFN-gamma, greatly impaired bacterial clearance, while removing cytokines such as IL-10 alongside bacterial defence proteins such as SapM greatly favour clearance. Simulations indicate a high propensity of the pathogen to persist under different conditions.

Temporal associations in fig-wasp-ant interactions: diel and phenological patterns

In a complex multitrophic plant-animal interaction system in which there are direct and indirect interactions between species, comprehending the dynamics of these multiple partners is very important for an understanding of how the system is structured. We investigated the plant Ficus racemosa L. (Moraceae) and its community of obligatory mutualistic and parasitic fig wasps (Hymenoptera: Chalcidoidea) that develop within the fig inflorescence or syconium, as well as their interaction with opportunistic ants. We focused on temporal resource partitioning among members of the fig wasp community over the development cycle of the fig syconia during which wasp oviposition and development occur and we studied the activity rhythm of the ants associated with this community. We found that the seven members of the wasp community partitioned their oviposition across fig syconium development phenology and showed interspecific variation in activity across the day-night cycle. The wasps presented a distinct sequence in their arrival at fig syconia for oviposition, with the parasitoid wasps following the galling wasps. Although fig wasps are known to be largely diurnal, we documented night oviposition in several fig wasp species for the first time. Ant activity on the fig syconia was correlated with wasp activity and was dependent on whether the ants were predatory or trophobiont-tending species; only numbers of predatory ants increased during peak arrivals of the wasps.

Taking Multiple Infections of Cells and Recombination into Account Leads to Small Within-Host Effective-Population-Size Estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N-e, are widely varying. Models assuming HIV-1 evolution to be neutral estimate N-e similar to 10(2)-10(4), smaller than the inverse mutation rate of HIV-1 (similar to 10(5)), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N-e>10(5), suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N-e similar to 10(3)-10(4), implying predominantly stochastic evolution. Interestingly, we find that N-e and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N-e>10(5) reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N-e similar to 10(3)-10(4) may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.

Interactions of interfacial arc cracks

The interaction of two interfacial arc cracks around a circular elastic inclusion embedded in an elastic matrix is examined. New results for stress intensity factors for a pair of interacting cracks are derived for a concentrated force acting in the matrix. For verifying the point load solutions, stress intensity factors under uniform loading are obtained by superposing point force results. For achieving this objective, a general method for generating desired stress fields inside a test region using point loads is described. The energetics of two interacting interfacial arc cracks is discussed in order to shed more light on the debonding of hard or soft inclusions from the matrix. The analysis based on complex variables is developed in a general way to handle the interactions of multiple interfacial arc cracks/straight cracks.

Modeling Architecture-OS interactions using Layered Queuing Network Models

The prevalent virtualization technologies provide QoS support within the software layers of the virtual machine monitor(VMM) or the operating system of the virtual machine(VM). The QoS features are mostly provided as extensions to the existing software used for accessing the I/O device because of which the applications sharing the I/O device experience loss of performance due to crosstalk effects or usable bandwidth. In this paper we examine the NIC sharing effects across VMs on a Xen virtualized server and present an alternate paradigm that improves the shared bandwidth and reduces the crosstalk effect on the VMs. We implement the proposed hardwaresoftware changes in a layered queuing network (LQN) model and use simulation techniques to evaluate the architecture. We find that simple changes in the device architecture and associated system software lead to application throughput improvement of up to 60%. The architecture also enables finer QoS controls at device level and increases the scalability of device sharing across multiple virtual machines. We find that the performance improvement derived using LQN model is comparable to that reported by similar but real implementations.

Computational Analysis of Flow through a Multiple Nozzle Driven Laser Cavity and Diffuser

Computations have been carried out for simulating supersonic flow through a set of converging-diverging nozzles with their expanding jets forming a laser cavity and flow patterns through diffusers, past the cavity. A thorough numerical investigation with 3-D RANS code is carried out to capture the flow distribution which comprises of shock patterns and multiple supersonic jet interactions. The analysis of pressure recovery characteristics during the flow through the diffusers is an important parameter of the simulation and is critical for the performance of the laser device. The results of the computation have shown a close agreement with the experimentally measured parameters as well as other established results indicating that the flow analysis done is found to be satisfactory.

A model investigation of vegetation-atmosphere interactions on a millennial timescale

A terrestrial biosphere model with dynamic vegetation capability, Integrated Biosphere Simulator (IBIS2), coupled to the NCAR Community Atmosphere Model (CAM2) is used to investigate the multiple climate-forest equilibrium states of the climate system. A 1000-year control simulation and another 1000-year land cover change simulation that consisted of global deforestation for 100 years followed by re-growth of forests for the subsequent 900 years were performed. After several centuries of interactive climate-vegetation dynamics, the land cover change simulation converged to essentially the same climate state as the control simulation. However, the climate system takes about a millennium to reach the control forest state. In the absence of deep ocean feedbacks in our model, the millennial time scale for converging to the original climate state is dictated by long time scales of the vegetation dynamics in the northern high latitudes. Our idealized modeling study suggests that the equilibrium state reached after complete global deforestation followed by re-growth of forests is unlikely to be distinguishable from the control climate. The real world, however, could have multiple climate-forest states since our modeling study is unlikely to have represented all the essential ecological processes (e. g. altered fire regimes, seed sources and seedling establishment dynamics) for the reestablishment of major biomes.

NMR Analysis of Cross Strand Aromatic Interactions in an 8 Residue Hairpin and a 14 Residue Three Stranded beta-Sheet Peptide

Cross strand aromatic interactions between a facing pair of phenylalanine residues in antiparallel beta-sheet structures have been probed using two structurally defined model peptides. The octapeptide Boc-(LFVPPLFV)-P-D-P-L-OMe (peptide 1) favors the beta-hairpin conformation nucleated by the type II' beta-turn formed by the (D)Pro-(L)Pro segment, placing Phe2 and Phe7 side chains in proximity. Two centrally positioned (D)Pro-(L)Pro segments facilitate the three stranded beta-sheet formation in the 14 residue peptide Boc-LFV(D)P(L)PLFVA(D)P(L)PLFV-OMe (peptide 2) in which the Phe2/Phe7 orientations are similar to that in the octapeptide. The anticipated folded conformations of peptides 1 and 2 are established by the delineation of intramolecularly hydrogen bonded NH groups and by the observation of specific cross strand NOEs. The observation of ring current shifted aromatic protons is a diagnostic of close approach of the Phe2 and Phe7 side chains. Specific assignment of aromatic proton resonances using HSQC and HSQC-TOCSY methods allow an analysis of interproton NOEs between the spatially proximate aromatic rings. This approach facilitates specific assignments in systems containing multiple aromatic rings in spectra at natural abundance. Evidence is presented for a dynamic process which invokes a correlated conformational change about the C-alpha-C-beta(chi(1)) bond for the pair of interacting Phe residues. NMR results suggest that aromatic ring orientations observed in crystals are maintained in solution. Anomalous temperature dependence of ring current induced proton chemical shifts suggests that solvophobic effects may facilitate aromatic ring clustering in apolar solvents.

The Role of Weak Interactions in the Phase Transition and Distinct Mechanical Behavior of Two Structurally Similar Caffeine Co-crystal Polymorphs Studied by Nanoindentation

Although weak interactions, such as C-H center dot center dot center dot O and pi-stacking, are generally considered to be insignificant, it is their reorganization that holds the key for many a solid-state phenomenon, such as phase transitions, plastic deformation, elastic flexibility, and mechanochromic luminescence in solid-state fluorophores. Despite this, the role of weak interactions in these dynamic phenomena is poorly understood. In this study, we investigate two co-crystal polymorphs of caffeine:4-chloro-3-nitrobenzoic acid, which have close structural similarity (2D layered structures), but surprisingly show distinct mechanical behavior. Form I is brittle, but shows shear-induced phase instability and, upon grinding, converts to Form II, which is soft and plastically shearable. This observation is in contrast to those reported in earlier studies on aspirin, wherein the metastable drug forms are softer and convert to stable and harder forms upon stressing To establish a molecular level understanding, have investigated the two co-crystal polymorphs I and II by single crystal X-ray diffraction, nanoindentation to quantify mechanical properties, and theoretical calculations. The lower hardness (from nanoindentation) and smooth potential surfaces (from theoretical studies) for shearing of layers in Form II allowed us to rationalize the role of stronger intralayer (sp(2))C-H center dot center dot center dot O and nonspecific interlayer pi-stacking interactions in the structure of II. Although the Form I also possesses the same type of interactions, its strength is clearly opposite, that is, weaker intralayer (sp(3))C-H center dot center dot center dot O and specific interlayer pi-stacking interactions. Hence, Form I is harder than Form IL Theoretical calculations and indentation on (111) of Form I suggested the low resistance of this face to mechanical stress; thus, Form I converts to II upon mechanical action. Hence, our approach demonstrates the usefulness of multiple techniques for establishing the role of weak noncovalent interactions in solid-state dynamic phenomena, such as stress induced phase transformation, and hence is important in the context of solid-state pharmaceutical chemistry and crystal engineering.

Conformational Diversity in Contryphans from Conus Venom: cis-trans Isomerisation and Aromatic/Proline Interactions in the 23-Membered Ring of a 7-Residue Peptide Disulfide Loop

Conformational diversity or shapeshifting in cyclic peptide natural products can, in principle, confer a single molecular entity with the property of binding to multiple receptors. Conformational equilibria have been probed in the contryphans, which are peptides derived from Conus venom possessing a 23-membered cyclic disulfide moiety. The natural sequences derived from Conus inscriptus, GCV(D)LYPWC* (In936) and Conus loroisii, GCP(D)WDPWC* (Lo959) differ in the number of proline residues within the macrocyclic ring. Structural characterisation of distinct conformational states arising from cis-trans equilibria about Xxx-Pro bonds is reported. Isomerisation about the C2-P3 bond is observed in the case of Lo959 and about the Y5-P6 bond in In936. Evidence is presented for as many as four distinct species in the case of the synthetic analogue V3P In936. The Tyr-Pro-Trp segment in In936 is characterised by distinct sidechain orientations as a consequence of aromatic/proline interactions as evidenced by specific sidechain-sidechain nuclear Overhauser effects and ring current shifted proton chemical shifts. Molecular dynamics simulations suggest that Tyr5 and Trp7 sidechain conformations are correlated and depend on the geometry of the Xxx-Pro bond. Thermodynamic parameters are derived for the cis trans equilibrium for In936. Studies on synthetic analogues provide insights into the role of sequence effects in modulating isomerisation about Xxx-Pro bonds.

Ovarian development in a primitively eusocial wasp: Social interactions affect behaviorally dominant and subordinate wasps in opposite directions relative to solitary females

In many primitively eusocial wasp species new nests are founded either by a single female or by a small group of females. In the single foundress nests, the lone female develops her ovaries, lays eggs as well as tends her brood. In multiple foundress nests social interactions, especially dominance-subordinate interactions, result in only one `dominant' female developing her ovaries and laying eggs. Ovaries of the remaining `subordinate' cofoundresses remain suppressed and these individuals function as workers and tend the dominant's brood. Using the tropical, primitively eusocial polistine wasp Ropalidia marginata and by comparing wasps held in isolation and those kept as pairs in the laboratory, we demonstrate that social interactions affect ovarian development of dominant and subordinate wasps among the pairs in opposite directions, suppressing the ovaries of the subordinate member of the pair below that of solitary wasps and boosting the ovaries of dominant member of the pair above that of solitary females. In addition to being of physiological interest, such mirror image effects of aggression on the ovaries of the aggressors and their victims, suggest yet another mechanism by which subordinates can enhance their indirect fitness and facilitate the evolution of worker behavior by kin selection. (C) 2014 Elsevier B.V. All rights reserved.

Optical properties of an inhomogeneously broadened Lambda V-system with multiple excited states

We present a theoretical model using a density matrix approach to show the influence of multiple excited states on the optical properties of an inhomogeneously broadened Lambda V-system of the Rb-87 D2 line. These closely spaced multiple excited states cause asymmetry in absorption and dispersion profiles. We observe the reduced absorption profiles, due to dressed state interactions of the applied electromagnetic fields, which results the Mollow sideband-like transparency windows. In a room temperature vapor, we obtain a narrow enhanced absorption and steep positive dispersion at the line center when the strengths of control and pump fields are equal. Here, we show how the probe transmittance varies when it passes through the atomic medium. We also discuss the transient behavior of our system which agrees well with the corresponding absorption and dispersion profiles. This study has potential applications in controllability of group velocity, and for optical and quantum information processing.

Quick re-introduction of selective scalar interactions in a pure-shift NMR spectrum

A new 1D NMR experiment cited as `Quick G-SERF', which re-introduces selective proton-proton scalar interactions in a pure shift spectrum during real time data acquisition, is reported. The method provides information on multiple proton-proton couplings from a single experiment, analogous to the 2D G-SERF technique, while significantly shortening the experimental time by 1-2 orders of magnitude due to reduced dimension and enhanced sensitivity.

The Self-Association of Graphane Is Driven by London Dispersion and Enhanced Orbital Interactions

We investigated the nature of the cohesive energy between graphane sheets via multiple CH center dot center dot center dot HC interactions, using density functional theory (DFT) including dispersion correction (Grimmes D3 approach) computations of n]graphane sigma dimers (n = 6-73). For comparison, we also evaluated the binding between graphene sheets that display prototypical pi/pi interactions. The results were analyzed using the block-localized wave function (BLW) method, which is a variant of ab initio valence bond (VB) theory. BLW interprets the intermolecular interactions in terms of frozen interaction energy (Delta E-F) composed of electrostatic and Pauli repulsion interactions, polarization (Delta E-pol), charge-transfer interaction (Delta E-CT), and dispersion effects (Delta E-disp). The BLW analysis reveals that the cohesive energy between graphane sheets is dominated by two stabilizing effects, namely intermolecular London dispersion and two-way charge transfer energy due to the sigma CH -> sigma*(HC) interactions. The shift of the electron density around the nonpolar covalent C-H bonds involved in the intermolecular interaction decreases the C-H bond lengths uniformly by 0.001 angstrom. The Delta E-CT term, which accounts for similar to 15% of the total binding energy, results in the accumulation of electron density in the interface area between two layers. This accumulated electron density thus acts as an electronic glue for the graphane layers and constitutes an important driving force in the self-association and stability of graphane under ambient conditions. Similarly, the double faced adhesive tape style of charge transfer interactions was also observed among graphene sheets in which it accounts for similar to 18% of the total binding energy. The binding energy between graphane sheets is additive and can be expressed as a sum of CH center dot center dot center dot HC interactions, or as a function of the number of C-H bonds.