Network approach for capturing ligand-induced subtle global changes in protein structures

Ligand-induced conformational changes in proteins are of immense functional relevance. It is a major challenge to elucidate the network of amino acids that are responsible for the percolation of ligand-induced conformational changes to distal regions in the protein from a global perspective. Functionally important subtle conformational changes (at the level of side-chain noncovalent interactions) upon ligand binding or as a result of environmental variations are also elusive in conventional studies such as those using root-mean-square deviations (r.m.s.d.s). In this article, the network representation of protein structures and their analyses provides an efficient tool to capture these variations (both drastic and subtle) in atomistic detail in a global milieu. A generalized graph theoretical metric, using network parameters such as cliques and/or communities, is used to determine similarities or differences between structures in a rigorous manner. The ligand-induced global rewiring in the protein structures is also quantified in terms of network parameters. Thus, a judicious use of graph theory in the context of protein structures can provide meaningful insights into global structural reorganizations upon perturbation and can also be helpful for rigorous structural comparison. Data sets for the present study include high-resolution crystal structures of serine proteases from the S1A family and are probed to quantify the ligand-induced subtle structural variations.

Unwinding of heterologous DNA by RecA protein during the search for homologous sequences

The search for homologous sequences promoted by RecA protein in vitro involves a presynaptic filament and naked duplex DNA, the multiple contacts of which produce nucleoprotein networks or coaggregates. The single-stranded DNA within the presynaptic filaments, however, is extended to an axial spacing 1.5 times that of B-form DNA. To investigate this paradoxical difference between the spacing of bases in the RecA presynaptic filament versus the target duplex DNA, we explored the effect of heterologous contacts on the conformation of DNA, and vice versa. In the presence of wheat germ topoisomerase I, RecA presynaptic filaments induced a rapid, limited reduction in the linking number of heterologous circular duplex DNA. This limited unwinding of heterologous duplex DNA, termed heterologous unwinding, was detected within 30 seconds and reached a steady state within a few minutes. Presynaptic filaments that were formed in the presence of ATP?S and separated from free RecA protein by gel filtration also generated a ladder of topoisomers upon incubation with relaxed duplex DNA and topoisomerase. The inhibition of heterologous contacts by 60 mImage -NaCl or 5 mImage -ADP resulted in a corresponding decrease in heterologous unwinding. In reciprocal fashion, the stability or number of heterologous contacts with presynaptic filaments was inversely related to the linking number of circular duplex DNA. These observations show that heterologous contacts with the presynaptic filament cause a limited unwinding of the duplex DNA, and conversely that the ability of the DNA to unwind stabilizes transient heterologous contacts.

DNA Condensation by the Rat Spermatidal Protein TP2 Shows GC-Rich Sequence Preference and Is Zinc Dependent

Transition protein-2 (TP2), isolated from rat testes, was recently shown to be a zinc metalloprotein. We have now carried out a detailed analysis of the DNA condensing properties of TP2 with various polynucleotides using circular dichroism spectroscopy. The condensation of the alternating copolymers by TP2 (incubated with 10 mu M ZnSO4), namely, poly(dG-dC). poly(dG-dC) and poly(dA-dT). poly(dA-dT), was severalfold higher than condensation of either of the homoduplexes poly(dG). poly-(dC) and poly(dA). poly(dT) or rat oligonucleosomal DNA. Between the two alternating copolymers, poly(dG-dC). poly(dG-dC) was condensed 3.2-fold more effectively than poly(dA-dT). poly(dA-dT). Preincubation of TP2 with 5 mM EDTA significantly reduced its DNA-condensing property. Interestingly, condensation of the alternating copolymer poly(dI-dC). poly(dI-dC) by TP2 was much less as compared to that of poly(dG-dC). poly(dG-dC). The V8 protease-derived N-terminal fragment (88 aa) condensed poly(dA-dT). poly(dA-dT) to a very small extent but did not have any effect on poly(dG-dC). poly-(dG-dC). The C-terminal fragment (28 aa) was able to condense poly(dA-dT) . poly(dA-dT) more effectively than poly(dG-dC). poly(dG-dC). These results suggest that TP2 in its zinc-coordinated form condenses GC-rich polynucleotides much more effectively than other types of polynucleotides. Neither the N-terminal two-thirds of TP2 which is the zinc-binding domain nor the C-terminal basic domain are as effective as intact TP2 in bringing about condensation of DNA.

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.

Primary structure of a cytotoxin-like basic protein from Naja naja naja (Indian Cobra) venom

The complete amino acid sequence of a cytotoxin-like basic protein (CLBP) from the venom of Naja naja naja (Indian Cobra) was determined by manual degradation using a 4-dimethylaminoazobenzene-4'-isothiocyanate double-coupling method. Peptide fragments obtained by chemical cleavage with cyanogen bromide and enzymic cleavages with trypsin and Staphylococcus aureus proteases for sequence analysis were purified by reversed-phase chromatography. The total number of amino acid residues was 61, with leucine as the C-terminal residue. (C) Munksgaard 1995.

Purification and partial characterization of acyl carrier proteins from developing oil seeds of pisa (Actinodaphne hookeri) and ground nut (Arachis hypogaea)

Acyl carrier proteins (ACP) were purified to homogeneity in the active form from developing seeds of pisa (Actinodaphne hookeri) which synthesizes exclusively trilaurin and from ground nut (Arachis hypogaea) which synthesizes triacylglycerols containing long chain fatty acids. Two major isoforms of ACPs were purified from developing pisa seeds using DEAE-cellulose, Superose-6 FPLC and C-4 reversed phase HPLC chromatographic methods. In contrast, only a single form of ACP was present in ground nut seeds which was purified by anion-exchange and activated thiol-Sepharose 4B affinity chromatography. The two isoforms of ACPs from pisa showed nearly the same specific activity of 6,706 and 7,175 pmol per min per mg protein while ground nut ACP showed a specific activity of 3,893 pmol per min per mg protein when assayed using E. coli acyl-ACP synthetase and [1-C-14]palmitic acid. When compared with E. coli ACP, the purified ACPs from both the seeds showed considerable difference in their mobility in native PAGE, but showed similar mobility in SDS-PAGE under reducing conditions. In the absence of reducing agents formation of dimers was quite prominent. The ACPs from both the seed sources were acid- and heat-stable. The major isoform of pisa seed ACP and the ground nut ACP contain 91 amino acids with M(r) 11,616 and 1,228 respectively. However, there is significant variation in their amino acid composition. A comparision of the amino acid sequence in the N-terminal region of pisa and ground nut seed ACPs showed considerable homology between themselves and with other plant ACPs but not with E. coli ACP.

Primary Sequence That Determines the Functional Overlap between Mitochondrial Heat Shock Protein 70 Ssc1 and Ssc3 of Saccharomyces cerevisiae

The evolutionary diversity of the HSP70 gene family at the genetic level has generated complex structural variations leading to altered functional specificity and mode of regulation in different cellular compartments. By utilizing Saccharomyces cerevisiae as a model system for better understanding the global functional cooperativity between Hsp70 paralogs, we have dissected the differences in functional properties at the biochemical level between mitochondrial heat shock protein 70 (mtHsp70) Ssc1 and an uncharacterized Ssc3 paralog. Based on the evolutionary origin of Ssc3 and a high degree of sequence homology with Ssc1, it has been proposed that both have a close functional overlap in the mitochondrial matrix. Surprisingly, our results demonstrate that there is no functional cross-talk between Ssc1 and Ssc3 paralogs. The lack of in vivo functional overlap is due to altered conformation and significant lower stability associated with Ssc3. The substrate-binding domain of Ssc3 showed poor affinity toward mitochondrial client proteins and Tim44 due to the open conformation in ADP-bound state. In addition to that, the nucleotide-binding domain of Ssc3 showed an altered regulation by the Mge1 co-chaperone due to a high degree of conformational plasticity, which strongly promotes aggregation. Besides, Ssc3 possesses a dysfunctional inter-domain interface thus rendering it unable to perform functions similar to generic Hsp70s. Moreover, we have identified the critical amino acid sequence of Ssc1 and Ssc3 that can ``make or break'' mtHsp70 chaperone function. Together, our analysis provides the first evidence to show that the nucleotide-binding domain of mtHsp70s plays a critical role in determining the functional specificity among paralogs and orthologs across kingdoms.

RPS: Repeats in Protein Sequences

Repeats are two or more contiguous segments of amino acid residues that are believed to have arisen as a result of intragenic duplication, recombination and mutation events. These repeats can be utilized for protein structure prediction and can provide insights into the protein evolution and phylogenetic relationship. Therefore, to aid structural biologists and phylogeneticists in their research, a computing resource (a web server and a database), Repeats in Protein Sequences (RPS), has been created. Using RPS, users can obtain useful information regarding identical, similar and distant repeats (of varying lengths) in protein sequences. In addition, users can check the frequency of occurrence of the repeats in sequence databases such as the Genome Database, PIR and SWISS-PROT and among the protein sequences available in the Protein Data Bank archive. Furthermore, users can view the three-dimensional structure of the repeats using the Java visualization plug-in Jmol. The proposed computing resource can be accessed over the World Wide Web at http://bioserver1.physics.iisc.ernet.in/rps/.

Angelman Syndrome Protein UBE3A Interacts with Primary Microcephaly Protein ASPM, Localizes to Centrosomes and Regulates Chromosome Segregation

Many proteins associated with the phenotype microcephaly have been localized to the centrosome or linked to it functionally. All the seven autosomal recessive primary microcephaly (MCPH) proteins localize at the centrosome. Microcephalic osteodysplastic primordial dwarfism type II protein PCNT and Seckel syndrome (also characterized by severe microcephaly) protein ATR are also centrosomal proteins. All of the above findings show the importance of centrosomal proteins as the key players in neurogenesis and brain development. However, the exact mechanism as to how the loss-of-function of these proteins leads to microcephaly remains to be elucidated. To gain insight into the function of the most commonly mutated MCPH gene ASPM, we used the yeast two-hybrid technique to screen a human fetal brain cDNA library with an ASPM bait. The analysis identified Angelman syndrome gene product UBE3A as an ASPM interactor. Like ASPM, UBE3A also localizes to the centrosome. The identification of UBE3A as an ASPM interactor is not surprising as more than 80% of Angelman syndrome patients have microcephaly. However, unlike in MCPH, microcephaly is postnatal in Angelman syndrome patients. Our results show that UBE3A is a cell cycle regulated protein and its level peaks in mitosis. The shRNA knockdown of UBE3A in HEK293 cells led to many mitotic abnormalities including chromosome missegregation, abnormal cytokinesis and apoptosis. Thus our study links Angelman syndrome protein UBE3A to ASPM, centrosome and mitosis for the first time. We suggest that a defective chromosome segregation mechanism is responsible for the development of microcephaly in Angelman syndrome.

Angelica archengelica extract induced perturbation of rat skin and tight junctional protein (ZO-1) of HaCaT cells

Background and purpose of the study: Herbal enhancers compared to the synthetic ones have shown less toxis effects. Coumarins have been shown at concentrations inhibiting phospoliphase C-Y (Phc-Y) are able to enhance tight junction (TJ) permeability due to hyperpoalation of Zonolous Occludense-1 (ZO-1) proteins. The purpose of this study was to evaluate the influence of ethanolic extract of Angelica archengelica (AA-E) which contain coumarin on permeation of repaglinide across rat epidermis and on the tight junction plaque protein ZO-1 in HaCaT cells. Methods: Transepidermal water loss (TEWL) from the rat skin treated with different concentrations of AA-E was assessed by Tewameter. Scanning and Transmission Electron Microscopy (TEM) on were performed on AA-E treated rat skin portions. The possibility of AA-E influence on the architecture of tight junctions by adverse effect on the cytoplasmic ZO-1 in HaCaT cells was investigated. Finally, the systemic delivery of repaglinide from the optimized transdermal formulation was investigated in rats. Results: The permeation of repaglinide across excised rat epidermis was 7-fold higher in the presence of AA-E (5% w/v) as compared to propylene glycol:ethanol (7:3) mixture. The extract was found to perturb the lipid microconstituents in both excised and viable rat skin, although, the effect was less intense in the later. The enhanced permeation of repaglinide across rat epidermis excised after treatment with AA-E (5% w/v) for different periods was in concordance with the high TEWL values of similarly treated viable rat skin. Further, the observed increase in intercellular space, disordering of lipid structure and corneocyte detachment indicated considerable effect on the ultrastructure of rat epidermis. Treatment of HaCaT cell line with AA-E (0.16% w/v) for 6 hrs influenced ZO-1 as evidenced by reduced immunofluorescence of anti-TJP1 (ZO-1) antibody in Confocal Laser Scanning Microscopy studies (CLSM) studies. The plasma concentration of repaglinide from transdermal formulation was maintained higher and for longer time as compared to oral administration of repaglinide. Major conclusion: Results suggest the overwhelming influence of Angelica archengelica in enhancing the percutaneous permeation of repaglinide to be mediated through perturbation of skin lipids and tight junction protein (ZO-1).

Structural studies of protein unfolding

Methods for macromolecular structure determination (NMR and crystallography) are now being used to get structural information on partially folded and unfolded states of proteins. These techniques, in combination with proton hydrogen exchange studies are powerful tools to extract information on non-native states of proteins. This review discusses progress In this area of protein folding.

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.

Water-mediated ionic interactions in protein structures

It is well known that water molecules play an indispensable role in the structure and function of biological macromolecules. The water-mediated ionic interactions between the charged residues provide stability and plasticity and in turn address the function of the protein structures. Thus, this study specifically addresses the number of possible water-mediated ionic interactions, their occurrence, distribution and nature found in 90% non-redundant protein chains. Further, it provides a statistical report of different charged residue pairs that are mediated by surface or buried water molecules to form the interactions. Also, it discusses its contributions in stabilizing various secondary structural elements of the protein. Thus, the present study shows the ubiquitous nature of the interactions that imparts plasticity and flexibility to a protein molecule.

Stochastic dynamics modeling of the protein sequence length distribution in genomes: implications for microbial evolution

In this paper, we report an analysis of the protein sequence length distribution for 13 bacteria, four archaea and one eukaryote whose genomes have been completely sequenced, The frequency distribution of protein sequence length for all the 18 organisms are remarkably similar, independent of genome size and can be described in terms of a lognormal probability distribution function. A simple stochastic model based on multiplicative processes has been proposed to explain the sequence length distribution. The stochastic model supports the random-origin hypothesis of protein sequences in genomes. Distributions of large proteins deviate from the overall lognormal behavior. Their cumulative distribution follows a power-law analogous to Pareto's law used to describe the income distribution of the wealthy. The protein sequence length distribution in genomes of organisms has important implications for microbial evolution and applications. (C) 1999 Elsevier Science B.V. All rights reserved.