Progressive structure-based alignment of homologous proteins: Adopting sequence comparison strategies

Comparison of multiple protein structures has a broad range of applications in the analysis of protein structure, function and evolution. Multiple structure alignment tools (MSTAs) are necessary to obtain a simultaneous comparison of a family of related folds. In this study, we have developed a method for multiple structure comparison largely based on sequence alignment techniques. A widely used Structural Alphabet named Protein Blocks (PBs) was used to transform the information on 3D protein backbone conformation as a ID sequence string. A progressive alignment strategy similar to CLUSTALW was adopted for multiple PB sequence alignment (mulPBA). Highly similar stretches identified by the pairwise alignments are given higher weights during the alignment. The residue equivalences from PB based alignments are used to obtain a three dimensional fit of the structures followed by an iterative refinement of the structural superposition. Systematic comparisons using benchmark datasets of MSTAs underlines that the alignment quality is better than MULTIPROT, MUSTANG and the alignments in HOMSTRAD, in more than 85% of the cases. Comparison with other rigid-body and flexible MSTAs also indicate that mulPBA alignments are superior to most of the rigid-body MSTAs and highly comparable to the flexible alignment methods. (C) 2012 Elsevier Masson SAS. All rights reserved.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. (c) 2012 Wiley Periodicals, Inc.

Theory and Algorithms for Hop-Count-Based Localization with Random Geometric Graph Models of Dense Sensor Networks

Wireless sensor networks can often be viewed in terms of a uniform deployment of a large number of nodes in a region of Euclidean space. Following deployment, the nodes self-organize into a mesh topology with a key aspect being self-localization. Having obtained a mesh topology in a dense, homogeneous deployment, a frequently used approximation is to take the hop distance between nodes to be proportional to the Euclidean distance between them. In this work, we analyze this approximation through two complementary analyses. We assume that the mesh topology is a random geometric graph on the nodes; and that some nodes are designated as anchors with known locations. First, we obtain high probability bounds on the Euclidean distances of all nodes that are h hops away from a fixed anchor node. In the second analysis, we provide a heuristic argument that leads to a direct approximation for the density function of the Euclidean distance between two nodes that are separated by a hop distance h. This approximation is shown, through simulation, to very closely match the true density function. Localization algorithms that draw upon the preceding analyses are then proposed and shown to perform better than some of the well-known algorithms present in the literature. Belief-propagation-based message-passing is then used to further enhance the performance of the proposed localization algorithms. To our knowledge, this is the first usage of message-passing for hop-count-based self-localization.

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background: Bacteria such as Escherichia coli and Salmonella typhimurium can utilize acetate as the sole source of carbon and energy. Acetate kinase (AckA) and phosphotransacetylase (Pta), key enzymes of acetate utilization pathway, regulate flux of metabolites in glycolysis, gluconeogenesis, TCA cycle, glyoxylate bypass and fatty acid metabolism. Results: Here we report kinetic characterization of S. typhimurium AckA (StAckA) and structures of its unliganded (Form-I, 2.70 angstrom resolution) and citrate-bound (Form-II, 1.90 angstrom resolution) forms. The enzyme showed broad substrate specificity with k(cat)/K-m in the order of acetate > propionate > formate. Further, the K-m for acetyl-phosphate was significantly lower than for acetate and the enzyme could catalyze the reverse reaction (i.e. ATP synthesis) more efficiently. ATP and Mg2+ could be substituted by other nucleoside 5'-triphosphates (GTP, UTP and CTP) and divalent cations (Mn2+ and Co2+), respectively. Form-I StAckA represents the first structural report of an unliganded AckA. StAckA protomer consists of two domains with characteristic beta beta beta alpha beta alpha beta alpha topology of ASKHA superfamily of proteins. These domains adopt an intermediate conformation compared to that of open and closed forms of ligand-bound Methanosarcina thermophila AckA (MtAckA). Spectroscopic and structural analyses of StAckA further suggested occurrence of inter-domain motion upon ligand-binding. Unexpectedly, Form-II StAckA structure showed a drastic change in the conformation of residues 230-300 compared to that of Form-I. Further investigation revealed electron density corresponding to a citrate molecule in a pocket located at the dimeric interface of Form-II StAckA. Interestingly, a similar dimeric interface pocket lined with largely conserved residues could be identified in Form-I StAckA as well as in other enzymes homologous to AckA suggesting that ligand binding at this pocket may influence the function of these enzymes. Conclusions: The biochemical and structural characterization of StAckA reported here provides insights into the biochemical specificity, overall fold, thermal stability, molecular basis of ligand binding and inter-domain motion in AckA family of enzymes. Dramatic conformational differences observed between unliganded and citrate-bound forms of StAckA led to identification of a putative ligand-binding pocket at the dimeric interface of StAckA with implications for enzymatic function.

Multilevel dodecagonal space-vector generation for induction motor drives by cascading three-level and two-level inverters

This study proposes an inverter circuit topology capable of generating multilevel dodecagonal (12-sided polygon) voltage space vectors by the cascaded connection of two-level and three-level inverters. By the proper selection of DC-link voltages and resultant switching states for the inverters, voltage space vectors whose tips lie on three concentric dodecagons, are obtained. A rectifier circuit for the inverter is also proposed, which significantly improves the power factor. The topology offers advantages such as the complete elimination of the fifth and seventh harmonics in phase voltages and an extension of the linear modulation range. In this study, a simple method for the calculation of pulse width modulation timing was presented along with extensive simulation and experimental results in order to validate the proposed concept.

Fluctuations and correlations of conserved charges in the (2+1) Polyakov quark meson model

We consider the (2 + 1) flavor Polyakov quark meson model and study the fluctuations (correlations) of conserved charges up to sixth (fourth) order. A comparison is made with lattice data wherever available and overall good qualitative agreement is found, more so for the case of the normalized susceptibilities. The model predictions for the ratio of susceptibilities go to that of an ideal gas of hadrons as in hadron resonance gas model at low temperatures while at high temperature the values are close to that of an ideal gas of massless quarks. Our study provides a strong basis for the use of the Polyakov quark meson model as an effective model to understand the topology of the QCD phase diagram. DOI: 10.1103/PhysRevD.86.114021 PACS numbers: 12.39.-x, 05.40.-a, 12.38.Aw, 12.38.Mh

Computing Reeb Graphs as a Union of Contour Trees

The Reeb graph of a scalar function tracks the evolution of the topology of its level sets. This paper describes a fast algorithm to compute the Reeb graph of a piecewise-linear (PL) function defined over manifolds and non-manifolds. The key idea in the proposed approach is to maximally leverage the efficient contour tree algorithm to compute the Reeb graph. The algorithm proceeds by dividing the input into a set of subvolumes that have loop-free Reeb graphs using the join tree of the scalar function and computes the Reeb graph by combining the contour trees of all the subvolumes. Since the key ingredient of this method is a series of union-find operations, the algorithm is fast in practice. Experimental results demonstrate that it outperforms current generic algorithms by a factor of up to two orders of magnitude, and has a performance on par with algorithms that are catered to restricted classes of input. The algorithm also extends to handle large data that do not fit in memory.

The origins of the evolutionary signal used to predict protein-protein interactions

Background: The correlation of genetic distances between pairs of protein sequence alignments has been used to infer protein-protein interactions. It has been suggested that these correlations are based on the signal of co-evolution between interacting proteins. However, although mutations in different proteins associated with maintaining an interaction clearly occur (particularly in binding interfaces and neighbourhoods), many other factors contribute to correlated rates of sequence evolution. Proteins in the same genome are usually linked by shared evolutionary history and so it would be expected that there would be topological similarities in their phylogenetic trees, whether they are interacting or not. For this reason the underlying species tree is often corrected for. Moreover processes such as expression level, are known to effect evolutionary rates. However, it has been argued that the correlated rates of evolution used to predict protein interaction explicitly includes shared evolutionary history; here we test this hypothesis. Results: In order to identify the evolutionary mechanisms giving rise to the correlations between interaction proteins, we use phylogenetic methods to distinguish similarities in tree topologies from similarities in genetic distances. We use a range of datasets of interacting and non-interacting proteins from Saccharomyces cerevisiae. We find that the signal of correlated evolution between interacting proteins is predominantly a result of shared evolutionary rates, rather than similarities in tree topology, independent of evolutionary divergence. Conclusions: Since interacting proteins do not have tree topologies that are more similar than the control group of non-interacting proteins, it is likely that coevolution does not contribute much to, if any, of the observed correlations.

The Coil-to-Helix Transition in IlvN Regulates the Allosteric Control of Escherichia coli Acetohydroxyacid Synthase I

The solution structure of IlvN, the regulatory subunit of Escherichia coil acetohydroxyacid synthase I, in the valine-bound form has been determined using high-resolution multidimensional, multinuclear nuclear magnetic resonance (NMR) methods. IlvN in the presence or absence of the effector molecule is present as a 22.5 kDa dimeric molecule. The ensemble of 20 low-energy structures shows a backbone root-mean-square deviation of 0.73 +/- 0.13 angstrom and a root-mean-square deviation of 1.16 +/- 0.13 angstrom for all heavy atoms. Furthermore, more than 98% of the backbone phi and psi dihedral angles occupy the allowed and additionally allowed regions of the Ramachandran map, which is indicative of the fact that the structures are of high stereochemical quality. Each protomer exhibits a beta alpha beta beta alpha beta alpha topology that is a characteristic feature of the ACT domain seen in metabolic enzymes. In the valine-bound form, IlvN exists apparently as a single conformer. In the free form, IlvN exists as a mixture of conformational states that are in intermediate exchange on the NMR time scale. Thus, a large shift in the conformational equilibrium is observed upon going from the free form to the bound form. The structure of the valine-bound form of IlvN was found to be similar to that of the ACT domain of the unliganded form of IlvH. Comparisons of the structures of the unliganded forms of these proteins suggest significant differences. The structural and conformational properties of IlvN determined here have allowed a better understanding of the mechanism of regulation of branched chain amino acid biosynthesis.

A Hybrid Five-Level Inverter With Common-Mode Voltage Elimination Having Single Voltage Source for IM Drive Applications

A new hybrid five-level inverter topology with common-mode voltage (CMV) elimination for induction motor drive is proposed in this paper. This topology has only one dc source, and different voltage levels are generated by using this voltage source along with floating capacitors charged to asymmetrical voltage levels. The pulsewidth modulation (PWM) scheme employed in this topology balances the capacitor voltages at the required levels at any power factor and modulation index while eliminating the CMV. This inverter has good fault-tolerant capability as it can be operated in three-or two-level mode with CMV elimination, in case of any failure in the H-bridges. More voltage levels with CMV elimination can be realized from this topology but only in a limited range of modulation index and power factor. Extensive simulation is done to validate the PWM technique for CMV elimination and balancing of the capacitor voltages. The experimental verification of the proposed inverter-fed induction motor is carried out in the linear modulation and overmodulation regions. The steady-state and transient operations of the drive are verified. The dynamics of the capacitor voltage balancing is also tested. The experimental results demonstrate that the proposed topology can be considered for industrial drive applications.

Binding of two DNA molecules by type II topoisomerases for decatenation

Topoisomerases (topos) maintain DNA topology and influence DNA transaction processes by catalysing relaxation, supercoiling and decatenation reactions. In the cellular milieu, division of labour between different topos ensures topological homeostasis and control of central processes. In Escherichia coli, DNA gyrase is the principal enzyme that carries out negative supercoiling, while topo IV catalyses decatenation, relaxation and unknotting. DNA gyrase apparently has the daunting task of undertaking both the enzyme functions in mycobacteria, where topo IV is absent. We have shown previously that mycobacterial DNA gyrase is an efficient decatenase. Here, we demonstrate that the strong decatenation property of the enzyme is due to its ability to capture two DNA segments in trans. Topo IV, a strong dedicated decatenase of E. coli, also captures two distinct DNA molecules in a similar manner. In contrast, E. coli DNA gyrase, which is a poor decatenase, does not appear to be able to hold two different DNA molecules in a stable complex. The binding of a second DNA molecule to GyrB/ParE is inhibited by ATP and the non-hydrolysable analogue, AMPPNP, and by the substitution of a prominent positively charged residue in the GyrB N-terminal cavity, suggesting that this binding represents a potential T-segment positioned in the cavity. Thus, after the GyrA/ParC mediated initial DNA capture, GyrB/ParE would bind efficiently to a second DNA in trans to form a T-segment prior to nucleotide binding and closure of the gate during decatenation.

A Hybrid Multilevel Inverter System Based on Dodecagonal Space Vectors for Medium Voltage IM Drives

Dodecagonal (12-sided) space vector pulsewidth modulation (PWM) schemes are characterized by the complete absence of (6n +/- 1)th-order harmonics (for odd n) in the phase voltages, within the linear modulation range and beyond, including over-modulation. This paper presents a new topology suitable for the realization of such multilevel inverter schemes for induction motor (IM) drives, by cascading two-level inverters with flying-capacitor-inverter fed floating H-bridge cells. Now, any standard IM may be used to get the dodecagonal operation which hitherto was possible only with open-end winding IM. To minimize the current total harmonic distortion (THD), a strategy for synchronous PWM is also proposed. It is shown that the proposed method is capable of obtaining better THD figures, compared to conventional dodecagonal schemes. The topology and the PWM strategy are validated through analysis and subsequently verified experimentally.

Medium Voltage Drive for Induction Motors Using Multilevel Octadecagonal Voltage Space Vectors

Multilevel inverters with hexagonal and dodecagonal voltage space vector structures have improved harmonic profile compared to two-level inverters. Further improvement in the quality of the waveform is possible using multilevel octadecagonal (18-sided polygon) voltage space vectors. This paper proposes an inverter circuit topology capable of generating multilevel octadecagonal voltage space vectors, by cascading two asymmetric three-level inverters. By the proper selection of dc-link voltages and the resultant switching states for the inverters, voltage space vectors, whose tips lie on three concentric octadecagons, are obtained. The advantages of octadecagonal voltage space vector-based pulsewidth modulation (PWM) techniques are the complete elimination of fifth, seventh, eleventh, and thirteenth harmonics in phase voltages and the extension of linear modulation range. In this paper, a simple PWM timing calculation method is also proposed. Experimental results have been presented in this paper to validate the proposed concept.

Analogy Between Rotating Euler-Bernoulli and Timoshenko Beams and Stiff Strings

The governing differential equation of a rotating beam becomes the stiff-string equation if we assume uniform tension. We find the tension in the stiff string which yields the same frequency as a rotating cantilever beam with a prescribed rotating speed and identical uniform mass and stiffness. This tension varies for different modes and are found by solving a transcendental equation using bisection method. We also find the location along the rotating beam where equivalent constant tension for the stiff string acts for a given mode. Both Euler-Bernoulli and Timoshenko beams are considered for numerical results. The results provide physical insight into relation between rotating beams and stiff string which are useful for creating basis functions for approximate methods in vibration analysis of rotating beams.

Reuse, recycle to de-bloat software

Most Java programmers would agree that Java is a language that promotes a philosophy of “create and go forth”. By design, temporary objects are meant to be created on the heap, possibly used and then abandoned to be collected by the garbage collector. Excessive generation of temporary objects is termed “object churn” and is a form of software bloat that often leads to performance and memory problems. To mitigate this problem, many compiler optimizations aim at identifying objects that may be allocated on the stack. However, most such optimizations miss large opportunities for memory reuse when dealing with objects inside loops or when dealing with container objects. In this paper, we describe a novel algorithm that detects bloat caused by the creation of temporary container and String objects within a loop. Our analysis determines which objects created within a loop can be reused. Then we describe a source-to-source transformation that efficiently reuses such objects. Empirical evaluation indicates that our solution can reduce upto 40% of temporary object allocations in large programs, resulting in a performance improvement that can be as high as a 20% reduction in the run time, specifically when a program has a high churn rate or when the program is memory intensive and needs to run the GC often.