Computer simulation of protein—carbohydrate complexes: application to arabinose-binding protein and pea lectin

The CCEM method (Contact Criteria and Energy Minimisation) has been developed and applied to study protein-carbohydrate interactions. The method uses available X-ray data even on the native protein at low resolution (above 2.4 Å) to generate realistic models of a variety of proteins with various ligands.The two examples discussed in this paper are arabinose-binding protein (ABP) and pea lectin. The X-ray crystal structure data reported on ABP-β-l-arabinose complex at 2.8, 2.4 and 1.7 Å resolution differ drastically in predicting the nature of the interactions between the protein and ligand. It is shown that, using the data at 2.4 Å resolution, the CCEM method generates complexes which are as good as the higher (1.7 Å) resolution data. The CCEM method predicts some of the important hydrogen bonds between the ligand and the protein which are missing in the interpretation of the X-ray data at 2.4 Å resolution. The theoretically predicted hydrogen bonds are in good agreement with those reported at 1.7 Å resolution. Pea lectin has been solved only in the native form at 3 Å resolution. Application of the CCEM method also enables us to generate complexes of pea lectin with methyl-α-d-glucopyranoside and methyl-2,3-dimethyl-α-d-glucopyranoside which explain well the available experimental data in solution.

An improved representation of vehicle incompatibility in frontal NCAP tests using a modified rigid barrier

With the objective of better understanding the significance of New Car Assessment Program (NCAP) tests conducted by the National Highway Traffic Safety Administration (NHTSA), head-on collisions between two identical cars of different sizes and between cars and a pickup truck are studied in the present paper using LS-DYNA models. Available finite element models of a compact car (Dodge Neon), midsize car (Dodge Intrepid), and pickup truck (Chevrolet C1500) are first improved and validated by comparing theanalysis-based vehicle deceleration pulses against corresponding NCAP crash test histories reported by NHTSA. In confirmation of prevalent perception, simulation-bascd results indicate that an NCAP test against a rigid barrier is a good representation of a collision between two similar cars approaching each other at a speed of 56.3 kmph (35 mph) both in terms of peak deceleration and intrusions. However, analyses carried out for collisions between two incompatible vehicles, such as an Intrepid or Neon against a C1500, point to the inability of the NCAP tests in representing the substantially higher intrusions in the front upper regions experienced by the cars, although peak decelerations in cars arc comparable to those observed in NCAP tests. In an attempt to improve the capability of a front NCAP test to better represent real-world crashes between incompatible vehicles, i.e., ones with contrasting ride height and lower body stiffness, two modified rigid barriers are studied. One of these barriers, which is of stepped geometry with a curved front face, leads to significantly improved correlation of intrusions in the upper regions of cars with respect to those yielded in the simulation of collisions between incompatible vehicles, together with the yielding of similar vehicle peak decelerations obtained in NCAP tests.

Spectroscopic, microscopic and first rheological investigations in charge-transfer interaction induced organogels

This article describes two-component charge-transfer interaction mediated organogels (CT-gels) derived from anthracene carboxamides obtained from 2-amino 2-hydroxymethyl-1,3-propanediol (TRIS), and 2,3-dialkoxyanthracenes as donors, with 2,4,7-trinitrofluorenone (TNF) as the common acceptor. We demonstrate the versatility of TNF as an electron acceptor in the formation of these gels. The effect of subtle changes in the donor structure on the gelation ability has been investigated by varying the alkyl chain length in the dialkoxyanthracene donors, and by varying the position of the TRIS substituent in the anthracene carboxamide donors. Distinct differences have been observed in the nature of the CT-gels based on these two kinds of anthracene donors. It has been reported in the literature that 2,3-dialkoxyanthracenes form gels on their own in various aliphatic hydrocarbons and alcohols for linear alkyl chains bearing at least 6mcarbon atoms (C-6). In the present study, it is shown that themCT-complex of these molecules with TNF is able to gel many alcoholic and a few hydrocarbon solvents. Also, in the presence of TNF, the 2,3-dialkoxyanthracenes (C-4-C-5) which were non-gelators on their own at ambient temperatures, form CT-gels in a number of alcohols. The other series of gelators discussed, the anthracene carboxamides, require the mandatory presence of TNF to form gels. This donor-acceptor complex forms gels in various aliphatic alcohols. Interestingly, the formation of these CT-gels requires rapid cooling in most of the cases. Thermal stability studies with both types of CT-gels indicate an optimum stoichiometry of 1 : 1 between the donor and the acceptor. Dynamic rheological experiments reveal these gels as viscoelastic soft materials, with the mechanical strength of these gels depending on the amount of TNF present. This provides a means to tune the strength of the gel by varying the doping concentration of the acceptor.

Sliding Modes for the Simulation of Mechanical and Electrical Systems Defined by Differential-Algebraic Equations

We offer a technique, motivated by feedback control and specifically sliding mode control, for the simulation of differential-algebraic equations (DAEs) that describe common engineering systems such as constrained multibody mechanical structures and electric networks. Our algorithm exploits the basic results from sliding mode control theory to establish a simulation environment that then requires only the most primitive of numerical solvers. We circumvent the most important requisite for the conventionalsimulation of DAEs: the calculation of a set of consistent initial conditions. Our algorithm, which relies on the enforcement and occurrence of sliding mode, will ensure that the algebraic equation is satisfied by the dynamic system even for inconsistent initial conditions and for all time thereafter. [DOI:10.1115/1.4001904]

Electron microscopic investigation of high-temperature oxide superconductors

Electron microscopic investigations have been carried out on superconducting YBa2Cu3 O7−δ, NdBa2Cu3 O7−δ and related oxides. All these orthorhombic oxides exhibit twin domains. Based on high resolution electron microscopy, it is shown that there is no significant change in the structure across the twins. Oxides of the La2−x Sr x (Ba x )CuO4 system do not show twins, but exhibit other types of defects. Twins appear to be characteristic of only the orthorhombic 123 structures.

Monte carlo simulation for molecular gas dynamics

The dynamics of low-density flows is governed by the Boltzmann equation of the kinetic theory of gases. This is a nonlinear integro-differential equation and, in general, numerical methods must be used to obtain its solution. The present paper, after a brief review of Direct Simulation Monte Carlo (DSMC) methods due to Bird, and Belotserkovskii and Yanitskii, studies the details of theDSMC method of Deshpande for mono as well as multicomponent gases. The present method is a statistical particle-in-cell method and is based upon the Kac-Prigogine master equation which reduces to the Boltzmann equation under the hypothesis of molecular chaos. The proposed Markoff model simulating the collisions uses a Poisson distribution for the number of collisions allowed in cells into which the physical space is divided. The model is then extended to a binary mixture of gases and it is shown that it is necessary to perform the collisions in a certain sequence to obtain unbiased simulation.

Study Of The Hvdc-Torsional Interactions Through Digital Dynamic Simulation

This paper presents the detailed dynamic digital simulation for the study of phenomenon of torsional interaction between HVDC-Turbine generator shaft, dynamics using the novel converter model presented in [ 1 ] The system model includes detailed representation of the synchronous generator and the shaft dynamics, the ac and dc network transients. The results of a case study indicate the various factors that influence the torsional interaction.

Microscopic expression for frequency and wave vector dependent dielectric constant of a dipolar liquid

A microscopic expression for the frequency and wave vector dependent dielectric constant of a dense dipolar liquid is derived starting from the linear response theory. The new expression properly takes into account the effects of the translational modes in the polarization relaxation. The longitudinal and the transverse components of the dielectric constant show vastly different behavior at the intermediate values of the wave vector k. We find that the microscopic structure of the dense liquid plays an important role at intermediate wave vectors. The continuum model description of the dielectric constant, although appropriate at very small values of wave vector, breaks down completely at the intermediate values of k. Numerical results for the longitudinal and the transverse dielectric constants are obtained by using the direct correlation function from the mean‐spherical approximation for dipolar hard spheres. We show that our results are consistent with all the limiting expressions known for the dielectric function of matter.

Numerical simulation of solidification of liquid aluminum alloy flowing on cooling slope

Preparation of semisolid slurry using a cooling slope is increasingly becoming popular, primarily because of the simplicity in design and ease control of the process. In this process, liquid alloy is poured down an inclined surface which is cooled from underneath. The cooling enables partial solidification and the incline provides the necessary shear for producing semisolid slurry. However, the final microstructure of the ingot depends on several process parameters such as cooling rate, incline angle of the cooling slope, length of the slope and initial melt superheat. In this work, a CFD model using volume of fluid (VOF) method for simulating flow along the cooling slope was presented. Equations for conservation of mass, momentum, energy and species were solved to predict hydrodynamic and thermal behavior, in addition to predicting solid fraction distribution and macrosegregation. Solidification was modeled using an enthalpy approach and a volume averaged technique for the different phases. The mushy region was modeled as a multi-layered porous medium consisting of fixed columnar dendrites and mobile equiaxed/fragmented grains. The alloy chosen for the study was aluminum alloy A356, for which adequate experimental data were available in the literature. The effects of two key process parameters, namely the slope angle and the pouring temperature, on temperature distribution, velocity distribution and macrosegregation were also studied.

Complexation of siRNA with Dendrimer: A Molecular Modeling Approach

This paper reports the structural behavior and thermodynamics of the complexation of siRNA with poly(amidoamine) (PAMAM) dendrimers of generation 3 (G3) and 4 (G4) through fully atomistic molecular dynamics (MD) simulations accompanied by free energy calculations and inherent structure determination. We have also done simulation with one siRNA and two dendrimers (2 x G3 or 2xG4) to get the microscopic picture of various binding modes. Our simulation results reveal the formation of stable siRNA-dendrimer complex over nanosecond time scale. With the increase in dendrimcr generation, the charge ratio increases and hence the binding energy between siRNA and dendrimer also increases in accordance with available experimental measurements. Calculated radial distribution functions of amines groups of various subgenerations in a given generation of dendrimer and phosphate in backbone of siRNA reveals that one dendrimer of generation 4 shows better binding with siRNA almost wrapping the dendrimer when compared to the binding with lower generation dendrimer like G3. In contrast, two dendrimers of generation 4 show binding without siRNA wrapping the den-rimer because of repulsion between two dendrimers. The counterion distribution around the complex and the water molecules in the hydration shell of siRNA give microscopic picture of the binding dynamics. We see a clear correlation between water. counterions motions and the complexation i.e. the water molecules and counterions which condensed around siRNA are moved away from the siRNA backbone when dendrimer start binding to the siRNA back hone. As siRNA wraps/bind to the dendrimer counterions originally condensed onto siRNA (Na-1) and dendrimer (Cl-) get released. We give a quantitative estimate of the entropy of counterions and show that there is gain in entropy due to counterions release during the complexation. Furthermore, the free energy of complexation of IG3 and IG4 at two different salt concentrations shows that increase in salt concentration leads to the weakening of the binding affinity of siRNA and dendrimer.

Simulation and graph mining tools for improving gene mapping efficiency

Gene mapping is a systematic search for genes that affect observable characteristics of an organism. In this thesis we offer computational tools to improve the efficiency of (disease) gene-mapping efforts. In the first part of the thesis we propose an efficient simulation procedure for generating realistic genetical data from isolated populations. Simulated data is useful for evaluating hypothesised gene-mapping study designs and computational analysis tools. As an example of such evaluation, we demonstrate how a population-based study design can be a powerful alternative to traditional family-based designs in association-based gene-mapping projects. In the second part of the thesis we consider a prioritisation of a (typically large) set of putative disease-associated genes acquired from an initial gene-mapping analysis. Prioritisation is necessary to be able to focus on the most promising candidates. We show how to harness the current biomedical knowledge for the prioritisation task by integrating various publicly available biological databases into a weighted biological graph. We then demonstrate how to find and evaluate connections between entities, such as genes and diseases, from this unified schema by graph mining techniques. Finally, in the last part of the thesis, we define the concept of reliable subgraph and the corresponding subgraph extraction problem. Reliable subgraphs concisely describe strong and independent connections between two given vertices in a random graph, and hence they are especially useful for visualising such connections. We propose novel algorithms for extracting reliable subgraphs from large random graphs. The efficiency and scalability of the proposed graph mining methods are backed by extensive experiments on real data. While our application focus is in genetics, the concepts and algorithms can be applied to other domains as well. We demonstrate this generality by considering coauthor graphs in addition to biological graphs in the experiments.

Systems analysis for planning of air fleets and maintenance facilities

The high cost and extraordinary demands made on sophisticated air defence systems, pose hard challenges to the managers and engineers who plan the operation and maintenance of such systems. This paper presents a study aimed at developing simulation and systems analysis techniques for the effective planning and efficient operation of small fleets of aircraft, typical of the air force of a developing country. We consider an important aspect of fleet management: the problem of resource allocation for achieving prescribed operational effectiveness of the fleet. At this stage, we consider a single flying-base, where the operationally ready aircraft are stationed, and a repair-depot, where the planes are overhauled. An important measure of operational effectiveness is ‘ availability ’, which may be defined as the expected fraction of the fleet fit for use at a given instant. The tour of aircraft in a flying-base, repair-depot system through a cycle of ‘ operationally ready ’ and ‘ scheduled overhaul ’ phases is represented first by a deterministic flow process and then by a cyclic queuing process. Initially the steady-state availability at the flying-base is computed under the assumptions of Poisson arrivals, exponential service times and an equivalent singleserver repair-depot. This analysis also brings out the effect of fleet size on availability. It defines a ‘ small ’ fleet essentially in terms of the important ‘ traffic ’ parameter of service rate/maximum arrival rate.A simulation model of the system has been developed using GPSS to study sensitivity to distributional assumptions, to validate the principal assumptions of the analytical model such as the single-server assumption and to obtain confidence intervals for the statistical parameters of interest.

Towards Understanding Powder Behavior Via Simulation

The aim of this study was to investigate powder and tablet behavior at the level of mechanical interactions between single particles. Various aspects of powder packing, mixing, compression, and bond formation were examined with the aid of computer simulations. The packing and mixing simulations were based on spring forces interacting between particles. Packing and breakage simulations included systems in which permanent bonds were formed and broken between particles, based on their interaction strengths. During the process, a new simulation environment based on Newtonian mechanics and elementary interactions between the particles was created, and a new method for evaluating mixing was developed. Powder behavior is a complicated process, and many of its aspects are still unclear. Powders as a whole exhibit some aspects of solids and others of liquids. Therefore, their physics is far from clear. However, using relatively simple models based on particle-particle interaction, many powder properties could be replicated during this work. Simulated packing densities were similar to values reported in the literature. The method developed for describing powder mixing correlated well with previous methods. The new method can be applied to determine mixing in completely homogeneous materials, without dividing them into different components. As such, it can describe the efficiency of the mixing method, regardless of the powder's initial setup. The mixing efficiency at different vibrations was examined, and we found that certain combinations of amplitude, direction, and frequencies resulted in better mixing while using less energy. Simulations using exponential force potentials between particles were able to explain the elementary compression behavior of tablets, and create force distributions that were similar to the pressure distributions reported in the literature. Tablet-breaking simulations resulted in breaking strengths that were similar to measured tablet breaking strengths. In general, many aspects of powder behavior can be explained with mechanical interactions at the particle level, and single particle properties can be reliably linked to powder behavior with accurate simulations.