A new parallel overlapped domain decomposition method for nonlinear dynamic finite element analysis

In this paper a new parallel algorithm for nonlinear transient dynamic analysis of large structures has been presented. An unconditionally stable Newmark-beta method (constant average acceleration technique) has been employed for time integration. The proposed parallel algorithm has been devised within the broad framework of domain decomposition techniques. However, unlike most of the existing parallel algorithms (devised for structural dynamic applications) which are basically derived using nonoverlapped domains, the proposed algorithm uses overlapped domains. The parallel overlapped domain decomposition algorithm proposed in this paper has been formulated by splitting the mass, damping and stiffness matrices arises out of finite element discretisation of a given structure. A predictor-corrector scheme has been formulated for iteratively improving the solution in each step. A computer program based on the proposed algorithm has been developed and implemented with message passing interface as software development environment. PARAM-10000 MIMD parallel computer has been used to evaluate the performances. Numerical experiments have been conducted to validate as well as to evaluate the performance of the proposed parallel algorithm. Comparisons have been made with the conventional nonoverlapped domain decomposition algorithms. Numerical studies indicate that the proposed algorithm is superior in performance to the conventional domain decomposition algorithms. (C) 2003 Elsevier Ltd. All rights reserved.

AlignHUSH: Alignment of HMMs using structure and hydrophobicity information

Background: Sensitive remote homology detection and accurate alignments especially in the midnight zone of sequence similarity are needed for better function annotation and structural modeling of proteins. An algorithm, AlignHUSH for HMM-HMM alignment has been developed which is capable of recognizing distantly related domain families The method uses structural information, in the form of predicted secondary structure probabilities, and hydrophobicity of amino acids to align HMMs of two sets of aligned sequences. The effect of using adjoining column(s) information has also been investigated and is found to increase the sensitivity of HMM-HMM alignments and remote homology detection. Results: We have assessed the performance of AlignHUSH using known evolutionary relationships available in SCOP. AlignHUSH performs better than the best HMM-HMM alignment methods and is observed to be even more sensitive at higher error rates. Accuracy of the alignments obtained using AlignHUSH has been assessed using the structure-based alignments available in BaliBASE. The alignment length and the alignment quality are found to be appropriate for homology modeling and function annotation. The alignment accuracy is found to be comparable to existing methods for profile-profile alignments. Conclusions: A new method to align HMMs has been developed and is shown to have better sensitivity at error rates of 10% and above when compared to other available programs. The proposed method could effectively aid obtaining clues to functions of proteins of yet unknown function. A web-server incorporating the AlignHUSH method is available at http://crick.mbu.iisc.ernet.in/similar to alignhush/

The carboxy terminal WD domain of the pre-mRNA splicing factor Prp17p is critical for function

In Saccharomyces cerevisiae, Prp17p is required for the efficient completion of the second step of pre-mRNA splicing. The function and interacting factors for this protein have not been elucidated. We have performed a mutational analysis of yPrp17p to identify protein domains critical for function. A series of deletions were made throughout the region spanning the N-terminal 158 amino acids of the protein, which do not contain any identified structural motifs. The C-terminal portion (amino acids 160–455) contains a WD domain containing seven WD repeats. We determined that a minimal functional Prp17p consists of the WD domain and 40 amino acids N-terminal to it. We generated a three-dimensional model of the WD repeats in Prp17p based on the crystal structure of the [beta]-transducin WD domain. This model was used to identify potentially important amino acids for in vivo functional characterization. Through analysis of mutations in four different loops of Prp17p that lie between [beta] strands in the WD repeats, we have identified four amino acids, 235TETG238, that are critical for function. These amino acids are predicted to be surface exposed and may be involved in interactions that are important for splicing. Temperature-sensitive prp17 alleles with mutations of these four amino acids are defective for the second step of splicing and are synthetically lethal with a U5 snRNA loop I mutation, which is also required for the second step of splicing. These data reinforce the functional significance of this region within the WD domain of Prp17p in the second step of splicing.

Improvement of protein structure comparison using a structural alphabet

The three dimensional structure of a protein provides major insights into its function. Protein structure comparison has implications in functional and evolutionary studies. A structural alphabet (SA) is a library of local protein structure prototypes that can abstract every part of protein main chain conformation. Protein Blocks (PBS) is a widely used SA, composed of 16 prototypes, each representing a pentapeptide backbone conformation defined in terms of dihedral angles. Through this description, the 3D structural information can be translated into a 1D sequence of PBs. In a previous study, we have used this approach to compare protein structures encoded in terms of PBs. A classical sequence alignment procedure based on dynamic programming was used, with a dedicated PB Substitution Matrix (SM). PB-based pairwise structural alignment method gave an excellent performance, when compared to other established methods for mining. In this study, we have (i) refined the SMs and (ii) improved the Protein Block Alignment methodology (named as iPBA). The SM was normalized in regards to sequence and structural similarity. Alignment of protein structures often involves similar structural regions separated by dissimilar stretches. A dynamic programming algorithm that weighs these local similar stretches has been designed. Amino acid substitutions scores were also coupled linearly with the PB substitutions. iPBA improves (i) the mining efficiency rate by 6.8% and (ii) more than 82% of the alignments have a better quality. A higher efficiency in aligning multi-domain proteins could be also demonstrated. The quality of alignment is better than DALI and MUSTANG in 81.3% of the cases. Thus our study has resulted in an impressive improvement in the quality of protein structural alignment. (C) 2011 Elsevier Masson SAS. All rights reserved.

Novel Pentameric Structure of the Diarrhea-Inducing Region of the Rotavirus Enterotoxigenic Protein NSP4

A novel pentameric structure which differs from the previously reported tetrameric form of the diarrhea-inducing region of the rotavirus enterotoxin NSP4 is reported here. A significant feature of this pentameric form is the absence of the calcium ion located in the core region of the tetrameric structures. The lysis of cells, the crystallization of the region spanning residues 95 to 146 of NSP4 (NSP4(95-146)) of strain ST3 (ST3: NSP4(95-146)) at acidic pH, and comparative studies of the recombinant purified peptide under different conditions by size-exclusion chromatography (SEC) and of the crystal structures suggested pH-, Ca(2+)-, and protein concentration-dependent oligomeric transitions in the peptide. Since the NSP4(95-146) mutant lacks the N-terminal amphipathic domain (AD) and most of the C-terminal flexible region (FR), to demonstrate that the pentameric transition is not a consequence of the lack of the N- and C-terminal regions, glutaraldehyde cross-linking of the Delta N72 and Delta N94 mutant proteins, which contain or lack the AD, respectively, but possess the complete C-terminal FR, was carried out. The results indicate the presence of pentamers in preparations of these longer mutants. Detailed SEC analyses of Delta N94 prepared under different conditions, however, revealed protein concentration-dependent but metal ion-and pH-independent pentamer accumulation at high concentrations which dissociated into tetramers and lower oligomers at low protein concentrations. While calcium appeared to stabilize the tetramer, magnesium in particular stabilized the dimer. Delta N72 existed primarily in the multimeric form under all conditions. These findings of a calcium-free NSP4 pentamer and its concentration-dependent and largely calcium-independent oligomeric transitions open up a new dimension in an understanding of the structural basis of its multitude of functions.

Crystal structure of a beta-prism II lectin from Remusatia vivipara

The crystal structure of a beta-prism II (BP2) fold lectin from Remusatia vivipara, a plant of traditional medicinal value, has been determined at a resolution of 2.4 A. This lectin (RVL, Remusatia vivipara lectin) is a dimer with each protomer having two distinct BP2 domains without a linker between them. It belongs to the ``monocot mannose-binding'' lectin family, which consists of proteins of high sequence and structural similarity. Though the overall tertiary structure is similar to that of lectins from snowdrop bulbs and garlic, crucial differences in the mannose-binding regions and oligomerization were observed. Unlike most of the other structurally known proteins in this family, only one of the three carbohydrate recognition sites (CRSs) per BP2 domain is found to be conserved. RVL does not recognize simple mannose moieties. RVL binds to only N-linked complex glycans like those present on the gp120 envelope glycoprotein of HIV and mannosylated blood proteins like fetuin, but not to simple mannose moieties. The molecular basis for these features and their possible functional implications to understand the different levels of carbohydrate affinities in this structural family have been investigated through structure analysis, modeling and binding studies. Apart from being the first structure of a lectin to be reported from the Araceae/Arum family, this protein also displays a novel mode of oligomerization among BP2 lectins.

A new pentameric structure of rotavirus NSP4 revealed by molecular replacement

The region spanning residues 95-146 of the rotavirus nonstructural protein NSP4 from the asymptomatic human strain ST3 has been purified and crystallized and diffraction data have been collected to a resolution of 2.6 angstrom. Several attempts to solve the structure by the molecular-replacement method using the available tetrameric structures of this domain were unsuccessful despite a sequence identity of 73% to the already known structures. A more systematic approach with a dimer as the search model led to an unexpected pentameric structure using the program Phaser. The various steps involved in arriving at this molecular-replacement solution, which unravelled a case of subtle variation between different oligomeric states unknown at the time of solving the structure, are presented in this paper.

Structure and Mechanistic Insights into Novel Iron-mediated Moonlighting Functions of Human J-protein Cochaperone, Dph4

J-proteins are obligate cochaperones of Hsp70s and stimulate their ATPase activity via the J-domain. Although the functions of J-proteins have been well understood in the context of Hsp70s, their additional co-evolved ``physiological functions'' are still elusive. We report here the solution structure and mechanism of novel iron-mediated functional roles of human Dph4, a type III J-protein playing a vital role in diphthamide biosynthesis and normal development. The NMR structure of Dph4 reveals two domains: a conserved J-domain and a CSL-domain connected via a flexible linker-helix. The linker-helix modulates the conformational flexibility between the two domains, regulating thereby the protein function. Dph4 exhibits a unique ability to bind iron in tetrahedral coordination geometry through cysteines of its CSL-domain. The oxidized Fe-Dph4 shows characteristic UV-visible and electron paramagnetic resonance spectral properties similar to rubredoxins. Iron-bound Dph4 (Fe-Dph4) also undergoes oligomerization, thus potentially functioning as a transient ``iron storage protein,'' thereby regulating the intracellular iron homeostasis. Remarkably, Fe-Dph4 exhibits vital redox and electron carrier activity, which is critical for important metabolic reactions, including diphthamide biosynthesis. Further, we observed that Fe-Dph4 is conformationally better poised to perform Hsp70-dependent functions, thus underlining the significance of iron binding in Dph4. Yeast Jjj3, a functional ortholog of human Dph4 also shows a similar iron-binding property, indicating the conserved nature of iron sequestration across species. Taken together, our findings provide invaluable evidence in favor of additional co-evolved specialized functions of J-proteins, previously not well appreciated.

Crystal Structure of Escherichia coli Diaminopropionate Ammonia-lyase Reveals Mechanism of Enzyme Activation and Catalysis

Pyridoxal 5'-phosphate (PLP)-dependent enzymes utilize the unique chemistry of a pyridine ring to carry out diverse reactions involving amino acids. Diaminopropionate (DAP) ammonia-lyase (DAPAL) is a prokaryotic PLP-dependent enzyme that catalyzes the degradation of D-and L-forms of DAP to pyruvate and ammonia. Here, we report the first crystal structure of DAPAL from Escherichia coli (EcDAPAL) in tetragonal and monoclinic forms at 2.0 and 2.2 angstrom resolutions, respectively. Structures of EcDAPAL soaked with substrates were also determined. EcDAPAL has a typical fold type II PLP-dependent enzyme topology consisting of a large and a small domain with the active site at the interface of the two domains. The enzyme is a homodimer with a unique biological interface not observed earlier. Structure of the enzyme in the tetragonal form had PLP bound at the active site, whereas the monoclinic structure was in the apo-form. Analysis of the apo and holo structures revealed that the region around the active site undergoes transition from a disordered to ordered state and assumes a conformation suitable for catalysis only upon PLP binding. A novel disulfide was found to occur near a channel that is likely to regulate entry of ligands to the active site. EcDAPAL soaked with DL-DAP revealed density at the active site appropriate for the reaction intermediate aminoacrylate, which is consistent with the observation that EcDAPAL has low activity under crystallization conditions. Based on the analysis of the structure and results of site-directed mutagenesis, a two-base mechanism of catalysis involving Asp(120) and Lys(77) is suggested.

Crystal Structure of a Monomeric Thiolase-Like Protein Type 1 (TLP1) from Mycobacterium smegmatis

An analysis of the Mycobacterium smegmatis genome suggests that it codes for several thiolases and thiolase-like proteins. Thiolases are an important family of enzymes that are involved in fatty acid metabolism. They occur as either dimers or tetramers. Thiolases catalyze the Claisen condensation of two acetyl-Coenzyme A molecules in the synthetic direction and the thiolytic cleavage of 3-ketoacyl-Coenzyme A molecules in the degradative direction. Some of the M. smegmatis genes have been annotated as thiolases of the poorly characterized SCP2-thiolase subfamily. The mammalian SCP2-thiolase consists of an N-terminal thiolase domain followed by an additional C-terminal domain called sterol carrier protein-2 or SCP2. The M. smegmatis protein selected in the present study, referred to here as the thiolase-like protein type 1 (MsTLP1), has been biochemically and structurally characterized. Unlike classical thiolases, MsTLP1 is a monomer in solution. Its structure has been determined at 2.7 angstrom resolution by the single wavelength anomalous dispersion method. The structure of the protomer confirms that the N-terminal domain has the thiolase fold. An extra C-terminal domain is indeed observed. Interestingly, it consists of six beta-strands forming an anti-parallel beta-barrel which is completely different from the expected SCP2-fold. Detailed sequence and structural comparisons with thiolases show that the residues known to be essential for catalysis are not conserved in MsTLP1. Consistent with this observation, activity measurements show that MsTLP1 does not catalyze the thiolase reaction. This is the first structural report of a monomeric thiolase-like protein from any organism. These studies show that MsTLP1 belongs to a new group of thiolase related proteins of unknown function.

Interaction Signatures Stabilizing the NAD(P)-Binding Rossmann Fold: A Structure Network Approach

The fidelity of the folding pathways being encoded in the amino acid sequence is met with challenge in instances where proteins with no sequence homology, performing different functions and no apparent evolutionary linkage, adopt a similar fold. The problem stated otherwise is that a limited fold space is available to a repertoire of diverse sequences. The key question is what factors lead to the formation of a fold from diverse sequences. Here, with the NAD(P)-binding Rossmann fold domains as a case study and using the concepts of network theory, we have unveiled the consensus structural features that drive the formation of this fold. We have proposed a graph theoretic formalism to capture the structural details in terms of the conserved atomic interactions in global milieu, and hence extract the essential topological features from diverse sequences. A unified mathematical representation of the different structures together with a judicious concoction of several network parameters enabled us to probe into the structural features driving the adoption of the NAD(P)-binding Rossmann fold. The atomic interactions at key positions seem to be better conserved in proteins, as compared to the residues participating in these interactions. We propose a ``spatial motif'' and several ``fold specific hot spots'' that form the signature structural blueprints of the NAD(P)-binding Rossmann fold domain. Excellent agreement of our data with previous experimental and theoretical studies validates the robustness and validity of the approach. Additionally, comparison of our results with statistical coupling analysis (SCA) provides further support. The methodology proposed here is general and can be applied to similar problems of interest.

N-Terminal disordered domain of saccharomyces cerevisiae Hop1 protein Is dispensable for DNA binding, bridging, and synapsis of double-stranded DNA molecules but is ecessary for spore formation

The cytological architecture of the synaptonemal complex (SC), a meiosis-specific proteinaceous structure, is evolutionarily conserved among eukaryotes. However, little is known about the biochemical properties of SC components or the mechanisms underlying their roles in meiotic chromosome synapsis and recombination. Functional analysis of Saccharomyces cerevisiae Hop1, a key structural component of SC, has begun to reveal important insights into its function in interhomolog recombination. Previously, we showed that Hop1 is a structure-specific DNA-binding protein, exhibits higher binding affinity for the Holliday junction, and induces structural distortion at the core of the junction. Furthermore, Hop1 promotes DNA condensation and intra- and intermolecular synapsis between duplex DNA molecules. Here, we show that Hop1 possesses a modular domain organization, consisting of an intrinsically disordered N-terminal domain and a protease-resistant C-terminal domain (Hop1CTD). Furthermore, we found that Hop1CTD exhibits strong homotypic as well as heterotypic protein protein interactions, and its biochemical activities were similar to those of the full-length Hop1 protein. However, Hop1CTD failed to complement the meiotic recombination defects of the Delta hop1 strain, indicating that both N- and C-terminal domains of Hop1 are essential for meiosis and spore formation. Altogether, our findings reveal novel insights into the structure-function relationships of Hop1 and help to further our understanding of its role in meiotic chromosome synapsis and recombination.

A Guided Wave approach for the detection of damage in a structure having elements with periodic damage

A wave propagation based approach for the detection of damage in components of structures having periodic damage has been proposed. Periodic damage pattern may arise in a structure due to periodicity in geometry and in loading. The method exploits the Block-Floquet band formation mechanism, a feature specific to structures with periodicity, to identify propagation bands (pass bands) and attenuation bands (stop bands) at different frequency ranges. The presence of damage modifies the wave propagation behaviour forming these bands. With proper positioning of sensors a damage force indicator (DFI) method can be used to locate the defect at an accuracy level of sensor to sensor distance. A wide range of transducer frequency may be used to obtain further information about the shape and size of the damage. The methodology is demonstrated using a few 1-D structures with different kinds of periodicity and damage. For this purpose, dynamic stiffness matrix is formed for the periodic elements to obtain the dispersion relationship using frequency domain spectral element and spectral super element method. The sensitivity of the damage force indicator for different types of periodic damages is also analysed.

Understanding the role of domain-domain linkers in the spatial orientation of domains in multi-domain proteins

Inter-domain linkers (IDLs)' bridge flanking domains and support inter-domain communication in multi-domain proteins. Their sequence and conformational preferences enable them to carry out varied functions. They also provide sufficient flexibility to facilitate domain motions and, in conjunction with the interacting interfaces, they also regulate the inter-domain geometry (IDG). In spite of the basic intuitive understanding of the inter-domain orientations with respect to linker conformations and interfaces, we still do not entirely understand the precise relationship among the three. We show that IDG is evolutionarily well conserved and is constrained by the domain-domain interface interactions. The IDLs modulate the interactions by varying their lengths, conformations and local structure, thereby affecting the overall IDG. Results of our analysis provide guidelines in modelling of multi-domain proteins from the tertiary structures of constituent domain components.

Investigations of Ramachandran disallowed conformations in protein domain families

In peptide and protein structures, occurrence of (phi,psi.) angles in the disallowed region of the Ramachandran map almost always suggests local regions of error or poor accuracy. However, very rarely genuine disallowed conformations occur as noted in the current study in proteins of known structure available at ultra-high resolution (<= 1.2 (A) over circle). In the current work, extent of conservation of genuine disallowed conformations in homologous proteins of known structures has been analyzed. From a dataset of 124 protein domain families, with structure of at least one constituent member in each family available at a resolution of 1.2 (A) over circle or better, we have analyzed the conservation of 221 disallowed conformations. It is observed that the disallowed conformation is only moderately conservedin protein domain families. In the gross dataset no particular residue type adopting disallowed conformation elicit high conservation of residue type though there are alignment positions in the dataset with complete conservation of both the residue type and the disallowed conformation. Conserved disallowed conformation in protein domain families play biologically significant role in roughly 50% of the cases. The residues with the disallowed conformation or its flanking residues are often located within or around the functional site of the protein. (C) 2013 Elsevier B.V. All rights reserved.