Family-based association analysis with ordered categorical phenotypes, covariates and interactions

Genetic association analyses of family-based studies with ordered categorical phenotypes are often conducted using methods either for quantitative or for binary traits, which can lead to suboptimal analyses. Here we present an alternative likelihood-based method of analysis for single nucleotide polymorphism (SNP) genotypes and ordered categorical phenotypes in nuclear families of any size. Our approach, which extends our previous work for binary phenotypes, permits straightforward inclusion of covariate, gene-gene and gene-covariate interaction terms in the likelihood, incorporates a simple model for ascertainment and allows for family-specific effects in the hypothesis test. Additionally, our method produces interpretable parameter estimates and valid confidence intervals. We assess the proposed method using simulated data, and apply it to a polymorphism in the c-reactive protein (CRP) gene typed in families collected to investigate human systemic lupus erythematosus. By including sex interactions in the analysis, we show that the polymorphism is associated with anti-nuclear autoantibody (ANA) production in females, while there appears to be no effect in males.

A likelihood ratio appropach to family-based association studies with covariates

We introduce a procedure for association based analysis of nuclear families that allows for dichotomous and more general measurements of phenotype and inclusion of covariate information. Standard generalized linear models are used to relate phenotype and its predictors. Our test procedure, based on the likelihood ratio, unifies the estimation of all parameters through the likelihood itself and yields maximum likelihood estimates of the genetic relative risk and interaction parameters. Our method has advantages in modelling the covariate and gene-covariate interaction terms over recently proposed conditional score tests that include covariate information via a two-stage modelling approach. We apply our method in a study of human systemic lupus erythematosus and the C-reactive protein that includes sex as a covariate.

Hepatitis C virus NS5A protein binds the SH3 domain of the Fyn tyrosine kinase with high affinity: mutagenic analysis of residues within the SH3 domain that contribute to the interaction

Background: The hepatitis C virus (HCV) non-structural 5A protein (NS5A) contains a highly conserved C-terminal polyproline motif with the consensus sequence Pro-X-X- Pro-X-Arg that is able to interact with the Src-homology 3 (SH3) domains of a variety of cellular proteins. Results: To understand this interaction in more detail we have expressed two N-terminally truncated forms of NS5A in E. coli and examined their interactions with the SH3 domain of the Src-family tyrosine kinase, Fyn. Surface plasmon resonance analysis revealed that NS5A binds to the Fyn SH3 domain with what can be considered a high affinity SH3 domain-ligand interaction (629 nM), and this binding did not require the presence of domain I of NS5A (amino acid residues 32-250). Mutagenic analysis of the Fyn SH3 domain demonstrated the requirement for an acidic cluster at the C-terminus of the RT-Src loop of the SH3 domain, as well as several highly conserved residues previously shown to participate in SH3 domain peptide binding. Conclusion: We conclude that the NS5A: Fyn SH3 domain interaction occurs via a canonical SH3 domain binding site and the high affinity of the interaction suggests that NS5A would be able to compete with cognate Fyn ligands within the infected cell.

The structure of echovirus type 12 bound to a two-domain fragment of its cellular attachment protein decay-accelerating factor (CD 55)

Echovirus type 12 (EV12), an enterovirus of the Picornaviridae family, uses the complement regulator, decay-accelerating factor (DAF, CD55) as a cellular receptor. We have calculated a three-dimensional reconstruction of EV12 bound to a fragment of DAF, consisting of short consensus repeat domains 3 and 4, from cryo-negative stain electron microscopy data (EMD #1057). This shows that, as for an earlier reconstruction of the related echovirus type 7 bound to DAF, attachment is not within the viral canyon but occurs close to the two-fold symmetry axes. Despite this general similarity, our reconstruction reveals a receptor interaction that is quite different from that observed for EV7. Fitting of the crystallographic co-ordinates for DAF34 and EV11 into the reconstruction shows a close agreement between the crystal structure of the receptor fragment and the density for the virus-bound receptor, allowing unambiguous positioning of the receptor with respect to the virion (PDB #1UPN). Our finding that the mode of virus-receptor interaction in EV12 is distinct from that seen for EV7 raises interesting questions regarding the evolution and biological significance of the DAF-binding phenotype in these viruses.

Mapping of domains on HIV envelope protein mediating association with calnexin and protein-disulfide isomerase.

The cell catalysts calnexin (CNX) and protein-disulfide isomerase (PDI) cooperate in establishing the disulfide bonding of the HIV envelope (Env) glycoprotein. Following HIV binding to lymphocytes, cell-surface PDI also reduces Env to induce the fusogenic conformation. We sought to define the contact points between Env and these catalysts to illustrate their potential as therapeutic targets. In lysates of Env-expressing cells, 15% of the gp160 precursor, but not gp120, coprecipitated with CNX, whereas only 0.25% of gp160 and gp120 coprecipitated with PDI. Under in vitro conditions, which mimic the Env/PDI interaction during virus/cell contact, PDI readily associated with Env. The domains of Env interacting in cellulo with CNX or in vitro with PDI were then determined using anti-Env antibodies whose binding site was occluded by CNX or PDI. Antibodies against domains V1/V2, C2, and the C terminus of V3 did not bind CNX-associated Env, whereas those against C1, V1/V2, and the CD4-binding domain did not react with PDI-associated Env. In addition, a mixture of the latter antibodies interfered with PDI-mediated Env reduction. Thus, Env interacts with intracellular CNX and extracellular PDI via discrete, largely nonoverlapping, regions. The sites of interaction explain the mode of action of compounds that target these two catalysts and may enable the design of further new competitive agents.

GTP binding controls complex formation by the human ROCO protein MASL1

The human ROCO proteins are a family of multi-domain proteins sharing a conserved ROC-COR supra-domain. The family has four members: leu- cine-rich repeat kinase 1 (LRRK1), leucine-rich repeat kinase 2 (LRRK2), death-associated protein kinase 1 (DAPK1) and malignant fibrous histiocy- toma amplified sequences with leucine-rich tandem repeats 1 (MASL1). Previous studies of LRRK1/2 and DAPK1 have shown that the ROC (Ras of complex proteins) domain can bind and hydrolyse GTP, but the cellular consequences of this activity are still unclear. Here, the first biochemical characterization of MASL1 and the impact of GTP binding on MASL1 complex formation are reported. The results demonstrate that MASL1, similar to other ROCO proteins, can bind guanosine nucleotides via its ROC domain. Furthermore, MASL1 exists in two distinct cellular com- plexes associated with heat shock protein 60, and the formation of a low molecular weight pool of MASL1 is modulated by GTP binding. Finally, loss of GTP enhances MASL1 toxicity in cells. Taken together, these data point to a central role for the ROC/GTPase domain of MASL1 in the reg- ulation of its cellular function.

Regulation of mitogen-activated protein kinases in cardiac myocytes through the small G protein Rac1.

Small guanine nucleotide-binding proteins of the Ras and Rho (Rac, Cdc42, and Rho) families have been implicated in cardiac myocyte hypertrophy, and this may involve the extracellular signal-related kinase (ERK), c-Jun N-terminal kinase (JNK), and/or p38 mitogen-activated protein kinase (MAPK) cascades. In other systems, Rac and Cdc42 have been particularly implicated in the activation of JNKs and p38-MAPKs. We examined the activation of Rho family small G proteins and the regulation of MAPKs through Rac1 in cardiac myocytes. Endothelin 1 and phenylephrine (both hypertrophic agonists) induced rapid activation of endogenous Rac1, and endothelin 1 also promoted significant activation of RhoA. Toxin B (which inactivates Rho family proteins) attenuated the activation of JNKs by hyperosmotic shock or endothelin 1 but had no effect on p38-MAPK activation. Toxin B also inhibited the activation of the ERK cascade by these stimuli. In transfection experiments, dominant-negative N17Rac1 inhibited activation of ERK by endothelin 1, whereas activated V12Rac1 cooperated with c-Raf to activate ERK. Rac1 may stimulate the ERK cascade either by promoting the phosphorylation of c-Raf or by increasing MEK1 and/or -2 association with c-Raf to facilitate MEK1 and/or -2 activation. In cardiac myocytes, toxin B attenuated c-Raf(Ser-338) phosphorylation (50 to 70% inhibition), but this had no effect on c-Raf activity. However, toxin B decreased both the association of MEK1 and/or -2 with c-Raf and c-Raf-associated ERK-activating activity. V12Rac1 cooperated with c-Raf to increase expression of atrial natriuretic factor (ANF), whereas N17Rac1 inhibited endothelin 1-stimulated ANF expression, indicating that the synergy between Rac1 and c-Raf is potentially physiologically important. We conclude that activation of Rac1 by hypertrophic stimuli contributes to the hypertrophic response by modulating the ERK and/or possibly the JNK (but not the p38-MAPK) cascades.

The Xanthomonas citri effector protein PthA interacts with citrus proteins involved in nuclear transport, protein folding and ubiquitination associated with DNA repair

P>Xanthomonas axonopodis pv. citri utilizes the type III effector protein PthA to modulate host transcription to promote citrus canker. PthA proteins belong to the AvrBs3/PthA family and carry a domain comprising tandem repeats of 34 amino acids that mediates protein-protein and protein-DNA interactions. We show here that variants of PthAs from a single bacterial strain localize to the nucleus of plant cells and form homo- and heterodimers through the association of their repeat regions. We hypothesize that the PthA variants might also interact with distinct host targets. Here, in addition to the interaction with alpha-importin, known to mediate the nuclear import of AvrBs3, we describe new interactions of PthAs with citrus proteins involved in protein folding and K63-linked ubiquitination. PthAs 2 and 3 preferentially interact with a citrus cyclophilin (Cyp) and with TDX, a tetratricopeptide domain-containing thioredoxin. In addition, PthAs 2 and 3, but not 1 and 4, interact with the ubiquitin-conjugating enzyme complex formed by Ubc13 and ubiquitin-conjugating enzyme variant (Uev), required for K63-linked ubiquitination and DNA repair. We show that Cyp, TDX and Uev interact with each other, and that Cyp and Uev localize to the nucleus of plant cells. Furthermore, the citrus Ubc13 and Uev proteins complement the DNA repair phenotype of the yeast Delta ubc13 and Delta mms2/uev1a mutants, strongly indicating that they are also involved in K63-linked ubiquitination and DNA repair. Notably, PthA 2 affects the growth of yeast cells in the presence of a DNA damage agent, suggesting that it inhibits K63-linked ubiquitination required for DNA repair.

PILZ Protein Structure and Interactions with PILB and the FIMX EAL Domain: Implications for Control of Type IV Pilus Biogenesis

The PilZ protein was originally identified as necessary for type IV pilus (T4P) biogenesis. Since then, a large and diverse family of bacterial PilZ homology domains have been identified, some of which have been implicated in signaling pathways that control important processes, including motility, virulence and biofilm formation. Furthermore, many PilZ homology domains, though not PilZ itself, have been shown to bind the important bacterial second messenger bis(3`-> 5`)cyclic diGMP (c-diGMP). The crystal structures of the PilZ orthologs from Xanthomonas axonopodis pv Citri (PilZ(XAC1133), this work) and from Xanthomonas campestris pv campestris (XC1028) present significant structural differences to other PilZ homologs that explain its failure to bind c-diGMP. NMR analysis of PilZ(XAC1133) shows that these structural differences are maintained in solution. In spite of their emerging importance in bacterial signaling, the means by which NZ proteins regulate specific processes is not clear. In this study, we show that PilZ(XAC1133) binds to PilB, an ATPase required for TV polymerization, and to the EAL domain of FiMX(XAC2398), which regulates TV biogenesis and localization in other bacterial species. These interactions were confirmed in NMR, two-hybrid and far-Western blot assays and are the first interactions observed between any PilZ domain and a target protein. While we were unable to detect phosphodiesterase activity for FimXX(AC2398) in vitro, we show that it binds c-diGMP both in the presence and in the absence of PilZ(XAC1133). Site-directed mutagenesis studies for conserved and exposed residues suggest that PilZ(XAC1133) interactions with FimX(XAC2398) and PilB(XAC3239) are mediated through a hydrophobic surface and an unstructured C-terminal extension conserved only in PilZ orthologs. The FimX-PilZ-PilB interactions involve a full set of ""degenerate"" GGDEF, EAL and PilZ domains and provide the first evidence of the means by which PilZ orthologs and FimX interact directly with the TP4 machinery. (C) 2009 Elsevier Ltd. All rights reserved.

The Saccharomyces cerevisiae COQ10 gene encodes a START domain protein required for function of coenzyme Q in respiration

Deletion of the Saccharomyces cerevisiae gene YOL008W, here referred to as COQ10, elicits a respiratory defect as a result of the inability of the mutant to oxidize NADH and succinate. Both activities are restored by exogenous coenzyme Q(2). Respiration is also partially rescued by COQ2, COQ7, or COQ8/ABC1, when these genes are present in high copy. Unlike other coq mutants, all of which lack Q(6), the coq10 mutant has near normal amounts of Q(6) in mitochondria. Coq10p is widely distributed in bacteria and eukaryotes and is homologous to proteins of the aromatic-rich protein family Pfam03654 and to members of the START domain superfamily that have a hydrophobic tunnel implicated in binding lipophilic molecules such as cholesterol and polyketides. Analysis of coenzyme Q in polyhistidine-tagged Coq10p purified from mitochondria indicates the presence 0.032-0.034 mol of Q(6)/mol of protein. We propose that Coq10p is a Q(6)-binding protein and that in the coq10 mutant Q(6) it is not able to act as an electron carrier, possibly because of improper localization.

KRAS gene mutations in Noonan syndrome familial cases cluster in the vicinity of the switch II region of the G-domain: Report of another family with metopic craniosynostosis

Noonan syndrome (NS) and Noonan-related disorders [cardio-facio-cutaneous (CFC), Costello, Noonan syndrome with multiple lentigines (NS-ML), and neurofibromatosis-Noonan syndromes (NFNS)] are a group of developmental disorders caused by mutations in genes of the RAS/MAPK pathway. Mutations in the KRAS gene account for only a small proportion of affected Noonan and CFC syndrome patients that present an intermediate phenotype between these two syndromes, with more frequent and severe intellectual disability in NS and less ectodermal involvement in CFC syndrome, as well as atypical clinical findings such as craniosynostosis. Recently, the first familial case with a novel KRAS mutation was described. We report on a second vertical transmission (a mother and two siblings) with a novel mutation (p.M72L), in which the proband has trigonocephaly and the affected mother and sister, prominent ectodermal involvement. Metopic suture involvement has not been described before, expanding the main different cranial sutures which can be affected in NS and KRAS gene mutations. The gene alteration found in the studied family is in close proximity to the one reported in the other familial case (close to the switch II region of the G-domain), suggesting that this specific region of the gene could have less severe effects on intellectual ability than the other KRAS gene mutations found in NS patients and be less likely to hamper reproductive fitness. (c) 2012 Wiley Periodicals, Inc.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

BPAG1 isoform-b: Complex distribution pattern in striated and heart muscle and association with plectin and α-actinin

BPAG1-b is the major muscle-specific isoform encoded by the dystonin gene, which expresses various protein isoforms belonging to the plakin protein family with complex, tissue-specific expression profiles. Recent observations in mice with either engineered or spontaneous mutations in the dystonin gene indicate that BPAG1-b serves as a cytolinker important for the establishment and maintenance of the cytoarchitecture and integrity of striated muscle. Here, we studied in detail its distribution in skeletal and cardiac muscles and assessed potential binding partners. BPAG1-b was detectable in vitro and in vivo as a high molecular mass protein in striated and heart muscle cells, co-localizing with the sarcomeric Z-disc protein alpha-actinin-2 and partially with the cytolinker plectin as well as with the intermediate filament protein desmin. Ultrastructurally, like alpha-actinin-2, BPAG1-b was predominantly localized at the Z-discs, adjacent to desmin-containing structures. BPAG1-b was able to form complexes with both plectin and alpha-actinin-2, and its NH(2)-terminus, which contains an actin-binding domain, directly interacted with that of plectin and alpha-actinin. Moreover, the protein level of BPAG1-b was reduced in muscle tissues from plectin-null mutant mice versus wild-type mice. These studies provide new insights into the role of BPAG1-b in the cytoskeletal organization of striated muscle.

Alba-domain proteins of Trypanosoma brucei are cytoplasmic RNA-binding proteins that interact with the translation machinery

Trypanosoma brucei and related pathogens transcribe most genes as polycistronic arrays that are subsequently processed into monocistronic mRNAs. Expression is frequently regulated post-transcriptionally by cis-acting elements in the untranslated regions (UTRs). GPEET and EP procyclins are the major surface proteins of procyclic (insect midgut) forms of T. brucei. Three regulatory elements common to the 3' UTRs of both mRNAs regulate mRNA turnover and translation. The glycerol-responsive element (GRE) is unique to the GPEET 3' UTR and regulates its expression independently from EP. A synthetic RNA encompassing the GRE showed robust sequence-specific interactions with cytoplasmic proteins in electromobility shift assays. This, combined with column chromatography, led to the identification of 3 Alba-domain proteins. RNAi against Alba3 caused a growth phenotype and reduced the levels of Alba1 and Alba2 proteins, indicative of interactions between family members. Tandem-affinity purification and co-immunoprecipitation verified these interactions and also identified Alba4 in sub-stoichiometric amounts. Alba proteins are cytoplasmic and are recruited to starvation granules together with poly(A) RNA. Concomitant depletion of all four Alba proteins by RNAi specifically reduced translation of a reporter transcript flanked by the GPEET 3' UTR. Pulldown of tagged Alba proteins confirmed interactions with poly(A) binding proteins, ribosomal protein P0 and, in the case of Alba3, the cap-binding protein eIF4E4. In addition, Alba2 and Alba3 partially cosediment with polyribosomes in sucrose gradients. Alba-domain proteins seem to have exhibited great functional plasticity in the course of evolution. First identified as DNA-binding proteins in Archaea, then in association with nuclear RNase MRP/P in yeast and mammalian cells, they were recently described as components of a translationally silent complex containing stage-regulated mRNAs in Plasmodium. Our results are also consistent with stage-specific regulation of translation in trypanosomes, but most likely in the context of initiation.

Plectin interacts with the rod domain of type III intermediate filament proteins desmin and vimentin

Plectin is a versatile cytolinker protein critically involved in the organization of the cytoskeletal filamentous system. The muscle-specific intermediate filament (IF) protein desmin, which progressively replaces vimentin during differentiation of myoblasts, is one of the important binding partners of plectin in mature muscle. Defects of either plectin or desmin cause muscular dystrophies. By cell transfection studies, yeast two-hybrid, overlay and pull-down assays for binding analysis, we have characterized the functionally important sequences for the interaction of plectin with desmin and vimentin. The association of plectin with both desmin and vimentin predominantly depended on its fifth plakin repeat domain and downstream linker region. Conversely, the interaction of desmin and vimentin with plectin required sequences contained within the segments 1A-2A of their central coiled-coil rod domain. This study furthers our knowledge of the interaction between plectin and IF proteins important for maintenance of cytoarchitecture in skeletal muscle. Moreover, binding of plectin to the conserved rod domain of IF proteins could well explain its broad interaction with most types of IFs.