The Critical Role of N- and C-Terminal Contact in Protein Stability and Folding of a Family 10 Xylanase under Extreme Conditions

Background: Stabilization strategies adopted by proteins under extreme conditions are very complex and involve various kinds of interactions. Recent studies have shown that a large proportion of proteins have their N- and C-terminal elements in close contact and suggested they play a role in protein folding and stability. However, the biological significance of this contact remains elusive. Methodology: In the present study, we investigate the role of N- and C-terminal residue interaction using a family 10 xylanase (BSX) with a TIM-barrel structure that shows stability under high temperature,alkali pH, and protease and SDS treatment. Based on crystal structure,an aromatic cluster was identified that involves Phe4, Trp6 and Tyr343 holding the Nand C-terminus together; this is a unique and important feature of this protein that might be crucial for folding and stabilityunder poly-extreme conditions. Conclusion: A series of mutants was created to disrupt this aromatic cluster formation and study the loss of stability and function under given conditions. While the deletions of Phe4 resulted in loss of stability, removal of Trp6 and Tyr343 affected in vivo folding and activity. Alanine substitution with Phe4, Trp6 and Tyr343 drastically decreased stability under all parameters studied. Importantly,substitution of Phe4 with Trp increased stability in SDS treatment.Mass spectrometry results of limited proteolysis further demonstrated that the Arg344 residue is highly susceptible to trypsin digestion in sensitive mutants such as DF4, W6A and Y343A, suggesting again that disruption of the Phe4-Trp6-Tyr343 (F-W-Y) cluster destabilizes the N-and C-terminal interaction. Our results underscore the importance of N- and C-terminal contact through aromatic interactions in protein folding and stability under extreme conditions, and these results may be useful to improve the stability of other proteins under suboptimal conditions.

Protein Structure and Function: Looking through the Network of Side-Chain Interactions

Network theory has become an excellent method of choice through which biological data are smoothly integrated to gain insights into complex biological problems. Understanding protein structure, folding, and function has been an important problem, which is being extensively investigated by the network approach. Since the sequence uniquely determines the structure, this review focuses on the networks of non-covalently connected amino acid side chains in proteins. Questions in structural biology are addressed within the framework of such a formalism. While general applications are mentioned in this review, challenging problems which have demanded the attention of scientific community for a long time, such as allostery and protein folding, are considered in greater detail. Our aim has been to explore these important problems through the eyes of networks. Various methods of constructing protein structure networks (PSN) are consolidated. They include the methods based on geometry, edges weighted by different schemes, and also bipartite network of protein-nucleic acid complexes. A number of network metrics that elegantly capture the general features as well as specific features related to phenomena, such as allostery and protein model validation, are described. Additionally, an integration of network theory with ensembles of equilibrium structures of a single protein or that of a large number of structures from the data bank has been presented to perceive complex phenomena from network perspective. Finally, we discuss briefly the capabilities, limitations, and the scope for further explorations of protein structure networks.

Intra and Inter-Molecular Communications Through Protein Structure Network

Communication within and across proteins is crucial for the biological functioning of proteins. Experiments such as mutational studies on proteins provide important information on the amino acids, which are crucial for their function. However, the protein structures are complex and it is unlikely that the entire responsibility of the function rests on only a few amino acids. A large fraction of the protein is expected to participate in its function at some level or other. Thus, it is relevant to consider the protein structures as a completely connected network and then deduce the properties, which are related to the global network features. In this direction, our laboratory has been engaged in representing the protein structure as a network of non-covalent connections and we have investigated a variety of problems in structural biology, such as the identification of functional and folding clusters, determinants of quaternary association and characterization of the network properties of protein structures. We have also addressed a few important issues related to protein dynamics, such as the process of oligomerization in multimers, mechanism on protein folding, and ligand induced communications (allosteric effect). In this review we highlight some of the investigations which we have carried out in the recent past. A review on protein structure graphs was presented earlier, in which the focus was on the graphs and graph spectral properties and their implementation in the study of protein structure graphs/networks (PSN). In this article, we briefly summarize the relevant parts of the methodology and the focus is on the advancement brought out in the understanding of protein structure-function relationships through structure networks. The investigations of structural/biological problems are divided into two parts, in which the first part deals with the analysis of PSNs based on static structures obtained from x-ray crystallography. The second part highlights the changes in the network, associated with biological functions, which are deduced from the network analysis on the structures obtained from molecular dynamics simulations.

Similarities in the Genomic Sequence and Coat Protein-Structure of Plant Viruses

The genomic sequences of several RNA plant viruses including cucumber mosaic virus, brome mosaic virus, alfalfa mosaic virus and tobacco mosaic virus have become available recently. The former two viruses are icosahedral while the latter two are bullet and rod shaped, respectively in particle morphology. The non-structural 3a proteins of cucumber mosaic virus and brome mosaic virus have an amino acid sequence homology of 35% and hence are evolutionarily related. In contrast, the coat proteins exhibit little homology, although the circular dichroism spectrum of these viruses are similar. The non-coding regions of the genome also exhibit variable but extensive homology. Comparison of the brome mosaic virus and alfalfa mosaic virus sequences reveals that they are probably related although with a much larger evolutionary distance. The polypeptide folds of the coat protein of three biologically distinct isometric plant viruses, tomato bushy stunt virus, southern bean mosaic virus and satellite tobacco necrosis virus have been shown to display a striking resemblance. All of them consist of a topologically similar 8-standard β-barrel. The implications of these studies to the understanding of the evolution of plant viruses will be discussed.

Allostery in tRNA Synthetases Elucidated from MD Simulations and Protein Structure Networks

tRNA synthetases (aaRS) are enzymes crucial in the translation of genetic code. The enzyme accylates the acceptor stem of tRNA by the congnate amino acid bound at the active site, when the anti-codon is recognized by the anti-codon site of aaRS. In a typical aaRS, the distance between the anti-codon region and the amino accylation site is approximately 70 Å. We have investigated this allosteric phenomenon at molecular level by MD simulations followed by the analysis of protein structure networks (PSN) of non-covalent interactions. Specifically, we have generated conformational ensembles by performing MD simulations on different liganded states of methionyl tRNA synthetase (MetRS) from Escherichia coli and tryptophenyl tRNA synthetase (TrpRS) from Human. The correlated residues during the MD simulations are identified by cross correlation maps. We have identified the amino acids connecting the correlated residues by the shortest path between the two selected members of the PSN. The frequencies of paths have been evaluated from the MD snapshots[1]. The conformational populations in different liganded states of the protein have been beautifully captured in terms of network parameters such as hubs, cliques and communities[2]. These parameters have been associated with the rigidity and plasticity of the protein conformations and can be associated with free energy landscape. A comparison of allosteric communication in MetRS and TrpRS [3] elucidated in this study highlights diverse means adopted by different enzymes to perform a similar function. The computational method described for these two enzymes can be applied to the investigation of allostery in other systems.

Empirical torsional potential functions from protein structure data. Phi- and psi-potentials for non-glycyl amino acid residues.

The torsional potential functions Vt(phi) and Vt(psi) around single bonds N--C alpha and C alpha--C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (phi, psi)-plane with the value of Vtot(phi, psi), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in L-configuration, are Vt(phi) = 1.0 cos (phi + 60 degrees); Vt(psi) = 0.5 cos (psi + 60 degrees) - 1.0 cos (2 psi + 30 degrees) - 0.5 cos (3 psi + 30 degrees). The dipeptide energy maps Vtot(phi, psi) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line psi = 0 degrees. These functions derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

Protein Hydration and Water Structure: X-ray Analysis of a Closely Packed Protein Crystal with Very Low Solvent Content

Low-humidity monoclinic lysozyme, resulting from a water-mediated transformation, has one of the lowest solvent contents (22% by volume) observed in a protein crystal. Its structure has been solved by the molecular replacement method and refined to an R value of 0.175 for 7684 observed reflections in the 10–1.75 Å resolution shell. 90% of the solvent in the well ordered crystals could be located. Favourable sites of hydration on the protein surface include side chains with multiple hydrogen-bonding centres, and regions between short hydrophilic side chains and the main-chain CO or NH groups of the same or nearby residues. Major secondary structural features are not disrupted by hydration. However, the free CO groups at the C terminii and, to a lesser extent, the NH groups at the N terminii of helices provide favourable sites for water interactions, as do reverse turns and regions which connect β-structure and helices. The hydration shell consists of discontinuous networks of water molecules, the maximum number of molecules in a network being ten. The substrate-binding cleft is heavily hydrated, as is the main loop region which is stabilized by water interactions. The protein molecules are close packed in the crystals with a molecular coordination number of 14. Arginyl residues are extensively involved in intermolecular hydrogen bonds and water bridges. The water molecules in the crystal are organized into discrete clusters. A distinctive feature of the clusters is the frequent occurrence of three-membered rings. The protein molecules undergo substantial rearrangement during the transformation from the native to the low-humidity form. The main-chain conformations in the two forms are nearly the same, but differences exist in the side-chain conformation. The differences are particularly pronounced in relation to Trp 62 and Trp 63. The shift in Trp 62 is especially interesting as it is also known to move during inhibitor binding.

Structure and dynamics of protein excited states with millisecond lifetimes

An in-depth understanding of biological processes often requires detailed atomic resolution structures of the molecules involved. However in solution where most of these processes occur the conformation of biomolecules like RNA, DNA and proteins is not static but fluctuates. Routinely used structural techniques like X-ray crystallography, NMR spectroscopy and cryo-electron microscopy have almost always been used to determine the structure of the dominant conformation or obtain an average structure of the biomolecule in solution with very little detailed information regarding the dynamics of these molecules in solution. Over the last few years, NMR based methods have been developed to study the dynamics of these biomolecules in solution in a site-specific manner with the aim of generating structures of the different conformations that these molecules can adopt in solution. One powerful technique is the Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiment, which can be used to detect and characterize protein excited states that are populated for as less as 0.5% of the time with ∼0.5–10 millisecond lifetimes. Due to recent advances in NMR pulse sequences and labeling methodology, it is now possible to determine the structures of these transiently populated excited states with millisecond lifetimes by obtaining accurate chemical shifts, residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) of these excited states. In these excited states the dynamics of some methyl containing residues can also be studied.

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.

Single-Molecule DNA Analysis Reveals That Yeast Hop1 Protein Promotes DNA Folding and Synapsis: Implications for Condensation of Meiotic Chromosomes

During meiosis, long-range interaction between homologous chromosomes is thought to be crucial for homology recognition, exchange of DNA strands, and production of normal haploid gametes. However, little is known about the identity of the proteins involved and the actual molecular mechanism(s) by which chromosomes recognize and recombine with their appropriate homologous partners. Single-molecule analyses have the potential to provide insights into our understanding of this fascinating and long-standing question. Using atomic force microscopy and magnetic tweezers techniques, we discovered that Hop1 protein, a key structural component of Saccharomyces cerevisiae synaptonemal complex, exhibits the ability to bridge noncontiguous DNA segments into intramolecular stem-loop structures in which the DNA segments appear to be fully synapsed within the filamentous protein stems. Additional evidence suggests that Hop1 folds DNA into rigid protein DNA filaments and higher-order nucleoprotein structures. Importantly, Hop1 promotes robust intra- and intermolecular synapsis between double-stranded DNA molecules, suggesting that juxtaposition of DNA sequences may assist in strand exchange between homologues by recombination-associated proteins. Finally, the evidence from ensemble experiments is consistent with the notion that Hop1 causes rigidification of DNA molecules. These results provide the first direct evidence for long-range protein-mediated DNA DNA synapsis, independent of crossover recombination, which is presumed to occur during meiotic recombination.

Effect of signal peptide on stability and folding of escherichia coli thioredoxin

The signal peptide plays a key role in targeting and membrane insertion of secretory and membrane proteins in both prokaryotes and eukaryotes. In E. coli, recombinant proteins can be targeted to the periplasmic space by fusing naturally occurring signal sequences to their N-terminus. The model protein thioredoxin was fused at its N-terminus with malE and pelB signal sequences. While WT and the pelB fusion are soluble when expressed, the malE fusion was targeted to inclusion bodies and was refolded in vitro to yield a monomeric product with identical secondary structure to WT thioredoxin. The purified recombinant proteins were studied with respect to their thermodynamic stability, aggregation propensity and activity, and compared with wild type thioredoxin, without a signal sequence. The presence of signal sequences leads to thermodynamic destabilization, reduces the activity and increases the aggregation propensity, with malE having much larger effects than pelB. These studies show that besides acting as address labels, signal sequences can modulate protein stability and aggregation in a sequence dependent manner.

Dissecting pi-helices: sequence, structure and function

A new procedure for the identification of regular secondary structures using a C-alpha trace has identified 659 pi-helices in 3582 protein chains, solved at high resolution. Taking advantage of this significantly expanded database of pi-helices, we have analysed the functional and structural roles of helices and determined the position-wise amino acid propensity within and around them. These helices range from 5 to 18 residues in length with the average twist and rise being 85.2 +/- 7.2 and 1.28 +/- 0.31 angstrom, respectively. A total of 546 (similar to 83%) out of 659 pi-helices occur in conjunction with alpha-helices, with 101 pi-helices being interspersed between two alpha-helices. The majority of interspersed pi-helices were found to be conserved across a large number of structures within a protein family and produce a significant bend in the overall helical segment as well as local distortions in the neighbouring a-helices. The presence of a pi-helical fragment leads to appropriate orientation of the constituent residues, so as to facilitate favourable interactions and also help in proper folding of the protein chain. In addition to intra helical 6 -> 1 N H center dot center dot center dot O hydrogen bonds, pi-helices are also stabilized by several other non-bonded interactions. pi-Helices show distinct positional residue preferences, which are different from those of a-helices.

Understanding Protein Structure from a Percolation Perspective

Underlying the unique structures and diverse functions of proteins area vast range of amino-acid sequences and a highly limited number of folds taken up by the polypeptide backbone. By investigating the role of noncovalent connections at the backbone level and at the detailed side-chain level, we show that these unique structures emerge from interplay between random and selected features. Primarily, the protein structure network formed by these connections shows simple (bond) and higher order (clique) percolation behavior distinctly reminiscent of random network models. However, the clique percolation specific to the side-chain interaction network bears signatures unique to proteins characterized by a larger degree of connectivity than in random networks. These studies reflect some salient features of the manner in which amino acid sequences select the unique structure of proteins from the pool of a limited number of available folds.

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.