Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

Free and ATP-bound structures of Ap(4)A hydrolase from Aquifex aeolicus V5

Asymmetric diadenosine tetraphosphate (Ap(4)A) hydrolases degrade the metabolite Ap(4)A back into ATP and AMP. The three-dimensional crystal structure of Ap(4)A hydrolase (16 kDa) from Aquifex aeolicus has been determined in free and ATP-bound forms at 1.8 and 1.95 angstrom resolution, respectively. The overall three-dimensional crystal structure of the enzyme shows an alpha beta alpha-sandwich architecture with a characteristic loop adjacent to the catalytic site of the protein molecule. The ATP molecule is bound in the primary active site and the adenine moiety of the nucleotide binds in a ring-stacking arrangement equivalent to that observed in the X-ray structure of Ap(4)A hydrolase from Caenorhabditis elegans. Binding of ATP in the active site induces local conformational changes which may have important implications in the mechanism of substrate recognition in this class of enzymes. Furthermore, two invariant water molecules have been identified and their possible structural and/or functional roles are discussed. In addition, modelling of the substrate molecule at the primary active site of the enzyme suggests a possible path for entry and/or exit of the substrate and/or product molecule.

Force-free equilibria of magnetized jets

Force-free equilibrium configurations of magnetic-pressure-dominated magnetized supersonic jets confined by slowly varying external pressure are investigated analytically. For the case where internal dissipation mechanisms are active, the lowest-energy field configuration is found to be the superposition of an axisymmetric mode and a helical mode with a wavelength equal to 5 times the jet radius, and the pressure below which the nonaxisymmetric mode becomes energetically favorable is given as 2700 times the product of the 4th power of the magnetic helicity per unit length and the -6th power of the magnetic flux. A model of the total and polarized emission of such a configuration is developed and applied to the extended well-collimated astronomically resolved jet NGC 6251. The model is shown to reproduce significant features such as transverse oscillations of the ridge line, width oscillations and emission knots, the projected magnetic-field configuration, oscillations of the degree of polarization, and the distribution of the Faraday rotation measure.

Magnetic energy dissipation in force-free jets

It is shown that a magnetic-pressure-dominated, supersonic jet which expands (or contracts) in response to variations in the confining external pressure can dissipate magnetic energy through field-line reconnection as it relaxes to a minimum-energy configuration. In order for a continuous dissipation to take place, the effective reconnection time must be a fraction ɛ ⪉ 1 of the expansion time. The amount of energy dissipation is calculated, and it is concluded that magnetic energy dissipation could, in principle, power the observed synchrotron emission in extragalactic radio jets such as NGC 6251. However, this mechanism is only viable if the reconnection time is substantially shorter than the nominal resistive tearing time in the jet.

Mechanism of failure in a free-standing Pt-aluminide bond coat during tensile testing at room temperature

The room temperature (RT) tensile behaviour of a free-standing high activity Pt-aluminide bond coat has been evaluated by microtensile testing technique. The coating had a typical three-layer microstructure. The stress-strain plot for the free-standing coating was linear, indicating the coating to be brittle at RT. Different fracture features were observed across the coating layers, namely quasi-cleavage in the outer layer and inner interdiffusion zone, and cleavage in the intermediate layer. By employing interrupted tensile test and observing the cross-sectional microstructure of the tested specimens, it was determined that failure of the microtensile samples occurred by the initiation of a single crack in the intermediate layer of the coating and its subsequent inside-out propagation. Such a mechanism of failure has been explained in terms of the fracture features observed across the sample thickness. This mechanism of failure is consistent with fracture toughness values of the individual coating layers. (C) 2009 Elsevier B.V. All rights reserved.

Computation of Thermal Stress Intensity Factors for Bimaterial Interface Cracks Using Domain Integral Method

The concept of domain integral used extensively for J integral has been applied in this work for the formulation of J(2) integral for linear elastic bimaterial body containing a crack at the interface and subjected to thermal loading. It is shown that, in the presence of thermal stresses, the J(k) domain integral over a closed path, which does not enclose singularities, is a function of temperature and body force. A method is proposed to compute the stress intensity factors for bimaterial interface crack subjected to thermal loading by combining this domain integral with the J(k) integral. The proposed method is validated by solving standard problems with known solutions.

Behaviour of Bi-material interface cracks in the presence of material nonlinearity of adherends

The well known features of crack face interpenetration/contact at the tip of an interface crack is re-examined using finite element analysis and assuming material nonlinear properties for the adherends. It was assumed in literature that the crack tips are fully open at all load levels in the presence of material nonlinearity of the adherends. Analysis for the case of remote tension shows that even in the presence of material nonlinearity, crack tip closes at small load levels and opens above a certain load level. Mixed-mode fracture parameters are evaluated for the situation when the crack tips are fully open. Due to the presence of nonlinearity, the mixed-mode fracture parameters are measured with the symmetric and anti-symmetric components of J-integral. The present analysis explains the sequence of events at the interface crack tip with progressively increasing remote tension load for the case of adherends with material nonlinear behaviour.

Free-volume dependent pressure sensitivity of Zr-based bulk metallic glass

Instrumented indentation experiments on a Zr-based bulk metallic glass (BMG) in as-cast, shot-peened and structurally relaxed conditions were conducted to examine the dependence of plastic deformation on its structural state. Results show significant differences in hardness, H, with structural relaxation increasing it and shot peening markedly reducing it, and slightly changed morphology of shear bands around the indents. This is in contrast to uniaxial compressive yield strength, sigma(y), which remains invariant with the change in the structural state of the alloys investigated. The plastic constraint factor, C = H/sigma(y), of the relaxed BMG increases compared with that of the as-cast glass, indicating enhanced pressure sensitivity upon annealing. In contrast, C of the shot-peened layer was found to be similar to that observed in crystalline metals, indicating that severe plastic deformation could eliminate pressure sensitivity. Microscopic origins for this result, in terms of shear transformation zones and free volume, are discussed.

On the Cubicity of AT-Free Graphs and Circular-Arc Graphs

A unit cube in k dimensions (k-cube) is defined as the Cartesian product R-1 x R-2 x ... x R-k where R-i (for 1 <= i <= k) is a closed interval of the form [a(i), a(i) + 1] on the real line. A graph G on n nodes is said to be representable as the intersection of k-cubes (cube representation in k dimensions) if each vertex of C can be mapped to a k-cube such that two vertices are adjacent in G if and only if their corresponding k-cubes have a non-empty intersection. The cubicity of G denoted as cub(G) is the minimum k for which G can be represented as the intersection of k-cubes. An interesting aspect about cubicity is that many problems known to be NP-complete for general graphs have polynomial time deterministic algorithms or have good approximation ratios in graphs of low cubicity. In most of these algorithms, computing a low dimensional cube representation of the given graph is usually the first step. We give an O(bw . n) algorithm to compute the cube representation of a general graph G in bw + 1 dimensions given a bandwidth ordering of the vertices of G, where bw is the bandwidth of G. As a consequence, we get O(Delta) upper bounds on the cubicity of many well-known graph classes such as AT-free graphs, circular-arc graphs and cocomparability graphs which have O(Delta) bandwidth. Thus we have: 1. cub(G) <= 3 Delta - 1, if G is an AT-free graph. 2. cub(G) <= 2 Delta + 1, if G is a circular-arc graph. 3. cub(G) <= 2 Delta, if G is a cocomparability graph. Also for these graph classes, there axe constant factor approximation algorithms for bandwidth computation that generate orderings of vertices with O(Delta) width. We can thus generate the cube representation of such graphs in O(Delta) dimensions in polynomial time.

Free vibration and wave propagation analysis of uniform and tapered rotating beams using spectrally formulated finite elements

This paper presents a formulation of an approximate spectral element for uniform and tapered rotating Euler-Bernoulli beams. The formulation takes into account the varying centrifugal force, mass and bending stiffness. The dynamic stiffness matrix is constructed using the weak form of the governing differential equation in the frequency domain, where two different interpolating functions for the transverse displacement are used for the element formulation. Both free vibration and wave propagation analysis is performed using the formulated elements. The studies show that the formulated element predicts results, that compare well with the solution available in the literature, at a fraction of the computational effort. In addition, for wave propagation analysis, the element shows superior convergence. (C) 2007 Elsevier Ltd. All rights reserved.

Energy identities in water wave theory for free-surface boundary condition with higher-order derivatives

A modified form of Green's integral theorem is employed to derive the energy identity in any water wave diffraction problem in a single-layer fluid for free-surface boundary condition with higher-order derivatives. For a two-layer fluid with free-surface boundary condition involving higher-order derivatives, two forms of energy identities involving transmission and reflection coefficients for any wave diffraction problem are also derived here by the same method. Based on this modified Green's theorem, hydrodynamic relations such as the energy-conservation principle and modified Haskind–Hanaoka relation are derived for radiation and diffraction problems in a single as well as two-layer fluid.

'Wet N2O oxidation' process and interface state density characterization of nanoscale nitrided SiO2 for flash memory application

In this paper we first present the 'wet N2O' furnace oxidation process to grow nitrided tunnel oxides in the thickness range 6 to 8 nm on silicon at a temperature of 800 degrees C. Electrical characteristics of MOS capacitors and MOSFETs fabricated using this oxide as gate oxide have been evaluated and the superior features of this oxide are ascertained The frequency response of the interface states, before and after subjecting the MOSFET gate oxide to constant current stress, is studied using a simple analytical model developed in this work.

Phase-Field Simulation of Fusion Interface Events during Solidification of Dissimilar Welds: Effect of Composition Inhomogeneity

We investigate the events near the fusion interfaces of dissimilar welds using a phase-field model developed for single-phase solidification of binary alloys. The parameters used here correspond to the dissimilar welding of a Ni/Cu couple. The events at the Ni and the Cu interface are very different, which illustrate the importance of the phase diagram through the slope of the liquidus curves. In the Ni side, where the liquidus temperature decreases with increasing alloying, solutal melting of the base metal takes place; the resolidification, with continuously increasing solid composition, is very sluggish until the interface encounters a homogeneous melt composition. The growth difficulty of the base metal increases with increasing initial melt composition, which is equivalent to a steeper slope of the liquidus curve. In the Cu side, the initial conditions result in a deeply undercooled melt and contributions from both constrained and unconstrained modes of growth are observed. The simulations bring out the possibility of nucleation of a concentrated solid phase from the melt, and a secondary melting of the substrate due to the associated recalescence event. The results for the Ni and Cu interfaces can be used to understand more complex dissimilar weld interfaces involving multiphase solidification.

Specific Small-Molecule Activator of Aurora Kinase A Induces Autophosphorylation in a Cell-Free System

Aurora kinases are essential for chromosomal segregation and cell division and thereby important for maintaining the proper genomic integrity. There are three classes of aurora kinases in humans: A, B, and C. Aurora kinase A is frequently overexpressed in various cancers. The link of the overexpression and tumorigenesis is yet to be understood. By employing virtual screening, we have found that anacardic acid, a pentadecane aliphatic chain containing hydroxylcarboxylic acid, from cashew nut shell liquid could be docked in Aurora kinases A and B. Remarkably, we found that anacardic acid could potently activate the Aurora kinase A mediated phosphorylation of histone H3, but at a similar concentration the activity of aurora kinase B remained unaffected in vitro. Mechanistically, anacardic acid induces the structural changes and also the autophosphorylation of the aurora kinase A to enhance the enzyme activity. This data thus indicate anacardic acid as the first small-molecule activator of Aurora kinase, which could be highly useful for probing the function of hyperactive (overexpressed) Aurora kinase A.