Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety-six of these alignment positions correspond to conserved Gly residue with (phi, ) values allowed for non-glycyl residues. Reasons for this observation were investigated by in-silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (phi, ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (phi, ) values which are disallowed for Ala. In-silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.

Accurate prediction of interfacial residues in two-domain proteins using evolutionary information: Implications for three-dimensional modeling

With the preponderance of multidomain proteins in eukaryotic genomes, it is essential to recognize the constituent domains and their functions. Often function involves communications across the domain interfaces, and the knowledge of the interacting sites is essential to our understanding of the structure-function relationship. Using evolutionary information extracted from homologous domains in at least two diverse domain architectures (single and multidomain), we predict the interface residues corresponding to domains from the two-domain proteins. We also use information from the three-dimensional structures of individual domains of two-domain proteins to train naive Bayes classifier model to predict the interfacial residues. Our predictions are highly accurate (approximate to 85%) and specific (approximate to 95%) to the domain-domain interfaces. This method is specific to multidomain proteins which contain domains in at least more than one protein architectural context. Using predicted residues to constrain domain-domain interaction, rigid-body docking was able to provide us with accurate full-length protein structures with correct orientation of domains. We believe that these results can be of considerable interest toward rational protein and interaction design, apart from providing us with valuable information on the nature of interactions. Proteins 2014; 82:1219-1234. (c) 2013 Wiley Periodicals, Inc.

The HORMA domain: an evolutionarily conserved domain discovered in chromatin-associated proteins, has unanticipated diverse functions

The HORMA domain (for Hop1p, Rev7p and MAD2) was discovered in three chromatin-associated proteins in the budding yeast Saccharomyces cerevisiae. This domain has also been found in proteins with similar functions in organisms including plants, animals and nematodes. The HORMA domain containing proteins are thought to function as adaptors for meiotic checkpoint protein signaling and in the regulation of meiotic recombination. Surprisingly, new work has disclosed completely unanticipated and diverse functions for the HORMA domain containing proteins. A. M. Villeneuve and colleagues (Schvarzstein et al., 2013) show that meiosis-specific HORMA domain containing proteins plays a vital role in preventing centriole disengagement during Caenorhabditis elegans spermatocyte meiosis. Another recent study reveals that S. cerevisiae Atg13 HORMA domain acts as a phosphorylation-dependent conformational switch in the cellular autophagic process. (C) 2014 Elsevier B.V. All rights reserved.

Sequence and conformational preferences at termini of alpha-helices in membrane proteins: Role of the helix environment

-helices are amongst the most common secondary structural elements seen in membrane proteins and are packed in the form of helix bundles. These -helices encounter varying external environments (hydrophobic, hydrophilic) that may influence the sequence preferences at their N and C-termini. The role of the external environment in stabilization of the helix termini in membrane proteins is still unknown. Here we analyze -helices in a high-resolution dataset of integral -helical membrane proteins and establish that their sequence and conformational preferences differ from those in globular proteins. We specifically examine these preferences at the N and C-termini in helices initiating/terminating inside the membrane core as well as in linkers connecting these transmembrane helices. We find that the sequence preferences and structural motifs at capping (Ncap and Ccap) and near-helical (N' and C') positions are influenced by a combination of features including the membrane environment and the innate helix initiation and termination property of residues forming structural motifs. We also find that a large number of helix termini which do not form any particular capping motif are stabilized by formation of hydrogen bonds and hydrophobic interactions contributed from the neighboring helices in the membrane protein. We further validate the sequence preferences obtained from our analysis with data from an ultradeep sequencing study that identifies evolutionarily conserved amino acids in the rat neurotensin receptor. The results from our analysis provide insights for the secondary structure prediction, modeling and design of membrane proteins. Proteins 2014; 82:3420-3436. (c) 2014 Wiley Periodicals, Inc.

A carboxy terminal domain of the L protein of rinderpest virus possesses RNA triphosphatase activity - The first enzyme in the viral mRNA capping pathway

The large protein L of negative-sense RNA viruses is a multifunctional protein involved in transcription and replication of genomic RNA. It also possesses enzymatic activities involved in capping and methylation of viral mRNAs. The pathway for mRNA capping followed by the L protein of the viruses in the Morbillivirus genus has not been established, although it has been speculated that these viruses may follow the unconventional capping pathway as has been shown for some viruses of Rhabdoviridae family. We had earlier shown that the large protein L of Rinderpest virus expressed as recombinant L-P complex in insect cells as well as the ribonucleoprotein complex from purified virus possesses RNA triphosphatase (RTPase) and guanylyltransferase activities, in addition to RNA dependent RNA polymerase activity. In the present work, we demonstrate that RTPase as well as nucleoside triphosphatase (NTPase) activities are exhibited by a subdomain of the L protein in the C terminal region (a.a. 1640 1840). The RTPase activity depends absolutely on a divalent cation, either magnesium or manganese. Both the RTPase and NTPase activities of the protein show dual metal specificity. Two mutant proteins having alanine mutations in the glutamic acid residues in motif-A of the RTPase domain did not show RTPase activity, while exhibiting reduced NTPase activity suggesting overlapping active sites for the two enzymatic functions. The RTPase and NTPase activities of the L subdomain resemble those of the Vaccinia capping enzyme D1 and the baculovirus LEF4 proteins. (C) 2015 Elsevier Inc. All rights reserved.

BRCA1 and BRCA2 germline mutation analysis among Indian women from south India: identification of four novel mutations and high-frequency occurrence of 185delAG mutation

Mutations in the BRCA1 and BRCA2 genes profoundly increase the risk of developing breast and/or ovarian cancer among women. To explore the contribution of BRCA1 and BRCA2 mutations in the development of hereditary breast cancer among Indian women, we carried out mutation analysis of the BRCA1 and BRCA2 genes in 61 breast or ovarian cancer patients from south India with a positive family history of breast and/or ovarian cancer. Mutation analysis was carried out using conformation-sensitive gel electrophoresis (CSGE) followed by sequencing. Mutations were identified in 17 patients (28.0%); 15 (24.6%) had BRCA1 mutations and two (3.28%) had BRCA2 mutations. While no specific association between BRCA1 or BRCA2 mutations with cancer type was seen, mutations were more often seen in families with ovarian cancer. While 40% (4/10) and 30.8% (4/12) of families with ovarian or breast and ovarian cancer had mutations, only 23.1% (9/39) of families with breast cancer carried mutations in the BRCA1 and BRCA2 genes. In addition, while BRCA1 mutations were found in all age groups, BRCA2 mutations were found only in the age group of <= 40 years. Of the BRCA1 mutations, there were three novel mutations (295delCA; 4213T -> A; 5267T -> G) G) and three mutations that have been reported earlier. Interestingly, 185delAG, a BRCA1 mutation which occurs at a very high frequency in Ashkenazi Jews, was found at a frequency of 16.4% (10/61). There was one novel mutation (4866insT) and one reported mutation in BRCA2. Thus, our study emphasizes the importance of mutation screening in familial breast and/or ovarian cancers, and the potential implications of these findings in genetic counselling and preventive therapy.

Structure, dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X-ray studies

Molecular dynamics simulations have been carried out on all the jacalin-carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X-ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin-carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate-binding residues are consistent with the known thermodynamic parameters of jacalin-carbohydrate interactions. The simulations, along with X-ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin-carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit. Proteins 2009; 77:760-777.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

MIPS: metal interactions in protein structures

MIPS (metal interactions in protein structures) is a database of metals in the three-dimensional acromolecular structures available in the Protein Data Bank. Bound metal ions in proteins have both catalytic and structural functions. The proposed database serves as an open resource for the analysis and visualization of all metals and their interactions with macromolecular (protein and nucleic acid) structures. MIPS can be searched via a user-friendly interface, and the interactions between metals and protein molecules, and the geometric parameters, can be viewed in both textual and graphical format using the freely available graphics plug-in Jmol. MIPS is updated regularly, by means of programmed scripts to find metal-containing proteins from newly released protein structures. The database is useful for studying the properties of coordination between metals and protein molecules. It also helps to improve understanding of the relationship between macromolecular structure and function. This database is intended to serve the scientific community working in the areas of chemical and structural biology, and is freely available to all users, around the clock, at http://dicsoft2.physics.iisc.ernet.in/mips/.

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.

A distinct physiological role of MutY in mutation prevention in mycobacteria

Oxidative damage to DNA results in the occurrence of 7,8-dihydro-B-oxoguanine (8-oxoG) in the genome. In eubacteria, repair of such damage is initiated by two major base-excision repair enzymes, MutM and MutY. We generated a MutY-deficient strain of Mycobacterium smegmatis to investigate the role of this enzyme in DNA repair. The MutY deficiency in M. smegmatis did not result in either a noteworthy susceptibility to oxidative stress or an increase in the mutation rate. However, rifampicin resistant isolates of the MutY-deficient strain showed distinct mutations in the rifampicin-resistance-determining region of rpoB. Besides the expected C to A (or G to T) mutations, an increase in A to C (or T to G) mutations was also observed. Biochemical characterization of mycobacterial MutY (M. smegmatis and M. tuberculosis) revealed an expected excision of A opposite 8-oxoG in DNA. Additionally, excision of G and T opposite 8-oxoG was detected. MutY formed complexes with DNA containing 8-oxoG: A, 8-oxoG: G or 8-oxoG: T but not 8-oxoG : C pairs. Primer extension reactions in cell-free extracts of M. smegmatis suggested error-prone incorporation of nucleotides into the DNA. Based on these observations, we discuss the physiological role of MutY in specific mutation prevention in mycobacteria.

A Method to Find Sequentially Separated Motifs in Biological Sequences (SSMBS)

Sequence motifs occurring in a particular order in proteins or DNA have been proved to be of biological interest. In this paper, a new method to locate the occurrences of up to five user-defined motifs in a specified order in large proteins and in nucleotide sequence databases is proposed. It has been designed using the concept of quantifiers in regular expressions and linked lists for data storage. The application of this method includes the extraction of relevant consensus regions from biological sequences. This might be useful in clustering of protein families as well as to study the correlation between positions of motifs and their functional sites in DNA sequences.

Role of 16S ribosomal RNA methylations in translation initiation in Escherichia coli

Translation initiation from the ribosomal P-site is the specialty of the initiator tRNAs (tRNA(fMet)). Presence of the three consecutive G-C base pairs (G29-C41, G30-C40 and G31-C39) in their anticodon stems, a highly conserved feature of the initiator tRNAs across the three kingdoms of life, has been implicated in their preferential binding to the P-site. How this feature is exploited by ribosomes has remained unclear. Using a genetic screen, we have isolated an Escherichia coli strain, carrying a G122D mutation in folD, which allows initiation with the tRNA(fMet) containing mutations in one, two or all the three G-C base pairs. The strain shows a severe deficiency of methionine and S-adenosylmethionine, and lacks nucleoside methylations in rRNA. Targeted mutations in the methyltransferase genes have revealed a connection between the rRNA modifications and the fundamental process of the initiator tRNA selection by the ribosome.

Analysis of the impact of a uracil DNA glycosylase attenuated in AP-DNA binding in maintenance of the genomic integrity in Escherichia coli

Uracil DNA glycosylase (Ung)initiates the uracil excision repair pathway. We have earlier characterized the Y66W and Y66H mutants of Ung and shown that they are compromised by similar to 7- and similar to 170-fold, respectively in their uracil excision activities. In this study, fluorescence anisotropy measurements show that compared with the wild-type, the Y66W protein is moderately compromised and attenuated in binding to AP-DNA. Allelic exchange of ung in Escherichia coli with ung::kan, ungY66H:amp or ungY66W:amp alleles showed similar to 5-, similar to 3.0- and similar to 2.0-fold, respectively increase in mutation frequencies. Analysis of mutations in the rifampicin resistance determining region of rpoB revealed that the Y66W allele resulted in an increase in A to G (or T to C) mutations. However, the increase in A to G mutations was mitigated upon expression of wild-type Ung from a plasmid borne gene. Biochemical and computational analyses showed that the Y66W mutant maintains strict specificity for uracil excision from DNA. Interestingly, a strain deficient in AP-endonucleases also showed an increase in A to G mutations. We discuss these findings in the context of a proposal that the residency of DNA glycosylase(s) onto the AP-sites they generate shields them until recruitment of AP-endonucleases for further repair.

Novel insertion and deletion mutants of RpoB that render Mycobacterium smegmatis RNA polymerase resistant to rifampicin-mediated inhibition of transcription

The startling increase in the occurrence of rifampicin (Rif) resistance in the clinical isolates of Mycobacterium tuberculosis worldwide is posing a serious concern to tuberculosis management. The majority of Rif resistance in bacteria arises from mutations in the RpoB subunit of the RNA polymerase. We isolated M. smegmatis strains harbouring either an insertion (6 aa) or a deletion (10 aa) in their RpoB proteins. Although these strains showed a compromised fitness for growth in 7H9 Middlebrook medium, their resistance to Rif was remarkably high. The attenuated growth of the strains correlated with decreased specific activities of the RNA polymerases from the mutants. While the RNA polymerases from the parent or a mutant strain (harbouring a frequently occurring mutation, H442Y, in RpoB) were susceptible to Rif-mediated inhibition of transcription from calf thymus DNA, those from the insertion and deletion mutants were essentially refractory to such inhibition. Three-dimensional structure modelling revealed that the RpoB amino acids that interact with Rif are either deleted or unable to interact with Rif due to their unsuitable spatial positioning in these mutants. We discuss possible uses of the RpoB mutants in studying transcriptional regulation in mycobacteria and as potential targets for drug design.