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 Block Expert (PBE): a web-based protein structure analysis server using a structural alphabet

Encoding protein 3D structures into 1D string using short structural prototypes or structural alphabets opens a new front for structure comparison and analysis. Using the well-documented 16 motifs of Protein Blocks (PBs) as structural alphabet, we have developed a methodology to compare protein structures that are encoded as sequences of PBs by aligning them using dynamic programming which uses a substitution matrix for PBs. This methodology is implemented in the applications available in Protein Block Expert (PBE) server. PBE addresses common issues in the field of protein structure analysis such as comparison of proteins structures and identification of protein structures in structural databanks that resemble a given structure. PBE-T provides facility to transform any PDB file into sequences of PBs. PBE-ALIGNc performs comparison of two protein structures based on the alignment of their corresponding PB sequences. PBE-ALIGNm is a facility for mining SCOP database for similar structures based on the alignment of PBs. Besides, PBE provides an interface to a database (PBE-SAdb) of preprocessed PB sequences from SCOP culled at 95% and of all-against-all pairwise PB alignments at family and superfamily levels. PBE server is freely available at http://bioinformatics.univ-reunion.fr/ PBE/.

Cascade PSI-BLAST web server: a remote homology search tool for relating protein domains

Owing to high evolutionary divergence, it is not always possible to identify distantly related protein domains by sequence search techniques. Intermediate sequences possess sequence features of more than one protein and facilitate detection of remotely related proteins. We have demonstrated recently the employment of Cascade PSI-BLAST where we perform PSI-BLAST for many 'generations', initiating searches from new homologues as well. Such a rigorous propagation through generations of PSI-BLAST employs effectively the role of intermediates in detecting distant similarities between proteins. This approach has been tested on a large number of folds and its performance in detecting superfamily level relationships is similar to 35% better than simple PSI-BLAST searches. We present a web server for this search method that permits users to perform Cascade PSI-BLAST searches against the Pfam, SCOP and SwissProt databases. The URL for this server is http://crick.mbu.iisc.ernet.in/similar to CASCADE/CascadeBlast.html.

HARMONY: a server for the assessment of protein structures

Protein structure validation is an important step in computational modeling and structure determination. Stereochemical assessment of protein structures examine internal parameters such as bond lengths and Ramachandran (phi, psi) angles. Gross structure prediction methods such as inverse folding procedure and structure determination especially at low resolution can sometimes give rise to models that are incorrect due to assignment of misfolds or mistracing of electron density maps. Such errors are not reflected as strain in internal parameters. HARMONY is a procedure that examines the compatibility between the sequence and the structure of a protein by assigning scores to individual residues and their amino acid exchange patterns after considering their local environments. Local environments are described by the backbone conformation, solvent accessibility and hydrogen bonding patterns. We are now providing HARMONY through a web server such that users can submit their protein structure files and, if required, the alignment of homologous sequences. Scores are mapped on the structure for subsequent examination that is useful to also recognize regions of possible local errors in protein structures. HARMONY server is located at http://caps.ncbs.res.in/harmony/

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Expression of cyclin E in endomitotic silk-gland cells from mulberry silkworm

The silk glands of mulberry silkworm Bombyx mori are endoreplicating tissues in which the genomic DNA undergoes multiple rounds of replication without mitosis and nuclear division. In the absence of normal mitotic division, the cell cycle essentially alternates between the G1 and S phases. Cyclin E is crucial for the G1/S transition in both mitotic and endoreplicating cycles. We have cloned and characterized cyclin E (cyclin box) from B. mori, which is nearly identical to the Drosophila cyclin E box except for an insertion of 21 amino acids. Two distinct cyclin E transcripts (1.7 and 2.1 kb) were detected in the silk-gland cells of B. mori and in the B. mori-derived embryonic cell line, BmN. Using anti Cyclin E antibodies two protein bands of 52 and 44 kDa were detected in silk glands and BmN cells at Comparable levels. Both BmN- and the silk-gland cells showed the presence of the interacting kinase Cdk2. Transcripts of the mitotic cyclin, cyclin B, were barely detectable in the endoreplicating silk-gland cells and amounted to only 4-7% of that seen in the mitotically dividing BmN cells. The near absence of cyclin B transcripts and the abundant expression of cyclin E in the silk glands correlate well with the alternation of only G1 and S phases without the intervening mitosis in these cells. (C) 2000 Elsevier Science B.V. All rights reserved.

Computer modelling approach to study the modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein

The modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein (ABP) have been studied by energy minimization using the low resolution (2.4 A) X-ray data of the protein. These studies suggest that these sugars preferentially bind in the alpha-form to ABP, unlike L-arabinose where both alpha- and beta-anomers bind almost equally. The best modes of binding of alpha- and beta-anomers of D-galactose and D-fucose differ slightly in the nature of the possible hydrogen bonds with the protein. The residues Arg 151 and Asn 232 of ABP from bidentate hydrogen bonds with both L-arabinose and D-galactose, but not with D-fucose or D-glucose. However in the case of L-arabinose, Arg 151 forms hydrogen bonds with the hydroxyl group at the C-4 atom and the ring oxygen, whereas in case of D-galactose it forms bonds with the hydroxyl groups at the C-4 and C-6 atoms of the pyranose ring. The calculated conformational energies also predict that D-galactose is a better inhibitor than D-fucose and D-glucose, in agreement with kinetic studies. The weak inhibitor D-glucose binds preferentially to one domain of ABP leading to the formation of a weaker complex. Thus these studies provide information about the most probable binding modes of these sugars and also provide a theoretical explanation for the observed differences in their binding affinities.

Preparation and characterization of monoclonal antibodies to nucleocapsid protein N and H glycoprotein of rinderpest virus

Fifteen stable mouse spleen cell myeloma hybrids (hybridomas) producing monoclonal antibodies to rinderpest virus proteins were produced. The specificity of these monoclonal antibodies was established by radioimmunoprecipitation followed by polyacrylamide gel analysis and immunofluorescence. Nine antibodies were specific for the surface glycoprotein H. All the nine clones showed inhibition of haemagglutination by measles virus. The antibodies from two clones (A7D2 and B2F6) neutralise infectious virus. Six clones produce antibodies reacting with the nucleocapsid protein N. Three antigenic sites designated I–III, with sites I and II partially overlapping, were topographically mapped on the H molecule by competitive binding assay. Similarly, two antigenic sites I and II were delineated on the N protein. The monoclonal antibodies were used to study the antigenic relationships of H and N proteins of rinderpest virus, measles virus and canine distemper virus.

Computer simulation of protein—carbohydrate complexes: application to arabinose-binding protein and pea lectin

The CCEM method (Contact Criteria and Energy Minimisation) has been developed and applied to study protein-carbohydrate interactions. The method uses available X-ray data even on the native protein at low resolution (above 2.4 Å) to generate realistic models of a variety of proteins with various ligands.The two examples discussed in this paper are arabinose-binding protein (ABP) and pea lectin. The X-ray crystal structure data reported on ABP-β-l-arabinose complex at 2.8, 2.4 and 1.7 Å resolution differ drastically in predicting the nature of the interactions between the protein and ligand. It is shown that, using the data at 2.4 Å resolution, the CCEM method generates complexes which are as good as the higher (1.7 Å) resolution data. The CCEM method predicts some of the important hydrogen bonds between the ligand and the protein which are missing in the interpretation of the X-ray data at 2.4 Å resolution. The theoretically predicted hydrogen bonds are in good agreement with those reported at 1.7 Å resolution. Pea lectin has been solved only in the native form at 3 Å resolution. Application of the CCEM method also enables us to generate complexes of pea lectin with methyl-α-d-glucopyranoside and methyl-2,3-dimethyl-α-d-glucopyranoside which explain well the available experimental data in solution.

cAMP-regulated Protein Lysine Acetylases in Mycobacteria

Cyclic AMP synthesized by Mycobacterium tuberculosis has been shown to play a role in pathogenesis. However, the high levels of intracellularcAMP found in both pathogenic and nonpathogenic mycobacteria suggest that additional and important biological processes are regulated by characterization of novel cAMP-binding proteins in M. smegmatis and M. tuberculosis (MSMEG_5458 and Rv0998, respectively) that contain a cyclic nucleotide binding domain fused to a domain that shows similarity to the GNAT family of acetyltransferases. We detect protein lysine acetylation in mycobacteria and identify a universal stress protein (USP) as a substrate of MSMEG_5458. Acetylation of a lysine residue in USP is regulated by cAMP, and using a strain deleted for MSMEG_5458, we show that USP is indeed an in vivo substrate for MSMEG_5458. The Rv0998 protein shows a strict cAMP-dependent acetylation of USP, despite a lower affinity for cAMP than MSMEG_5458. Thus, this report not only represents the first demonstration of protein lysine acetylation in mycobacteria but also describes a unique functional interplay between a cyclic nucleotide binding domain and a protein acetyltransferase.

Primary structure of belladonna mottle virus coat protein.

The coat protein of belladonna mottle virus (a tymovirus) was cleaved by trypsin and chymotrypsin, and the peptides were separated by high performance liquid chromatography using a combination of gel permeation, reverse phase, and ion pair chromatography. The peptides were sequenced manually using the 4-N, N-dimethylaminoazobenzene-4'-isothiocyanate/phenyl isothiocyanate double-coupling method. The chymotryptic peptides were aligned by overlapping sequences of tryptic peptides and by homology with another tymovirus, eggplant mosaic virus. The belladonna mottle virus is more closely related to eggplant mosaic virus than to turnip yellow mosaic virus, the type member of this group, as evident from the sequence homologies of 57 and 32%, respectively. The accumulation of basic residues at the amino terminus implicated in RNA-protein interactions in many spherical plant viruses was absent in all the three sequences. Interestingly, the amino-terminal region is the least conserved among the tymoviruses. The longest stretch of conserved sequence between belladonna mottle virus and eggplant mosaic virus was residues 34-44, whereas it was residues 96-102 in the case of belladonna mottle virus and turnip yellow mosaic virus. A tetrapeptide in the region (residues 154-157) was found to be common for all the three sequences. It is possible that these conserved regions (residues 34-44, 96-102, 154-157) are involved in either intersubunit or RNA-protein interactions.

Dynamic changes in mitogen-activated protein kinase (MAPK) activities in the corpus luteum of the bonnet monkey (Macaca radiata) during development, induced luteolysis, and simulated early pregnancy: A role for p38 MAPK in the regulation of luteal function

Changes in MAPK activities were examined in the corpus luteum (CL) during luteolysis and pregnancy, employing GnRH antagonist (Cetrorelix)-induced luteolysis, stages of CL, and hCG treatment to mimic early pregnancy as model systems in the bonnet monkey. We hypothesized that MAPKs could serve to phosphorylate critical phosphoproteins to regulate luteal function. Analysis of several indices for structural (caspase-3 activity and DNA fragmentation) and functional (progesterone and steroidogenic acute regulatory protein expression) changes in the CL revealed that the decreased luteal function observed during Cetrorelix treatment and late luteal phase was associated with increased caspase-3 activity and DNA fragmentation. As expected, human chorionic gonadotropin treatment dramatically increased luteal function, but the indices for structural changes were only partially attenuated. All three MAPKs appeared to be constitutively active in the mid-luteal-phase CL, and activities of ERK-1/2 and p38-MAPK (p38), but not Jun N-terminal kinase (JNK)-1/2, decreased significantly (P < 0.05) within 12 - 24 h after Cetrorelix treatment. During the late luteal phase, in contrast to decreased ERK-1/2 and p38 activities, JNK-1/2 activities increased significantly (P < 0.05). Although human chorionic gonadotropin treatment increased ERK-1/2 and p38 activities, it decreased JNK-1/2 activities. The activation status of p38 was correlated with the phosphorylation status of an upstream activator, MAPK kinase-3/6 and the expression of MAPK activated protein kinase-3, a downstream target. Intraluteal administration of p38 kinase inhibitor (SB203580), but not MAPK kinase-1/2 inhibitor (PD98059), decreased the luteal function. Together, these data suggest an important role for p38 in the regulation of CL function in primates.