Protein-protein interactions are important in carrying out many biological processes and functions. These interactions may be either permanent or of temporary nature. Several studies have employed tools like solvent accessibility and graph theory to identify these interactions, but still more studies need to be performed to quantify and validate them. Although we now have many databases available with predicted and experimental results on protein-protein interactions, we still do not have many databases which focus on providing structural details of the interacting complexes, their oligomerisation state and homologues. In this work, protein-protein interactions have been thoroughly investigated within the structural regime and quantified for their strength using calculated pseudoenergies. The PPCheck server, an in-house webserver, has been used for calculating the pseudoenergies like van der Waals, hydrogen bonds and electrostatic energy based on distances between atoms of amino acids from two interacting proteins. PPCheck can be visited at . Based on statistical data, as obtained by studying established protein-protein interacting complexes from earlier studies, we came to a conclusion that an average protein-protein interface consisted of about 51 to 150 amino acid residues and the generalized energy per residue ranged from -2 kJ mol(-1) to -6 kJ mol(-1). We found that some of the proteins have an exceptionally higher number of amino acids at the interface and it was purely because of their elaborate interface or extended topology i.e. some of their secondary structure regions or loops were either inter-mixing or running parallel to one another or they were taking part in domain swapping. Residue networks were prepared for all the amino acids of the interacting proteins involved in different types of interactions (like van der Waals, hydrogen-bonding, electrostatic or intramolecular interactions) and were analysed between the query domain-interacting partner pair and its remote homologue-interacting partner pair. We found that, in exceptional cases, homologous proteins belonging to the same superfamily, but with remote sequence similarity, can share similar interfaces.

BACKGROUNDModeling protein–protein interactions (PPIs) using docking algorithms is useful for understanding biomolecular interactions and mechanisms. Typically, a docking algorithm generates a large number of docking poses, and it is often challenging to select the best native-like pose. A further challenge is to recognize key residues, termed as hotspots, at protein–protein interfaces, which contribute more in stabilizing a protein–protein interface.RESULTSWe had earlier developed a computer algorithm, called PPCheck, which ascribes pseudoenergies to measure the strength of PPIs. Native-like poses could be successfully identified in 27 out of 30 test cases, when applied on a separate set of decoys that were generated using FRODOCK. PPCheck, along with conservation and accessibility scores, was able to differentiate ‘native-like and non-native-like poses from 1883 decoys of Critical Assessment of Prediction of Interactions (CAPRI) targets with an accuracy of 60%. PPCheck was trained on a 10-fold mixed dataset and tested on a 10-fold mixed test set for hotspot prediction. We obtain an accuracy of 72%, which is in par with other methods, and a sensitivity of 59%, which is better than most existing methods available for hotspot prediction that uses similar datasets. Other relevant tests suggest that PPCheck can also be reliably used to identify conserved residues in a protein and to perform computational alanine scanning.CONCLUSIONSPPCheck webserver can be successfully used to differentiate native-like and non-native-like docking poses, as generated by docking algorithms. The webserver can also be a convenient platform for calculating residue conservation, for performing computational alanine scanning, and for predicting protein–protein interface hotspots. While PPCheck can differentiate the generated decoys into native-like and non-native-like decoys with a fairly good accuracy, the results improve dramatically when features like conservation and accessibility are included. The method can be successfully used in ranking/scoring the decoys, as obtained from docking algorithms.

BackgroundKrishna Tulsi, a member of Lamiaceae family, is a herb well known for its spiritual, religious and medicinal importance in India. The common name of this plant is ‘Tulsi’ (or ‘Tulasi’ or ‘Thulasi’) and is considered sacred by Hindus. We present the draft genome of Ocimum tenuiflurum L (subtype Krishna Tulsi) in this report. The paired-end and mate-pair sequence libraries were generated for the whole genome sequenced with the Illumina Hiseq 1000, resulting in an assembled genome of 374 Mb, with a genome coverage of 61 % (612 Mb estimated genome size). We have also studied transcriptomes (RNA-Seq) of two subtypes of O. tenuiflorum, Krishna and Rama Tulsi and report the relative expression of genes in both the varieties.ResultsThe pathways leading to the production of medicinally-important specialized metabolites have been studied in detail, in relation to similar pathways in Arabidopsis thaliana and other plants. Expression levels of anthocyanin biosynthesis-related genes in leaf samples of Krishna Tulsi were observed to be relatively high, explaining the purple colouration of Krishna Tulsi leaves. The expression of six important genes identified from genome data were validated by performing q-RT-PCR in different tissues of five different species, which shows the high extent of urosolic acid-producing genes in young leaves of the Rama subtype. In addition, the presence of eugenol and ursolic acid, implied as potential drugs in the cure of many diseases including cancer was confirmed using mass spectrometry.ConclusionsThe availability of the whole genome of O.tenuiflorum and our sequence analysis suggests that small amino acid changes at the functional sites of genes involved in metabolite synthesis pathways confer special medicinal properties to this herb.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-015-0562-x) contains supplementary material, which is available to authorized users.

No abstract

Several catalysis, cellular regulation, immune function, cell wall assembly, transport, signaling and inhibition occur through Protein- Protein Interactions (PPI). This is possible with the formation of specific yet stable protein-protein interfaces. Therefore, it is of interest to understand its molecular principles using structural data in relation to known function. Several interface features have been documented using known X-ray structures of protein complexes since 1975. This has improved our understanding of the interface using structural features such as interface area, binding energy, hydrophobicity, relative hydrophobicity, salt bridges and hydrogen bonds. The strength of binding between two proteins is dependent on interface size (number of residues at the interface) and thus its corresponding interface area. It is known that large interfaces have high binding energy (sum of (van der Waals) vdW, H-bonds, electrostatics). However, the selective role played by each of these energy components and more especially that of vdW is not explicitly known. Therefore, it is important to document their individual role in known protein-protein structural complexes. It is of interest to relate interface size with vdW, H-bonds and electrostatic interactions at the interfaces of protein structural complexes with known function using statistical and multiple linear regression analysis methods to identify the prominent force. We used the manually curated non-redundant dataset of 278 hetero-dimeric protein structural complexes grouped using known functions by Sowmya et al. (2015) to gain additional insight to this phenomenon using a robust inter-atomic non-covalent interaction analyzing tool PPCheck (Anshul and Sowdhamini, 2015). This dataset consists of obligatory (enzymes, regulator, biological assembly), immune and nonobligatory (enzyme and regulator inhibitors) complexes. Results show that the total binding energy is more for large interfaces. However, this is not true for its individual energy factors. Analysis shows that vdW energies contribute to about 75% ± 11% on average among all complexes and it also increases with interface size (r2 ranging from 0.67 to 0.89 with p<0.01) at 95% confidence limit irrespective of molecular function. Thus, vdW is both dominant and proportional at the interface independent of molecular function. Nevertheless, H bond energy contributes to 15% ± 6.5% on average in these complexes. It also moderately increases with interface size (r2 ranging from 0.43 to 0.61 with p<0.01) only among obligatory and immune complexes. Moreover, there is about 11.3% ± 8.7% contribution by electrostatic energy. It increases with interface size specifically among non-obligatory regulator-inhibitors (r2 = 0.44). It is implied that both H-bonds and electrostatics are neither dominant nor proportional at the interface. Nonetheless, their presence cannot be ignored in binding. Therefore, H-bonds and (or) electrostatic energy having specific role for improved stability in complexes is implied. Thus, vdW is commo...

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.

No abstract

This protocol describes a stepwise process to identify proteins of interest from a query proteome derived from NGS data. We implemented this protocol on Moringa oleifera transcriptome to identify proteins involved in secondary metabolite and vitamin biosynthesis and ion transport. This knowledge-driven protocol identifies proteins using an integrated approach involving sensitive sequence search and evolutionary relationships. We make use of functionally important residues (FIR) specific for the query protein family identified through its homologous sequences and literature. We screen protein hits based on the clustering with true homologues through phylogenetic tree reconstruction complemented with the FIR mapping. The protocol was validated for the protein hits through qRT-PCR and transcriptome quantification. Our protocol demonstrated a higher specificity as compared to other methods, particularly in distinguishing cross-family hits. This protocol was effective in transcriptome data analysis of M. oleifera as described in Pasha et al. Knowledge-driven protocol to identify secondary metabolite synthesizing protein in a highly specific manner. Use of functionally important residues for screening of true hits. Beneficial for metabolite pathway reconstruction in any (species, metagenomics) NGS data.