A generic individual-based spatially explicit model as a novel tool for investigating insect-plant interactions: A case study of the behavioural ecology of frugivorous Tephritidae

Computational modelling of mechanisms underlying processes in the real world can be of great value in understanding complex biological behaviours. Uptake in general biology and ecology has been rapid. However, it often requires specific data sets that are overly costly in time and resources to collect. The aim of the current study was to test whether a generic behavioural ecology model constructed using published data could give realistic outputs for individual species. An individual-based model was developed using the Pattern-Oriented Modelling (POM) strategy and protocol, based on behavioural rules associated with insect movement choices. Frugivorous Tephritidae (fruit flies) were chosen because of economic significance in global agriculture and the multiple published data sets available for a range of species. The Queensland fruit fly (Qfly), Bactrocera tryoni, was identified as a suitable individual species for testing. Plant canopies with modified architecture were used to run predictive simulations. A field study was then conducted to validate our model predictions on how plant architecture affects fruit flies’ behaviours. Characteristics of plant architecture such as different shapes, e.g., closed-canopy and vase-shaped, affected fly movement patterns and time spent on host fruit. The number of visits to host fruit also differed between the edge and centre in closed-canopy plants. Compared to plant architecture, host fruit has less contribution to effects on flies’ movement patterns. The results from this model, combined with our field study and published empirical data suggest that placing fly traps in the upper canopy at the edge should work best. Such a modelling approach allows rapid testing of ideas about organismal interactions with environmental substrates in silico rather than in vivo, to generate new perspectives. Using published data provides a saving in time and resources. Adjustments for specific questions can be achieved by refinement of parameters based on targeted experiments.

sRNAscanner: A Computational Tool for Intergenic Small RNA Detection in Bacterial Genomes

Background:Bacterial non-coding small RNAs (sRNAs) have attracted considerable attention due to their ubiquitous nature and contribution to numerous cellular processes including survival, adaptation and pathogenesis. Existing computational approaches for identifying bacterial sRNAs demonstrate varying levels of success and there remains considerable room for improvement. Methodology/Principal Findings: Here we have proposed a transcriptional signal-based computational method to identify intergenic sRNA transcriptional units (TUs) in completely sequenced bacterial genomes. Our sRNAscanner tool uses position weight matrices derived from experimentally defined E. coli K-12 MG1655 sRNA promoter and rho-independent terminator signals to identify intergenic sRNA TUs through sliding window based genome scans. Analysis of genomes representative of twelve species suggested that sRNAscanner demonstrated equivalent sensitivity to sRNAPredict2, the best performing bioinformatics tool available presently. However, each algorithm yielded substantial numbers of known and uncharacterized hits that were unique to one or the other tool only. sRNAscanner identified 118 novel putative intergenic sRNA genes in Salmonella enterica Typhimurium LT2, none of which were flagged by sRNAPredict2. Candidate sRNA locations were compared with available deep sequencing libraries derived from Hfq-co-immunoprecipitated RNA purified from a second Typhimurium strain (Sittka et al. (2008) PLoS Genetics 4: e1000163). Sixteen potential novel sRNAs computationally predicted and detected in deep sequencing libraries were selected for experimental validation by Northern analysis using total RNA isolated from bacteria grown under eleven different growth conditions. RNA bands of expected sizes were detected in Northern blots for six of the examined candidates. Furthermore, the 5'-ends of these six Northern-supported sRNA candidates were successfully mapped using 5'-RACE analysis. Conclusions/Significance: We have developed, computationally examined and experimentally validated the sRNAscanner algorithm. Data derived from this study has successfully identified six novel S. Typhimurium sRNA genes. In addition, the computational specificity analysis we have undertaken suggests that similar to 40% of sRNAscanner hits with high cumulative sum of scores represent genuine, undiscovered sRNA genes. Collectively, these data strongly support the utility of sRNAscanner and offer a glimpse of its potential to reveal large numbers of sRNA genes that have to date defied identification. sRNAscanner is available from: http://bicmku.in:8081/sRNAscanner or http://cluster.physics.iisc.ernet.in/sRNAscanner/.

Computational Methods for Studies of Quantum Dots and Boson Condensates

This thesis presents ab initio studies of two kinds of physical systems, quantum dots and bosons, using two program packages of which the bosonic one has mainly been developed by the author. The implemented models, \emph{i.e.}, configuration interaction (CI) and coupled cluster (CC) take the correlated motion of the particles into account, and provide a hierarchy of computational schemes, on top of which the exact solution, within the limit of the single-particle basis set, is obtained. The theory underlying the models is presented in some detail, in order to provide insight into the approximations made and the circumstances under which they hold. Some of the computational methods are also highlighted. In the final sections the results are summarized. The CI and CC calculations on multiexciton complexes in self-assembled semiconductor quantum dots are presented and compared, along with radiative and non-radiative transition rates. Full CI calculations on quantum rings and double quantum rings are also presented. In the latter case, experimental and theoretical results from the literature are re-examined and an alternative explanation for the reported photoluminescence spectra is found. The boson program is first applied on a fictitious model system consisting of bosonic electrons in a central Coulomb field for which CI at the singles and doubles level is found to account for almost all of the correlation energy. Finally, the boson program is employed to study Bose-Einstein condensates confined in different anisotropic trap potentials. The effects of the anisotropy on the relative correlation energy is examined, as well as the effect of varying the interaction potential.}

Computational Fluid Dynamics in Combustion

Computational fluid dynamics has reached a stage where flow field in practical situation can be predicted to aid the design and to probe into the fundamental flow physics to understand and resolve the issues in fundamental fluid mechanics The study examines the computation of reacting flows After exploring the conservation equations for species and energy, the methods of closing the reaction rate terms in turbulent flow have been examined briefly Two cases of computation where combustion-flow interaction plays important role, have been discussed to illustrate the computational aspects and the physical insight that can be gained by the reacting flow computation

Theoretical and Computational Analysis of Static and Dynamic Anomalies in Water-DMSO Binary Mixture at Low DMSO Concentrations

Experiments have repeatedly observed both thermodynamic and dynamic anomalies in aqueous binary mixtures, surprisingly at low solute concentration. Examples of such binary mixtures include water-DMSO, water-ethanol, water-tertiary butyl alcohol (TBA), and water-dioxane, to name a few. The anomalies have often been attributed to the onset of a structural transition, whose nature, however, has been left rather unclear. Here we study the origin of such anomalies using large scale computer simulations and theoretical analysis in water-DMSO binary mixture. At very low DMSO concentration (below 10%), small aggregates of DMSO are solvated by water through the formation of DMSO-(H2O)(2) moieties. As the concentration is increased beyond 10-12% of DMSO, spanning clusters comprising the same moieties appear in the system. Those clusters are formed and stabilized not only through H-bonding but also through the association of CH3 groups of DMSO. We attribute the experimentally observed anomalies to a continuum percolation-like transition at DMSO concentration X-DMSO approximate to 12-15%. The largest cluster size of CH3-CH3 aggregation clearly indicates the formation of such percolating clusters. As a result, a significant slowing down is observed in the decay of associated rotational auto time correlation functions (of the S = O bond vector of DMSO and O-H bond vector of water). Markedly unusual behavior in the mean square fluctuation of total dipole moment again suggests a structural transition around the same concentration range. Furthermore, we map our findings to an interacting lattice model which substantiates the continuum percolation model as the reason for low concentration anomalies in binary mixtures where the solutes involved have both hydrophilic and hydrophobic moieties.

A numerical study of division of flow in open channels

A two-dimensional numerical model which employs the depth-averaged forms of continuity and momentum equations along with k-e turbulence closure scheme is used to simulate the flow at the open channel divisions. The model is generalised to flows of arbitrary geometries and MacCormack finite volume method is used for solving governing equations. Application of cartesian version of the model to analyse the flow at right-angled junction is presented. The numerical predictions are compared with experimental data of earlier investigators and measurements made as part of the present study. Performance of the model in predicting discharge distribution, surface profiles, separation zone parameters and energy losses is evaluated and discussed in detail. To illustrate the application of the numerical model to analyse the flow in acute angled offtakes and streamlined branch entries, a few computational results are presented.

Role of indenter angle on the plastic deformation underneath a sharp indenter and on representative strains: An experimental and numerical study

The subsurface microhardness mapping technique of Chaudhri was utilized to determine the shape, size and distribution of plastic strain underneath conical indenters of varying semi-apex angles, alpha (55 degrees, 65 degrees and 75 degrees). Results show that the elastic-plastic boundary under the indenters is elliptical in nature, contradicting the expanding cavity model, and the ellipticity increases with alpha. The maximum plastic strain immediately under the indenter was found to decrease with increasing alpha. Complementary finite-element analysis was conducted to examine the ability of simulations to capture the experimental observations. A comparison of computational and experimental results indicates that the plastic strain distributions as well as the maximum strains immediately beneath the indenter do not match, suggesting that simulation of sharp indentation requires further detailed studies for complete comprehension. Representative strains, epsilon(r), evaluated as the volume-average strains within the elastic-plastic boundary, decrease with increasing alpha and are in agreement with those estimated by using the dimensional analysis. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Computational and experimental studies of fluid flow and heat transfer in a Calandria based reactor

CFD investigations are carried out to study the heat flux and temperature distribution in the calandria using a 3–Dimensional RANS code. Internal flow computations and experimental studies are carried out for a calandria embedded with a matrix of tubes working together as a reactor. Numerical investigations are carried on the Calandria reactor vessel with horizontal inlets and outlets located on top and the bottom to study the flow pattern and the associated temperature distribution. The computations have been carried out to simulate fluid flow and convective heat transfer for assigned near–to working conditions with different moderator injection rates and reacting heat fluxes. The results of computations provide an estimate of the tolerance bands for safe working limits for the heat dissipation at different working conditions by virtue of prediction of the hot spots in the calandria. The isothermal CFD results are validated by a set of experiments on a specially designed scaled model conducted over a range of flows and simulation parameters. The comparison of CFD results with experiments show good agreement.

Surface effects on stacking fault and twin formation in fcc nanofilms: A first-principles study

The free surface effects on stacking fault and twin formation in fcc metals (Al, Cu, and Ni) were examined by first-principles calculations based on density functional theory (DFT). It is found that the generalized planar fault (GPF) energies of Ni are much larger than bulk Ni with respect to Al and Cu. The discrepancy is attributed to the localized relaxation of Ni nanofilm to accommodate the large expansion of the inter-planar separation induced at the fault plane. The localized relaxation can be coupled to the electronic structure of Ni nanofilms. (C) 2011 Elsevier B.V. All rights reserved.

Study of non-local wave properties of nanotubes with surface effects

In macroscopic and even microscopic structural elements, surface effects can be neglected and classical theories are sufficient. As the structural size decreases towards the nanoscale regime, the surface-to-bulk energy ratio increases and surface effects must be taken into account. In the present work, the terahertz wave dispersion characteristics of a nanotube are studied with consideration of the surface effects as well as the non-local small scale effects. Non-local elasticity theory is used to derive the general governing differential equation based on equilibrium approach to include those scale effects. Scale and surface property dependent wave characteristic equations are obtained via spectral analysis. For the present study the material properties of an anodic alumina nanotube with crystallographic of < 111 > direction are considered. The present analysis shows that the effect of surface properties (surface integrated residual stress and surface integrated modulus) on the flexural wave characteristics of anodic nanotubes are more significant. It has been found that the flexural wavenumbers with surface effects are high as compared to that without surface effects. It has also been shown that, with consideration of surface effects the flexural wavenumbers are under compressive nature. The effect of the small scale and the size of the nanotube on wave dispersion properties are also captured in the present work. (C) 2012 Elsevier B.V. All rights reserved.

Study of Fluid Flow and Mixing Behaviour of a Vacuum Degasser

In recent times the demand of ultra-low carbon steel (ULCS) with improved mechanical properties such as good ductility and good workability has been increased as it is used to produce cold-rolled steel sheets for automobiles. For producing ULCS efficiently, it is necessary to improve the productivity of the vacuum degassers such as RH, DH and tank degasser. Recently, it has been claimed that using a new process, called REDA (revolutionary degassing activator), one can achieve the carbon content below 10 ppm in less time. As such, REDA process has not been studied thoroughly in terms of fluid flow and mass transfer which is a necessary precursor to understand and design this process. Therefore, momentum and mass transfer of the process has been studied by solving momentum and species balance equations along with k-epsilon turbulent model in two-dimension (2D) for REDA process. Similarly, computational fluid dynamic studies have been made in 2D for tank and RH degassers to compare them with REDA process. Computational results have been validated with published experimental and theoretical data. It is found that REDA process is the most efficient among all these processes in terms of mixing efficiency. Fluid flow phenomena have been studied in details for REDA process by varying gas flow rate, depth of immersed snorkel in the steel, diameter of the snorkel and change in vacuum pressure. It is found that design of snorkel affects the melt circulation in the bath significantly.

Mesh Simplification Based on Edge Collapsing Could Improve Computational Efficiency in Near Infrared Optical Tomographic Imaging

The diffusion equation-based modeling of near infrared light propagation in tissue is achieved by using finite-element mesh for imaging real-tissue types, such as breast and brain. The finite-element mesh size (number of nodes) dictates the parameter space in the optical tomographic imaging. Most commonly used finite-element meshing algorithms do not provide the flexibility of distinct nodal spacing in different regions of imaging domain to take the sensitivity of the problem into consideration. This study aims to present a computationally efficient mesh simplification method that can be used as a preprocessing step to iterative image reconstruction, where the finite-element mesh is simplified by using an edge collapsing algorithm to reduce the parameter space at regions where the sensitivity of the problem is relatively low. It is shown, using simulations and experimental phantom data for simple meshes/domains, that a significant reduction in parameter space could be achieved without compromising on the reconstructed image quality. The maximum errors observed by using the simplified meshes were less than 0.27% in the forward problem and 5% for inverse problem.

Bromenium-Catalysed Tandem Ring Opening/Cyclisation of Vinylcyclopropanes and Vinylcyclobutanes: A Metal-Free 3+2+1]/4+2+1] Cascade for the Synthesis of Chiral Amidines and Computational Investigation

We present a detailed study of a 3+2+1] cascade cyclisation of vinylcyclopropanes (VCP) catalysed by a bromenium species (Brd+?Xd-) generated in situ, which results in the synthesis of chiral bicyclic amidines in a tandem one-pot operation. The formation of amidines involves the ring-opening of VCPs with Br?X, followed by a Ritter-type reaction with chloramine-T and a tandem cyclisation. The reaction has been further extended to vinylcyclobutane systems and involves a 4+2+1] cascade cyclisation with the same reagents. The versatility of the methodology has been demonstrated by careful choice of VCPs and VCBs to yield bicyclo4.3.0]-, -4.3.1]- and -4.4.0]amidines in enantiomerically pure form. On the basis of the experimental observations and DFT calculations, a reasonable mechanism has been put forth to account for the formation of the products and the observed stereoselectivity. We propose the existence of a p-stabilised homoallylic carbocation at the cyclopropane carbon as the reason for high stereoselectivity. DFT studies at B3LYP/6-311+G** and M06-2X/6-31+G* levels of theory in gas-phase calculations suggest the ring-opening of VCP is initiated at the p-complex stage (between the double bond and Br?X). This can be clearly perceived from the solution-phase (acetonitrile) calculations using the polarisable continuum model (PCM) solvation model, from which the extent of the ring opening of VCP was found to be noticeably high. Studies also show that the formation of zero-bridge bicyclic amidines is favoured over other bridged bicyclic amidines. The energetics of competing reaction pathways is compared to explain the product selectivity.

Conformational dynamics of HIV-1 protease: a comparative molecular dynamics simulation study with multiple amber force fields

Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the er effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.

Single nucleotide polymorphisms in genes that are common targets of luteotropin and luteolysin in primate corpus luteum: Computational exploration

Luteal insufficiency affects fertility and hence study of mechanisms that regulate corpus luteum (CL) function is of prime importance to overcome infertility problems. Exploration of human genome sequence has helped to study the frequency of single nucleotide polymorphisms (SNPs). Clinical benefits of screening SNPs in infertility are being recognized well in recent times. Examining SNPs in genes associated with maintenance and regression of CL may help to understand unexplained luteal insufficiency and related infertility. Publicly available microarray gene expression databases reveal the global gene expression patterns in primate CL during the different functional state. We intend to explore computationally the deleterious SNPs of human genes reported to be common targets of luteolysin and luteotropin in primate CL Different computational algorithms were used to dissect out the functional significance of SNPs in the luteinizing hormone sensitive genes. The results raise the possibility that screening for SNPs might be integrated to evaluate luteal insufficiency associated with human female infertility for future studies. (C) 2012 Elsevier B.V. All rights reserved,