Observation of a power-law memory kernel for fluctuations within a single protein molecule

The fluctuation of the distance between a fluorescein-tyrosine pair within a single protein complex was directly monitored in real time by photoinduced electron transfer and found to be a stationary, time-reversible, and non-Markovian Gaussian process. Within the generalized Langevin equation formalism, we experimentally determine the memory kernel K(t), which is proportional to the autocorrelation function of the random fluctuating force. K(t) is a power-law decay, t(-0.51 +/- 0.07) in a broad range of time scales (10(-3)-10 s). Such a long-time memory effect could have implications for protein functions.

Glioblastoma-Specific Protein Interaction Network Identifies PP1A and CSK21 as Connecting Molecules between Cell Cycle-Associated Genes

Glioblastoma (GBM; grade IV astrocytoma) is a very aggressive form of brain cancer with a poor survival and few qualified predictive markers. This study integrates experimentally validated genes that showed specific upregulation in GBM along with their protein-protein interaction information. A system level analysis was used to construct GBM-specific network. Computation of topological parameters of networks showed scale-free pattern and hierarchical organization. From the large network involving 1,447 proteins, we synthesized subnetworks and annotated them with highly enriched biological processes. A careful dissection of the functional modules, important nodes, and their connections identified two novel intermediary molecules CSK21 and protein phosphatase 1 alpha (PP1A) connecting the two subnetworks CDC2-PTEN-TOP2A-CAV1-P53 and CDC2-CAV1-RB-P53-PTEN, respectively. Real-time quantitative reverse transcription-PCR analysis revealed CSK21 to be moderately upregulated and PP1A to be overexpressed by 20-fold in GBM tumor samples. Immunohistochemical staining revealed nuclear expression of PP1A only in GBM samples. Thus, CSK21 and PP1A, whose functions are intimately associated with cell cycle regulation, might play key role in gliomagenesis. Cancer Res; 70(16); 6437-47. (C)2010 AACR.

Interaction of rat testis protein, TP, with nucleosome core particle

Circular dichroism studies have revealed that addition of testis specific protein, TP in vitro, to rat testes nucleosome core particle resulted in a decrease in the compaction of the core particle DNA. This was also corroborated by thermal denaturation analysis. Addition of TP to nucleosome core particle resulted in the conversion of a biphasic transition towards a single phase. However, at the same time there was a 20% reduction in the overall hyperchromicity of core particle DNA at core particle to TP molar ratios of 1:2 and 1:3. These observations along with our earlier report, showing the DNA melting properties of TP, suggest that TP may play an important role in the disassembly process of nucleosome core particle during spermiogenesis.

Functional linkages can reveal protein complexes for structure determination

The study of molecular machines, and protein complexes in general, is a growth area of biology. Is there a computational method for inferring which combinations of proteins in an organism are likely to form a crystallizable complex? We use the Protein Data Bank (PDB) to assess the usefulness of inferred functional protein linkages for this task. We find that of 242 nonredundant prokaryotic protein complexes (complexes excluding structural variants of the same protein) from organisms that are shared between the current PDB and the Prolinks functional linkage database, 44% (107/242) contain proteins that are linked at high-confidence by one or more methods of computed functional linkages. This suggests that computing functional linkages will be useful in defining protein complexes for structural studies. We offer a database of such inferred linkages corresponding to likely protein complexes for some 629,952 pairs of proteins in 154 prokaryotes and archea.

Interaction of water molecules in non-identical protein structures

The interaction of the protein atoms with the surrounding water oxygen atoms has been computed for 392 protein chains from 369 protein structures belonging to 90% non-homologous high resolution (<= 1.5 angstrom) protein Structures with a crystallographic R-factor <= 20%. The percentage composition of the polar atoms is found to be 36.3%. An average of 82.55% of water oxygen atoms are found to be in the primary hydration shell and 15.12% in the secondary hydration shell. The average Percentage of interactions of water oxygen atoms with the polar atoms of the main chain and side chain are 54% and 46%. respectively. The interaction of the acidic residues, aspartate and glutamate, with the water oxygen atoms is more when compared to that of the other residues.

Isolation of a phosphoryl choline-binding protein from the hemolymph of the snail, Achatina fulica

A phosphorylcholine-binding protein from the hemolymph of the snail Achatina fulica was purified to near homogeneity using a Sepharose phenylphosphorylcholine affinity column. The protein bound to the affinity column was eluted with 5 mM phosphorylcholine as a single symmetrical peak. The purified protein (400 Kda) contained 35–40% carbohydrate. On SDS-PAGE the protein separated into two bands of 20 and 24 Kda, and had a pI of 5.9. On immunodiffusion, antiserum to the snail phosphorylcholine binding protein did not cross-react against other phosphorylcholine binding proteins, like rat serum phosphorylcholine-binding protein (PCBP), limulus C-reactive protein (CRP), or human CRP. On pretreatment of the snail hemolymph with this antiserum, the hemagglutination titer of the hemolymph was markedly decreased. The purified snail phosphorylcholine binding protein agglutinated rabbit erythrocytes in the absence of divalent cation (Ca+2) but trace amount of Ca+2 increased its binding. The strongest inhibitor of the agglutination reaction was lactose, followed by melibiose and 2-deoxygalactose. The relationships of the snail phosphorylcholine binding protein to other hemolymph agglutinins and to CRPs are discussed in light of common phylogeny.

Characterization and immune response of a protein produced by a cDNA clone of foot and mouth disease virus, type Asia 1 63/72

A 0.9 kb double stranded cDNA of foot and mouth disease virus (FMDV) Type Asia 1, 63/72 was cloned in an expression vector, pUR222. A protein of 38 kd was produced by the clone which reacted with the antibodies raised against the virus. A 20 kd protein which may be derived from the 38 kd protein contained the antigenic epitopes of the protein VP1 of the virus. Injection of 10-20 micrograms of the partially purified 38 and 20 kd proteins or a lysate of cells containing 240 micrograms of the proteins elicited high titers of FMDV specific antibodies in guinea pigs and cattle respectively. Also, at these concentrations, the proteins protected 5 of 8 guinea pigs and 3 of 8 cattle when challenged with a virulent virus.

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.