Intra and Inter-Molecular Communications Through Protein Structure Network

Communication within and across proteins is crucial for the biological functioning of proteins. Experiments such as mutational studies on proteins provide important information on the amino acids, which are crucial for their function. However, the protein structures are complex and it is unlikely that the entire responsibility of the function rests on only a few amino acids. A large fraction of the protein is expected to participate in its function at some level or other. Thus, it is relevant to consider the protein structures as a completely connected network and then deduce the properties, which are related to the global network features. In this direction, our laboratory has been engaged in representing the protein structure as a network of non-covalent connections and we have investigated a variety of problems in structural biology, such as the identification of functional and folding clusters, determinants of quaternary association and characterization of the network properties of protein structures. We have also addressed a few important issues related to protein dynamics, such as the process of oligomerization in multimers, mechanism on protein folding, and ligand induced communications (allosteric effect). In this review we highlight some of the investigations which we have carried out in the recent past. A review on protein structure graphs was presented earlier, in which the focus was on the graphs and graph spectral properties and their implementation in the study of protein structure graphs/networks (PSN). In this article, we briefly summarize the relevant parts of the methodology and the focus is on the advancement brought out in the understanding of protein structure-function relationships through structure networks. The investigations of structural/biological problems are divided into two parts, in which the first part deals with the analysis of PSNs based on static structures obtained from x-ray crystallography. The second part highlights the changes in the network, associated with biological functions, which are deduced from the network analysis on the structures obtained from molecular dynamics simulations.

On the relationship of thermodynamic parameters with the buried surface area in protein-ligand complex formation

Prediction of thermodynamic parameters of protein-protein and antigen-antibody complex formation from high resolution structural parameters has recently received much attention, since an understanding of the contributions of different fundamental processes like hydrophobic interactions, hydrogen bonding, salt bridge formation, solvent reorganization etc. to the overall thermodynamic parameters and their relations with the structural parameters would lead to rational drug design. Using the results of the dissolution of hydrocarbons and other model compounds the changes in heat capacity (DeltaCp), enthalpy (DeltaH) and entropy (DeltaS) have been empirically correlated with the polar and apolar surface areas buried during the process of protein folding/unfolding and protein-ligand complex formation. In this regard, the polar and apolar surfaces removed from the solvent in a protein-ligand complex have been calculated from the experimentally observed values of changes in heat capacity (DeltaCp) and enthalpy (DeltaH) for protein-ligand complexes for which accurate thermodynamic and high resolution structural data are available, and the results have been compared with the x-ray crystallographic observations. Analyses of the available results show poor correlation between the thermodynamic and structural parameters. Probable reasons for this discrepancy are mostly related with the reorganization of water accompanying the reaction which is indeed proven by the analyses of the energetics of the binding of the wheat germ agglutinin to oligosaccharides.

Novel CDNF/MANF protein family : Molecular structure, expression and neurotrophic activity

Neurotrophic factors (NTFs) are secreted proteins which promote the survival of neurons, formation and maintenance of neuronal contacts and regulate synaptic plasticity. NTFs are also potential drug candidates for the treatment of neurodegenerative diseases. Parkinson’s disease (PD) is mainly caused by the degeneration of midbrain dopaminergic neurons. Current therapies for PD do not stop the neurodegeneration or repair the affected neurons. Thus, search of novel neurotrophic factors for midbrain dopaminergic neurons, which could also be used as therapeutic proteins, is highly warranted. In the present study, we identified and characterized a novel protein named conserved dopamine neurotrophic factor (CDNF), a homologous protein to mesencephalic astrocyte-derived neurotrophic factor (MANF). Others have shown that MANF supports the survival of embryonic midbrain dopaminergic neurons in vitro, and protects cultured cells against endoplasmic reticulum (ER) stress. CDNF and MANF form a novel evolutionary conserved protein family with characteristic eight conserved cysteine residues in their primary structure. The vertebrates have CDNF and MANF encoding genes, whereas the invertebrates, including Drosophila and Caenorhabditis have a single homologous CDNF/MANF gene. In this study we show that CDNF and MANF are secreted proteins. They are widely expressed in the mammalian brain, including the midbrain and striatum, and in several non-neuronal tissues. We expressed and purified recombinant human CDNF and MANF proteins, and tested the neurotrophic activity of CDNF on midbrain dopaminergic neurons using a 6-hydroxydopamine (6-OHDA) rat model of PD. In this model, a single intrastriatal injection of CDNF protected midbrain dopaminergic neurons and striatal dopaminergic fibers from the 6-OHDA toxicity. Importantly, an intrastriatal injection of CDNF also restored the functional activity of the nigrostriatal dopaminergic system when given after the striatal 6-OHDA lesion. Thus, our study shows that CDNF is a potential novel therapeutic protein for the treatment of PD. In order to elucidate the molecular mechanisms of CDNF and MANF activity, we resolved their crystal structure. CDNF and MANF proteins have two domains; an amino (N)-terminal saposin-like domain and a presumably unfolded carboxy (C)-terminal domain. The saposin-like domain, which is formed by five α-helices and stabilized by three intradomain disulphide bridges, may bind to lipids or membranes. The C-terminal domain contains an internal cysteine bridge in a CXXC motif similar to that of thiol/disulphide oxidoreductases and isomerases, and may thus facilitate protein folding in the ER. Our studies suggest that CDNF and MANF are novel potential therapeutic proteins for the treatment of neurodegenerative diseases. Future studies will reveal the neurotrophic and cytoprotective mechanisms of CDNF and MANF in more detail.

Neurotrophic factors and their receptors

Four GDNF ligands (GDNF, neurturin, artemin and persephin), and mesencephalic astrocyte-derived neurotrophic factor (MANF) and conserved dopamine neurotrophic factor (CDNF) protect midbrain dopaminergic neurons that degenerate in Parkinson's disease. Each GDNF ligand binds a specific coreceptor GDNF family receptor α (GFRα), leading to the formation of a heterotetramer complex, which then interacts with receptor tyrosine kinase RET, the signalling receptor. The present thesis describes the structural and biochemical characterization of the GDNF2-GFRα12 complex and the MANF and CDNF proteins. Previous and current mutation data and comparison between GDNF-GFRα1 and artemin-GFRα3 binding interfaces show that N162GFRα1, I175GFRα1, V230GFRα1, Y120GDNF and L114GDNF are the specificity determinants among different ligand-coreceptor pairs. The structure suggests that sucrose octasulphate, a heparin mimic, interacts with a region R190-K202 within domain 2 of GFRα1. Mutating these residues on the GFRα1 surface, which are not in the GDNF binding region, affected RET phosphorylation, which provides a putative RET binding region in domain 2 and 3 of GFRα1. The structural comparison of the GDNF-GFRα1 and artemin-GFRα3 complexes shows a difference in bend angle between the ligand monomers. This variation in bend angle of the ligand may affect the kinetics of RET phosphorylation. To confirm that the difference is not due to crystallization artefacts, I crystallized the GDNF-GFRα1 complex without SOS in different cell dimensions. The structure of the second GDNF-GFRα1 complex is very similar to the previous one, suggesting that the difference between the artemin-GFRα3 and GDNF-GFRα1 complexes are intrinsic, not due to crystal packing. Finally, MANF and CDNF are bifunctional proteins with extracellular neurotrophic activity and ER resident cytoprotective role. The crystal structures of MANF and CDNF are presented here. Intriguingly, the structures of both the neurotrophic factors do not show structural similarity to any of previously known growth factor superfamilies; instead they are similar to saposins, the lipid-binding proteins. The N-terminal domain of MANF and CDNF contain conserved lysines and arginines on its surface, which may interact with negatively charged head groups of phospholipids, as saposins do. Thus MANF and CDNF may provide neurotrophic activities by interacting with a lipo-receptor. The structure of MANF shows a CXXC motif forming internal disulphide bridge in the natively unfolded C-terminus. This motif is common to reductases and disulphide isomerases. It is thus tempting to speculate that the CXXC motif of MANF and CDNF may be involved in oxidative protein folding, which may explain its cytoprotective role in the ER.

An incipient 310 helix in Piv-Pro-Pro-Ala-NHMe as a model for peptide folding

The molecular mechanism of helix nucleation in peptides and proteins is not yet understood and the question of whether sharp turns in the polypeptide backbone serve as nuclei for protein folding has evoked controversy1,2. A recent study of the conformation of a tetrapeptide containing the stereochemically constrained residue alpha-aminoisobutyric acid, both in solution and the solid state, yielded a structure consisting of two consecutive beta-turns, leading to an incipient 310 helical conformation3,4. This led us to speculate that specific tri- and tetra-peptide sequences may indeed provide a helical twist to the amino-terminal segment of helical regions in proteins and provide a nucleation site for further propagation. The transformation from a 310 helical structure to an alpha-helix should be facile and requires only small changes in the phi and psi conformational angles and a rearrangement of the hydrogen bonding pattern5. If such a mechanism is involved then it should be possible to isolate an incipient 310 helical conformation in a tripeptide amide or tetrapeptide sequence, based purely on the driving force derived from short-range interactions. We have synthesised and studied the model peptide pivaloyl-Pro-Pro-Ala-NHMe (compound I) and provide here spectroscopic evidence for a 310 helical conformation in compound I.

Post-transcriptional Repair of a Split Heat Shock Protein 90 Gene by mRNA trans-Splicing

Heat shock protein 90 participates in diverse biological processes ranging from protein folding, cell cycle, signal transduction and development to evolution in all eukaryotes. It is also critically involved in regulating growth of protozoa such as Dictyostelium discoideum, Leishmania donovani, Plasmodium falciparum, Trypanosoma cruzi, and Trypanosoma evansi. Selective inhibition of Hsp90 has also been explored as an intervention strategy against important human diseases such as cancer, malaria, or trypanosomiasis. Giardia lamblia, a simple protozoan parasite of humans and animals, is an important cause of diarrheal disease with significant morbidity and some mortality in tropical countries. Here we show that the G. lamblia cytosolic hsp90 ( glhsp90) is split in two similar sized fragments located 777 kb apart on the same scaffold. Intrigued by this unique arrangement, which appears to be specific for the Giardiinae, we have investigated the biosynthesis of GlHsp90. We used genome sequencing to confirm the split nature of the giardial hsp90. However, a specific antibody raised against the peptide detected a product with a mass of about 80 kDa, suggesting a post-transcriptional rescue of the genomic defect. We show evidence for the joining of the two independent Hsp90 transcripts in-trans to one long mature mRNA presumably by RNA splicing. The splicing junction carries hallmarks of classical cis-spliced introns, suggesting that the regular cis-splicing machinery may be sufficient for repair of the open reading frame. A complementary 26-nt sequence in the ``intron'' regions adjacent to the splice sites may assist in positioning the two pre-mRNAs for processing. This is the first example of post-transcriptional rescue of a split gene by trans-splicing.

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.

Network properties of protein-decoy structures

Convergence of the vast sequence space of proteins into a highly restricted fold/conformational space suggests a simple yet unique underlying mechanism of protein folding that has been the subject of much debate in the last several decades. One of the major challenges related to the understanding of protein folding or in silico protein structure prediction is the discrimination of non-native structures/decoys from the native structure. Applications of knowledge-based potentials to attain this goal have been extensively reported in the literature. Also, scoring functions based on accessible surface area and amino acid neighbourhood considerations were used in discriminating the decoys from native structures. In this article, we have explored the potential of protein structure network (PSN) parameters to validate the native proteins against a large number of decoy structures generated by diverse methods. We are guided by two principles: (a) the PSNs capture the local properties from a global perspective and (b) inclusion of non-covalent interactions, at all-atom level, including the side-chain atoms, in the network construction accommodates the sequence dependent features. Several network parameters such as the size of the largest cluster, community size, clustering coefficient are evaluated and scored on the basis of the rank of the native structures and the Z-scores. The network analysis of decoy structures highlights the importance of the global properties contributing to the uniqueness of native structures. The analysis also exhibits that the network parameters can be used as metrics to identify the native structures and filter out non-native structures/decoys in a large number of data-sets; thus also has a potential to be used in the protein `structure prediction' problem.

Identification of an exported heat shock protein 70 in Plasmodium falciparum

Host cell remodelling is a hallmark of malaria pathogenesis. It involves protein folding, unfolding and trafficking events and thus participation of chaperones such as Hsp70s and Hsp40s is well speculated. Until recently, only Hsp40s were thought to be the sole representative of the parasite chaperones in the exportome. However, based on the re-annotated Plasmodium falciparum genome sequence, a putative candidate for exported Hsp70 has been reported, which otherwise was known to be a pseudogene. We raised a specific antiserum against a C-terminal peptide uniquely present in PfHsp70-x. Immunoblotting and immunofluorescence-based approaches in combination with sub-cellular fractionation by saponin and streptolysin-O have been taken to determine the expression and localization of PfHsp70-x in infected erythrocyte. The re-annotated sequence of PfHsp70-x reveals it to be a functional protein with an endoplasmic reticulum signal peptide. It gets maximally expressed at the schizont stage of intra-erythrocytic life cycle. Majority of the protein localizes to the parasitophorous vacuole and some of it gets exported to the erythrocyte compartment where it associates with Maurer's clefts. The identification of an exported parasite Hsp70 chaperone presents us with the fact that the parasite has evolved customized chaperones which might be playing crucial roles in aspects of trafficking and host cell remodelling.

Protein Structure and Function: Looking through the Network of Side-Chain Interactions

Network theory has become an excellent method of choice through which biological data are smoothly integrated to gain insights into complex biological problems. Understanding protein structure, folding, and function has been an important problem, which is being extensively investigated by the network approach. Since the sequence uniquely determines the structure, this review focuses on the networks of non-covalently connected amino acid side chains in proteins. Questions in structural biology are addressed within the framework of such a formalism. While general applications are mentioned in this review, challenging problems which have demanded the attention of scientific community for a long time, such as allostery and protein folding, are considered in greater detail. Our aim has been to explore these important problems through the eyes of networks. Various methods of constructing protein structure networks (PSN) are consolidated. They include the methods based on geometry, edges weighted by different schemes, and also bipartite network of protein-nucleic acid complexes. A number of network metrics that elegantly capture the general features as well as specific features related to phenomena, such as allostery and protein model validation, are described. Additionally, an integration of network theory with ensembles of equilibrium structures of a single protein or that of a large number of structures from the data bank has been presented to perceive complex phenomena from network perspective. Finally, we discuss briefly the capabilities, limitations, and the scope for further explorations of protein structure networks.

Cloning and expression of glucose regulated protein 78 (GRP78) in Fenneropenaeus chinensis

GRP78 (78 kDa glucose-regulated protein), also known as BiP (immunoglobulin heavy-chain-binding protein), is an essential regulator of endoplasmic reticulum (ER) homeostasis because of its multiple functions in protein folding, ER calcium binding, and controlling of the activation of transmembrane ER stress sensors. In this report, we cloned the full length cDNA of GRP78 (FcGRP78) from Chinese shrimp Fenneropenaeus chinensis. This cDNA revealed a 2,325 bp with 1,968 bp open reading frame encoding 655 amino acids. This is the first reported GRP78 gene in Crustacea. The deduced amino acid sequence of FcGRP78 shared high identity with previously reported insect GRP78s: 86, 87 and 85% identity with GRP78s of Drosophila melanogaster, Aedes aegypti and Bombyx mori, respectively. Northern blot analysis shows that FcGRP78 is ubiquitously expressed in tissues of shrimp. Heat shock at 35A degrees C significantly enhanced the expression of FcGRP78 at the first hour, reached the maximum at 4 h post heat shock, dropped after that and resumed to the normal level until 48 h of post recovery at 25A degrees C. Additionally, differential expression of FcGRP78 was detected in haemocytes, hepatopancreas and lymphoid organ when shrimp were challenged by white spot syndrome virus (WSSV). We inferred that FcGRP78 may play important roles in chaperoning, protein folding and immune function of shrimp.

Protein Docking by the Underestimation of Free Energy Funnels in the Space of Encounter Complexes

Similarly to protein folding, the association of two proteins is driven by a free energy funnel, determined by favorable interactions in some neighborhood of the native state. We describe a docking method based on stochastic global minimization of funnel-shaped energy functions in the space of rigid body motions (SE(3)) while accounting for flexibility of the interface side chains. The method, called semi-definite programming-based underestimation (SDU), employs a general quadratic function to underestimate a set of local energy minima and uses the resulting underestimator to bias further sampling. While SDU effectively minimizes functions with funnel-shaped basins, its application to docking in the rotational and translational space SE(3) is not straightforward due to the geometry of that space. We introduce a strategy that uses separate independent variables for side-chain optimization, center-to-center distance of the two proteins, and five angular descriptors of the relative orientations of the molecules. The removal of the center-to-center distance turns out to vastly improve the efficiency of the search, because the five-dimensional space now exhibits a well-behaved energy surface suitable for underestimation. This algorithm explores the free energy surface spanned by encounter complexes that correspond to local free energy minima and shows similarity to the model of macromolecular association that proceeds through a series of collisions. Results for standard protein docking benchmarks establish that in this space the free energy landscape is a funnel in a reasonably broad neighborhood of the native state and that the SDU strategy can generate docking predictions with less than 5 � ligand interface Ca root-mean-square deviation while achieving an approximately 20-fold efficiency gain compared to Monte Carlo methods.

Unfolded protein response genes regulated by CED-1 are required for Caenorhabditis elegans innate immunity.

The endoplasmic reticulum stress response, also known as the unfolded protein response (UPR), has been implicated in the normal physiology of immune defense and in several disorders, including diabetes, cancer, and neurodegenerative disease. Here, we show that the apoptotic receptor CED-1 and a network of PQN/ABU proteins involved in a noncanonical UPR response are required for proper defense to pathogen infection in Caenorhabditis elegans. A full-genome microarray analysis indicates that CED-1 functions to activate the expression of pqn/abu genes. We also show that ced-1 and pqn/abu genes are required for the survival of C. elegans exposed to live Salmonella enterica, and that overexpression of pqn/abu genes confers protection against pathogen-mediated killing. The results indicate that unfolded protein response genes, regulated in a CED-1-dependent manner, are involved in the C. elegans immune response to live bacteria.