Efficient sequential assignments in proteins with reduced dimensionality 3D HN(CA)NH

We present reduced dimensionality (RD) 3D HN(CA)NH for efficient sequential assignment in proteins. The experiment correlates the N-15 and H-1 chemical shift of a residue ('i') with those of its immediate N-terminal (i - 1) and C-terminal (i + 1) neighbors and provides four-dimensional chemical shift correlations rapidly with high resolution. An assignment strategy is presented which combines the correlations observed in this experiment with amino acid type information obtained from 3D CBCA(CO)NH. By classifying the 20 amino acid types into seven distinct categories based on C-13(beta) chemical shifts, it is observed that a stretch of five sequentially connected residues is sufficient to map uniquely on to the polypeptide for sequence specific resonance assignments. This method is exemplified by application to three different systems: maltose binding protein (42 kDa), intrinsically disordered domain of insulin-like growth factor binding protein-2 and Ubiquitin. Fast data acquisition is demonstrated using longitudinal H-1 relaxation optimization. Overall, 3D HN(CA)NH is a powerful tool for high throughput resonance assignment, in particular for unfolded or intrinsically disordered polypeptides.

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins

Short range side chain-backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution = 1.5 angstrom). Hydrogen bonds involving residues i - 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (Cn, where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C10 (COdi .NHi + 2), C13 (COdi .NHi + 3) and C17 (NdHi .COi - 4) structures. In contrast, Gln predominantly forms C16 (COei .NHi - 3), C12 (NeHi .COi - 2), C15 (NeHi .COi - 3) and C18 (NeHi .COi - 4) motifs, with only the C18motif being analogous to the Asn C17structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone beta-turns and a-turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain-backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered. Proteins 2012; (c) 2011 Wiley Periodicals, Inc.

Reduced dimensionality 3D HNCAN for unambiguous HN, CA and N assignment in proteins

We present here an improvisation of HNN (Panchal, Bhavesh et al., 2001) called RD 3D HNCAN for backbone (HN, CA and N-15) assignment in both folded and unfolded proteins. This is a reduced dimensionality experiment which employs CA chemical shifts to improve dispersion. Distinct positive and negative peak patterns of various triplet segments along the polypeptide chain observed in HNN are retained and these provide start and check points for the sequential walk. Because of co-incrementing of CA and N-15, peaks along one of the dimensions appear at sums and differences of the CA and N-15 chemical shifts. This changes the backbone assignment protocol slightly and we present this in explicit detail. The performance of the experiment has been demonstrated using Ubiquitin and Plasmodium falciparum P2 proteins. The experiment is particularly valuable when two neighboring amino acid residues have nearly identical backbone N-15 chemical shifts. (C) 2012 Elsevier Inc. All rights reserved.

Mutations in Drosophila Myosin Rod Cause Defects in Myofibril Assembly

The roles of myosin during muscle contraction are well studied, but how different domains of this protein are involved in myofibril assembly in vivo is far less understood. The indirect flight muscles (IFMs) of Drosophila melanogaster provide a good model for understanding muscle development and function in vivo. We show that two missense mutations in the rod region of the myosin heavy-chain gene, Mhc, give rise to IFM defects and abnormal myofibrils. These defects likely result from thick filament abnormalities that manifest during early sarcomere development or later by hypercontraction. The thick filament defects are accompanied by marked reduction in accumulation of flightin, a myosin binding protein, and its phosphorylated forms, which are required to stabilise thick filaments. We investigated with purified rod fragments whether the mutations affect the coiled-coil structure, rod aggregate size or rod stability. No significant changes in these parameters were detected, except for rod thermodynamic stability in one mutation. Molecular dynamics simulations suggest that these mutations may produce localised rod instabilities. We conclude that the aberrant myofibrils are a result of thick filament defects, but that these in vivo effects cannot be detected in vitro using the biophysical techniques employed. The in vivo investigation of these mutant phenotypes in IFM development and function provides a useful platform for studying myosin rod and thick filament formation generically, with application to the aetiology of human myosin rod myopathies. (C) 2012 Elsevier Ltd. All rights reserved.

IDH1 Mutations in Diffusely Infiltrating Astrocytomas Grade Specificity, Association With Protein Expression, and Clinical Relevance

IDH1 mutations are frequent genetic alterations in low-grade diffuse gliomas and secondary glioblastoma (GBM). To validate mutation frequency, IDH1 gene at codon 132 was sequenced in 74 diffusely infiltrating astrocytomas: diffuse astrocytoma (DA; World Health Organization WHO] grade II), anaplastic astrocytoma (AA; WHO grade III), and GBM (WHO grade IV). All cases were immunostained with IDH1-R132H monoclonal antibody. Mutational status was correlated with mutant protein expression, patient age, duration of symptoms, and prognosis of patients with GBM. We detected 31 (41.9%) heterozygous IDH1 mutations resulting in arginine-to-histidine substitution (R132H;CGT-CAT). All 12 DAs (100%), 13 of 14 AAs (92.9%), and 6 of 48 GBMs (12.5%) (5/6 83.3%] secondary, and 1/42 2.4%] primary) harbored IDH1 mutations. The correlation between mutational status and protein expression was significant (P < .001). IDH1 mutation status, though not associated with prognosis of patients with GBM, showed significant association with younger age and longer duration of symptoms in the whole cohort (P < .001). Our study validates IDH1 mutant protein expression across various grades of astrocytoma, and demonstrates a high incidence of IDH1 mutations in DA, AA, and secondary GBM.

An algorithm to find all palindromic sequences in proteins

A palindrome is a set of characters that reads the same forwards and backwards. Since the discovery of palindromic peptide sequences two decades ago, little effort has been made to understand its structural, functional and evolutionary significance. Therefore, in view of this, an algorithm has been developed to identify all perfect palindromes (excluding the palindromic subset and tandem repeats) in a single protein sequence. The proposed algorithm does not impose any restriction on the number of residues to be given in the input sequence. This avant-garde algorithm will aid in the identification of palindromic peptide sequences of varying lengths in a single protein sequence.

Clinical phenotype and the lack of mutations in the CHRNG, CHRND, and CHRNA1 genes in two Indian families with Escobar syndrome

The objective of this study was to report the clinical phenotype and genetic analysis of two Indian families with Escobar syndrome (ES). The diagnosis of ES in both families was made on the basis of published clinical features. Blood samples were collected from members of both families and used in genomic DNA isolation. The entire coding regions and intron-exon junctions of the ES gene CHRNG (cholinergic receptor, nicotinic, gamma), and two other related genes, CHRND and CHRNA1, were amplified and sequenced to search for mutations in both families. Both families show a typical form of ES. Sequencing of the entire coding regions including the intron-exon junctions of the three genes did not yield any mutations in these families. In conclusion, it is possible that the mutations in these genes are located in the promoter or deep intronic regions that we failed to identify or the ES in these families is caused by mutations in a different gene. The lack of mutations in CHRNG has also been reported in several families, suggesting the possibility of at least one more gene for this syndrome. Clin Dysmorphol 22:54-58 (C) 2013 Wolters Kluwer Health vertical bar Lippincott Williams & Wilkins.

Investigating proline puckering states in diproline segments in proteins

In the current study, the puckering states of the Proline ring occurring in diproline segments (LPro-LPro) in proteins has been investigated with a segregation made on the basis of cis and trans states for the Pro-Pro peptide bond and the conformational states for the diproline segment to investigate the effects of conformation of the diproline segment on the corresponding puckering state of the Proline ring in the segment if any. The value of the endocyclic ring torsional angles of the pyrrolidine ring has been used for calculating and visualizing various puckering states using a proposed new sign convention (+/-) nomenclature. The results have been compared to that obtained in a previous study on peptides from this group. In this study, quite interestingly, the Planar (G) conformation that was present in 14.3% of the cases in peptides, appears to be nearly a rare conformation in the case of proteins (1.9%). The present study indicates that the (C-exo/C-exo), (C-exo/Twisted C-exo-C-endo) and (Twisted C-endo-C-exo/Twisted C-endo-C-exo) categories are the most preferred combinations. For Proline rings in proteins, the states C-exo, Twisted C-exo-C-endo and Twisted C-endo-C-exo are the most preferred states. Within diproline segments, the pyrrolidine ring conformations do not show a strong co-relation to the backbone conformation in which they are observed. It is likely that five-membered rings have a considerable plasticity of structure and are readily deformed to accommodate a variety of energetically preferred backbone conformations.

Depth: a web server to compute depth, cavity sizes, detect potential small-molecule ligand-binding cavities and predict the pK(a) of ionizable residues in proteins

Residue depth accurately measures burial and parameterizes local protein environment. Depth is the distance of any atom/residue to the closest bulk water. We consider the non-bulk waters to occupy cavities, whose volumes are determined using a Voronoi procedure. Our estimation of cavity sizes is statistically superior to estimates made by CASTp and VOIDOO, and on par with McVol over a data set of 40 cavities. Our calculated cavity volumes correlated best with the experimentally determined destabilization of 34 mutants from five proteins. Some of the cavities identified are capable of binding small molecule ligands. In this study, we have enhanced our depth-based predictions of binding sites by including evolutionary information. We have demonstrated that on a database (LigASite) of similar to 200 proteins, we perform on par with ConCavity and better than MetaPocket 2.0. Our predictions, while less sensitive, are more specific and precise. Finally, we use depth (and other features) to predict pK(a)s of GLU, ASP, LYS and HIS residues. Our results produce an average error of just <1 pH unit over 60 predictions. Our simple empirical method is statistically on par with two and superior to three other methods while inferior to only one. The DEPTH server (http://mspc.bii.a-star.edu.sg/depth/) is an ideal tool for rapid yet accurate structural analyses of protein structures.

Non-Syndromic Hearing Impairment in India: High Allelic Heterogeneity among Mutations in TMPRSS3, TMC1, USHIC, CDH23 and TMIE

Mutations in the autosomal genes TMPRSS3, TMC1, USHIC, CDH23 and TMIE are known to cause hereditary hearing loss. To study the contribution of these genes to autosomal recessive, non-syndromic hearing loss (ARNSHL) in India, we examined 374 families with the disorder to identify potential mutations. We found four mutations in TMPRSS3, eight in TMC1, ten in USHIC, eight in CDH23 and three in TMIE. Of the 33 potentially pathogenic variants identified in these genes, 23 were new and the remaining have been previously reported. Collectively, mutations in these five genes contribute to about one-tenth of ARNSHL among the families examined. New mutations detected in this study extend the allelic heterogeneity of the genes and provide several additional variants for structure-function correlation studies. These findings have implications for early DNA-based detection of deafness and genetic counseling of affected families in the Indian subcontinent.

PROPENSITIES OF AMINO ACID RESIDUES IN PROTEINS FOR DIFFERENT REGIONS OF THE RAMACHANDRAN MAP

Using a dataset of 1164 crystal structures of largely non-homologous proteins defined at a resolution of 1.5 angstrom or better, we have investigated the (phi,psi) preferences of 20 residue types by considering the residues which occur in loops. Propensities of residue types to occur in the loops with (phi,psi) values in the aa region of the Ramachandran map has a poor correlation coefficient of 0.48 to the Chou-Fasman propensities of the residue types to occur in the a-helical segments. However the correlation coefficient between propensities of residues in loops to adopt beta conformations and those in beta-sheet is much higher (0.95). These observations suggest that a-helix formation is well influenced by the local amino acid sequence while intrinsic preference of residue types for beta-sheet plays a major role in the formation of beta-sheet. The main chain polar groups of residues in loops, that can affect the (phi,psi) values, can be involved in intra-molecular hydrogen bonding. Therefore we investigated further by considering subset of residues in loops with low (0 to 2) number of intra-molecular hydrogen bonds per residue involving main chain polar atoms. For this subset, the correlation coefficients between propensities for alpha-helix and alpha(R) region and between beta-sheet and beta-region are 0.26 and 0.64 respectively. This reiterates higher intrinsic tendency of beta-region favouring residues to adopt beta-sheet than alpha(R) region favouring residues to adopt alpha-helical structure.

Identifying functionally important cis-peptide containing segments in proteins and their utility in molecular function annotation

Cis-peptide embedded segments are rare in proteins but often highlight their important role in molecular function when they do occur. The high evolutionary conservation of these segments illustrates this observation almost universally, although no attempt has been made to systematically use this information for the purpose of function annotation. In the present study, we demonstrate how geometric clustering and level-specific Gene Ontology molecular-function terms (also known as annotations) can be used in a statistically significant manner to identify cis-embedded segments in a protein linked to its molecular function. The present study identifies novel cis-peptide fragments, which are subsequently used for fragment-based function annotation. Annotation recall benchmarks interpreted using the receiver-operator characteristic plot returned an area-under-curve >0.9, corroborating the utility of the annotation method. In addition, we identified cis-peptide fragments occurring in conjunction with functionally important trans-peptide fragments, providing additional insights into molecular function. We further illustrate the applicability of our method in function annotation where homology-based annotation transfer is not possible. The findings of the present study add to the repertoire of function annotation approaches and also facilitate engineering, design and allied studies around the cis-peptide neighborhood of proteins.

Isotopic fractionation in proteins as a measure of hydrogen bond length

If a deuterated molecule containing strong intramolecular hydrogen bonds is placed in a hydrogenated solvent, it may preferentially exchange deuterium for hydrogen. This preference is due to the difference between the vibrational zero-point energy for hydrogen and deuterium. It is found that the associated fractionation factor (I) is correlated with the strength of the intramolecular hydrogen bonds. This correlation has been used to determine the length of the H-bonds (donor-acceptor separation) in a diverse range of enzymes and has been argued to support the existence of short low-barrier H-bonds. Starting with a potential energy surface based on a simple diabatic state model for H-bonds, we calculate (I) as a function of the proton donor-acceptor distance R. For numerical results, we use a parameterization of the model for symmetric 0-H. ``.0 bonds R. H. McKenzie, Chem. Phys. Lett. 535, 196 (2012)]. We consider the relative contributions of the 0-H stretch vibration, O-H bend vibrations (both in plane and out of plane), tunneling splitting effects at finite temperature, and the secondary geometric isotope effect. We compare our total (I) as a function of R with NMR experimental results for enzymes, and in particular with an earlier model parametrization (D(R), used previously to determine bond lengths. (C) 2015 AIP Publishing LLC.

Mutations in the nucleotide binding and hydrolysis domains of Helicobacter pylori MutS2 lead to altered biochemical activities and inactivation of its in vivo function

Background: Helicobacter pylori MutS2 (HpMutS2), an inhibitor of recombination during transformation is a non-specific nuclease with two catalytic sites, both of which are essential for its anti-recombinase activity. Although HpMutS2 belongs to a highly conserved family of ABC transporter ATPases, the role of its ATP binding and hydrolysis activities remains elusive. Results: To explore the putative role of ATP binding and hydrolysis activities of HpMutS2 we specifically generated point mutations in the nucleotide-binding Walker-A (HpMutS2-G338R) and hydrolysis Walker-B (HpMutS2-E413A) domains of the protein. Compared to wild-type protein, HpMutS2-G338R exhibited similar to 2.5-fold lower affinity for both ATP and ADP while ATP hydrolysis was reduced by similar to 3-fold. Nucleotide binding efficiencies of HpMutS2-E413A were not significantly altered; however the ATP hydrolysis was reduced by similar to 10-fold. Although mutations in the Walker-A and Walker-B motifs of HpMutS2 only partially reduced its ability to bind and hydrolyze ATP, we demonstrate that these mutants not only exhibited alterations in the conformation, DNA binding and nuclease activities of the protein but failed to complement the hyper-recombinant phenotype displayed by mutS2-disrupted strain of H. pylori. In addition, we show that the nucleotide cofactor modulates the conformation, DNA binding and nuclease activities of HpMutS2. Conclusions: These data describe a strong crosstalk between the ATPase, DNA binding, and nuclease activities of HpMutS2. Furthermore these data show that both, ATP binding and hydrolysis activities of HpMutS2 are essential for the in vivo anti-recombinase function of the protein.