mulPBA: an efficient multiple protein structure alignment method based on a structural alphabet

The increasing number of available protein structures requires efficient tools for multiple structure comparison. Indeed, multiple structural alignments are essential for the analysis of function, evolution and architecture of protein structures. For this purpose, we proposed a new web server called multiple Protein Block Alignment (mulPBA). This server implements a method based on a structural alphabet to describe the backbone conformation of a protein chain in terms of dihedral angles. This sequence-like' representation enables the use of powerful sequence alignment methods for primary structure comparison, followed by an iterative refinement of the structural superposition. This approach yields alignments superior to most of the rigid-body alignment methods and highly comparable with the flexible structure comparison approaches. We implement this method in a web server designed to do multiple structure superimpositions from a set of structures given by the user. Outputs are given as both sequence alignment and superposed 3D structures visualized directly by static images generated by PyMol or through a Jmol applet allowing dynamic interaction. Multiple global quality measures are given. Relatedness between structures is indicated by a distance dendogram. Superimposed structures in PDB format can be also downloaded, and the results are quickly obtained. mulPBA server can be accessed at www.dsimb.inserm.fr/dsimb_tools/mulpba/.

Fourier transform of the collagen triple-helical structure and its significance

The Fourier transforms of the collagen molecular structure have been calculated taking into consideration various side chain atoms, as well as the presence of bound water molecules. There is no significant change in the calculated intensity distribution on including the side chain atoms of non-imino-acid residues. Taking into account the presence of about two bound water molecules per tripeptide unit, the agreement with the observed x-ray pattern is slightly improved. Fourier transforms have also been calculated for the detailed molecular geometries proposed from other laboratories. It is found that there are no major differences between them, as compared to our structure, either in the positions of peak intensity or in the intensity distribution. Hence it is not possible to judge the relative merits of the various molecular geometries for the collagen triple helix from a comparison of the calculated transforms with the meagre data available from its x-ray fibre pattern. It is also concluded that the collagen molecular structure should be regarded as a somewhat flexible chain structure, capable of adapting itself to the requirements of the different side groups which occur in each local region.

The influence of esterification of carboxyl groups of ribonuclease-a on its structure and immunological activity

The effect of modification of carboxyl groups of Ribonuclease-Aa on the enzymatic activity and the antigenic structure of the protein has been studied. Modification of four of the eleven free carboxyl groups of the protein by esterification in anhydrous methanol/0.1 M hydrochloric acid resulted in nearly 80% loss in enzymatic activity but had very little influence on the antigenic structure of the protein. Further increases in the modification of the carboxyl groups caused a progressive loss in immunological activity, and the fully methylated RNase-A exhibited nearly 30% immunological activity. Concomitant with this change in the antigenic structure of the protein, the ability of the molecule to complement with RNase-S-protein increased, clearly indicating the unfolding of the peptide "tail" from the remainder of the molecule. The susceptibility to proteolysis, accessibility of methionine residues for orthobenzoquinone reaction and the loss in immunological activity of the more extensively esterified derivatives of RNase-A are suggestive of the more flexible conformation of these derivatives as compared with the compact native conformation. The fact that even the fully methylated RNase-A retains nearly 30% of its immunological activity suggested that the modified protein contained antibody recognizable residual native structure, which presumably accommodates some antigenic determinants.

Structure, dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X-ray studies

Molecular dynamics simulations have been carried out on all the jacalin-carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X-ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin-carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate-binding residues are consistent with the known thermodynamic parameters of jacalin-carbohydrate interactions. The simulations, along with X-ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin-carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit. Proteins 2009; 77:760-777.

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.

Asymptotic analysis for the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell: The beam mode

Using asymptotics, the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell vibrating in the beam mode (viz. circumferential wave order n = 1) are studied. Initially, the uncoupled wavenumbers of the acoustic fluid and the cylindrical shell structure are discussed. Simple closed form expressions for the structural wavenumbers (longitudinal, torsional and bending) are derived using asymptotic methods for low- and high-frequencies. It is found that at low frequencies the cylinder in the beam mode behaves like a Timoshenko beam. Next, the coupled dispersion equation of the system is rewritten in the form of the uncoupled dispersion equation of the structure and the acoustic fluid, with an added fluid-loading term involving a parameter mu due to the coupling. An asymptotic expansion involving mu is substituted in this equation. Analytical expressions are derived for the coupled wavenumbers (as modifications to the uncoupled wavenumbers) separately for low- and high-frequency ranges and further, within each frequency range, for large and small values of mu. Only the flexural wavenumber, the first rigid duct acoustic cut-on wavenumber and the first pressure-release acoustic cut-on wavenumber are considered. The general trend found is that for small mu, the coupled wavenumbers are close to the in vacuo structural wavenumber and the wavenumbers of the rigid-acoustic duct. With increasing mu, the perturbations increase, until the coupled wavenumbers are better identified as perturbations to the pressure-release wavenumbers. The systematic derivation for the separate cases of small and large mu gives more insight into the physics and helps to continuously track the wavenumber solutions as the fluid-loading parameter is varied from small to large values. Also, it is found that at any frequency where two wavenumbers intersect in the uncoupled analysis, there is no more an intersection in the coupled case, but a gap is created at that frequency. This method of asymptotics is simple to implement using a symbolic computation package (like Maple). (C) 2008 Elsevier Ltd. All rights reserved.

Asymptotic analysis for the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell: The axisymmetric mode

The coupled wavenumbers of a fluid-filled flexible cylindrical shell vibrating in the axisymmetric mode are studied. The coupled dispersion equation of the system is rewritten in the form of the uncoupled dispersion equation of the structure and the acoustic fluid, with an added fluid-loading term involving a parameter e due to the coupling. Using the smallness of Poisson's ratio (v), a double-asymptotic expansion involving e and v 2 is substituted in this equation. Analytical expressions are derived for the coupled wavenumbers (for large and small values of E). Different asymptotic expansions are used for different frequency ranges with continuous transitions occurring between them. The wavenumber solutions are continuously tracked as e varies from small to large values. A general trend observed is that a given wavenumber branch transits from a rigidwalled solution to a pressure-release solution with increasing E. Also, it is found that at any frequency where two wavenumbers intersect in the uncoupled analysis, there is no more an intersection in the coupled case, but a gap is created at that frequency. Only the axisymmetric mode is considered. However, the method can be extended to the higher order modes.

Three-dimensional structure of heat shock protein 90 from Plasmodium falciparum: molecular modelling approach to rational drug design against malaria

We have recently implicated heat shock protein 90 from Plasmodium falciparum (PfHsp90) as a potential drug target against malaria. Using inhibitors specific to the nucleotide binding domain of Hsp90, we have shown potent growth inhibitory effects on development of malarial parasite in human erythrocytes. To gain better understanding of the vital role played by PfHsp90 in parasite growth, we have modeled its three dimensional structure using recently described full length structure of yeast Hsp90. Sequence similarity found between PfHsp90 and yeast Hsp90 allowed us to model the core structure with high confidence. The superimposition of the predicted structure with that of the template yeast Hsp90 structure reveals an RMSD of 3.31 angstrom. The N-terminal and middle domains showed the least RMSD (1.76 angstrom) while the more divergent C-terminus showed a greater RMSD (2.84 angstrom) with respect to the template. The structure shows overall conservation of domains involved in nucleotide binding, ATPase activity, co-chaperone binding as well as inter-subunit interactions. Important co-chaperones known to modulate Hsp90 function in other eukaryotes are conserved in malarial parasite as well. An acidic stretch of amino acids found in the linker region, which is uniquely extended in PfHsp90 could not be modeled in this structure suggesting a flexible conformation. Our results provide a basis to compare the overall structure and functional pathways dependent on PfHsp90 in malarial parasite. Further analysis of differences found between human and parasite Hsp90 may make it possible to design inhibitors targeted specifically against malaria.

Conformation of the flexible bent helix of Leu1-zervamicin in crystal C and a possible gating action for ion passage

The membrane channel-forming polypeptide, Leu(1)-zervamicin, Ac-Leu-Ile-Gln-Iva-Ile(5)-Thr-Aib-Leu-Aib-Hyp(10) -Gln-Aib-Hyp-Aib-Pro(15)-Phol (Aib: alpha-aminoisobutyric acid; Iva: isovaline; Hyp: 4-hydroxyproline; Phol: phenylalininol) has been analyzed by x-ray diffraction in a third crystal form. Although the bent helix is quite similar to the conformations found in crystals A and B, the amount of bending is more severe with a bending angle approximate to 47 degrees, The water channel formed by the convex polar faces of neighboring helices is larger at the mouth than in crystals A and B, and the water sites have become disordered. The channel is interrupted in the middle by a hydrogen bond between the OH of Hyp(10) and the NH2 of the Gln(11) of a neighboring molecule. The side chain of Gln(11) is wrapped around the helix backbone in an unusual fashion in order that it can augment the polar side of the helix. In the present crystal C there appears to be an additional conformation for the Gln(11) side chain (with approximate to 20% occupancy) that opens the channel for possible ion passage. Structure parameters for C85H140N18O22.xH(2)O.C2H5OH are space group P2(1)2(1)2(1), a = 10.337 (2) Angstrom, b = 28.387 (7) Angstrom, c = 39.864 (11) Angstrom, Z = 4, agreement factor R = 12.99% for 3250 data observed > 3 sigma(F), resolution = 1.2 Angstrom. (C) 1994 John Wiley & Sons, Inc.

Multiple support excitations of open-plane frames by a filtered white noise and soil-structure interaction

Seismic structural design is essentially the estimation of structural response to a forced motion, which may be deterministic or stochastic, imposed on the ground. The assumption that the same ground motion acts at every point of the base of the structure (or at every support) is not always justifiable; particularly in case of very large structures when considerable spatial variability in ground motion can exist over significant distances example long span bridges. This variability is partly due to the delay in arrival of the excitation at different supports (which is called the wave passage effect) and due to heterogeneity in the ground medium which results in incoherency and local effects. The current study examines the influence of the wave passage effect (in terms of delay in arrival of horizontal ground excitation at different supports and neglecting transmission through the structure) on the response of a few open-plane frame building structures with soil-structure interaction. The ground acceleration has been modeled by a suitably filtered white noise. As a special case, the ground excitation at different supports has also been treated as statistically independent to model the extreme case of incoherence due to local effects and due to modifications to the ground motion resulting from wave reflections and refractions in heterogeneous soil media. The results indicate that, even for relatively short spanned building frames, wave passage effect can be significant. In the absence of soil-structure interaction, it can significantly increase the root mean square (rms) value of the shear in extreme end columns for the stiffer frames but has negligible effect on the flexible frames when total displacements are considered. It is seen that pseudo-static displacements increasingly contribute to the rms value of column shear as the time delay increases both for the stiffer and for the more flexible frames. When soil-structure interaction is considered, wave passage effect (in terms of total displacements) is significant only for low soil shear modulus, G. values (where soil-structure interaction significantly lowers the fundamental frequency) and for stiff frames. The contribution of pseudo-static displacement to these rms values is found to decrease with increase in G. In general, wave passage effect for most interactive frames is insignificant compared to the attenuating effect a decrease in G, has on the response of the interactive structure to uniform support excitation. When the excitations at different supports are statistically independent, it is seen that for both the stiff and flexible frames, the rms value of the column shear in extreme end columns is several times larger (more for the stiffer frames) than the value corresponding to uniform base excitation with the pseudo-static displacements contributing over 99% of the rms value of column shear. Soil-structure interaction has an attenuating effect on the rms value of the column shear, the effect decreasing with increase in G,. Here too, the pseudo-static displacements contribute very largely to the column shear. The influence of the wave passage effect on the response of three 2-bay frames with and without soil-structure interaction to a recorded horizontal accelerogram is also examined. (C) 2010 Elsevier Ltd. All rights reserved.

Rigid and flexible region in lysozyme and the invariant features in its hydration shell

Water-mediated transformations provide a useful handle for exploring the flexibility in protein molecules and the invariant features in their hydration shells. Low-humidity monoclinic hen egg white lysozyme, resulting from such a transformation, has perhaps the lowest solvent content observed in any protein crystal so far and has a well-ordered structure. A detailed comparison involving this structure, low-humidity tetragonal lysozyme, and the other available refined crystal structures of the enzyme permits the delineation of the relatively rigid, moderately flexible and highly flexible regions of the molecule. The relatively rigid region forms a contiguous structural unit close to the molecular centroid and encompasses parts of of the main beta-structure and three alpha-helices. The hydration shell of the protein contains 30 invariant water molecules. Many of them are involved in holding different parts of the molecule together or in stabilizing local structure. Five of the six invariant water molecules attached to the substrate-binding region form part of a water cluster contiguous with the side-chains of the catalytic residues Glu-35 and Asp-52.

Synthesis of path generating flexible-link mechanisms

Flexible-link mechanisms are those linkage mechanisms (or structures) which are capable of motion by virtue of elastic deformation of one or more;links. In such mechanisms a single flexible link; can replace several rigid links and joints resulting in fewer links, fewer pin joints, reduced overall weight and reduced mechanical error. In spite of such clear advantages, contributions toward flexible-link mechanisms remain very scarce. The area of flexible-link mechanisms offers much scope for further exploration. This paper attempts to show the potential of flexible-link mechanisms in accomplishing a kinematic task like path generation. Synthesis of a four-bar mechanism with a flexible rocker for circular and straight line path generation is carried out. Displacement analysis of the structure is carried out using finite element method (FEM) and synthesis is formulated and solved as an optimization problem. Several numerical examples are presented for illustration. Based on the results obtained with these examples, the flexible-link mechanism considered shows good promise for-path generation.

Solid state structural and solution studies on the formation of a flexible cavity for anions by copper(I) and 1,2-bis(diphenylphosphino)ethane

Copper(l) complexes of 1,2-bis(diphenylphosphino)ethane (dppe) with a stoichiometry Cu-2(dppe)(3)(X)(2) [X- = CN- (1), SCN- (2), NO3- (3)] are obtained from direct reactions of CuX and dppe. The complexes are structurally and spectroscopically (NMR and IR) characterized. The structure of the [Cu-2(dPPe)(3)](2+) dication is similar to the structural motif observed in many other complexes with a chelating dppe and a bridging dppe connecting two copper centers. In complexes 1 -3, the anions are confined to the cavity formed by the phosphines which force a monodentate coordination mode despite the predominant bidentate/bridging character of the anions. The coordination angles rather than the thermochemical radii dictate the steric requirement of anions. While the solution behavior of 3, with nitrate, is similar to complexes studied earlier, complexes with pseudohalides exhibit new solution behavior. (C) 2002 Elsevier Science Ltd. 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.

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.