The origins of the evolutionary signal used to predict protein-protein interactions

Background: The correlation of genetic distances between pairs of protein sequence alignments has been used to infer protein-protein interactions. It has been suggested that these correlations are based on the signal of co-evolution between interacting proteins. However, although mutations in different proteins associated with maintaining an interaction clearly occur (particularly in binding interfaces and neighbourhoods), many other factors contribute to correlated rates of sequence evolution. Proteins in the same genome are usually linked by shared evolutionary history and so it would be expected that there would be topological similarities in their phylogenetic trees, whether they are interacting or not. For this reason the underlying species tree is often corrected for. Moreover processes such as expression level, are known to effect evolutionary rates. However, it has been argued that the correlated rates of evolution used to predict protein interaction explicitly includes shared evolutionary history; here we test this hypothesis. Results: In order to identify the evolutionary mechanisms giving rise to the correlations between interaction proteins, we use phylogenetic methods to distinguish similarities in tree topologies from similarities in genetic distances. We use a range of datasets of interacting and non-interacting proteins from Saccharomyces cerevisiae. We find that the signal of correlated evolution between interacting proteins is predominantly a result of shared evolutionary rates, rather than similarities in tree topology, independent of evolutionary divergence. Conclusions: Since interacting proteins do not have tree topologies that are more similar than the control group of non-interacting proteins, it is likely that coevolution does not contribute much to, if any, of the observed correlations.

Solution nuclear magnetic resonance structure of the GATase subunit and structural Basis of the interaction between GATase and ATPPase subunits in a two-subunit-type GMPS from methanocaldococcus jannaschii

The solution structure of the monomeric glutamine amidotransferase (GATase) subunit of the Methanocaldococcus janaschii (Mj) guanosine monophosphate synthetase (GMPS) has been determined using high-resolution nuclear magnetic resonance methods. Gel filtration chromatography and N-15 backbone relaxation studies have shown that the Mj GATase subunit is present in solution as a 21 kDa (188-residue) monomer. The ensemble of 20 lowest-energy structures showed root-mean-square deviations of 0.35 +/- 0.06 angstrom for backbone atoms and 0.8 +/- 0.06 angstrom for all heavy atoms. Furthermore, 99.4% of the backbone dihedral angles are present in the allowed region of the Ramachandran map, indicating the stereochemical quality of the structure. The core of the tertiary structure of the GATase is composed of a seven-stranded mixed beta-sheet that is fenced by five alpha-helices. The Mj GATase is similar in structure to the Pyrococcus horikoshi (Ph) GATase subunit. Nuclear magnetic resonance (NMR) chemical shift perturbations and changes in line width were monitored to identify residues on GATase that were responsible for interaction with magnesium and the ATPPase subunit, respectively. These interaction studies showed that a common surface exists for the metal ion binding as well as for the protein-protein interaction. The dissociation constant for the GATase-Mg2+ interaction has been found to be similar to 1 mM, which implies that interaction is very weak and falls in the fast chemical exchange regime. The GATase-ATPPase interaction, on the other hand, falls in the intermediate chemical exchange regime on the NMR time scale. The implication of this interaction in terms of the regulation of the GATase activity of holo GMPS is discussed.

The beta hairpin structure within ribosomal protein S5 mediates interplay between domains II and IV and regulates HCV IRES function

Translation initiation in Hepatitis C Virus (HCV) is mediated by Internal Ribosome Entry Site (IRES), which is independent of cap-structure and uses a limited number of canonical initiation factors. During translation initiation IRES-40S complex formation depends on high affinity interaction of IRES with ribosomal proteins. Earlier, it has been shown that ribosomal protein S5 (RPS5) interacts with HCV IRES. Here, we have extensively characterized the HCV IRES-RPS5 interaction and demonstrated its role in IRES function. Computational modelling and RNA-protein interaction studies demonstrated that the beta hairpin structure within RPS5 is critically required for the binding with domains II and IV. Mutations disrupting IRES-RPS5 interaction drastically reduced the 80S complex formation and the corresponding IRES activity. Computational analysis and UV cross-linking experiments using various IRES-mutants revealed interplay between domains II and IV mediated by RPS5. In addition, present study demonstrated that RPS5 interaction is unique to HCV IRES and is not involved in 40S-3 ` UTR interaction. Further, partial silencing of RPS5 resulted in preferential inhibition of HCV RNA translation. However, global translation was marginally affected by partial silencing of RPS5. Taken together, results provide novel molecular insights into IRES-RPS5 interaction and unravel its functional significance in mediating internal initiation of translation.

The Structure of the Holo-Acyl Carrier Protein of Leishmania major Displays a Remarkably Different Phosphopantetheinyl Transferase Binding Interface

The genome of Leishmania major encodes a type II fatty acid biosynthesis pathway for which no structural or biochemical information exists. Here, for the first time, we have characterized the central player of the pathway, the acyl carrier protein (LmACP), using nuclear magnetic resonance (NMR). Structurally, the LmACP molecule is similar to other type II ACPs, comprising a four-helix bundle, enclosing a hydrophobic core. Dissimilarities in sequence, however, exist in helix II (recognition helix) of the protein. The enzymatic conversion of apo-LmACP into the holo form using type I (Escherichia coli AcpS) and type II (Sfp type) phosphopantetheinyl transferases (PPTs) is relatively slow. Mutagenesis studies underscore the importance of the residues present at the protein protein interaction interface of LmACP in modulating the activity of PPTs. Interestingly, the cognate PPT for this ACP, the L. major 4'-phosphopantetheinyl transferase (LmPPT), does not show any enzymatic activity toward it, though it readily converts other type I and type II ACPs into their holo forms. NMR chemical shift perturbation studies suggest a moderately tight complex between LmACP and its cognate PPT, suggesting inhibition. We surmise that the unique surface of LmACP might have evolved to complement its cognate enzyme (LmPPT), possibly for the purpose of regulation.

Influence of deamidation(s) in the 67-74 region of ribonuclease on its refolding

The influence of chemical mutation featuring the selective conversion of asparagine or glutamine to aspartic or glutamic acid, respectively, on the kinetics of refolding of reduced RNase has been studied. The monodeamidated derivatives of RNase A, viz. RNase Aa1a, Aa1b, and Aa1c having their deamidations in the region 67-74, were found to regain nearly their original enzymatic activity. However, a marked difference in the kinetics of refolding is seen, the order of regain of enzymic activity being RNase A greater than Aa1c congruent to Aa1a greater than Aa1b. The similarities in the distinct elution positions on Amberlite XE-64, gel electrophoretic mobilities, and u.v. spectra of reoxidized and native derivatives indicated that the native structures are formed. The slower rate of reappearance of enzymic activity in the case of the monodeamidated derivatives appears to result from altered interactions in the early stages of refolding. The roles of some amino acid residues of the 67-74 region in the pathway of refolding of RNase A are discussed.

On the reactivity of carboxyl groups of ribonuclease-A in anhydrous methanol

The esterification of Ribonuclease-A in methanol/0.1 M hydrochloric acid has been studied by measuring the decrease in the number of titratable groups of the protein and estimating the amount of methanol incorporated. Esterification of nearly five of the 11 free carboxyl groups of the protein resulted in almost complete inactivation of the enzyme. The initial products of esterification have been chromatographed on Amberlite columns, and five partially active methyl ester derivatives of Ribonuclease-A have been isolated. The dimethyl ester, the initial product of esterification with reduced catalytic activity, has the carboxyl groups of Glu-49 and Asp-53 modified. Even in the non-aqueous solvent, as in the native structure of the protein in aqueous solution, these carboxyl groups are the fast reacting ones. Subsquently, the esterification reaction appears to proceed preferentially at the C-terminal region of the molecule. Comparison of the reactivities of carboxyl groups of Ribonuclease-A in acidic methanol to that known in aqueous solutions (with carbodiimides) suggests that the structure of Ribonuclease-A in the non-aqueous solvent resembles, at least in part, the structure in aqueous environment.

Observation of a mixed antiparallel and parallel beta-sheet motif in the crystal structure of Boc-Ala-Ile-Aib-OMe

A diastereomeric mixture of the tripeptide Boc-Ala-Ile-Aib-OMe crystallized in the space group P1 from CH3OH/H2O. The unit cell parameters are a = 10.593(2) A, b = 14.377(3) A, c = 17.872(4) A, alpha = 104.41(2) degrees, beta = 90.55(2) degrees, gamma = 106.91(2) degrees, V = 2512.4 A3, Z = 4. X-Ray crystallographic studies show the presence of four molecules in the asymmetric unit consisting of two pairs of diastereomeric peptides, Boc-L-Ala-L-Ile-Aib-OMe and Boc-L-Ala-D-Ile-Aib-OMe. The four molecules in the asymmetric unit form a rarely found mixed antiparallel and parallel beta-sheet hydrogen bond motif. The Ala and (L,D)-Ile residues in all the four molecules adopt the extended conformations, while the phi, psi values of the Aib residues are in the right-handed helical region. In one of the molecules the Ile sidechain adopts the unusual gauche conformation about the C beta-C gamma bond.

NMR and X-ray crystallographic studies on cyclic tetrapeptide, cyclo (D-Phe-Pro-Sar-Gly)

The conformation of the synthetic cyclic tetrapeptide cyclo(D-Phe-Pro-Sar-Gly) has been determined in solution using the nuclear magnetic resonance technique and in the crystal state by X-ray crystallography. Results showed that the peptide exhibited two different conformations in solution, conformer 1 having cis-trans-cis-trans peptide bonds and conformer 2 having trans-cis-trans-cis peptide bonds. No intramolecular hydrogen bonds were observed in the structures. The X-ray diffraction studies showed the crystals to be orthorhombic with space group P2(1)2(1)2(1) with unit-cell dimensions, a = 5.790, b = 10.344, c = 31.446 A, Z = 4, R = 0.104 for 2301 observed reflections. The crystal structure showed only one type of conformer having cis-trans-cis-trans peptide bonds similar to the conformer 1 in solution.

Prediction of Candidate Primary Immunodeficiency Disease Genes Using a Support Vector Machine Learning Approach

Screening and early identification of primary immunodeficiency disease (PID) genes is a major challenge for physicians. Many resources have catalogued molecular alterations in known PID genes along with their associated clinical and immunological phenotypes. However, these resources do not assist in identifying candidate PID genes. We have recently developed a platform designated Resource of Asian PDIs, which hosts information pertaining to molecular alterations, protein-protein interaction networks, mouse studies and microarray gene expression profiling of all known PID genes. Using this resource as a discovery tool, we describe the development of an algorithm for prediction of candidate PID genes. Using a support vector machine learning approach, we have predicted 1442 candidate PID genes using 69 binary features of 148 known PID genes and 3162 non-PID genes as a training data set. The power of this approach is illustrated by the fact that six of the predicted genes have recently been experimentally confirmed to be PID genes. The remaining genes in this predicted data set represent attractive candidates for testing in patients where the etiology cannot be ascribed to any of the known PID genes.

Stereochemistry of peptides containing 1-aminocycloheptane-1-carboxylic acid (Ac7c)

The crystal structures of four peptides incorporating 1-aminocycloheptane-1-carboxylic acid (Ac7c) are described. Boc-Aib-Ac7c-NHMe and Boc-Pro-Ac7c-Ala-OMe adopt beta-turn conformations stabilized by an intramolecular 4----1 hydrogen bond, the former folding into a type-I/III beta-turn and the latter into a type-II beta-turn. In the dipeptide esters, Boc-Aib-Ac7c-OMe and Boc-Pro-Ac7c-OMe, the Ac7c and Aib residues adopt helical conformations, while the Pro residue remains semi-extended in both the molecules of Boc-Pro-Ac7c-OMe found in the asymmetric unit. The cycloheptane ring of Ac7c residues adopts a twist-chair conformation in all the peptides studied. 1H-NMR studies in CDCl3 and (CD3)2SO and IR studies in CDCl3 suggest that Boc-Aib-Ac7c-NHMe and Boc-Pro-Ac7c-Ala-OMe maintain the beta-turn conformations in solution.

Monoclinic polymorph of Boc-Trp-Ile-Ala-Aib-Ile-Val-Aib-Leu-Aib-Pro-OMe(anhydrous). Parallel packing of 3(10)-/alpha-helices and a transition of helix type

The structures of two crystal forms of Boc-Trp-Ile-Ala-Aib-Ile-Val-Aib-Leu-Aib-Pro-OMe have been determined. The triclinic form (P1, Z = 1) from DMSO/H2O crystallizes as a dihydrate (Karle, Sukumar & Balaram (1986) Proc, Natl, Acad. Sci. USA 83, 9284-9288). The monoclinic form (P2(1), Z = 2) crystallized from dioxane is anhydrous. The conformation of the peptide is essentially the same in both crystal system, but small changes in conformational angles are associated with a shift of the helix from a predominantly alpha-type to a predominantly 3(10)-type. The r.m.s. deviation of 33 atoms in the backbone and C beta positions of residues 2-8 is only 0.29 A between molecules in the two polymorphs. In both space groups, the helical molecules pack in a parallel fashion, rather than antiparallel. The only intermolecular hydrogen bonding is head-to-tail between helices. There are no lateral hydrogen bonds. In the P2(1) cell, a = 9.422(2) A, b = 36.392(11) A, c = 10.548(2) A, beta = 111.31(2) degrees and V = 3369.3 A for 2 molecules of C60H97N11O13 per cell.

Aggregation studies in crystals of apolar helical peptides: Boc-Aib-Val-Ala-Leu-Aib-Val-Ala-Leu-Aib-OMe

In the crystal, the backbone of Boc-(Aib-Val-Ala-Leu)2-Aib-OMe adopts a helical form with four alpha-type hydrogen bonds in the middle, flanked by 3(10)-type hydrogen bonds at either end. The helical molecules stack in columns with head-to-tail hydrogen bonds, either directly between NH and CO, or bridged by solvent molecules. The packing of the helices is parallel, even in space group P2(1). Cell parameters are a = 9.837(2) A, b = 15.565(3) A, c = 20.087(5) A, beta = 96.42(2) degrees, dcalc = 1.091 g/cm3 for C46H83N9O12.1.5H2O.0.67CH3OH. There appears to be some hydration of the backbone in this apolar helix.

Parallel and antiparallel aggregation of alpha-helices. Crystal structures of two apolar decapeptides X-Trp-Ile-Ala-Aib-Ile-Val-Aib-Leu-Aib-Pro-OMe (X = Boc, Ac)

An apolar helical decapeptide with different end groups, Boc- or Ac-, crystallizes in a completely parallel fashion for the Boc-analog and in an antiparallel fashion for the Ac-analog. In both crystals, the packing motif consists of rows of parallel molecules. In the Boc-crystals, adjacent rows assemble with the helix axes pointed in the same direction. In the Ac-crystals, adjacent rows assemble with the helix axes pointed in opposite directions. The conformations of the molecules in both crystals are quite similar, predominantly alpha-helical, except for the tryptophanyl side chain where chi 1 congruent to 60 degrees in the Boc- analog and congruent to 180 degrees in the Ac-analog. As a result, there is one lateral hydrogen bond between helices, N(1 epsilon)...O(7), in the Ac-analog. The structures do not provide a ready rationalization of packing preference in terms of side-chain interactions and do not support a major role for helix dipole interactions in determining helix orientation in crystals. The crystal parameters are as follow. Boc-analog: C60H97N11O13.C3H7OH, space group Pl with a = 10.250(3) A, b = 12.451(4) A, c = 15.077(6) A, alpha = 96.55(3) degrees, beta = 92.31(3) degrees, gamma = 106.37(3) degrees, Z = 1, R = 5.5% for 5581 data ([F] greater than 3.0 sigma(F)), resolution 0.89 A. Ac-analog: C57H91N11O12, space group P2(1) with a = 9.965(1) A, b = 19.707(3) A, c = 16.648(3) A, beta = 94.08(1), Z = 2, R = 7.2% for 2530 data ([F] greater than 3.0 sigma(F)), resolution 1.00 A.

Stereochemistry of peptides containing 1-aminocycloheptane-1-carboxylic acid (Ac-c)

The crystal structures of four peptides incorporating 1-aminocycloheptane-1-carboxylic acid (Ac7c) are described. Boc-Aib-Ac7c-NHMe and Boc-Pro-Ac7c-Ala-OMe adopt beta-turn conformations stabilized by an intramolecular 4----1 hydrogen bond, the former folding into a type-I/III beta-turn and the latter into a type-II beta-turn. In the dipeptide esters, Boc-Aib-Ac7c-OMe and Boc-Pro-Ac7c-OMe, the Ac7c and Aib residues adopt helical conformations, while the Pro residue remains semi-extended in both the molecules of Boc-Pro-Ac7c-OMe found in the asymmetric unit. The cycloheptane ring of Ac7c residues adopts a twist-chair conformation in all the peptides studied. 1H-NMR studies in CDCl3 and (CD3)2SO and IR studies in CDCl3 suggest that Boc-Aib-Ac7c-NHMe and Boc-Pro-Ac7c-Ala-OMe maintain the beta-turn conformations in solution.