Moonlighting Proteins of Lactobacillus crispatus : Extracellular Localization, Cell Wall Anchoring and Interactions with the Host

Moonlighting functions have been described for several proteins previously thought to localize exclusively in the cytoplasm of bacterial or eukaryotic cells. Moonlighting proteins usually perform conserved functions, e. g. in glycolysis or as chaperonins, and their traditional and moonlighting function(s) usually localize to different cell compartments. The most characterized moonlighting proteins in Grampositive bacteria are the glycolytic enzymes enolase and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which function in bacteria-host interactions, e. g. as adhesins or plasminogen receptors. Research on bacterial moonlighting proteins has focused on Gram-positive bacterial pathogens, where many of their functions have been associated with bacterial virulence. In this thesis work I show that also species of the genus Lactobacillus have moonlighting proteins that carry out functions earlier associated with bacterial virulence only. I identified enolase, GAPDH, glutamine synthetase (GS), and glucose-6-phosphate isomerase (GPI) as moonlighting proteins of Lactobacillus crispatus strain ST1 and demonstrated that they are associated with cell surface and easily released from the cell surface into incubation buffer. I also showed that these lactobacillar proteins moonlight either as adhesins with affinity for basement membrane and extracellular matrix proteins or as plasminogen receptors. The mechanisms of surface translocation and anchoring of bacterial moonlighting proteins have remained enigmatic. In this work, the surface localization of enolase, GAPDH, GS and GPI was shown to depend on environmental factors. The members of the genus Lactobacillus are fermentative organisms that lower the ambient pH by producing lactic acid. At acidic pH enolase, GAPDH, GS and GPI were associated with the cell surface, whereas at neutral pH they were released into the buffer. The release did not involve de novo protein synthesis. I showed that purified recombinant His6-enolase, His6-GAPDH, His6-GS and His6-GPI reassociate with cell wall and bind in vitro to lipoteichoic acids at acidic pH. The in-vitro binding of these proteins localizes to cell division septa and cell poles. I also show that the release of moonlighting proteins is enhanced in the presence of cathelicidin LL- 37, which is an antimicrobial peptide and a central part of the innate immunity defence. I found that the LL-37-induced detachment of moonlighting proteins from cell surface is associated with cell wall permeabilization by LL-37. The results in this thesis work are compatible with the hypothesis that the moonlighting proteins of L. crispatus associate to the cell wall via electrostatic or ionic interactions and that they are released into surroundings in stress conditions. Their surface translocation is, at least in part, a result from their release from dead or permeabilized cells and subsequent reassociation onto the cell wall. The results of this thesis show that lactobacillar cells rapidly change their surface architecture in response to environmental factors and that these changes influence bacterial interactions with the host.

The function and regulation of the U11-48K protein in U12-dependent splicing

The removal of noncoding sequences, or introns, from the eukaryotic messenger RNA precursors is catalyzed by a ribonucleoprotein complex known as the spliceosome. In most eukaryotes, two distinct classes of introns exist, each removed by a specific type of spliceosome. The major, U2-type introns account for over 99 % of all introns, and are almost ubiquitous. The minor, U12-type introns are found in most but not all eukaryotes, and reside in conserved locations in a specific set of genes. Due to their slow excision rates, the U12-type introns are expected to be involved in the regulation of the genes containing them by inhibiting the maturation of the messenger RNAs. However, little information is currently available on how the activity of the U12-dependent spliceosome itself is regulated. The levels of many known splicing factors are regulated through unproductive alternative splicing events, which lead to inclusion of premature STOP codons, targeting the transcripts for destruction by the nonsense-mediated decay pathway. These alternative splice sites are typically found in highly conserved sequence elements, which also contain binding sites for factors regulating the activation of the splice sites. Often, the activation is achieved by binding of products of the gene in question, resulting in negative feedback loops. In this study, I show that U11-48K, a protein factor specific to the minor spliceosome, specifically recognizes the U12-type 5' splice site sequence, and is essential for proper function of the minor spliceosome. Furthermore, the expression of U11-48K is regulated through a feedback mechanism, which functions through conserved sequence elements that activate alternative splicing and nonsense-mediated decay. This mechanism is conserved from plants to animals, highlighting both the importance and early origin of this mechanism in regulating splicing factors. I also show that the feedback regulation of U11-48K is counteracted by a component of the major spliceosome, the U1 small nuclear ribonucleoprotein particle, as well as members of the hnRNP F/H protein family. These results thus suggest that the feedback mechanism is finely tuned by multiple factors to achieve precise control of the activity of the U12-dependent spliceosome.

Structure of Sesbania Mosaic Virus at 4·7 Å Resolution and Partial Sequence of the Coat Protein

Sesbania mosaic virus (SMV) is a plant virus infecting Sesbania grandiflora plants in Andhra Pradesh, India. Amino acid sequence of the tryptic peptides of SMV coat protein were determined using a gas phase sequenator. These sequences showed identical amino acids at 69% of the positions when aligned with the corresponding residues of southern bean mosaic virus (SBMV).Crystals diffracting to better than 3 Å resolution were obtained by precipitating the virus with ammonium sulphate. The crystals belonged to rhombohedral space group R3 with α = 291·4 Å and α = 61·9°. Three-dimensional X-ray diffraction data on these crystals were collected to a resolution of 4·7 Å, using a Siemens-Nicolet area detector system. Self-rotation function studies revealed the icosahedral symmetry of the virus particles, as well as their precise orientation in the unit cell. Cross-rotation function and modelling studies with SBMV showed that it is a valid starting model for SMV structure determination. Low resolution phases computed using a polyalanine model of SBMV were subjected to refinement and extension by real-space electron density averaging and solvent flattening. The final electron density map revealed a polypeptide fold similar to SBMV. The single disulphide bridge of SBMV coat protein is retained in SMV. Four icosahedrally independent cation binding sites have been tentatively identified. Three of these sites, related by a quasi threefold axis, are also found in SBMV. The fourth site is situated on the quasi threefold axis. Aspartic acid residues, which replace Ile218 of SBMV from the quasi threefold-related subunits are suitable ligands to the cation at this site

Diversity in Functional Organization of Class I and Class II Biotin Protein Ligase

The cell envelope of Mycobacterium tuberculosis (M. tuberculosis) is composed of a variety of lipids including mycolic acids, sulpholipids, lipoarabinomannans, etc., which impart rigidity crucial for its survival and pathogenesis. Acyl CoA carboxylase (ACC) provides malonyl-CoA and methylmalonyl-CoA, committed precursors for fatty acid and essential for mycolic acid synthesis respectively. Biotin Protein Ligase (BPL/BirA) activates apo-biotin carboxyl carrier protein (BCCP) by biotinylating it to an active holo-BCCP. A minimal peptide (Schatz), an efficient substrate for Escherichia coli BirA, failed to serve as substrate for M. tuberculosis Biotin Protein Ligase (MtBPL). MtBPL specifically biotinylates homologous BCCP domain, MtBCCP87, but not EcBCCP87. This is a unique feature of MtBPL as EcBirA lacks such a stringent substrate specificity. This feature is also reflected in the lack of self/promiscuous biotinylation by MtBPL. The N-terminus/HTH domain of EcBirA has the selfbiotinable lysine residue that is inhibited in the presence of Schatz peptide, a peptide designed to act as a universal acceptor for EcBirA. This suggests that when biotin is limiting, EcBirA preferentially catalyzes, biotinylation of BCCP over selfbiotinylation. R118G mutant of EcBirA showed enhanced self and promiscuous biotinylation but its homologue, R69A MtBPL did not exhibit these properties. The catalytic domain of MtBPL was characterized further by limited proteolysis. Holo-MtBPL is protected from proteolysis by biotinyl-59 AMP, an intermediate of MtBPL catalyzed reaction. In contrast, apo-MtBPL is completely digested by trypsin within 20 min of co-incubation. Substrate selectivity and inability to promote self biotinylation are exquisite features of MtBPL and are a consequence of the unique molecular mechanism of an enzyme adapted for the high turnover of fatty acid biosynthesis.

Cloning and sequence of a cDNA encoding a novel hybrid prolme-rich protein associated with cytokinin-induced haustoria formation in Cuscuta reflexa

A complete cDNA encoding a novel hybrid Pro-rich protein (HyPRP) was identified by differentially screening 3x10(4) recombinant plaques of a Cuscuta reflexa cytokinin-induced haustorial cDNA library constructed in lambda gt10. The nucleotide (nt) sequence consists of: (i) a 424-bp 5'-non coding region having five start codons (ATGs) and three upstream open reading frames (uORFs); (ii) an ORF of 987 bp with coding potential for a 329-amino-acid (aa) protein of M(r), 35203 with a hydrophobic N-terminal region including a stretch of nine consecutive Phe followed by a Pro-rich sequence and a Cys-rich hydrophobic C terminus; and (iii) a 178-bp 3'-UTR (untranslated region). Comparison of the predicted aa sequence with the NBRF and SWISSPROT databases and with a recent report of an embryo-specific protein of maize [Jose-Estanyol et al., Plant Cell 4 (1992) 413-423] showed it to be similar to the class of HyPRPs encoded by genes preferentially expressed in young tomato fruits, maize embryos and in vitro-cultured carrot embryos. Northern analysis revealed an approx. 1.8-kb mRNA of this gene expressed in the subapical region of the C. reflexa vine which exhibited maximum sensitivity to cytokinin in haustorial induction.

Characterization of the DNA-binding domain of beta protein, a component of phage lambda Red-pathway by UV catalyzed cross-linking

beta protein, a key component of Red-pathway of phage lambda is necessary for its growth and general genetic recombination in recombination-deficient mutants of Escherichia coli. To facilitate studies on structure-function relationships, we overexpressed beta protein and purified it to homogeneity. A chemical cross-linking reagent, glutaraldehyde, was used to stabilize the physical association of beta protein in solution. A 67-kDa band, corresponding to homodimer, was identified after separation by SDS-polyacrylamide gel electrophoresis. Stoichiometric measurements indicated a site-size of 1 monomer of beta protein/5 nucleotide residues. Electrophoretic gel mobility shift assays suggested that beta protein formed stable nucleoprotein complexes with 36-mer, but not with 27- or 17-mer DNA. Interestingly, the interaction of beta protein with DNA and the stability of nucleoprotein complexes was dependent on the presence of MgCl2, and the binding was abolished by 250 mM NaCl. The K-d of beta protein binding to 36-mer DNA was on the order of 1.8 x 10(-6) M. Photochemical cross-linking of native beta protein or its fragments, generated by chymotrypsin, to 36-mer DNA was performed to identify its DNA-binding domain. Characterization of the cross-linked peptide disclosed that amino acids required for DNA-binding specificity resided within a 20-kDa peptide at the N-terminal end. These findings provide a basis for further understanding oi the structure and function of beta protein.

Comparative Nucleotide and Amino Acid Sequence Analysis of the Sequence-Specific RNA-Binding Rotavirus Nonstructural Protein NSP3

NSP3, an acidic nonstructural protein, encoded by gene 7 has been implicated as the key player in the assembly of the 11 viral plus-strand RNAs into the early replication intermediates during rotavirus morphogenesis. To date, the sequence or NSP3 from only three animal rotaviruses (SA11, SA114F, and bovine UK) has been determined and that from a human strain has not been reported. To determine the genetic diversity among gene 7 alleles from group A rotaviruses, the nucleotide sequence of the NSP3 gene from 13 strains belonging to nine different G serotypes, from both humans and animals, has been determined. Based on the amino acid sequence identity as well as phylogenetic analysis, NSP3 from group A rotaviruses falls into three evolutionarily related groups, i.e., the SA11 group, the Wa group, and the S2 group. The SA 11/SA114F gene appears to have a distant ancestral origin from that of the others and codes for a polypeptide of 315 amino acids (aa) in length. NSP3 from all other group A rotaviruses is only 313 aa in length because of a 2-amino-acid deletion near the carboxy-terminus, While the SA114F gene has the longest 3' untranslated region (UTR) of 132 nucleotides, that from other strains suffered deletions of varying lengths at two positions downstream of the translational termination codon. In spite of the divergence of the nucleotide (nt) sequence in the protein coding region, a stretch of about 80 nt in the 3' UTR is highly conserved in the NSP3 gene from all the strains. This conserved sequence in the 3' UTR might play an important role in the regulation of expression of the NSP3 gene. (C) 1995 Academic Press, Inc.

Lattice model for rapidly folding protein-like heteropolymers

Protein folding is a relatively fast process considering the astronomical number of conformations in which a protein could find itself. Within the framework of a lattice model, we show that one can design rapidly folding sequences by assigning the strongest attractive couplings to the contacts present in a target native state, Our protein design can be extended to situations with both attractive and repulsive contacts. Frustration is minimized by ensuring that all the native contacts are again strongly attractive. Strikingly, this ensures the inevitability of folding and accelerates the folding process by an order of magnitude, The evolutionary implications of our findings are discussed.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

Essential dynamics and sidechain hydrogen bond cluster studies on eosinophil cationic protein

Eosinophil Cationic Protein (ECP) is a member of RNase A superfamily which carries out the obligatory catalytic role of cleaving RNA. It is involved in a variety of biological functions. Molecular dynamics simulations followed by essential dynamics analysis on this protein are carried out with the goal of gaining insights into the dynamical properties at atomic level. The top essential modes contribute to subspaces and to the transition phase. Further, the sidechain-sidechain/sidechain-mainchain hydrogen bond clusters are analyzed in the top modes, and compared with those of crystal structure. The role of residues identified by these methods is discussed in the context of concerted motion, structure and stability of the protein.

Finding a Box Representation for a Graph in $O(n^2\Delta^2\ln n)$ time

An axis-parallel box in $b$-dimensional space is a Cartesian product $R_1 \times R_2 \times \cdots \times R_b$ where $R_i$ (for $1 \leq i \leq b$) is a closed interval of the form $[a_i, b_i]$ on the real line. For a graph $G$, its boxicity is the minimum dimension $b$, such that $G$ is representable as the intersection graph of (axis-parallel) boxes in $b$-dimensional space. The concept of boxicity finds application in various areas of research like ecology, operation research etc. Chandran, Francis and Sivadasan gave an $O(\Delta n^2 \ln^2 n)$ randomized algorithm to construct a box representation for any graph $G$ on $n$ vertices in $\lceil (\Delta + 2)\ln n \rceil$ dimensions, where $\Delta$ is the maximum degree of the graph. They also came up with a deterministic algorithm that runs in $O(n^4 \Delta )$ time. Here, we present an $O(n^2 \Delta^2 \ln n)$ deterministic algorithm that constructs the box representation for any graph in $\lceil (\Delta + 2)\ln n \rceil$ dimensions.

Functional regulation of PVBV Nuclear Inclusion protein-a protease activity upon interaction with Viral Protein genome-linked and phosphorylation

Regulation of NIa-Pro is crucial for polyprotein processing and hence, for successful infection of potyviruses. We have examined two novel mechanisms that could regulate NIa-Pro activity. Firstly, the influence of VPg domain on the proteolytic activity of NIa-Pro was investigated. It was shown that the turnover number of the protease increases when these two domains interact (as: two-fold; trans: seven-fold) with each other. Secondly, the protease activity of NIa-Pro could also be modulated by phosphorylation at Ser129. A mutation of this residue either to aspartate (phosphorylation-mimic) or alanine (phosphorylation-deficient) drastically reduces the protease activity. Based on these observations and molecular modeling studies, we propose that interaction with VPg as well as phosphorylation of Ser129 could relay a signal through Trp143 present at the protein surface to the active site pocket by subtle conformational changes, thus modulating protease activity of NIa-Pro. (C) 2011 Elsevier Inc. All rights reserved.

Protein Model Discrimination Using Mutational Sensitivity Derived from Deep Sequencing

A major bottleneck in protein structure prediction is the selection of correct models from a pool of decoys. Relative activities of similar to 1,200 individual single-site mutants in a saturation library of the bacterial toxin CcdB were estimated by determining their relative populations using deep sequencing. This phenotypic information was used to define an empirical score for each residue (Rank Score), which correlated with the residue depth, and identify active-site residues. Using these correlations, similar to 98% of correct models of CcdB (RMSD <= 4 angstrom) were identified from a large set of decoys. The model-discrimination methodology was further validated on eleven different monomeric proteins using simulated RankScore values. The methodology is also a rapid, accurate way to obtain relative activities of each mutant in a large pool and derive sequence-structure-function relationships without protein isolation or characterization. It can be applied to any system in which mutational effects can be monitored by a phenotypic readout.