HARMONY: a server for the assessment of protein structures

Protein structure validation is an important step in computational modeling and structure determination. Stereochemical assessment of protein structures examine internal parameters such as bond lengths and Ramachandran (phi, psi) angles. Gross structure prediction methods such as inverse folding procedure and structure determination especially at low resolution can sometimes give rise to models that are incorrect due to assignment of misfolds or mistracing of electron density maps. Such errors are not reflected as strain in internal parameters. HARMONY is a procedure that examines the compatibility between the sequence and the structure of a protein by assigning scores to individual residues and their amino acid exchange patterns after considering their local environments. Local environments are described by the backbone conformation, solvent accessibility and hydrogen bonding patterns. We are now providing HARMONY through a web server such that users can submit their protein structure files and, if required, the alignment of homologous sequences. Scores are mapped on the structure for subsequent examination that is useful to also recognize regions of possible local errors in protein structures. HARMONY server is located at http://caps.ncbs.res.in/harmony/

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.

Elimination of cross-talk in diffuse optical tomographic reconstructions through a weighted update procedure: results of simulation study - art. no. 61640S

Reconstructions in optical tomography involve obtaining the images of absorption and reduced scattering coefficients. The integrated intensity data has greater sensitivity to absorption coefficient variations than scattering coefficient. However, the sensitivity of intensity data to scattering coefficient is not zero. We considered an object with two inhomogeneities (one in absorption and the other in scattering coefficient). The standard iterative reconstruction techniques produced results, which were plagued by cross talk, i.e., the absorption coefficient reconstruction has a false positive corresponding to the location of scattering inhomogeneity, and vice-versa. We present a method to remove cross talk in the reconstruction, by generating a weight matrix and weighting the update vector during the iteration. The weight matrix is created by the following method: we first perform a simple backprojection of the difference between the experimental and corresponding homogeneous intensity data. The built up image has greater weightage towards absorption inhomogeneity than the scattering inhomogeneity and its appropriate inverse is weighted towards the scattering inhomogeneity. These two weight matrices are used as multiplication factors in the update vectors, normalized backprojected image of difference intensity for absorption inhomogeneity and the inverse of the above for the scattering inhomogeneity, during the image reconstruction procedure. We demonstrate through numerical simulations, that cross-talk is fully eliminated through this modified reconstruction procedure.

A simple classical approach for the melting temperature of inert-gas nanoparticles

Like the metal and semiconductor nanoparticles, the melting temperature of free inert-gas nanoparticles decreases with decreasing size. The variation is linear with the inverse of the particle size for large nanoparticles and deviates from the linearity for small nanoparticles. The decrease in the melting temperature is slower for free nanoparticles with non-wetting surfaces, while the decrease is faster for nanoparticles with wetting surfaces. Though the depression of the melting temperature has been reported for inert-gas nanoparticles in porous glasses, superheating has also been observed when the nanoparticles are embedded in some matrices. By using a simple classical approach, the influence of size, geometry and the matrix on the melting temperature of nanoparticles is understood quantitatively and shown to be applicable for other materials. It is also shown that the classical approach can be applied to understand the size-dependent freezing temperature of nanoparticles.

Distributed nonlinear optimization algorithms for general network problems

There are a number of large networks which occur in many problems dealing with the flow of power, communication signals, water, gas, transportable goods, etc. Both design and planning of these networks involve optimization problems. The first part of this paper introduces the common characteristics of a nonlinear network (the network may be linear, the objective function may be non linear, or both may be nonlinear). The second part develops a mathematical model trying to put together some important constraints based on the abstraction for a general network. The third part deals with solution procedures; it converts the network to a matrix based system of equations, gives the characteristics of the matrix and suggests two solution procedures, one of them being a new one. The fourth part handles spatially distributed networks and evolves a number of decomposition techniques so that we can solve the problem with the help of a distributed computer system. Algorithms for parallel processors and spatially distributed systems have been described.There are a number of common features that pertain to networks. A network consists of a set of nodes and arcs. In addition at every node, there is a possibility of an input (like power, water, message, goods etc) or an output or none. Normally, the network equations describe the flows amoungst nodes through the arcs. These network equations couple variables associated with nodes. Invariably, variables pertaining to arcs are constants; the result required will be flows through the arcs. To solve the normal base problem, we are given input flows at nodes, output flows at nodes and certain physical constraints on other variables at nodes and we should find out the flows through the network (variables at nodes will be referred to as across variables).The optimization problem involves in selecting inputs at nodes so as to optimise an objective function; the objective may be a cost function based on the inputs to be minimised or a loss function or an efficiency function. The above mathematical model can be solved using Lagrange Multiplier technique since the equalities are strong compared to inequalities. The Lagrange multiplier technique divides the solution procedure into two stages per iteration. Stage one calculates the problem variables % and stage two the multipliers lambda. It is shown that the Jacobian matrix used in stage one (for solving a nonlinear system of necessary conditions) occurs in the stage two also.A second solution procedure has also been imbedded into the first one. This is called total residue approach. It changes the equality constraints so that we can get faster convergence of the iterations.Both solution procedures are found to coverge in 3 to 7 iterations for a sample network.The availability of distributed computer systems — both LAN and WAN — suggest the need for algorithms to solve the optimization problems. Two types of algorithms have been proposed — one based on the physics of the network and the other on the property of the Jacobian matrix. Three algorithms have been deviced, one of them for the local area case. These algorithms are called as regional distributed algorithm, hierarchical regional distributed algorithm (both using the physics properties of the network), and locally distributed algorithm (a multiprocessor based approach with a local area network configuration). The approach used was to define an algorithm that is faster and uses minimum communications. These algorithms are found to converge at the same rate as the non distributed (unitary) case.

Optimum atomic spacing for AlAs etching in GaAs epitaxial lift-off technology

With respect to GaAs epitaxial lift-off technology, we report here the optimum atomic spacing (5-10 nm) needed to etch off the AlAs release layer that is sandwiched between two GaAs epitaxial layers. The AlAs etching rate in hydrofluoric acid based solutions was monitored as a function of release layer thickness. We found a sudden quenching in the etching rate, approximately 20 times that of the peak value, at lower dimensions (similar to2.5 nm) of the AlAs epitaxial layer. Since this cannot be explained on the basis of a previous theory (inverse square root of release layer thickness), we propose a diffusion-limited mechanism to explain this reaction process. With the diffusion constant being a mean-free-path-dependent parameter, a relation between the mean free path and the width of the channel is considered. This relation is in reasonable agreement with the experimental results and gives a good physical insight to the reaction kinetics.

Influence of crossed electric and quantizing magnetic fields on the Einstein relation in nonlinear optical, optoelectronic and related materials: Simplified theory, relative comparison and suggestion for experimental determination

An attempt is made to study the Einstein relation for the diffusivity-to-mobility ratio (DMR) under crossed fields' configuration in nonlinear optical materials on the basis of a newly formulated electron dispersion law by incorporating the crystal field in the Hamiltonian and including the anisotropies of the effective electron mass and the spin-orbit splitting constants within the framework of kp formalisms. The corresponding results for III-V, ternary and quaternary compounds form a special case of our generalized analysis. The DMR has also been investigated for II-VI and stressed materials on the basis of various appropriate dispersion relations. We have considered n-CdGeAs2, n-Hg1-xCdxTe, n-In1-xGaxAsyP1-y lattice matched to InP, p-CdS and stressed n-InSb materials as examples. The DMR also increases with increasing electric field and the natures of oscillations are totally band structure dependent with different numerical values. It has been observed that the DMR exhibits oscillatory dependences with inverse quantizing magnetic field and carrier degeneracy due to the Subhnikov-de Haas effect. An experimental method of determining the DMR for degenerate materials in the present case has been suggested. (C) 2010 Elsevier B.V. All rights reserved.

Analysis of clearance fit pin joints

The analysis of clearance fit joints falls within the realm of mixed boundary problems with moving boundaries. In this paper, this problem is solved by a simple continuum method of analysis applying an inverse technique; the region of contact is specified and the corresponding causative load is evaluated. Illustrations are given for a rigid clearance fit pin in a large elastic plate with smooth zero-shear interface between pin and plate, under biaxial plate stress at infinity and due to load transfer through pin.

Segal-Bargmann transform and Paley-Wiener theorems on M(2)

We study the Segal-Bargmann transform on M(2). The range of this transform is characterized as a weighted Bergman space. In a similar fashion Poisson integrals are investigated. Using a Gutzmer's type formula we characterize the range as a class of functions extending holomorphically to an appropriate domain in the complexification of M(2). We also prove a Paley-Wiener theorem for the inverse Fourier transform.

A matrix identity and perturbation expressions

An identity expressing formally the diagonal and off-diagonal elements of an inverse of a matrix is deduced employing operator techniques. Several well-known perturbation expressions for the self-energy are deduced as special cases. A new approximation and other applications following from the above formalism are briefly indicated through illustrations from a perturbed harmonic oscillator, chemisorption approximations and Kelly's result in the problem of electron correlation.

Thyroidal control of hepatic release and metabolism of vitamin A

The inverse relationship that exists between thyroxine and the vitamin A level of plasma has been examined in chicken. Thyroxine treatment leads to a decrease in the level of vitamin A carrier proteins, retinol-binding protein and prealbumin-2 in plasma and liver. There is an accumulation of vitamin A in the liver, with a greater proportion of vitamin A alcohol being present compared to that of control birds. In thyroxine treatment there is enhanced plasma turnover of retinol-binding protein and prealbumin-2, while their rates of synthesis are marginally increased. Amino acid supplementation partially counteracts effects of thyroxine treatment. Amino acid supplementation of thyroxine-treated birds does not alter the plasma turnover rates of retinol-binding protein and prealbumin-2 but increases substentially their rates of synthesis. The release of vitamin A into circulation is interfered with in hyperthyroidism due to inadequate availability of retinol-binding protein being caused by enhanced plasma turnover rate not compensated for by synthesis.

Philosophy as a Path to Happiness : Attainment of Happiness in Arabic Peripatetic and Ismaili Philosophy

The aim of the dissertation is to explore the idea of philosophy as a path to happiness in classical Arabic philosophy. The starting point is in comparison of two distinct currents between the 10th and early 11th centuries, Peripatetic philosophy, represented by al-Fārābī and Ibn Sīnā, and Ismaili philosophy represented by al-Kirmānī and the Brethren of Purity. They initially offer two contrasting views about philosophy in that the attitude of the Peripatetics is rationalistic and secular in spirit, whereas for the Ismailis philosophy represents the esoteric truth behind revelation. Still, they converge in their view that the ultimate purpose of philosophy lies in its ability to lead man towards happiness. Moreover, they share a common concept of happiness as a contemplative ideal of human perfection, which refers primarily to an otherworldly state of the soul s ascent to the spiritual world. For both the way to happiness consists of two parts: theory and practice. The practical part manifests itself in the idea of the purification of the rational soul from its bodily attachments in order for it to direct its attention fully to the contemplative life. Hence, there appears an ideal of philosophical life with the goal of relative detachment from the worldly life. The regulations of the religious law in this context appear as the primary means for the soul s purification, but for all but al-Kirmānī they are complemented by auxiliary philosophical practices. The ascent to happiness, however, takes place primarily through the acquisition of theoretical knowledge. The saving knowledge consists primarily of the conception of the hierarchy of physical and metaphysical reality, but all of philosophy forms a curriculum through which the soul gradually ascends towards a spiritual state of being along an order that is inverse to the Neoplatonic emanationist hierarchy of creation. For Ismaili philosophy the ascent takes place from the exoteric religious sciences towards the esoteric philosophical knowledge. For Peripatetic philosophers logic performs the function of an instrument enabling the ascent, mathematics is treated either as propaedeutic to philosophy or as a mediator between physical and metaphysical knowledge, whereas physics and metaphysics provide the core of knowledge necessary for the attainment of happiness.

Irradiation effects in graphene and related materials

Nanomaterials with a hexagonally ordered atomic structure, e.g., graphene, carbon and boron nitride nanotubes, and white graphene (a monolayer of hexagonal boron nitride) possess many impressive properties. For example, the mechanical stiffness and strength of these materials are unprecedented. Also, the extraordinary electronic properties of graphene and carbon nanotubes suggest that these materials may serve as building blocks of next generation electronics. However, the properties of pristine materials are not always what is needed in applications, but careful manipulation of their atomic structure, e.g., via particle irradiation can be used to tailor the properties. On the other hand, inadvertently introduced defects can deteriorate the useful properties of these materials in radiation hostile environments, such as outer space. In this thesis, defect production via energetic particle bombardment in the aforementioned materials is investigated. The effects of ion irradiation on multi-walled carbon and boron nitride nanotubes are studied experimentally by first conducting controlled irradiation treatments of the samples using an ion accelerator and subsequently characterizing the induced changes by transmission electron microscopy and Raman spectroscopy. The usefulness of the characterization methods is critically evaluated and a damage grading scale is proposed, based on transmission electron microscopy images. Theoretical predictions are made on defect production in graphene and white graphene under particle bombardment. A stochastic model based on first-principles molecular dynamics simulations is used together with electron irradiation experiments for understanding the formation of peculiar triangular defect structures in white graphene. An extensive set of classical molecular dynamics simulations is conducted, in order to study defect production under ion irradiation in graphene and white graphene. In the experimental studies the response of carbon and boron nitride multi-walled nanotubes to irradiation with a wide range of ion types, energies and fluences is explored. The stabilities of these structures under ion irradiation are investigated, as well as the issue of how the mechanism of energy transfer affects the irradiation-induced damage. An irradiation fluence of 5.5x10^15 ions/cm^2 with 40 keV Ar+ ions is established to be sufficient to amorphize a multi-walled nanotube. In the case of 350 keV He+ ion irradiation, where most of the energy transfer happens through inelastic collisions between the ion and the target electrons, an irradiation fluence of 1.4x10^17 ions/cm^2 heavily damages carbon nanotubes, whereas a larger irradiation fluence of 1.2x10^18 ions/cm^2 leaves a boron nitride nanotube in much better condition, indicating that carbon nanotubes might be more susceptible to damage via electronic excitations than their boron nitride counterparts. An elevated temperature was discovered to considerably reduce the accumulated damage created by energetic ions in both carbon and boron nitride nanotubes, attributed to enhanced defect mobility and efficient recombination at high temperatures. Additionally, cobalt nanorods encapsulated inside multi-walled carbon nanotubes were observed to transform into spherical nanoparticles after ion irradiation at an elevated temperature, which can be explained by the inverse Ostwald ripening effect. The simulation studies on ion irradiation of the hexagonal monolayers yielded quantitative estimates on types and abundances of defects produced within a large range of irradiation parameters. He, Ne, Ar, Kr, Xe, and Ga ions were considered in the simulations with kinetic energies ranging from 35 eV to 10 MeV, and the role of the angle of incidence of the ions was studied in detail. A stochastic model was developed for utilizing the large amount of data produced by the molecular dynamics simulations. It was discovered that a high degree of selectivity over the types and abundances of defects can be achieved by carefully selecting the irradiation parameters, which can be of great use when precise pattering of graphene or white graphene using focused ion beams is planned.

The Alfvén Wave Equation for Magnetospheric Plasmas

It is shown that besides the continuous spectrum which damps away as inverse power of time, the coupled Alfvén wave equation, which gives coupling between a shear Alfvén wave and a surface wave, can also admit a well behaved harmonic solution in the closed form for a set of initial conditions. This solution, though valid for finite time intervals, points out that the Alfvén surface waves can have a band of frequency (instead of a monochromatic frequency for a nonsheared magnetic field) within which the local field line resonance frequency can lie, and thus can excite magnetic pulsations with latitude-dependent frequency. By considering magnetic fields not only varying in magnitude but also in direction, it is shown that the time interval for the validity of the harmonic solution depend upon the angle between the magnetic field directions on either side of the magnetopause. For small values of the angle the time interval can become appreciably large.

Thermodynamics of Cobalt (II, III) Oxide (Co3O4): Evidence of Phase Transition

The Gibbs' energy change for the reaction, 3CoO (r.s.)+1/2O2(g)→Co3O4(sp), has been measured between 730 and 1250 K using a solid state galvanic cell: Pt, CuO+Cu2O|(CaO)ZrO2|CoO+Co3O4,Pt. The emf of this cell varies nonlinearly with temperature between 1075 and 1150 K, indicating a second or higher order phase transition in Co3O4around 1120 (±20) K, associated with an entropy change of ∼43 Jmol-1K-1. The phase transition is accompanied by an anomalous increase in lattice parameter and electrical conductivity. The cubic spinel structure is retained during the transition, which is caused by the change in CO+3 ions from low spin to high spin state. The octahedral site preference energy of CO+3 ion in the high spin state has been evaluated as -24.8 kJ mol-1. This is more positive than the value for CO+2 ion (-32.9 kJ mol-1). The cation distribution therefore changes from normal to inverse side during the phase transition. The transformation is unique, coupling spin unpairing in CO+3 ion with cation rearrangement on the spinel lattice, DTA in pure oxygen revealed a small peak corresponding to the transition, which could be differentiated from the large peak due to decomposition. TGA showed that the stoichiometry of oxide is not significantly altered during the transition. The Gibbs' energy of formation of Co3O4 from CoO and O2 below and above phase transition can be represented by the equations:ΔG0=-205,685+170.79T(±200) J mol-1(730-1080 K) and ΔG0=-157,235+127.53T(±200) J mol-1(1150-1250 K).