A Search for Energy Minimized Sequences of Proteins

In this paper, we present numerical evidence that supports the notion of minimization in the sequence space of proteins for a target conformation. We use the conformations of the real proteins in the Protein Data Bank (PDB) and present computationally efficient methods to identify the sequences with minimum energy. We use edge-weighted connectivity graph for ranking the residue sites with reduced amino acid alphabet and then use continuous optimization to obtain the energy-minimizing sequences. Our methods enable the computation of a lower bound as well as a tight upper bound for the energy of a given conformation. We validate our results by using three different inter-residue energy matrices for five proteins from protein data bank (PDB), and by comparing our energy-minimizing sequences with 80 million diverse sequences that are generated based on different considerations in each case. When we submitted some of our chosen energy-minimizing sequences to Basic Local Alignment Search Tool (BLAST), we obtained some sequences from non-redundant protein sequence database that are similar to ours with an E-value of the order of 10(-7). In summary, we conclude that proteins show a trend towards minimizing energy in the sequence space but do not seem to adopt the global energy-minimizing sequence. The reason for this could be either that the existing energy matrices are not able to accurately represent the inter-residue interactions in the context of the protein environment or that Nature does not push the optimization in the sequence space, once it is able to perform the function.

Inhibition of Tyrosine Phosphorylation of Sperm Flagellar Proteins, Outer Dense Fiber Protein-2 and Tektin-2, Is Associated With Impaired Motility During Capacitation of Hamster Spermatozoa

In mammals, acquisition of fertilization competence of spermatozoa is dependent on the phenomenon of sperm capacitation. One of the critical molecular events of sperm capacitation is protein tyrosine phosphorylation. In a previous study, we demonstrated that a specific epidermal growth factor receptor (EGFR)-tyrosine kinase inhibitor, tyrphostin-A47, inhibited hamster sperm capacitation, accompanied by a reduced sperm protein tyrosine phosphorylation. Interestingly, a high percentage of tyrphostin-A47-treated spermatozoa exhibited circular motility, which was associated with a distinct hypo-tyrosine phosphorylation of flagellar proteins, predominantly of Mr 45,000-60,000. In this study, we provide evidence on the localization of capacitation-associated tyrosine-phosphorylated proteins to the nonmembranous, structural components of the sperm flagellum. Consistent with this, we show their ultrastructural localization in the outer dense fiber, axoneme, and fibrous sheath of spermatozoa. Among hypo-tyrosine phosphorylated major proteins of tyrphostin-A47-treated spermatozoa, we identified the 45 kDa protein as outer dense fiber protein-2 and the 51 kDa protein as tektin-2, components of the sperm outer dense fiber and axoneme, respectively. This study shows functional association of hypo-tyrosine-phosphorylation status of outer dense fiber protein-2 and tektin-2 with impaired flagellar bending of spermatozoa, following inhibition of EGFR-tyrosine kinase, thereby showing the critical importance of flagellar protein tyrosine phosphorylation during capacitation and hyperactivation of hamster spermatozoa.

Probiotics modify tight junction proteins in an animal model of non-alcoholic fatty liver disease

Objectives We have investigated the effects of a multi–species probiotic preparation containing a combination of probiotic bacterial genera that included Bifidobacteria, Lactobacilli and a Streptococcus in a mouse model of high fat diet/obesity induced liver steatosis. Methods Three groups of C57B1/6J mice were fed either a standard chow or a high fat diet for 20 weeks, while a third group was fed a high fat diet for 10 weeks and then concomitantly administered probiotics for a further 10 weeks. Serum, liver and large bowel samples were collected for analysis. Results The expression of the tight junction proteins ZO-1 and ZO-2 was reduced (p < 0.05) in high fat diet fed mice compared to chow fed mice. Probiotic supplementation helped to maintain tight ZO-1 and ZO-2 expression compared with the high fat diet group (p < 0.05), but did not restore ZO-1 or ZO-2 expression compared with chow fed mice. Mice fed a high fat diet ± probiotics had significant steatosis development compared to chow fed mice (p < 0.05); steatosis was less severe in the probiotics group compared to the high fat diet group. Hepatic triglycerides concentration was higher in mice fed a high fat diet ± probiotics compared to the chow group (p < 0.05), and was lower in the probiotics group compared to the high fat diet group (p < 0.05). Compared to chow fed mice, serum glucose and cholesterol concentrations, and the activity of alanine transaminase were higher (p < 0.05), whereas serum triglyceride concentration was lower (p < 0.05) in mice fed a high fat diet ± probiotics. Conclusions Supplementation with a multi-species probiotic formulation helped to maintain tight junction proteins ZO-1 and ZO-2, and reduced hepatic triglyceride concentrations compared with a HFD alone.

Molecular characterisation of ABC transporter type FtsE and FtsX proteins of Mycobacterium tuberculosis.

Elicitation of drug resistance and various survival strategies inside host macrophages have been the hallmarks of Mycobacterium tuberculosis as a successful pathogen. ATP Binding Cassette (ABC) transporter type proteins are known to be involved in the efflux of drugs in bacterial and mammalian systems. FtsE, an ABC transporter type protein, in association with the integral membrane protein FtsX, is involved in the assembly of potassium ion transport proteins and probably of cell division proteins as well, both of which being relevant to tubercle bacillus. In this study, we cloned ftsE gene of M. tuberculosis, overexpressed and purified. The recombinant MtFtsE-6xHis protein and the native MtFtsE protein were found localized on the membrane of E. coli and M. tuberculosis cells, respectively. MtFtsE-6xHis protein showed ATP binding in vitro, for which the K42 residue in the Walker A motif was found essential. While MtFtsE-6xHis protein could partially complement growth defect of E. coli ftsE temperature-sensitive strain MFT1181, co-expression of MtFtsE and MtFtsX efficiently complemented the growth defect, indicating that the MtFtsE and MtFtsX proteins might be performing an associated function. MtFtsE and MtFtsX-6xHis proteins were found to exist as a complex on the membrane of E. coli cells co-expressing the two proteins.

Structure of Mycobacterium tuberculosis RuvA, a protein involved in recombination

The process of recombinational repair is crucial for maintaining genomic integrity and generating biological diversity. In association with RuvB and RuvC, RuvA plays a central role in processing and resolving Holliday junctions, which are a critical intermediate in homologous recombination. Here, the cloning, purification and structure determination of the RuvA protein from Mycobacterium tuberculosis (MtRuvA) are reported. Analysis of the structure and comparison with other known RuvA proteins reveal an octameric state with conserved subunit-subunit interaction surfaces, indicating the requirement of octamer formation for biological activity. A detailed analysis of plasticity in the RuvA molecules has led to insights into the invariant and variable regions, thus providing a framework for understanding regional flexibility in various aspects of RuvA function.

Characterization of single-stranded DNA-binding proteins from Mycobacteria. The carboxyl-terminal of domain of SSB is essential for stable association with its cognate RecA protein

Single-stranded DNA-binding proteins (SSB) play an important role in most aspects of DNA metabolism including DNA replication, repair, and recombination. We report here the identification and characterization of SSB proteins of Mycobacterium smegmatis and Mycobacterium tuberculosis. Sequence comparison of M. smegmatis SSB revealed that it is homologous to M. tuberculosis SSB, except for a small spacer connecting the larger amino-terminal domain with the extreme carboxyl-terminal tail. The purified SSB proteins of mycobacteria bound single-stranded DNA with high affinity, and the association and dissociation constants were similar to that of the prototype SSB. The proteolytic signatures of free and bound forms of SSB proteins disclosed that DNA binding was associated with structural changes at the carboxyl-terminal domain. Significantly, SSB proteins from mycobacteria displayed high affinity for cognate RecA, whereas Escherichia coli SSB did not under comparable experimental conditions. Accordingly, SSB and RecA were coimmunoprecipitated from cell lysates, further supporting an interaction between these proteins in vivo. The carboxyl-terminal domain of M. smegmatis SSB, which is not essential for interaction with ssDNA, is the site of binding of its cognate RecA. These studies provide the first evidence for stable association of eubacterial SSB proteins with their cognate RecA, suggesting that these two proteins might function together during DNA repair and/or recombination.

Evolutionary Significance of Self-acylation Property in Acyl Carrier Proteins

Acyl carrier protein is an integral component of many cellular metabolic processes. A number of studies have reported self-acylation behavior in acyl carrier proteins. Although AM exhibit high levels of similarity in their primary and tertiary structures, self-acylation behavior is restricted to only some ACPs that can be classified into two major families based on their function. The first family of ACPs is involved in polyketide biosynthesis, whereas the second family participates in fatty acid synthesis. Facilitated by the growing number of genome sequences available for analyses, large-scale phylogenetic studies were used in these studies to uncover as to how self-acylation behavior of acyl carrier proteins is linked with the evolution of metabolic pathways in organisms. These studies show that self-acylation behavior in acyl carrier proteins was lost during the course of evolution, with certain organisms and organelles viz. plastids, retaining it for specified functions. (C) 2009 IUBMB IUBMB Life, 61(8): 853-859, 2009

Inter-helical Interactions in Membrane Proteins: Analysis Based on the Local Backbone Geometry and the Side Chain Interactions

The availability of a significant number of the Structures of helical membrane proteins has prompted us to investigate the mode of helix-helix packing. In the present study, we have considered a dataset of alpha-helical membrane proteins representing Structures solved from all the known superfamilies. We have described the geometry of all the helical residues in terms of local coordinate axis at the backbone level. Significant inter-helical interactions have been considered as contacts by weighing the number of atom-atom contacts, including all the side-chain atoms. Such a definition of local axis and the contact criterion has allowed us to investigate the inter-helical interaction in a systematic and quantitative manner. We show that a single parameter (designated as alpha), which is derived from the parameters representing the Mutual orientation of local axes, is able to accurately Capture the details of helix-helix interaction. The analysis has been carried Out by dividing the dataset into parallel, anti-parallel, and perpendicular orientation of helices. The study indicates that a specific range of alpha value is preferred for interactions among the anti-parallel helices. Such a preference is also seen among interacting residues of parallel helices, however to a lesser extent. No such preference is seen in the case of perpendicular helices, the contacts that arise mainly due to the interaction Of Surface helices with the end of the trans-membrane helices. The Study Supports the prevailing view that the anti-parallel helices are well packed. However, the interactions between helices of parallel orientation are non-trivial. The packing in alpha-helical membrane proteins, which is systematically and rigorously investigated in this study, may prove to be useful in modeling of helical membrane proteins.

Characterization of LRRTM and NGR gene families: Expression and functions

Some leucine-rich repeat (LRR) -containing membrane proteins are known regulators of neuronal growth and synapse formation. In this work I characterize two gene families encoding neuronal LRR membrane proteins, namely the LRRTM (leucine-rich repeat, transmembrane neuronal) and NGR (Nogo-66 receptor) families. I studied LRRTM and NGR family member's mRNA tissue distribution by RT-PCR and by in situ hybridization. Subcellular localization of LRRTM1 protein was studied in neurons and in non-neuronal cells. I discovered that LRRTM and NGR family mRNAs are predominantly expressed in the nervous system, and that each gene possesses a specific expression pattern. I also established that LRRTM and NGR family mRNAs are expressed by neurons, and not by glial cells. Within neurons, LRRTM1 protein is not transported to the plasma membrane; rather it localizes to endoplasmic reticulum. Nogo-A (RTN4), MAG, and OMgp are myelin-associated proteins that bind to NgR1 to limit axonal regeneration after central nervous system injury. To better understand the functions of NgR2 and NgR3, and to explore the possible redundancy in the signaling of myelin inhibitors of neurite growth, I mapped the interactions between NgR family and the known and candidate NgR1 ligands. I identified high-affinity interactions between RTN2-66, RTN3-66 and NgR1. I also demonstrate that Rtn3 mRNA is expressed in the same glial cell population of mouse spinal cord white matter as Nogo-A mRNA, and thus it could have a role in myelin inhibition of axonal growth. To understand how NgR1 interacts with multiple structurally divergent ligands, I aimed first to map in more detail the nature of Nogo-A:NgR1 interactions, and then to systematically map the binding sites of multiple myelin ligands in NgR1 by using a library of NgR1 expression constructs encoding proteins with one or multiple surface residues mutated to alanine. My analysis of the Nogo-A:NgR1 -interactions revealed a novel interaction site between the proteins, suggesting a trivalent Nogo-A:NgR1-interaction. Our analysis also defined a central binding region on the concave side of NgR1's LRR domain that is required for the binding of all known ligands, and a surrounding region critical for binding MAG and OMgp. To better understand the biological role of LRRTMs, I generated Lrrtm1 and Lrrtm3 knock out mice. I show here that reporter genes expressed from the targeted loci can be used for maping the neuronal connections of Lrrtm1 and Lrrtm3 expressing neurons in finer detail. With regard to LRRTM1's role in humans, we found a strong association between a 70 kb-spanning haplotype in the proposed promoter region of LRRTM1 gene and two possibly related phenotypes: left-handedness and schizophrenia. Interestingly, the responsible haplotype was linked to phenotypic variability only when paternally inherited. In summary, I identified two families of neuronal receptor-like proteins, and mapped their expression and certain protein-protein interactions. The identification of a central binding region in NgR1 shared by multiple ligands may facilitate the design and development of small molecule therapeutics blocking binding of all NgR1 ligands. Additionally, the genetic association data suggests that allelic variation upstream of LRRTM1 may play a role in the development of left-right brain asymmetry in humans. Lrrtm1 and Lrrtm3 knock out mice developed as a part of this study will likely be useful for schizophrenia and Alzheimer s disease research.

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of Delta C-p of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Structural analysis of the roles of influenza A virus membrane-associated proteins in assembly and morphology

The assembly of influenza A virus at the plasma membrane of infected cells leads to release of enveloped virions that are typically round in tissue culture-adapted strains but filamentous in strains isolated from patients. The viral proteins hemagglutinin (HA), neuraminidase (NA), matrix protein 1 (M1), and M2 ion channel all contribute to virus assembly. When expressed individually or in combination in cells, they can all, under certain conditions, mediate release of membrane-enveloped particles, but their relative roles in virus assembly, release, and morphology remain unclear. To investigate these roles, we produced membrane-enveloped particles by plasmid-derived expression of combinations of HA, NA, and M proteins (M1 and M2) or by infection with influenza A virus. We monitored particle release, particle morphology, and plasma membrane morphology by using biochemical methods, electron microscopy, electron tomography, and cryo-electron tomography. Our data suggest that HA, NA, or HANA (HA plus NA) expression leads to particle release through nonspecific induction of membrane curvature. In contrast, coexpression with the M proteins clusters the glycoproteins into filamentous membrane protrusions, which can be released as particles by formation of a constricted neck at the base. HA and NA are preferentially distributed to differently curved membranes within these particles. Both the budding intermediates and the released particles are morphologically similar to those produced during infection with influenza A virus. Together, our data provide new insights into influenza virus assembly and show that the M segment together with either of the glycoproteins is the minimal requirement to assemble and release membrane-enveloped particles that are truly virus-like.

Amino acid interaction preferences in proteins

Understanding the key factors that influence the interaction preferences of amino acids in the folding of proteins have remained a challenge. Here we present a knowledge-based approach for determining the effective interactions between amino acids based on amino acid type, their secondary structure, and the contact based environment that they find themselves in the native state structure as measured by their number of neighbors. We find that the optimal information is approximately encoded in a 60 x 60 matrix describing the 20 types of amino acids in three distinct secondary structures (helix, beta strand, and loop). We carry out a clustering scheme to understand the similarity between these interactions and to elucidate a nonredundant set. We demonstrate that the inferred energy parameters can be used for assessing the fit of a given sequence into a putative native state structure.

How Effective Is The Data On Co-Occurrence Of Domains In Multi-Domain Proteins In Prediction Of Protein-Protein Interactions?

We explore the fuse of information on co-occurrence of domains in multi-domain proteins in predicting protein-protein interactions. The basic premise of our work is the assumption that domains co-occurring in a polypeptide chain undergo either structural or functional interactions among themselves. In this study we use a template dataset of domains in multidomain proteins and predict protein-protein interactions in a target organism. We note that maximum number of correct predictions of interacting protein domain families (158) is made in S. cerevisiae when the dataset of closely related organisms is used as the template followed by the more diverse dataset of bacterial proteins (48) and a dataset of randomly chosen proteins (23). We conclude that use of multi-domain information from organisms closely-related to the target can aid prediction of interacting protein families.