Parallel implementation of AutoDock

Computational docking of ligands to protein structures is a key step in structure-based drug design. Currently, the time required for each docking run is high and thus limits the use of docking in a high-throughput manner, warranting parallelization of docking algorithms. AutoDock, a widely used tool, has been chosen for parallelization. Near-linear increases in speed were observed with 96 processors, reducing the time required for docking ligands to HIV-protease from 81 min, as an example, on a single IBM Power-5 processor ( 1.65 GHz), to about 1 min on an IBM cluster, with 96 such processors. This implementation would make it feasible to perform virtual ligand screening using AutoDock.

Time and Energy Complexity of Distributed Computation of a Class of Functions in Wireless Sensor Networks

We consider a scenario in which a wireless sensor network is formed by randomly deploying n sensors to measure some spatial function over a field, with the objective of computing a function of the measurements and communicating it to an operator station. We restrict ourselves to the class of type-threshold functions (as defined in the work of Giridhar and Kumar, 2005), of which max, min, and indicator functions are important examples: our discussions are couched in terms of the max function. We view the problem as one of message-passing distributed computation over a geometric random graph. The network is assumed to be synchronous, and the sensors synchronously measure values and then collaborate to compute and deliver the function computed with these values to the operator station. Computation algorithms differ in (1) the communication topology assumed and (2) the messages that the nodes need to exchange in order to carry out the computation. The focus of our paper is to establish (in probability) scaling laws for the time and energy complexity of the distributed function computation over random wireless networks, under the assumption of centralized contention-free scheduling of packet transmissions. First, without any constraint on the computation algorithm, we establish scaling laws for the computation time and energy expenditure for one-time maximum computation. We show that for an optimal algorithm, the computation time and energy expenditure scale, respectively, as Theta(radicn/log n) and Theta(n) asymptotically as the number of sensors n rarr infin. Second, we analyze the performance of three specific computation algorithms that may be used in specific practical situations, namely, the tree algorithm, multihop transmission, and the Ripple algorithm (a type of gossip algorithm), and obtain scaling laws for the computation time and energy expenditure as n rarr infin. In particular, we show that the computation time for these algorithms scales as Theta(radicn/lo- g n), Theta(n), and Theta(radicn log n), respectively, whereas the energy expended scales as , Theta(n), Theta(radicn/log n), and Theta(radicn log n), respectively. Finally, simulation results are provided to show that our analysis indeed captures the correct scaling. The simulations also yield estimates of the constant multipliers in the scaling laws. Our analyses throughout assume a centralized optimal scheduler, and hence, our results can be viewed as providing bounds for the performance with practical distributed schedulers.

Reliability assessment of internal stability of reinforced soil structures: A pseudo-dynamic approach

A methodology for reliability based optimum design of reinforced soil structures subjected to horizontal and vertical sinusoidal excitation based on pseudo-dynamic approach is presented. The tensile strength of reinforcement required to maintain the stability is computed using logarithmic spiral failure mechanism. The backfill soil properties, geometric and strength properties of reinforcement are treated as random variables. Effects of parameters like soil friction angle, horizontal and vertical seismic accelerations, shear and primary wave velocities, amplification factors for seismic acceleration on the component and system probability of failures in relation to tension and pullout capacities of reinforcement have been discussed. In order to evaluate the validity of the present formulation, static and seismic reinforcement force coefficients computed by the present method are compared with those given by other authors. The importance of the shear wave velocity in the estimation of the reliability of the structure is highlighted. The Ditlevsen's bounds of system probability of failure are also computed by taking into account the correlations between three failure modes, which is evaluated using the direction cosines of the tangent planes at the most probable points of failure. (c) 2009 Elsevier Ltd. All rights reserved.

Automatic Road Extraction using High Resolution Satellite Image Based on Texture Progressive Analysis and Normalized Cut Method

In this paper the approach for automatic road extraction for an urban region using structural, spectral and geometric characteristics of roads has been presented. Roads have been extracted based on two levels: Pre-processing and road extraction methods. Initially, the image is pre-processed to improve the tolerance by reducing the clutter (that mostly represents the buildings, parking lots, vegetation regions and other open spaces). The road segments are then extracted using Texture Progressive Analysis (TPA) and Normalized cut algorithm. The TPA technique uses binary segmentation based on three levels of texture statistical evaluation to extract road segments where as, Normalizedcut method for road extraction is a graph based method that generates optimal partition of road segments. The performance evaluation (quality measures) for road extraction using TPA and normalized cut method is compared. Thus the experimental result show that normalized cut method is efficient in extracting road segments in urban region from high resolution satellite image.

New measures for estimating surface complementarity and packing at protein-protein interfaces

A number of methods exist that use different approaches to assess geometric properties like the surface complementarity and atom packing at the protein-protein interface. We have developed two new and conceptually different measures using the Delaunay tessellation and interface slice selection to compute the surface complementarity and atom packing at the protein-protein interface in a straightforward manner. Our measures show a strong correlation among themselves and with other existing measures, and can be calculated in a highly time-efficient manner. The measures are discriminative for evaluating biological, as well as non-biological protein-protein contacts, especially from large protein complexes and large-scale structural studies(http://pallab.serc. iisc.ernet.in/nip_nsc). (C) 201 Federation of European Biochemical Societies. Published by Elsevier B. V. All rights reserved.

Handling sectional views in volume-based approach to automatically construct 3D solid from 2D views

This paper presents an algorithm for solid model reconstruction from 2D sectional views based on volume-based approach. None of the existing work in automatic reconstruction from 2D orthographic views have addressed sectional views in detail. It is believed that the volume-based approach is better suited to handle different types of sectional views. The volume-based approach constructs the 3D solid by a boolean combination of elementary solids. The elementary solids are formed by sweep operation on loops identified in the input views. The only adjustment to be made for the presence of sectional views is in the identification of loops that would form the elemental solids. In the algorithm, the conventions of engineering drawing for sectional views, are used to identify the loops correctly. The algorithm is simple and intuitive in nature. Results have been obtained for full sections, offset sections and half sections. Future work will address other types of sectional views such as removed and revolved sections and broken-out sections. (C) 2004 Elsevier Ltd. All rights reserved.

Adhesively-bonded Patch Repair with Composites

Adhesively-bonded composite patch repairs over cracked or corrosion-damaged metallic aircraft structures have shown great promise for extending life of ageing structures. This study presents the numerical investigation into the interface behaviour of adhesively-bonded cracked aluminum alloy substrate patched with fibre-reinforced composite material. The adhesive is modelled as an elasto-plastic bilinear material to characterise the debond behaviour, while the defective substrate is regarded as linear elastic continuum. Two typical patch shapes were selected based on information available in the literature. Geometric and material nonlinear analyses for square and octagonal patches were performed to capture peel and shear stresses developed between the substrate and the patch to examine the possibility of interface delamination/debonding. Parametric studies on adhesive thickness and patch thickness were carried out to predict their infuence on damage tolerance of repaired structures.

Long-range synthon Aufbau modules (LSAM) in crystal structures: systematic changes in C6H6-nFn (0 <= n <= 6) fluorobenzenes

We discuss the assembly of a three-dimensional molecular crystal in terms of short-range supramolecular synthons that spontaneously organize themselves according to Aufbau principles into long-range geometries characteristic of the molecules themselves. For this purpose we have examined the systematic changes in the known crystal structures of a family of fluorobenzenes, C6H6-nFn, where 0 <= n <= 6. Crystal assembly is initiated by forming long-range synthon Aufbau modules (LSAM) that carry the imprint of the synthons. For example, when 1 <= n <= 5 the short-range synthons use H center dot center dot center dot F interactions to form the LSAMs. In the n = 0 and n = 6 compounds, the synthons are H center dot center dot center dot C and F center dot center dot center dot C interactions, respectively. The LSAMs are usually one-dimensional. In this study we show that these 1D LSAMs assemble into 2D quasi-hexagonal close-packed layers. The 3D crystal structure is obtained from the various kinds of close-packing known for these 2D layers. The final stages of this 1D -> 2D -> 3D assembly seem to be more influenced by the packing of LSAMs than by any other factor. In these final stages, there may not be so much influence exerted by the stronger short-range synthons. We discuss the evolution of these fluorobenzene crystal structures in terms of putative LSAMs and the purely geometric relationships between the n and (6 - n) compounds that can thus be expected. Such particle-hole pairs show structural similarities. Our discussion is quantified by the interpretation of intermolecular distances in terms of atomic sizes and with qualitative predictions of magnetic model systems.

Long-range synthon Aufbau modules (LSAM) in crystal structures: systematic changes in C6H6-nFn (0 <= n <= 6) fluorobenzenes

We discuss the assembly of a three-dimensional molecular crystal in terms of short-range supramolecular synthons that spontaneously organize themselves according to Aufbau principles into long-range geometries characteristic of the molecules themselves. For this purpose we have examined the systematic changes in the known crystal structures of a family of fluorobenzenes, C6H6-nFn, where 0 <= n <= 6. Crystal assembly is initiated by forming long-range synthon Aufbau modules (LSAM) that carry the imprint of the synthons. For example, when 1 <= n <= 5 the short-range synthons use H center dot center dot center dot F interactions to form the LSAMs. In the n = 0 and n = 6 compounds, the synthons are H center dot center dot center dot C and F center dot center dot center dot C interactions, respectively. The LSAMs are usually one-dimensional. In this study we show that these 1D LSAMs assemble into 2D quasi-hexagonal close-packed layers. The 3D crystal structure is obtained from the various kinds of close-packing known for these 2D layers. The final stages of this 1D -> 2D -> 3D assembly seem to be more influenced by the packing of LSAMs than by any other factor. In these final stages, there may not be so much influence exerted by the stronger short-range synthons. We discuss the evolution of these fluorobenzene crystal structures in terms of putative LSAMs and the purely geometric relationships between the n and (6 - n) compounds that can thus be expected. Such particle-hole pairs show structural similarities. Our discussion is quantified by the interpretation of intermolecular distances in terms of atomic sizes and with qualitative predictions of magnetic model systems.

On Using Projective Relations for Calibration and Estimation in a Structured-Light Scanner

For structured-light scanners, the projective geometry between a projector-camera pair is identical to that of a camera-camera pair. Consequently, in conjunction with calibration, a variety of geometric relations are available for three-dimensional Euclidean reconstruction. In this paper, we use projector-camera epipolar properties and the projective invariance of the cross-ratio to solve for 3D geometry. A key contribution of our approach is the use of homographies induced by reference planes, along with a calibrated camera, resulting in a simple parametric representation for projector and system calibration. Compared to existing solutions that require an elaborate calibration process, our method is simple while ensuring geometric consistency. Our formulation using the invariance of the cross-ratio is also extensible to multiple estimates of 3D geometry that can be analysed in a statistical sense. The performance of our system is demonstrated on some cultural artifacts and geometric surfaces.

On the compactness of isometry groups in complex analysis

We prove that the group of continuous isometries for the Kobayashi or Caratheodory metrics of a strongly convex domain in C-n is compact unless the domain is biholomorphic to the ball. A key ingredient, proved using differential geometric ideas, is that a continuous isometry between a strongly convex domain and the ball has to be biholomorphic or anti-biholomorphic. Combining this with a metric version of Pinchuk's rescaling technique gives the main result.

Inverted V-Notch: Practical Proportional Weir

This paper presents a practical linear proportional weir of simple geometric shape in the form of an inverted V-notch or inward trapezium. The flow through this weir, of half-width w and altitude d, for depths above 0.22d is proportional to the depth of flow measured above a reference plane situated at 0.08d for all heads in the range 0.22d<=h<=0.94d, with a maximum percentage deviation of ±1.5 from the theoretical discharge. The linear relationship between head and discharge is based on numerical optimization procedures. Nearly 75% of the depth of inverted V-notch can be used effectively as the measuring range. Experiments with four weirs, with different vertex angles, show excellent agreement with the theory by giving an average coefficient of discharge for each weir varying from 0.61–0.62.

Behaviour of Bi-adhesive Joints

The use of relatively low modulus adhesive at the ends of overlap in a bi-adhesive bondline of a bonded joint can reduce the stress concentration significantly and, therefore, potentially lead to higher strength of the joint. This study presents the two-dimensional and three-dimensional nonlinear (geometric and material) finite element analyses of adhesively bonded single lap joints having modulus-graded bondline under monotonic loading conditions. The adhesives were modelled as an elasto-plastic multi-linear material, while the substrates were regarded as both linear elastic and bi-linear elasto-plastic material. The computational simulations have been performed to investigate the bondline behaviour by studying the stress and strain distributions both at the mid-plane as well as at the interface of the bondline. It has been observed that the static strength is higher for joints with bi-adhesive bondlines compared to those with single adhesives in bondline. Higher joint strength has also been observed for optimum bi-adhesive bondline ratio through parametric studies. Effects of load level, and bondline thickness on stress distribution in the bi-adhesive bondline have also been studied. 3D analysis results reveal the existence of complex multi-axial stress/strain state at the ends of the overlap in the bondline which cannot be observed in 2D plane strain analysis. About 1/3rd of the width of the joint from the free edge in the width direction has 3D stress state, especially in the compliant adhesive of the bondline. Magnitudes of longitudinal and lateral stress/strain components are comparable to peel stress/strain components. It has also been analytically shown that the in-plane global stiffness of the joint remains unaffected by modulus gradation of the bondline adhesive. (C) Koninklijke Brill NV, Leiden, 2010.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Adhesively-bonded Patch Repair with Composites

Adhesively-bonded composite patch repairs over cracked or corrosion-damaged metallic aircraft structures have shown great promise for extending life of ageing structures. This study presents the numerical investigation into the interface behaviour of adhesively-bonded cracked aluminum alloy substrate patched with fibre-reinforced composite material. The adhesive is modelled as an elasto-plastic bilinear material to characterise the debond behaviour, while the defective substrate is regarded as linear elastic continuum. Two typical patch shapes were selected based on information available in the literature. Geometric and material nonlinear analyses for square and octagonal patches were performed to capture peel and shear stresses developed between the substrate and the patch to examine the possibility of interface delamination/debonding. Parametric studies on adhesive thickness and patch thickness were carried out to predict their infuence on damage tolerance of repaired structures.