Backbones of Folded Proteins Reveal Novel Invariant Amino Acid Neighborhoods

Secretory Conservation in Insulin Producing Cells: Is There a System-Level Law of Mass Action in Biology?

Altered secretion of insulin from pancreatic β-cells can manifest into disorders. For example, a lack of endogenously produced and/or secreted insulin results in Type 1 diabetes (and other associated subtypes). Pancreatic β-cells are the endocrine secretory cells that promote insulin secretion in response to glucose stimulation. Secretion in response to extracellular triggers is an interplay among various signaling pathways, transcription factors, and molecular mechanisms. The Mouse Insulinoma 6 (MIN6) cell line serves as a model system for gaining mechanistic insights into pancreatic β-cell functions. It is obvious that higher glucose consumption and increased insulin secretion are correlated. However, it has been reported that intracellular ATP levels remain ∼ constant beyond the extracellular glucose (EG) concentration of 10 mM. Therefore, any cause-effect relationship between glucose consumption (GC) and enhanced insulin secretion (eIS) remains unclear. We also found that total cellular protein, as well as total protein content in the culture "supernatant," remains constant regardless of varying EG concentrations. This indicated that eIS may be at the cost of (a) intracellular synthesis of other proteins and (b) secretion of other secretory proteins, or both (a) and (b), somehow coupled with GC by cells. To gain insights into the above, we carried out a transcriptome study of MIN6 cells exposed to hypoglycemic (HoG = 2.8 mM EG) and hyperglycemic (HyG = 25 mM EG) conditions. Expression of transcripts was analyzed in terms of Fragments Per Kilobase of transcript per Million mapped reads and Transcripts Per Million (FPKM and TPM) as well as values obtained by normalizing w.r.t. "∑(FPKM)" and "∑(TPM)." We report that HyG extracellular conditions lead to an ∼2-fold increase in insulin secretion compared to HoG measured by the enzyme-linked immunosorbent assay (ELISA) and transcripts of secreted proteins as well as their isoforms decreased in HyG conditions compared to HoG. Our results show for the first time that eIS in HyG conditions is at the cost of reduced transcription of other secreted proteins and is coupled with higher GC. The higher GC at increased extracellular glucose also indicates a yet undiscovered role of glucose molecules enhancing insulin secretion, since ATP levels resulting from glucose metabolism have been reported to be constant above an EG concentration of 10 mM. While extrapolation of our results to clinical implications is ambitious at best, this work reports novel cellular level aspects that seem relevant in some clinical observations pertaining to Type 1 diabetes. In addition, the conservatory nature of cellular secretions in insulin-secreting cells, discovered here, may be a general feature in cell biology.

Pseudopeptosomes: non-lipidated vesicular assemblies from bispidine-appended pseudopeptides

We report a novel molecular topology-based approach for creating reproducible vesicular assemblies in different solvent environments (including aqueous) using specifically designed pseudopeptides. Deviating from the classical "polar head group and hydrophobic tail" model of amphiphiles, we showed (reversible) self-assembly of synthesized pseudopeptides into vesicles. Naming these new type/class of vesicles "pseudopetosomes", we characterized them by high-resolution microscopy (scanning electron, transmission electron, atomic force, epifluorescence and confocal) along with dynamic light scattering. While accounting for hydropathy index of the constituent amino acids (side chains) of pseudopeptides, we probed molecular interactions, resulting in assembly of pseudopeptosomes by spectroscopy (fourier-transform infrared and fluorescence). Molecular characterization by X-ray crystallography and circular dichroism revealed "tryptophan (Trp)-Zip" arrangements and/or hydrogen-bonded one-dimensional assembly depending on specific pseudopeptides and solvent environments. Our data indicated that pseudopeptosomes are formed in solutions by self-assembly of bispidine pseudopeptides (of Trp, leucine and alanine amino-acid constituents) into sheets that transform into vesicular structures. Thus, we showed that assembly of pseudopeptosomes utilizes the full spectrum of all four weak interactions essential in biological systems. Our findings have direct implications in chemical and synthetic biology, but may also provide a new avenue of investigations on origins of life via pseudopeptosome-like assemblies. We also showed that these designer peptides can act as carriers for cellular transport.

Aspects of Biological Replication and Evolution Independent of the Central Dogma: Insights from Protein-Free Vesicular Transformations and Protein-Mediated Membrane Remodeling

Biological membrane remodeling is central to living systems. In spite of serving as “containers” of whole-living systems and functioning as dynamic compartments within living systems, biological membranes still find a “blue collar” treatment compared to the “white collar” nucleic acids and proteins in biology. This may be attributable to the fact that scientific literature on biological membrane remodeling is only 50 years old compared to ~ 150 years of literature on proteins and a little less than 100 years on nucleic acids. However, recently, evidence for symbiotic origins of eukaryotic cells from data only on biological membranes was reported. This, coupled with appreciation of reproducible amphiphilic self-assemblies in aqueous environments (mimicking replication), has already initiated discussions on origins of life beyond nucleic acids and proteins. This work presents a comprehensive compilation and meta-analyses of data on self-assembly and vesicular transformations in biological membranes—starting from model membranes to establishment of Influenza Hemagglutinin-mediated membrane fusion as a prototypical remodeling system to a thorough comparison between enveloped mammalian viruses and cellular vesicles. We show that viral membrane fusion proteins, in addition to obeying “stoichiometry-driven protein folding”, have tighter compositional constraints on their amino acid occurrences than general-structured proteins, regardless of type/class. From the perspective of vesicular assemblies and biological membrane remodeling (with and without proteins) we find that cellular vesicles are quite different from viruses. Finally, we propose that in addition to pre-existing thermodynamic frameworks, kinetic considerations in de novo formation of metastable membrane structures with available “third-party” constituents (including proteins) were not only crucial for origins of life but also continue to offer morphological replication and/or functional mechanisms in modern life forms, independent of the central dogma.

O-H stretching frequency red shifts do not correlate with the dissociation energies in the dimethylether and dimethylsulfide complexes of phenol derivatives

In this perspective, we present a comprehensive report on the spectroscopic and computational investigations of the hydrogen bonded (H-bonded) complexes of Me2O and Me2S with seven para-substituted H-bond donor phenols. The salient finding was that although the dissociation energies, D0, of the Me2O complexes were consistently higher than those of the analogous Me2S complexes, the red-shifts in phenolic O-H frequencies, Δν(O-H), showed the exactly opposite trend. This is in contravention of the general perception that the red shift in the X-H stretching frequency in the X-HY hydrogen bonded complexes is a reliable indicator of H-bond strength (D0), a concept popularly known as the Badger-Bauer rule. This is also in contrast to the trend reported for the H-bonded complexes of H2S/H2O with several para substituted phenols of different pKa values wherein the oxygen centered hydrogen bonded (OCHB) complexes consistently showed higher Δν(O-H) and D0 compared to those of the analogous sulfur centered hydrogen bonded (SCHB) complexes. Our effort was to understand these intriguing observations based on the spectroscopic investigations of 1 : 1 complexes in combination with a variety of high level quantum chemical calculations. Ab initio calculations at the MP2 level and the DFT calculations using various dispersion corrected density functionals (including DFT-D3) were performed on counterpoise corrected surfaces to compute the dissociation energy, D0, of the H-bonded complexes. The importance of anharmonic frequency computations is underscored as they were able to correctly reproduce the observed trend in the relative OH frequency shifts unlike the harmonic frequency computations. We have attempted to find a unified correlation that would globally fit the observed red shifts in the O-H frequency with the H-bonding strength for the four bases, namely, H2S, H2O, Me2O, and Me2S, in this set of H-bond donors. It was found that the proton affinity normalized Δν(O-H) values scale very well with the H-bond strength.

Unique and exclusive peptide signatures directly identify intrinsically disordered proteins from sequences without structural information

Intrinsically disordered proteins are now widely accepted to play crucial roles in biological functions. Identification of signatures of intrinsic disorder is one of the key steps towards building a proper repertoire for their occurrence in proteomes. In this work, systematic computational synthesis of a library of all possible (3368400) dipeptides, tripeptides, tetrapeptides and pentapeptides using the natural 20 amino acids allowed us to identify 36 unique tetrapeptides present exclusively in intrinsically disordered proteins and absent in the complete primary sequence space of naturally occurring structured proteins. Further, out of more than 530000 known naturally occurring primary sequences without any structural information, 1349 sequences contain the above identified unique signatures of intrinsic disorder. These sequences, having cellular functions varying from housekeeping to metabolic to transport, more than double the number of the currently known intrinsically disordered proteins. On similar lines, we report that 26577 pentapeptide signatures exclusive to intrinsically disordered proteins, and absent in naturally occurring structured proteins, identify ∼50% of more than half-a-million curated protein sequences without structural information to be intrinsically disordered. The results reported are a major leap forward in exploring functional manifestations of intrinsically disordered proteins.Communicated by Ramaswamy H. Sarma.

Structural disorder originates beyond narrow stoichiometric margins of amino acids in naturally occurring folded proteins

Rigorous analyses of Euclidean distances between non-peptide bonded residues in structures of several thousand naturally occurring folded proteins yielded a surprising "margin of life" for percentage occurrence of individual amino acids in naturally occurring folded proteins. On one hand, the concept of "margin of life", referring to lower than expected variances in average stoichiometric occurrences of individual amino acids in folded proteins, remains unchallenged since its discovery a decade ago. On the other hand, within this past decade there has been a strong emergence of a gradual paradigm shift in biology, from sequence-structure-function in proteins to sequence-disorder-function, fuelled by discoveries on functional implications of intrinsically disordered proteins (primary sequences that do not form stable structures). Thus the applicability of "margin of life" to peptide-bonded residues in all known natural proteins, adopting stable structures vis-à-vis intrinsically disordered needs to be explored. Therefore in this work, we analyze compositions of the complete naturally occurring primary sequence space (over 560000 sequences) after dividing it into mutually exclusive subsets of structured and intrinsically disordered proteins along with a subset without any structural information. While finding that occurrence of different peptides (up to pentapeptides) is a direct consequence of the relative occurrences of their constituting residues in folded proteins, we report that structural disorder in natural proteins originates beyond the narrow stoichiometric margins of amino acids found in structured proteins.Communicated by Ramaswamy H. Sarma.

Comparing protein structures and inferring functions with a novel three-dimensional Yau-Hausdorff method

Structures and functions of proteins play various essential roles in biological processes. The functions of newly discovered proteins can be predicted by comparing their structures with that of known-functional proteins. Many approaches have been proposed for measuring the protein structure similarity, such as the template-modeling (TM)-score method, GRaphlet (GR)-Align method as well as the commonly used root-mean-square deviation (RMSD) measures. However, the alignment comparisons between the similarity of protein structure cost much time on large dataset, and the accuracy still have room to improve. In this study, we introduce a new three-dimensional (3D) Yau-Hausdorff distance between any two 3D objects. The (3D) Yau-Hausdorff distance can be used in particular to measure the similarity/dissimilarity of two proteins of any size and does not need aligning and superimposing two structures. We apply structural similarity to study function similarity and perform phylogenetic analysis on several datasets. The results show that (3D) Yau-Hausdorff distance could serve as a more precise and effective method to discover biological relationships between proteins than other methods on structure comparison. Communicated by Ramaswamy H. Sarma.

In silico approaches illustrate the evolutionary pattern and protein-small molecule interactions of quinolone synthase from Aegle marmelos Correa

Quinolone synthase from Aegle marmelos (AmQNS) is a Rutacean-specific plant type III polyketide synthase that synthesizes quinolone, acridone, and benzalacetone with therapeutic potential. Simple architecture and broad substrate affinity of AmQNS make it as one of the target enzymes to produce novel structural scaffolds. Another unique feature of AmQNS despite its high similarity to acridone forming type III polyketide synthase from Citrus microcarpa is the variation in the product formation. Hence, to explore the characteristic features of AmQNS, an in-depth sequence and structure-based bioinformatics analyses were performed. Our studies indicated that AmQNS and its nearest homologs have evolved by a series of gene duplication events and strong purifying selection pressure constrains them in the evolutionary process. Additionally, some amino acid alterations were identified in the functionally important region(s), which can contribute to the functional divergence of the enzyme. Prediction of favorable amino acid substitutions will be advantageous in the metabolic engineering of AmQNS for the production of novel compounds. Furthermore, comparative modeling and docking studies were utilized to investigate the structural behavior and small molecule interaction pattern of AmQNS. The observations and results reported here are crucial for advancing our understanding of AmQNS's phylogenetic position, selection pressure, evolvability, interaction pattern and thus providing the foundation for further studies on the structural and reaction mechanism.

Amino acid composition analysis of human secondary transport proteins and implications for reliable membrane topology prediction

Secondary transporters in humans are a large group of proteins that transport a wide range of ions, metals, organic and inorganic solutes involved in energy transduction, control of membrane potential and osmotic balance, metabolic processes and in the absorption or efflux of drugs and xenobiotics. They are also emerging as important targets for development of new drugs and as target sites for drug delivery to specific organs or tissues. We have performed amino acid composition (AAC) and phylogenetic analyses and membrane topology predictions for 336 human secondary transport proteins and used the results to confirm protein classification and to look for trends and correlations with structural domains and specific substrates and/or function. Some proteins showed statistically high contents of individual amino acids or of groups of amino acids with similar physicochemical properties. One recurring trend was a correlation between high contents of charged and/or polar residues with misleading results in predictions of membrane topology, which was especially prevalent in Mitochondrial Carrier family proteins. We demonstrate how charged or polar residues located in the middle of transmembrane helices can interfere with their identification by membrane topology tools resulting in missed helices in the prediction. Comparison of AAC in the human proteins with that in 235 secondary transport proteins from Escherichia coli revealed similar overall trends along with differences in average contents for some individual amino acids and groups of similar amino acids that are presumed to result from a greater number of functions and complexity in the higher organism.

Transmembrane Domain Lengths Serve as Signatures of Organismal Complexity and Viral Transport Mechanisms

It is known that membrane proteins are important in various secretory pathways, with a possible role of their transmembrane domains (TMDs) as sorting determinant factors. One key aspect of TMDs associated with various “checkposts” (i.e. organelles) of intracellular trafficking is their length. To explore possible linkages in organisms with varying “complexity” and differences in TMD lengths of membrane proteins associated with different organelles (such as Endoplasmic Reticulum, Golgi, Endosomes, Nucleus, Plasma Membrane), we analyzed ~70000 membrane protein sequences in over 300 genomes of fungi, plants, non-mammalian vertebrates and mammals. We report that as we move from simpler to complex organisms, variation in organellar TMD lengths decreases, especially compared to their respective plasma membranes, with increasing organismal complexity. This suggests an evolutionary pressure in modulating length of TMDs of membrane proteins with increasing complexity of communication between sub-cellular compartments. We also report functional applications of our findings by discovering remarkable distinctions in TMD lengths of membrane proteins associated with different intracellular transport pathways. Finally, we show that TMD lengths extracted from viral proteins can serve as somewhat weak indicators of viral replication sites in plant cells but very strong indicators of different entry pathways employed by animal viruses.

Distance-dependent hydrophobic-hydrophobic contacts in protein folding simulations

Successful prediction of protein folding from an amino acid sequence is a challenge in computational biology. In order to reveal the geometric constraints that drive protein folding, highlight those constraints kept or missed by distinct lattices and for establishing which class of intra- and inter-secondary structure element interactions is the most relevant for the correct folding of proteins, we have calculated inter-alpha carbon distances in a set of 42 crystal structures consisting of mainly helix, sheet or mixed conformations. The inter-alpha carbon distances were also calculated in several lattice "hydrophobic-polar" models built from the same protein set. We found that helix structures are more prone to form "hydrophobic-hydrophobic" contacts than beta-sheet structures. At a distance lower than or equal to 3.8 Å (very short-range interactions), "hydrophobic-hydrophobic" contacts are almost absent in the native structures, while they are frequent in all the analyzed lattice models. At distances in-between 3.8 and 9.5 Å (short-/medium-range interactions), the best performing lattice for reproducing mainly helix structures is the body-centered-cubic lattice. If protein structures contain sheet portions, lattice performances get worse, with few exceptions observed for double-tetrahedral and body-centered-cubic lattices. Finally, we can observe that ab initio protein folding algorithms, i.e. those based on the employment of lattices and Monte Carlo simulated annealings, can be improved simply and effectively by preventing the generation of "hydrophobic-hydrophobic" contacts shorter than 3.8 Å, by monitoring the "hydrophobic-hydrophobic/polar-polar" contact ratio in short-/medium distance ranges and by using preferentially a body-centered-cubic lattice.

[Effect of point substitutions of Asp-714 and Asp-720 residues on the structure and function of the H+ -ATPase of the yeast plasma membrane]

Membrane-spanning M5 and M6 segments, which play a role in the formation of cation transport sites in H(+)-, Ca2(+)-, K(+)-, Na(+)-, and other P2-ATPases, are connected by a short extracytoplasmic loop. In the yeast plasma membrane H(+)-ATPase, which belongs to a family of P2-ATPases, the loop is connected to M5 and M6 through the Asp-714 and Asp-720 residues. In this work, the effect of point amino, acidreplacements of Asp-714 and Asp-720 by Ala, Val, Asn, and Glu residues on the function of the enzyme was studied. The Asp714Asn point mutant possessed activities similar to those of the wild-type enzyme, whereas the replacement of Asp-714 by other amino acid residues disrupted biogenesis and led to a loss of activity. All mutants with substitution of Asp-720 were expressed and possessed relatively high activity. The D720V mutant displayed significantly reduced expression levels, activity, H+ transport, and ATP hydrolyzing activity. Thus, substitutions of Asp-714, except for the D714N mutant, led to significant defects in biogenesis and/or function of the enzyme. The results indicate the important role for the Asp-714 residue in biogenesis, structure stability, and enzyme function.

Role of loop L5-6 connecting transmembrane segments M5 and M6 in biogenesis and functioning of yeast Pma1 H+-ATPase

The L5-6 loop is a short extracytoplasmic stretch (714-DNSLDID) connecting transmembrane segments M5 and M6 and forming along with segments M4 and M8 the core through which cations are transported by H+-, Ca2+-, K+,Na+-, H+,K+-, and other P2-ATPases. To study structure-function relationships within this loop of the yeast plasma membrane Pma1 H+-ATPase, alanine- and cysteine-scanning mutagenesis has been employed. Ala and Cys substitutions for the most conserved residue (Leu717) led to complete block in biogenesis preventing the enzyme from reaching secretory vesicles. The Ala replacement at Asp714 led to five-fold decrease in the mutant expression and loss of its activity, while the Cys substitution blocked biogenesis completely. Replacements of other residues did not lead to loss of enzymatic activity. Additional replacements were made for Asp714 and Asp720 (Asp®Asn/Glu). Of the substitutions made at Asp714, only D714N partially restored the mutant enzyme biogenesis and functioning. However, all mutant enzymes with substituted Asp720 were active. The expressed mutants (34-95% of the wild-type level) showed activity high enough (35-108%) to be analyzed in detail. One of the mutants (I719A) had three-fold reduced coupling ratio between ATP hydrolysis and H+ transport; however, the I719C mutation was rather indistinguishable from the wild-type enzyme. Thus, substitutions at two of the seven positions seriously affected biogenesis and/or functioning of the enzyme. Taken together, these results suggest that the M5-M6 loop residues play an important role in protein stability and function, and they are probably responsible for proper arrangement of transmembrane segments M5 and M6 and other domains of the enzyme. This might also be important for the regulation of the enzyme.

A statistical anomaly indicates symbiotic origins of eukaryotic membranes

We report a remarkable statistical paradox from a compositional analysis of membrane lipids of archaea, eubacteria, and eukaryotes. The presence of Simpson's paradox in the analysis of lipids common to plasma membranes of all domains provides the first evidence based on the compositional analysis of lipidomics data for the symbiotic origins of eukaryotic cells.

Compositional analyses of nucleic acids and proteins have shed light on possible origins of living cells. In this work, rigorous compositional analyses of ∼5000 plasma membrane lipid constituents of 273 species in the three life domains (archaea, eubacteria, and eukaryotes) revealed a remarkable statistical paradox, indicating symbiotic origins of eukaryotic cells involving eubacteria. For lipids common to plasma membranes of the three domains, the number of carbon atoms in eubacteria was found to be similar to that in eukaryotes. However, mutually exclusive subsets of same data show exactly the opposite—the number of carbon atoms in lipids of eukaryotes was higher than in eubacteria. This statistical paradox, called Simpson's paradox, was absent for lipids in archaea and for lipids not common to plasma membranes of the three domains. This indicates the presence of interaction(s) and/or association(s) in lipids forming plasma membranes of eubacteria and eukaryotes but not for those in archaea. Further inspection of membrane lipid structures affecting physicochemical properties of plasma membranes provides the first evidence (to our knowledge) on the symbiotic origins of eukaryotic cells based on the “third front” (i.e., lipids) in addition to the growing compositional data from nucleic acids and proteins.

CARd-3D: Carbon Distribution in 3D Structure Program for Globular Proteins

Spatial arrangement of carbon in protein structure is analyzed here. Particularly, the carbon fractions around individual atoms are compared. It is hoped that it follows the principle of 31.45% carbon around individual atoms. The results reveal that globular protein's atoms follow this principle. A comparative study on monomer versus dimer reveal that carbon is better distributed in dimeric form than in its monomeric form. Similar study on solid versus liquid structures reveals that the liquid (NMR) structure has better carbon distribution over the corresponding solid (X-Ray) structure. The carbon fraction distributions in fiber and toxin protein are compared. Fiber proteins follow the principle of carbon fraction distribution. At the same time it has another broad spectrum of carbon distribution than in globular proteins. The toxin protein follows an abnormal carbon fraction distribution. The carbon fraction distribution plays an important role in deciding the structure and shape of proteins. It is hoped to help in understanding the protein folding and function.

Point mutations in the extracytosolic loop between transmembrane segments M5 and M6 of the yeast Pma1 H+-ATPase: alanine-scanning mutagenesis

Membrane-spanning segments M4, M5, M6, and M8 of the H(+)-, Ca(2+)-, and K(+), Na(+)-ATPases, which belong to the P2-type pumps are the core through which cations are transported. M5 and M6 loop is a short extracytoplasmic stretch of the seven amino acid residues (714-DNSLDID) connecting two of these segments, M5 and M6, where residues involved in the formation of the proton-binding site(s) are located. In the present study, we have used alanine-scanning mutagenesis to explore the structural and functional relationships within this loop of the yeast plasma membrane Pma1 H(+)-ATPase. Of the 7 Ala mutants made, substitution for the most conserved residue (Leu-717) has led to a severe misfolding and complete block in biogenesis of the mutant enzyme. The replacement of Asp-714 has also caused misfolding leading to significant decrease in the expression of the mutant and loss of activity. The remaining mutants were expressed in secretory vesicles at 21-119% of the wild-type level and were active enough to be analyzed in detail. One of these mutants (I719A) showed five- to threefold decrease in both expression and ATP hydrolyzing and H(+) pumping activities and also threefold reduction in the coupling ratio between ATP hydrolysis and H(+) transport. Thus, Ala substitutions at three positions of the seven seriously affected biogenesis, folding, stability and/or functioning of the enzyme. Taken together, these results lead to suggestion that M5 and M6 loop play an important role in the protein stability and function and is responsible for proper arrangement of transmembrane segments M5 and M6 and probably other domains of the enzyme. Results for additional conserved substitutions (Asn and Glu) at Asp-714 and Asp-720 confirmed this suggestion.

A homology/ab initio hybrid algorithm for sampling near-native protein conformations

One of the major challenges for protein tertiary structure prediction strategies is the quality of conformational sampling algorithms, which can effectively and readily search the protein fold space to generate near-native conformations. In an effort to advance the field by making the best use of available homology as well as fold recognition approaches along with ab initio folding methods, we have developed Bhageerath-H Strgen, a homology/ab initio hybrid algorithm for protein conformational sampling. The methodology is tested on the benchmark CASP9 dataset of 116 targets. In 93% of the cases, a structure with TM-score ≥ 0.5 is generated in the pool of decoys. Further, the performance of Bhageerath-H Strgen was seen to be efficient in comparison with different decoy generation methods. The algorithm is web enabled as Bhageerath-H Strgen web tool which is made freely accessible for protein decoy generation (http://www.scfbio-iitd.res.in/software/Bhageerath-HStrgen1.jsp).

Capturing native/native like structures with a physico-chemical metric (pcSM) in protein folding

Specification of the three dimensional structure of a protein from its amino acid sequence, also called a "Grand Challenge" problem, has eluded a solution for over six decades. A modestly successful strategy has evolved over the last couple of decades based on development of scoring functions (e.g. mimicking free energy) that can capture native or native-like structures from an ensemble of decoys generated as plausible candidates for the native structure. A scoring function must be fast enough in discriminating the native from unfolded/misfolded structures, and requires validation on a large data set(s) to generate sufficient confidence in the score. Here we develop a scoring function called pcSM that detects true native structure in the top 5 with 93% accuracy from an ensemble of candidate structures. If we eliminate the native from ensemble of decoys then pcSM is able to capture near native structure (RMSD<=5Ǻ) in top 10 with 86% accuracy. The parameters considered in pcSM are a C-alpha Euclidean metric, secondary structural propensity, surface areas and an intramolecular energy function. pcSM has been tested on 415 systems consisting 142,698 decoys (public and CASP-largest reported hitherto in literature). The average rank for the native is 2.38, a significant improvement over that existing in literature. In-silico protein structure prediction requires robust scoring technique(s). Therefore, pcSM is easily amenable to integration into a successful protein structure prediction strategy. The tool is freely available at http://www.scfbio-iitd.res.in/software/pcsm.jsp.

Genomes to hits in silico - a country path today, a highway tomorrow: a case study of chikungunya

These are exciting times for bioinformaticians, computational biologists and drug designers with the genome and proteome sequences and related structural databases growing at an accelerated pace. The post-genomic era has triggered high expectations for a rapid and successful treatment of diseases. However, in this biological information rich and functional knowledge poor scenario, the challenges are indeed grand, no less than the assembly of the genome of the whole organism. These include functional annotation of genes, identification of druggable targets, prediction of three-dimensional structures of protein targets from their amino acid sequences, arriving at lead compounds for these targets followed by a transition from bench to bedside. We propose here a "Genome to Hits In Silico" strategy (called Dhanvantari) and illustrate it on Chikungunya virus (CHIKV). "Genome to hits" is a novel pathway incorporating a series of steps such as gene prediction, protein tertiary structure determination, active site identification, hit molecule generation, docking and scoring of hits to arrive at lead compounds. The current state of the art for each of the steps in the pathway is high-lighted and the feasibility of creating an automated genome to hits assembly line is discussed.

In silico investigation of the disease-associated retinoschisin C110Y and C219G mutants

The juvenile X-linked retinoschisis (XLRS) is a retinal disease caused by mutations in the secretory protein, retinoschisin (RS1). Majority of the disease is resulted from single point mutations on the RS1 discoidin domain with cysteine mutations being related to some of the more severe cases of XLRS. Previous studies have indicated that two mutations (C110Y and C219G), which involve cysteines that form intramolecular disulfide bonds in the native discoidin domain, resulted in different oligomerization states of the proteins and did not correlate with the degree of protein stability as calculated by the change in folding free energy. Through homology modeling, bioinformatics predictions, molecular dynamics (MD) and docking simulations, we attempt to investigate the effects of these two mutations on the structure of the RS1 discoidin domain in relevance to the discrepancy found between structural stability and aggregation propensity. Based on our findings, this discrepancy can be explained by the ability of C110Y mutant to establish suitable modules for initiating amorphous aggregation and to expand the aggregating mass through predominantly hydrophobic interactions. The low capability of C219G mutant to oligomerize, on the other hand, may be due to its greater structural instability and lesser hydrophobic tendency, two properties that may be unsupportive of aggregation. The results, altogether, indicate that aggregation propensity in the RS1 C110Y mutant is dependent upon the formation of suitable aggregating substrates for propagation of aggregation and not directly related to or determined by overall structural instability. As for the wildtype protein, the binding specificity of the spikes for biological function and the formation of octameric structure are contributed by important loop interactions, as well as evolved structural and sequence-based properties that prevent aggregation.

LoCo: a novel main chain scoring function for protein structure prediction based on local coordinates

BACKGROUND

Successful protein structure prediction requires accurate low-resolution scoring functions so that protein main chain conformations that are close to the native can be identified. Once that is accomplished, a more detailed and time-consuming treatment to produce all-atom models can be undertaken. The earliest low-resolution scoring used simple distance-based "contact potentials," but more recently, the relative orientations of interacting amino acids have been taken into account to improve performance.

RESULTS

We developed a new knowledge-based scoring function, LoCo, that locates the interaction partners of each individual residue within a local coordinate system based only on the position of its main chain N, Cα and C atoms. LoCo was trained on a large set of experimentally determined structures and optimized using standard sets of modeled structures, or "decoys." No structure used to train or optimize the function was included among those used to test it. When tested against 29 other published main chain functions on a group of 77 commonly used decoy sets, our function outperformed all others in Cα RMSD rank of the best-scoring decoy, with statistically significant p-values < 0.05 for 26 out of the 29 other functions considered. LoCo is fast, requiring on average less than 6 microseconds per residue for interaction and scoring on commonly-used computer hardware.

CONCLUSIONS

Our function demonstrates an unmatched combination of accuracy, speed, and simplicity and shows excellent promise for protein structure prediction. Broader applications may include protein-protein interactions and protein design.