Improved Detection of Remote Homologues Using Cascade PSI-BLAST: Influence of Neighbouring Protein Families on Sequence Coverage

Background: Development of sensitive sequence search procedures for the detection of distant relationships between proteins at superfamily/fold level is still a big challenge. The intermediate sequence search approach is the most frequently employed manner of identifying remote homologues effectively. In this study, examination of serine proteases of prolyl oligopeptidase, rhomboid and subtilisin protein families were carried out using plant serine proteases as queries from two genomes including A. thaliana and O. sativa and 13 other families of unrelated folds to identify the distant homologues which could not be obtained using PSI-BLAST. Methodology/Principal Findings: We have proposed to start with multiple queries of classical serine protease members to identify remote homologues in families, using a rigorous approach like Cascade PSI-BLAST. We found that classical sequence based approaches, like PSI-BLAST, showed very low sequence coverage in identifying plant serine proteases. The algorithm was applied on enriched sequence database of homologous domains and we obtained overall average coverage of 88% at family, 77% at superfamily or fold level along with specificity of similar to 100% and Mathew's correlation coefficient of 0.91. Similar approach was also implemented on 13 other protein families representing every structural class in SCOP database. Further investigation with statistical tests, like jackknifing, helped us to better understand the influence of neighbouring protein families. Conclusions/Significance: Our study suggests that employment of multiple queries of a family for the Cascade PSI-BLAST searches is useful for predicting distant relationships effectively even at superfamily level. We have proposed a generalized strategy to cover all the distant members of a particular family using multiple query sequences. Our findings reveal that prior selection of sequences as query and the presence of neighbouring families can be important for covering the search space effectively in minimal computational time. This study also provides an understanding of the `bridging' role of related families.

Implicit and multigrid procedures for steady-state computations with upwind algorithms

A computational study for the convergence acceleration of Euler and Navier-Stokes computations with upwind schemes has been conducted in a unified framework. It involves the flux-vector splitting algorithms due to Steger-Warming and Van Leer, the flux-difference splitting algorithms due to Roe and Osher and the hybrid algorithms, AUSM (Advection Upstream Splitting Method) and HUS (Hybrid Upwind Splitting). Implicit time integration with line Gauss-Seidel relaxation and multigrid are among the procedures which have been systematically investigated on an individual as well as cumulative basis. The upwind schemes have been tested in various implicit-explicit operator combinations such that the optimal among them can be determined based on extensive computations for two-dimensional flows in subsonic, transonic, supersonic and hypersonic flow regimes. In this study, the performance of these implicit time-integration procedures has been systematically compared with those corresponding to a multigrid accelerated explicit Runge-Kutta method. It has been demonstrated that a multigrid method employed in conjunction with an implicit time-integration scheme yields distinctly superior convergence as compared to those associated with either of the acceleration procedures provided that effective smoothers, which have been identified in this investigation, are prescribed in the implicit operator.

Extrapolation and series acceleration procedures: Hard spheres and discs as illustrations

Three new procedures - in the context of estimation of virial coefficients and summation of the partial virial series for hard discs and hard spheres - are proposed. They are based on the parametrised Euler transformation, a novel resummation, identity and the ε-convergence methods respectively. A comparison with other estimates (molecular dynamics, graph theory and empirical methods) reveals satisfactory agreement.

High microwave susceptibility of NaH2PO4 center dot 2H(2)O: Rapid synthesis of crystalline and glassy phosphates with NASICON-type chemistry

High microwave susceptibility of NaH2PO4 . 2H(2)O has been discovered, This hydrated acid phosphate of sodium can be heated upto 1000 K or more when exposed to 2.45 GHz microwaves. Using this, a novel microwave-assisted preparation of a number of important crystalline and glassy materials with NASICON-type chemistry has been accomplished in less than 8 min which is only a fraction of the time required for conventional synthetic procedures, The present single-shot approach to the preparation of phosphates is attractive in terms of its simplicity, rapidity, and general applicability, A ''step-ladder'' heating mechanism has been proposed to account for the high microwave absorbing ability of NaH2PO4 . 2H(2)O.

In search of inorganic nonlinear optical materials for second harmonic generation

In a search for inorganic oxide materials showing second-order nonlinear optical (NLO) susceptibility, we investigated several berates, silicates, and a phosphate containing trans-connected MO6, octahedral chains or MO5 square pyramids, where, M = d(0): Ti(IV), Nb(V), or Ta(V), Our investigations identified two new NLO structures: batisite, Na2Ba(TiO)(2)Si4O12, containing trans-connected TiO5 octahedral chains, and fresnoite, Ba2TiOSi2O7, containing square-pyramidal TiO5. Investigation of two other materials containing square-pyramidal TiO5 viz,, Cs2TiOP2O7 and Na4Ti2Si8O22. 4H(2)O, revealed that isolated TiO5, square pyramids alone do not cause a second harmonic generation (SHG) response; rather, the orientation of TiO5 units to produce -Ti-O-Ti-O- chains with alternating long and short Ti-O distances in the fresnoite structure is most likely the origin of a strong SHG response in fresnoite,

Software simulation of the EEG

The literature contains many examples of digital procedures for the analytical treatment of electroencephalograms, but there is as yet no standard by which those techniques may be judged or compared. This paper proposes one method of generating an EEG, based on a computer program for Zetterberg's simulation. It is assumed that the statistical properties of an EEG may be represented by stationary processes having rational transfer functions and achieved by a system of software fillers and random number generators.The model represents neither the neurological mechanism response for generating the EEG, nor any particular type of EEG record; transient phenomena such as spikes, sharp waves and alpha bursts also are excluded. The basis of the program is a valid ‘partial’ statistical description of the EEG; that description is then used to produce a digital representation of a signal which if plotted sequentially, might or might not by chance resemble an EEG, that is unimportant. What is important is that the statistical properties of the series remain those of a real EEG; it is in this sense that the output is a simulation of the EEG. There is considerable flexibility in the form of the output, i.e. its alpha, beta and delta content, which may be selected by the user, the same selected parameters always producing the same statistical output. The filtered outputs from the random number sequences may be scaled to provide realistic power distributions in the accepted EEG frequency bands and then summed to create a digital output signal, the ‘stationary EEG’. It is suggested that the simulator might act as a test input to digital analytical techniques for the EEG, a simulator which would enable at least a substantial part of those techniques to be compared and assessed in an objective manner. The equations necessary to implement the model are given. The program has been run on a DEC1090 computer but is suitable for any microcomputer having more than 32 kBytes of memory; the execution time required to generate a 25 s simulated EEG is in the region of 15 s.

Low temperature fluorination of sulphur compounds with elemental fluorine

Low temperature fluorination technique is adopted for fluorination of the following sulphur compounds in freon-11 medium (1) Sulphur dioxide (2) Thionyl chloride (3) Sulphuryl chloride (4) Tetrasulphur tetra nitride and (5) Sulphur bromide. All the compounds undergo oxidative fluorination to give rise to sulphur-fluorine compounds except sulphuryl chloride which resists fluorination. Sulphuryl chloride thus behaves as a good solvent medium for fluorination of other reactive compounds like elemental sulphur. Details of the experimental procedures adopted and the identification of the products will be presented.

Use of α- and β-subunit specific antibodies in studying interaction of hCG with Leydig cell receptors

Antisera (a/s) raised to individual α- and β-subunits of human chorionic gonadotropin (hCG) have been characterized for specificity using immunoaffinity procedures and used to study the disposition of the two subunits when intact hCG is complexed with luteinizing hormone (LH) receptor of the Leydig cells. Three kinds of experiments were done. (a) The ability of the preformed hormone-antibody (H-Ab) complex to bind to receptor and stimulate a response; (b) the ability of the a/s to dissociate hCG from its complex with the receptor and thereby terminate response; and (c) the ability of the premixed antibody and receptor to compete for binding of labeled hCG. Although the subunit specific a/s used here were equipotent in binding hCG (capacity to bind and Ka being very similar), their behavior once the receptor preparation or Leydig cell is introduced into the system was drastically different. The β-subunit antibody relative to the α-subunit antibody, appeared to be poorly effective in preventing hCG from either binding to the receptor or inhibiting the continuation of response. The results suggest that hCG upon interaction with the receptor loses the determinants specific to the β-region more rapidly compared to those specific to the α-region suggesting thereby that the initial interaction of hCG with the receptor should be occurring through sites in the β-subunit. Although the α-subunit portion of the hCG molecule is available for binding to the antibody for a relatively longer time, the biological response of the cell seems very sensitive to such binding with the antibody as it invariably results in loss of response. In the Leydig cell system, the ability of the a/s to bind hCG that is already complexed to the receptor appears to be dependent upon the time of addition of the antibody to the incubation medium. The antisera were totally ineffective in inhibiting steroidogenic response to hCG if added 60 min after addition of hCG. This would suggest that the hormone-receptor complex once formed perhaps continues to change its orientation with the result that with time relatively less and less of antigenic determinants become available for antibody binding.

Identification of parameters in unconfined aquifers

A method is presented for identification of parameters in unconfined aquifers from pumping tests, based on the optimisation of the objective function using the least squares approach. Four parameters are to be evaluated, namely: The hydraulic conductivity in the radial and the vertical directions, the storage coefficient and the specific yield. The sensitivity analysis technique is used for solving the optimisation problem. Besides eliminating the subjectivity involved in the graphical procedure, the method takes into account the field data at all time intervals without classifying them into small and large time intervals and does not use the approximation that the ratio of the storage coefficient to the specific yield tends to zero. Two illustrative examples are presented and it is found that the parameter estimates from the computational and graphical procedures differ fairly significantly.

A computerized methodology for structural synthesis of kinematic chains: Part 1- Formulation

In an effort to develop a fully computerized approach for structural synthesis of kinematic chains the steps involved in the method of structural synthesis based on transformation of binary chains [38] have been recast in a format suitable for implementation on a digital computer. The methodology thus evolved has been combined with the algebraic procedures for structural analysis [44] to develop a unified computer program for structural synthesis and analysis of simple jointed kinematic chains with a degree of freedom 0. Applications of this program are presented in the succeeding parts of the paper.

Handpumps: problems and search for remedies

This paper gives a brief survey of research and development work done on hand pumps in India as well as elsewhere and sets out the approach adopted by ASTRA Working Group. Ten ways in which a hand pump breakdown in practice have been identified. The physical reasons behind each type of breakdown analysed. Remedial measures have been developed from this analysis. Laboratory test rigs fabricated to evaluate these measures have been described and some experimental results presented. The course of further work has been charted.

Infrared spectra of 1.3-dithiole-2-thione and its selenium analogues ând frequency assignment and molecular force constants

Infrared spectra of 1,3-dithiole-2-thione (DTT) and its four selenium analogues have been studied in the region 4000 to 20 cm�1. Assignment of all the fundamental frequencies was made by noting the band shifts on progressive selenation. Normal coordinate analysis procedures have been applied for both in-plane and out-of-plane vibrations to help the assignments. The Urey�Bradley force function supplemented with valence force constants for the out-of-plane vibrations was employed for coordinate calculations. A correlation of the infrared assignments of DTT with its different selenium analogues is accomplished. Further, the infrared assignments are compared with those of trithiocarbonate ion and its selenium analogues and other structurally related heterocyclic molecules.

Perturbation expansions and series acceleration procedures: Part-II. Extrapolation techniques

Three new procedures for the extrapolation of series coefficients from a given power series expansion are proposed. They are based on (i) a novel resummation identity, (ii) parametrised Euler transformation (pet) and (iii) a modifiedpet. Several examples taken from the Ising model series expansions, ferrimagnetic systems, etc., are illustrated. Apart from these applications, the higher order virial coefficients for hard spheres and hard discs have also been evaluated using the new techniques and these are compared with the estimates obtained by other methods. A satisfactory agreement is revealed between the two.

Computer-aided analysis of the structure of kinematic chains

Using the link-link incidence matrix to represent a simple-jointed kinematic chain algebraic procedures have been developed to determine its structural characteristics such as the type of freedom of the chain, the number of distinct mechanisms and driving mechanisms that can be derived from the chain. A computer program incorporating these graph theory based procedures has been applied successfully for the structural analysis of several typical chains.

DFL: A Data Flow Language

Many novel computer architectures like array and multiprocessors which achieve high performance through the use of concurrency exploit variations of the von Neumann model of computation. The effective utilization of the machines makes special demands on programmers and their programming languages, such as the structuring of data into vectors or the partitioning of programs into concurrent processes. In comparison, the data flow model of computation demands only that the principle of structured programming be followed. A data flow program, often represented as a data flow graph, is a program that expresses a computation by indicating the data dependencies among operators. A data flow computer is a machine designed to take advantage of concurrency in data flow graphs by executing data independent operations in parallel. In this paper, we discuss the design of a high level language (DFL: Data Flow Language) suitable for data flow computers. Some sample procedures in DFL are presented. The implementation aspects have not been discussed in detail since there are no new problems encountered. The language DFL embodies the concepts of functional programming, but in appearance closely resembles Pascal. The language is a better vehicle than the data flow graph for expressing a parallel algorithm. The compiler has been implemented on a DEC 1090 system in Pascal.