Comparative structural and energetic analysis of WW domain-peptide interactions

Systematic profiling and identification of the peptide-mediated interactions between human Yes-associated protein and its partners in esophageal cancer

Human Yes-associated protein (YAP) is involved in the Hippo signaling pathway and serves as a coactivator to modulate gene expression, which contains a transactivation domain (TD) responsible for binding to the downstream TEA domain family (TEAD) of transcription factors and two WW1/2 domains that recognize the proline-rich motifs (PRMs) present in a variety of upstream protein partners through peptide-mediated interactions (PMIs). The downstream YAP TD-TEAD interactions are closely associated with gastric cancer, and a number of therapeutic agents have been developed to target the interactions. In contrast, the upstream YAP WW1/2-partner interactions are thought to be involved in esophageal cancer but still remain largely unexplored. Here, we attempted to elucidate the complicated PMIs between the YAP WW1/2 domains and various PRMs of YAP-interacting proteins. A total of 106 peptide segments carrying the class I WW-binding motif [P/L]Px[Y/P] were extracted from 22 partner candidates, which are potential recognition sites of YAP WW1/2 domains. Structural and energetic analyses of the intermolecular interactions between the domains and peptides created a systematic domain-peptide binding profile, from which a number of biologically functional PMIs were identified and then substantiated in vitro using fluorescence spectroscopy assays. It is revealed that: (a) The sequence requirement for the partner recognition site binding to YAP WW1/2 domains is a decapeptide segment that contains a core PRM motif as well as two three-residue extensions from each side of the motif; the core motif and extended sections are responsible for the binding stability and recognition specificity of domain-peptide interaction, respectively. (b) There is an exquisite difference in the recognition specificity of the two domains; the LPxP and PPxP appear to more prefer WW1 than WW2, whereas the WW2 can bind more effectively to LPxY and PPxY than WW1. (c) WW2 generally exhibits a higher affinity to the panel of recognition site candidates than WW1. In addition, a number of partner peptides were found as promising recognition sites of the two domains and/or to have a good selectivity between the two domains. For example, the DVL1 peptide was determined to have moderate affinity to WW2 and strong selectivity for WW2 over WW1. Hydrogen bonds play a central role in selectivity.

Kibra Modulates Learning and Memory via Binding to Dendrin

Kibra is a synaptic scaffold protein regulating learning and memory. Alterations of Kibra-encoding gene WWC1 cause various neuronal disorders, including Alzheimer's disease and Tourette syndrome. However, the molecular mechanism underlying Kibra's function in neurons is poorly understood. Here we discover that Kibra, via its N-terminal WW12 tandem domains, binds to a postsynaptic density enriched protein, Dendrin, with a nanomolar dissociation constant. On the basis of the structure of Kibra WW12 in complex with Dendrin PY motifs, we developed a potent peptide inhibitor capable of specifically blocking the binding between Kibra and Dendrin in neurons. Systematic administration of the inhibitory peptide attenuated excitatory synaptic transmission, completely blocked long-term potentiation induction, and impaired spatial learning and memory. A Kibra mutation found in Tourette syndrome patients causes defects in binding to Dendrin. Thus, Kibra can modulate spatial learning and memory via binding to Dendrin.

Biophysical studies and modelling indicate the binding preference of TAZ WW domain for LATS1 PPxY motif

The Hippo tumor suppressor pathway is an important regulator of cell proliferation and apoptosis, and signal transduction occurs through phosphorylation of the effector protein TAZ by the serine/threonine kinase LATS1/2. Here, we report the biophysical and computational studies to characterize the interaction between TAZ and LATS1/2 through WW domain-PPxY motif binding. We show that the TAZ WW domain exhibits a binding preference for the second of the two PPxY motifs of LATS1 in vitro. We modelled the structure of the domain in complex with LATS1 PPxY2 peptide and, through molecular dynamics simulations, show that WW domain-PPxY2 complex is stable with some flexibility in the peptide region. Next, we predict and verify that L143 and T150 of the WW domain are important for TAZ binding with the PPxY2 peptide using mutational and isothermal titration calorimetric studies. Furthermore, we suggest that the electrostatic potential of charged residues within the binding pocket may influence the ligand affinity among otherwise highly similar WW domains.

Impact of the green tea ingredient epigallocatechin gallate and a short pentapeptide (Ile-Ile-Ala-Glu-Lys) on the structural organization of mixed micelles and the related uptake of cholesterol

BACKGROUND

High levels of blood cholesterol are conventionally linked to an increased risk of developing cardiovascular disease (Grundy, 1986). Here we examine the molecular mode of action of natural products with known cholesterol-lowering activity, such as for example the green tea ingredient epigallocatechin gallate and a short pentapeptide, Ile-Ile-Ala-Glu-Lys.

METHODS

Molecular Dynamics simulations are used to gain insight into the formation process of mixed micelles and, correspondingly, how active agents epigallocatechin gallate and Ile-Ile-Ala-Glu-Lys could possibly interfere with it.

RESULTS

Self-assembly of physiological micelles occurs on the order of 35-50 ns; most of the structural properties of mixed micelles are unaffected by epigallocatechin gallate or Ile-Ile-Ala-Glu-Lys which integrate into the micellar surface; the diffusive motion of constituting lipids palmitoyl-oleoyl-phosphatidylcholine and cholesterol is significantly down-regulated by both epigallocatechin gallate and Ile-Ile-Ala-Glu-Lys; CONCLUSIONS: The molecular mode of action of natural compounds epigallocatechin gallate and Ile-Ile-Ala-Glu-Lys is a significant down-regulation of the diffusive motion of micellar lipids.

GENERAL SIGNIFICANCE

Natural compounds like the green tea ingredient epigallocatechin gallate and a short pentapeptide, Ile-Ile-Ala-Glu-Lys, lead to a significant down-regulation of the diffusive motion of micellar lipids thereby modulating cholesterol absorption into physiological micelles.

WW domains of the yes-kinase-associated-protein (YAP) transcriptional regulator behave as independent units with different binding preferences for PPxY motif-containing ligands

YAP is a WW domain-containing effector of the Hippo tumor suppressor pathway, and the object of heightened interest as a potent oncogene and stemness factor. YAP has two major isoforms that differ in the number of WW domains they harbor. Elucidating the degree of co-operation between these WW domains is important for a full understanding of the molecular function of YAP. We present here a detailed biophysical study of the structural stability and binding properties of the two YAP WW domains aimed at investigating the relationship between both domains in terms of structural stability and partner recognition. We have carried out a calorimetric study of the structural stability of the two YAP WW domains, both isolated and in a tandem configuration, and their interaction with a set of functionally relevant ligands derived from PTCH1 and LATS kinases. We find that the two YAP WW domains behave as independent units with different binding preferences, suggesting that the presence of the second WW domain might contribute to modulate target recognition between the two YAP isoforms. Analysis of structural models and phage-display studies indicate that electrostatic interactions play a critical role in binding specificity. Together, these results are relevant to understand of YAP function and open the door to the design of highly specific ligands of interest to delineate the functional role of each WW domain in YAP signaling.

Invited review: Caseins and the casein micelle: their biological functions, structures, and behavior in foods

A typical casein micelle contains thousands of casein molecules, most of which form thermodynamically stable complexes with nanoclusters of amorphous calcium phosphate. Like many other unfolded proteins, caseins have an actual or potential tendency to assemble into toxic amyloid fibrils, particularly at the high concentrations found in milk. Fibrils do not form in milk because an alternative aggregation pathway is followed that results in formation of the casein micelle. As a result of forming micelles, nutritious milk can be secreted and stored without causing either pathological calcification or amyloidosis of the mother's mammary tissue. The ability to sequester nanoclusters of amorphous calcium phosphate in a stable complex is not unique to caseins. It has been demonstrated using a number of noncasein secreted phosphoproteins and may be of general physiological importance in preventing calcification of other biofluids and soft tissues. Thus, competent noncasein phosphoproteins have similar patterns of phosphorylation and the same type of flexible, unfolded conformation as caseins. The ability to suppress amyloid fibril formation by forming an alternative amorphous aggregate is also not unique to caseins and underlies the action of molecular chaperones such as the small heat-shock proteins. The open structure of the protein matrix of casein micelles is fragile and easily perturbed by changes in its environment. Perturbations can cause the polypeptide chains to segregate into regions of greater and lesser density. As a result, the reliable determination of the native structure of casein micelles continues to be extremely challenging. The biological functions of caseins, such as their chaperone activity, are determined by their composition and flexible conformation and by how the casein polypeptide chains interact with each other. These same properties determine how caseins behave in the manufacture of many dairy products and how they can be used as functional ingredients in other foods.

The interaction properties of the human Rab GTPase family--comparative analysis reveals determinants of molecular binding selectivity

Background

Rab GTPases constitute the largest subfamily of the Ras protein superfamily. Rab proteins regulate organelle biogenesis and transport, and display distinct binding preferences for effector and activator proteins, many of which have not been elucidated yet. The underlying molecular recognition motifs, binding partner preferences and selectivities are not well understood.

Methodology/Principal Findings

Comparative analysis of the amino acid sequences and the three-dimensional electrostatic and hydrophobic molecular interaction fields of 62 human Rab proteins revealed a wide range of binding properties with large differences between some Rab proteins. This analysis assists the functional annotation of Rab proteins 12, 14, 26, 37 and 41 and provided an explanation for the shared function of Rab3 and 27. Rab7a and 7b have very different electrostatic potentials, indicating that they may bind to different effector proteins and thus, exert different functions. The subfamily V Rab GTPases which are associated with endosome differ subtly in the interaction properties of their switch regions, and this may explain exchange factor specificity and exchange kinetics.

Conclusions/Significance

We have analysed conservation of sequence and of molecular interaction fields to cluster and annotate the human Rab proteins. The analysis of three dimensional molecular interaction fields provides detailed insight that is not available from a sequence-based approach alone. Based on our results, we predict novel functions for some Rab proteins and provide insights into their divergent functions and the determinants of their binding partner selectivity.

The emerging contribution of sequence context to the specificity of protein interactions mediated by PDZ domains

The canonical binding mode of PDZ domains to target motifs involves a small interface, unlikely to fully account for PDZ-target interaction specificities. Here, we review recent work on sequence context, defined as the regions surrounding not only the PDZ domains but also their target motifs. We also address the theoretical problem of defining the core of PDZ domains and the practical issue of designing PDZ constructs. Sequence context is found to introduce structural diversity, to impact the stability and solubility of constructs, and to deeply influence binding affinity and specificity, thereby increasing the difficulty of predicting PDZ-motif interactions. We expect that sequence context will have similar importance for other protein interactions mediated by globular domains binding to short linear motifs.

Clustering of HIV-1 Subtypes Based on gp120 V3 Loop electrostatic properties

BACKGROUND

The V3 loop of the glycoprotein gp120 of HIV-1 plays an important role in viral entry into cells by utilizing as coreceptor CCR5 or CXCR4, and is implicated in the phenotypic tropisms of HIV viruses. It has been hypothesized that the interaction between the V3 loop and CCR5 or CXCR4 is mediated by electrostatics. We have performed hierarchical clustering analysis of the spatial distributions of electrostatic potentials and charges of V3 loop structures containing consensus sequences of HIV-1 subtypes.

RESULTS

Although the majority of consensus sequences have a net charge of +3, the spatial distribution of their electrostatic potentials and charges may be a discriminating factor for binding and infectivity. This is demonstrated by the formation of several small subclusters, within major clusters, which indicates common origin but distinct spatial details of electrostatic properties. Some of this information may be present, in a coarse manner, in clustering of sequences, but the spatial details are largely lost. We show the effect of ionic strength on clustering of electrostatic potentials, information that is not present in clustering of charges or sequences. We also make correlations between clustering of electrostatic potentials and net charge, coreceptor selectivity, global prevalence, and geographic distribution. Finally, we interpret coreceptor selectivity based on the N6X7T8|S8X9 sequence glycosylation motif, the specific positive charge location according to the 11/24/25 rule, and the overall charge and electrostatic potential distribution.

CONCLUSIONS

We propose that in addition to the sequence and the net charge of the V3 loop of each subtype, the spatial distributions of electrostatic potentials and charges may also be important factors for receptor recognition and binding and subsequent viral entry into cells. This implies that the overall electrostatic potential is responsible for long-range recognition of the V3 loop with coreceptors CCR5/CXCR4, whereas the charge distribution contributes to the specific short-range interactions responsible for the formation of the bound complex. We also propose a scheme for coreceptor selectivity based on the sequence glycosylation motif, the 11/24/25 rule, and net charge.

Prediction of peptides binding to the PKA RIIalpha subunit using a hierarchical strategy

MOTIVATION

Favorable interaction between the regulatory subunit of the cAMP-dependent protein kinase (PKA) and a peptide in A-kinase anchoring proteins (AKAPs) is critical for translocating PKA to the subcellular sites where the enzyme phosphorylates its substrates. It is very hard to identify AKAPs peptides binding to PKA due to the high sequence diversity of AKAPs.

RESULTS

We propose a hierarchical and efficient approach, which combines molecular dynamics (MD) simulations, free energy calculations, virtual mutagenesis (VM) and bioinformatics analyses, to predict peptides binding to the PKA RIIα regulatory subunit in the human proteome systematically. Our approach successfully retrieved 15 out of 18 documented RIIα-binding peptides. Literature curation supported that many newly predicted peptides might be true AKAPs. Here, we present the first systematic search for AKAP peptides in the human proteome, which is useful to further experimental identification of AKAPs and functional analysis of their biological roles.

The effect of electrostatics on factor H function and related pathologies

Factor H (FH) contributes to the regulation of the complement system by binding to polyanionic surfaces and the proteins C3b/C3c/C3d. This implicates charge and electrostatic interactions in recognition and binding of FH. Despite the large amount of experimental and pathology data the exact mechanism at molecular level is not yet known. We have implemented a computational framework for comparative analysis of the charge and electrostatic diversity of FH modules and C3b domains to identify electrostatic hotspots and predict potential binding sites. Our electrostatic potential clustering analysis shows that charge distributions and electrostatic potential distributions are more useful in understanding C3b-FH interactions than net charges alone. We present a model of non-specific electrostatic interactions of FH with polyanion-rich surfaces and specific interactions with C3b, using our computational data and existing experimental data. We discuss the electrostatic contributions to the formation of the C3b-FH complex and the competition between FH and Factor Bb (Bb) for binding to C3b. We also discuss the significance of mutations of charged amino acids in the pathobiology of FH-mediated disease, such as age-related macular degeneration, atypical hemolytic uremic syndrome, and dense deposit disease. Our data can be used to guide future experimental studies.

A regression framework incorporating quantitative and negative interaction data improves quantitative prediction of PDZ domain-peptide interaction from primary sequence

Motivation: Predicting protein interactions involving peptide recognition domains is essential for understanding the many important biological processes they mediate. It is important to consider the binding strength of these interactions to help us construct more biologically relevant protein interaction networks that consider cellular context and competition between potential binders.

Results: We developed a novel regression framework that considers both positive (quantitative) and negative (qualitative) interaction data available for mouse PDZ domains to quantitatively predict interactions between PDZ domains, a large peptide recognition domain family, and their peptide ligands using primary sequence information. First, we show that it is possible to learn from existing quantitative and negative interaction data to infer the relative binding strength of interactions involving previously unseen PDZ domains and/or peptides given their primary sequence. Performance was measured using cross-validated hold out testing and testing with previously unseen PDZ domain–peptide interactions. Second, we find that incorporating negative data improves quantitative interaction prediction. Third, we show that sequence similarity is an important prediction performance determinant, which suggests that experimentally collecting additional quantitative interaction data for underrepresented PDZ domain subfamilies will improve prediction.

Availability and Implementation: The Matlab code for our SemiSVR predictor and all data used here are available at http://baderlab.org/Data/PDZAffinity.

Contact:gary.bader@utoronto.ca; dengnaiyang@cau.edu.cn

Supplementary information:Supplementary data are available at Bioinformatics online.

Using genome-wide measurements for computational prediction of SH2-peptide interactions

Peptide-recognition modules (PRMs) are used throughout biology to mediate protein–protein interactions, and many PRMs are members of large protein domain families. Recent genome-wide measurements describe networks of peptide–PRM interactions. In these networks, very similar PRMs recognize distinct sets of peptides, raising the question of how peptide-recognition specificity is achieved using similar protein domains. The analysis of individual protein complex structures often gives answers that are not easily applicable to other members of the same PRM family. Bioinformatics-based approaches, one the other hand, may be difficult to interpret physically. Here we integrate structural information with a large, quantitative data set of SH2 domain–peptide interactions to study the physical origin of domain–peptide specificity. We develop an energy model, inspired by protein folding, based on interactions between the amino-acid positions in the domain and peptide. We use this model to successfully predict which SH2 domains and peptides interact and uncover the positions in each that are important for specificity. The energy model is general enough that it can be applied to other members of the SH2 family or to new peptides, and the cross-validation results suggest that these energy calculations will be useful for predicting binding interactions. It can also be adapted to study other PRM families, predict optimal peptides for a given SH2 domain, or study other biological interactions, e.g. protein–DNA interactions.

Analysis of electrostatic contributions to the selectivity of interactions between RING-finger domains and ubiquitin-conjugating enzymes

The zinc-coordinated protein motifs known as RING-finger domains, present on a class of ubiquitin ligases (E3's), recruit ubiquitin-conjugating enzymes (E2s), tethering them to substrate proteins for covalent modification with ubiquitin. Each RING-finger domain can recruit different E2s, and these interactions are frequently selective, in that certain RING-finger domains associate preferentially with certain E2s. This selectivity acquires particular biological relevance when the recruited E2s exert specialized functions. We have explored the determinants that specify the presence or absence of experimentally detectable interaction between two RING-finger domains, those on RNF11 and RNF103, and two E2s, UBC13, a specialized E2 that catalyzes ubiquitin chain elongation through Lys63 of ubiquitin, and UbcH7, which mediates polyubiquitylation through Lys48. Through the iterative use of computational predictive tools and experimental validations, we have found that these interactions and their selectivity are partly governed by the combinations of electrostatic interactions linking specific residues of the contact interfaces. Our analysis also predicts that the main determinants of selectivity of these interactions reside on the RING-finger domains, rather than on the E2s. The application of some of these rules of interaction selectivity has permitted us to experimentally manipulate the selectivity of interaction of the RING-finger domain-E2 pairs under study.

Evolutionary and biophysical relationships among the papillomavirus E2 proteins

Infection by human papillomavirus (HPV) may result in clinical conditions ranging from benign warts to invasive cancer. The HPV E2 protein represses oncoprotein transcription and is required for viral replication. HPV E2 binds to palindromic DNA sequences of highly conserved four base pair sequences flanking an identical length variable 'spacer'. E2 proteins directly contact the conserved but not the spacer DNA. Variation in naturally occurring spacer sequences results in differential protein affinity that is dependent on their sensitivity to the spacer DNA's unique conformational and/or dynamic properties. This article explores the biophysical character of this core viral protein with the goal of identifying characteristics that associated with risk of virally caused malignancy. The amino acid sequence, 3d structure and electrostatic features of the E2 protein DNA binding domain are highly conserved; specific interactions with DNA binding sites have also been conserved. In contrast, the E2 protein's transactivation domain does not have extensive surfaces of highly conserved residues. Rather, regions of high conservation are localized to small surface patches. Implications to cancer biology are discussed.

Rapid comparison of properties on protein surface

The mapping of physicochemical characteristics onto the surface of a protein provides crucial insights into its function and evolution. This information can be further used in the characterization and identification of similarities within protein surface regions. We propose a novel method which quantitatively compares global and local properties on the protein surface. We have tested the method on comparison of electrostatic potentials and hydrophobicity. The method is based on 3D Zernike descriptors, which provides a compact representation of a given property defined on a protein surface. Compactness and rotational invariance of this descriptor enable fast comparison suitable for database searches. The usefulness of this method is exemplified by studying several protein families including globins, thermophilic and mesophilic proteins, and active sites of TIM beta/alpha barrel proteins. In all the cases studied, the descriptor is able to cluster proteins into functionally relevant groups. The proposed approach can also be easily extended to other surface properties. This protein surface-based approach will add a new way of viewing and comparing proteins to conventional methods, which compare proteins in terms of their primary sequence or tertiary structure.

Electrostatic similarity of proteins: application of three dimensional spherical harmonic decomposition

We present a method for describing and comparing global electrostatic properties of biomolecules based on the spherical harmonic decomposition of electrostatic potential data. Unlike other approaches our method does not require any prior three dimensional structural alignment. The electrostatic potential, given as a volumetric data set from a numerical solution of the Poisson or Poisson-Boltzmann equation, is represented with descriptors that are rotation invariant. The method can be applied to large and structurally diverse sets of biomolecules enabling to cluster them according to their electrostatic features.

webPIPSA: a web server for the comparison of protein interaction properties

Protein molecular interaction fields are key determinants of protein functionality. PIPSA (Protein Interaction Property Similarity Analysis) is a procedure to compare and analyze protein molecular interaction fields, such as the electrostatic potential. PIPSA may assist in protein functional assignment, classification of proteins, the comparison of binding properties and the estimation of enzyme kinetic parameters. webPIPSA is a web server that enables the use of PIPSA to compare and analyze protein electrostatic potentials. While PIPSA can be run with downloadable software (see http://projects.eml.org/mcm/software/pipsa), webPIPSA extends and simplifies a PIPSA run. This allows non-expert users to perform PIPSA for their protein datasets. With input protein coordinates, the superposition of protein structures, as well as the computation and analysis of electrostatic potentials, is automated. The results are provided as electrostatic similarity matrices from an all-pairwise comparison of the proteins which can be subjected to clustering and visualized as epograms (tree-like diagrams showing electrostatic potential differences) or heat maps. webPIPSA is freely available at: http://pipsa.eml.org.

qPIPSA: relating enzymatic kinetic parameters and interaction fields

Background

The simulation of metabolic networks in quantitative systems biology requires the assignment of enzymatic kinetic parameters. Experimentally determined values are often not available and therefore computational methods to estimate these parameters are needed. It is possible to use the three-dimensional structure of an enzyme to perform simulations of a reaction and derive kinetic parameters. However, this is computationally demanding and requires detailed knowledge of the enzyme mechanism. We have therefore sought to develop a general, simple and computationally efficient procedure to relate protein structural information to enzymatic kinetic parameters that allows consistency between the kinetic and structural information to be checked and estimation of kinetic constants for structurally and mechanistically similar enzymes.

Results

We describe qPIPSA: quantitative Protein Interaction Property Similarity Analysis. In this analysis, molecular interaction fields, for example, electrostatic potentials, are computed from the enzyme structures. Differences in molecular interaction fields between enzymes are then related to the ratios of their kinetic parameters. This procedure can be used to estimate unknown kinetic parameters when enzyme structural information is available and kinetic parameters have been measured for related enzymes or were obtained under different conditions. The detailed interaction of the enzyme with substrate or cofactors is not modeled and is assumed to be similar for all the proteins compared. The protein structure modeling protocol employed ensures that differences between models reflect genuine differences between the protein sequences, rather than random fluctuations in protein structure.

Conclusion

Provided that the experimental conditions and the protein structural models refer to the same protein state or conformation, correlations between interaction fields and kinetic parameters can be established for sets of related enzymes. Outliers may arise due to variation in the importance of different contributions to the kinetic parameters, such as protein stability and conformational changes. The qPIPSA approach can assist in the validation as well as estimation of kinetic parameters, and provide insights into enzyme mechanism.

ArhGAP9, a novel MAP kinase docking protein, inhibits Erk and p38 activation through WW domain binding

We have identified human ArhGAP9 as a novel MAP kinase docking protein that interacts with Erk2 and p38alpha through complementarily charged residues in the WW domain of ArhGAP9 and the CD domains of Erk2 and p38alpha. This interaction sequesters the MAP kinases in their inactive states through displacement of MAP kinase kinases targeting the same sites. While over-expression of wild type ArhGAP9 caused MAP kinase activation by the epidermal growth factor receptor (EGFR) to be suppressed and preserved the actin stress fibres in quiescent Swiss 3T3 fibroblasts, over-expression of an ArhGAP9 mutant defective in MAP kinase binding restored EGFR-induced MAP kinase activation and resulted in significant disruption of the stress fibres, consistent with the role of Erk activation in disassembly of actin stress fibres. The interaction between ArhGAP9 and the MAP kinases represents a novel mechanism of cross-talk between Rho GTPase and MAP kinase signaling.

2D-RNA-coupling numbers: A new computational chemistry approach to link secondary structure topology with biological function

Methods for prediction of proteins, DNA, or RNA function and mapping it onto sequence often rely on bioinformatics alignment approach instead of chemical structure. Consequently, it is interesting to develop computational chemistry approaches based on molecular descriptors. In this sense, many researchers used sequence-coupling numbers and our group extended them to 2D proteins representations. However, no coupling numbers have been reported for 2D-RNA topology graphs, which are highly branched and contain useful information. Here, we use a computational chemistry scheme: (a) transforming sequences into RNA secondary structures, (b) defining and calculating new 2D-RNA-coupling numbers, (c) seek a structure-function model, and (d) map biological function onto the folded RNA. We studied as example 1-aminocyclopropane-1-carboxylic acid (ACC) oxidases known as ACO, which control fruit ripening having importance for biotechnology industry. First, we calculated tau(k)(2D-RNA) values to a set of 90-folded RNAs, including 28 transcripts of ACO and control sequences. Afterwards, we compared the classification performance of 10 different classifiers implemented in the software WEKA. In particular, the logistic equation ACO = 23.8 . tau(1)(2D-RNA) + 41.4 predicts ACOs with 98.9%, 98.0%, and 97.8% of accuracy in training, leave-one-out and 10-fold cross-validation, respectively. Afterwards, with this equation we predict ACO function to a sequence isolated in this work from Coffea arabica (GenBank accession DQ218452). The tau(1)(2D-RNA) also favorably compare with other descriptors. This equation allows us to map the codification of ACO activity on different mRNA topology features. The present computational-chemistry approach is general and could be extended to connect RNA secondary structure topology to other functions.

An Expanded WW Domain Recognition Motif Revealed by the Interaction between Smad7 and the E3 Ubiquitin Ligase Smurf2

Smurf2 is an E3 ubiquitin ligase that drives degradation of the transforming growth factor-beta receptors and other targets. Recognition of the receptors by Smurf2 is accomplished through an intermediary protein, Smad7. Here we have demonstrated that the WW3 domain of Smurf2 can directly bind to the Smad7 polyproline-tyrosine (PY) motif. Of particular interest, the highly conserved WW domain binding site Trp, which interacts with target PY motifs, is a Phe in the Smurf2 WW3 domain. To examine this interaction, the solution structure of the complex between the Smad7 PY motif region (ELESPPPPYSRYPMD) and the Smurf2 WW3 domain was determined. The structure reveals that, in addition to binding the PY motif, the WW3 domain binds six residues C-terminal to the PY motif (PY-tail). Although the Phe in the WW3 domain binding site decreases affinity relative to the canonical Trp, this is balanced by additional interactions between the PY-tail and the beta1-strand and beta1-beta2 loop of the WW3 domain. The interaction between the Smurf2 WW3 domain and the Smad7 PY motif is the first example of PY motif recognition by a WW domain with a Phe substituted for the binding site Trp. This unusual interaction allows the Smurf2 WW3 domain to recognize a subset of PY motif-containing proteins utilizing an expanded surface to provide specificity.

Comparative analysis of enzyme activities and mRNA levels of peptidyl prolylcis/transisomerases in various organs of wild type andPin1−/−mice

We investigated the enzyme activity of peptidyl prolyl cis/trans isomerases (PPIases) in brain, testis, lung, liver, and mouse embryonic fibroblasts (MEF) of Pin1+/+ and Pin1-/- mice. The aim of this study is to determine if other PPIases can substitute for the loss of Pin1 activity in Pin1-/- mice and what influence Pin1 depletion has on the activities of other PPIases members. The results show that high PPIase activities of Pin1 are found in organs that have the tendency to develop Pin1 knockout phenotypes and, therefore, provide for the first time an enzymological basis for these observations. Furthermore we determined the specific activity (k(cat)/K(M)) of endogenous Pin1 and found that it is strongly reduced as compared with the recombinant protein in all investigated organs. These results suggest that posttranslational modifications may influence the PPIase activity in vivo. The activities originating from cyclophilin and FKBP are not influenced by the Pin1 knockout, but a basal enzymatic activity towards phosphorylated substrates could be found in Pin1-/- lysates. Real time PCR experiments of all PPIases in different mouse organs and MEF of Pin1+/+ and Pin1-/- mice support the finding and reveal the specific expression profiles of PPIases in mice.

Immunophysical properties and prediction of activities for vaccinia virus complement control protein and smallpox inhibitor of complement enzymes using molecular dynamics and electrostatics

We present immunophysical modeling for VCP, SPICE, and three mutants using MD simulations and Poisson-Boltzmann-type electrostatic calculations. VCP and SPICE are homologous viral proteins that control the complement system by imitating, structurally and functionally, natural regulators of complement activation. VCP and SPICE consist of four CCP modules connected with short flexible loops. MD simulations demonstrate that the rather complex modules of VCP/SPICE and their mutants exhibit a high degree of intermodular spatial mobility, which is affected by surface mutations. Electrostatic calculations using snapshots from the MD trajectories demonstrate variable spatial distribution of the electrostatic potentials, which suggests dynamic binding properties. We use covariance analysis to identify correlated modular oscillations. We also use electrostatic similarity indices to cluster proteins with common electrostatic properties. Our results are compared with experimental data to form correlations between the overall positive electrostatic potential of VCP/SPICE with binding and activity. We show how these correlations can be used to predict binding and activity properties. This work is expected to be useful for understanding the function of native CCP-containing regulators of complement activation and receptors and for the design of antiviral therapeutics and complement inhibitors.

Folding, misfolding, and amyloid protofibril formation of WW domain FBP28

We study the folding mechanism of a triple beta-strand WW domain from the Formin binding protein 28 (FBP28) at atomic resolution with explicit water model using replica exchange molecular dynamics computer simulations. Extended sampling over a wide range of temperatures to obtain the free energy, enthalpy, and entropy surfaces as a function of structural coordinates has been performed. Simulations were started from different configurations covering the folded and unfolded states. In the free energy landscape a transition state is identified and its structures and -values are compared with experimental data from a homologous protein, the prolyl-isomerase Pin1 WW domain. A stable intermediate state is found to accumulate during the simulation characterized by the carboxyl-terminal beta-strand 3 having misregistered hydrogen bonds and where the structural heterogeneity is due to nonnative turn II formation. Furthermore, the aggregation behavior of the FBP28 WW domain may be related to one such misfolded structure, which has a much lower free energy of dimer formation than that of the native dimer. Based on the misfolded dimer, aggregation to form protofibril structure is discussed.

Specificity and versatility of SH3 and other proline-recognition domains: structural basis and implications for cellular signal transduction

Protein-protein interactions occurring via the recognition of short peptide sequences by modular interaction domains play a central role in the assembly of signalling protein complexes and larger protein networks that regulate cellular behaviour. In addition to spatial and temporal factors, the specificity of signal transduction is intimately associated with the specificity of many co-operative, pairwise binding events upon which various pathways are built. Although protein interaction domains are usually identified via the recognition code, the consensus sequence motif, to which they selectively bind, they are highly versatile and play diverse roles in the cell. For example, a given interaction domain can bind to multiple sequences that exhibit no apparent identity, and, on the other hand, domains of the same class or different classes may favour a given consensus motif. This promiscuity in ligand selection is typified by the SH3 (Src homology 3) domain and several other interaction modules that commonly recognize proline-rich sequences. Furthermore, interaction domains are highly adaptable, a property that is essential for the evolution of novel pathways and modulation of signalling dynamics. The ability of certain interaction domains to perform multiple tasks, however, poses a challenge for the cell to control signalling specificity when cross-talk between pathways is undesired. Extensive structural and biochemical analysis of many interaction domains in recent years has started to shed light on the molecular basis underlying specific compared with diverse binding events that are mediated by interaction domains and the role affinity plays in affecting domain specificity and regulating cellular signal transduction.

Protein interactions within the TcZFP zinc finger family members of Trypanosoma cruzi: implications for their functions

The small zinc finger proteins tbZFP1 and tbZFP2 have been implicated in the control of Trypanosoma brucei differentiation to the procyclic form. Here, we report that the complete ZFP family in Trypanosoma cruzi is composed by four members, ZFP1A and B, and ZFP2A and B. ZFP1B is a paralog specific gene restricted to T. cruzi, while the ZFP2A and B paralogs diverged prior to the trypanosomatid lineage separation. Moreover, we demonstrate that TcZFP1 and TcZFP2 members interact with each other and that this interaction is mediated by a WW domain in TcZFP2. Also, TcZFP2B strongly homodimerizes by a glycine rich region absent in TcZFP2A. We propose a model to discuss the relevance of these protein-protein interactions in terms of the functions of these proteins.

Using 3D protein structures to derive 3D-QSARs

The three-dimensional structures of proteins are being solved apace, yet this information is often underused in quantitative structure-activity relationship (QSAR) studies. Here, we describe and compare methods for exploiting protein structures to derive 3D-QSARs. These methods can facilitate molecular design and lead optimization and should increasingly become a standard component of the drug designer's repertoire.: