Unveiling Protein Functions through the Dynamics of the Interaction Network

Protein interaction networks have become a tool to study biological processes, either for predicting molecular functions or for designing proper new drugs to regulate the main biological interactions. Furthermore, such networks are known to be organized in sub-networks of proteins contributing to the same cellular function. However, the protein function prediction is not accurate and each protein has traditionally been assigned to only one function by the network formalism. By considering the network of the physical interactions between proteins of the yeast together with a manual and single functional classification scheme, we introduce a method able to reveal important information on protein function, at both micro- and macro-scale. In particular, the inspection of the properties of oscillatory dynamics on top of the protein interaction network leads to the identification of misclassification problems in protein function assignments, as well as to unveil correct identification of protein functions. We also demonstrate that our approach can give a network representation of the meta-organization of biological processes by unraveling the interactions between different functional classes

Structural homologies with ATP- and folate-binding enzymes in the crystal structure of folylpolyglutamate synthetase

Folylpolyglutamate synthetase, which is responsible for the addition of a polyglutamate tail to folate and folate derivatives, is an ATP-dependent enzyme isolated from eukaryotic and bacterial sources, where it plays a key role in the retention of the intracellular folate pool. Here, we report the 2.4-Å resolution crystal structure of the MgATP complex of the enzyme from Lactobacillus casei. The structural analysis reveals that folylpolyglutamate synthetase is a modular protein consisting of two domains, one with a typical mononucleotide-binding fold and the other strikingly similar to the folate-binding enzyme dihydrofolate reductase. We have located the active site of the enzyme in a large interdomain cleft adjacent to an ATP-binding P-loop motif. Opposite this site, in the C domain, a cavity likely to be the folate binding site has been identified, and inspection of this cavity and the surrounding protein structure suggests that the glutamate tail of the substrate may project into the active site. A further feature of the structure is a well defined Ω loop, which contributes both to the active site and to interdomain interactions. The determination of the structure of this enzyme represents the first step toward the elucidation of the molecular mechanism of polyglutamylation of folates and antifolates.

Complete resolution of the solid-state NMR spectrum of a uniformly 15N-labeled membrane protein in phospholipid bilayers

Complete resolution of the amide resonances in a three-dimensional solid-state NMR correlation spectrum of a uniformly 15N-labeled membrane protein in oriented phospholipid bilayers is demonstrated. The three orientationally dependent frequencies, 1H chemical shift, 1H–15N dipolar coupling, and 15N chemical shift, associated with each amide resonance are responsible for resolution among resonances and provide sufficient angular restrictions for protein structure determination. Because the protein is completely immobilized by the phospholipids on the relevant NMR time scales (10 kHz), the linewidths will not degrade in the spectra of larger proteins. Therefore, these results demonstrate that solid-state NMR experiments can overcome the correlation time problem and extend the range of proteins that can have their structures determined by NMR spectroscopy to include uniformly 15N-labeled membrane proteins in phospholipid bilayers.

Native protein sequences are close to optimal for their structures

How large is the volume of sequence space that is compatible with a given protein structure? Starting from random sequences, low free energy sequences were generated for 108 protein backbone structures by using a Monte Carlo optimization procedure and a free energy function based primarily on Lennard–Jones packing interactions and the Lazaridis–Karplus implicit solvation model. Remarkably, in the designed sequences 51% of the core residues and 27% of all residues were identical to the amino acids in the corresponding positions in the native sequences. The lowest free energy sequences obtained for ensembles of native-like backbone structures were also similar to the native sequence. Furthermore, both the individual residue frequencies and the covariances between pairs of positions observed in the very large SH3 domain family were recapitulated in core sequences designed for SH3 domain structures. Taken together, these results suggest that the volume of sequence space optimal for a protein structure is surprisingly restricted to a region around the native sequence.

Protein fluctuations are sensed by stimulated infrared echoes of the vibrations of carbon monoxide and azide probes

The correlation functions of the fluctuations of vibrational frequencies of azide ions and carbon monoxide in proteins are determined directly from stimulated photon echoes generated with femtosecond infrared pulses. The asymmetric stretching vibration of azide bound to carbonic anhydrase II exhibits a pronounced evolution of its vibrational frequency distribution on the time scale of a few picoseconds, which is attributed to modifications of the ligand structure through interactions with the nearby Thr-199. When azide is bound in hemoglobin, a more complex evolution of the protein structure is required to interchange the different ligand configurations, as evidenced by the much slower relaxation of the frequency distribution in this case. The time evolution of the distribution of frequencies of carbon monoxide bound in hemoglobin occurs on the ≈10-ps time scale and is very nonexponential. The correlation functions of the frequency fluctuations determine the evolution of the protein structure local to the probe and the extent to which the probe can navigate those parts of the energy landscape where the structural configurations are able to modify the local potential energy function of the probe.

Diffractive optics-based heterodyne-detected four-wave mixing signals of protein motion: From “protein quakes” to ligand escape for myoglobin

Ligand transport through myoglobin (Mb) has been observed by using optically heterodyne-detected transient grating spectroscopy. Experimental implementation using diffractive optics has provided unprecedented sensitivity for the study of protein motions by enabling the passive phase locking of the four beams that constitute the experiment, and an unambiguous separation of the Real and Imaginary parts of the signal. Ligand photodissociation of carboxymyoglobin (MbCO) induces a sequence of events involving the relaxation of the protein structure to accommodate ligand escape. These motions show up in the Real part of the signal. The ligand (CO) transport process involves an initial, small amplitude, change in volume, reflecting the transit time of the ligand through the protein, followed by a significantly larger volume change with ligand escape to the surrounding water. The latter process is well described by a single exponential process of 725 ± 15 ns at room temperature. The overall dynamics provide a distinctive signature that can be understood in the context of segmental protein fluctuations that aid ligand escape via a few specific cavities, and they suggest the existence of discrete escape pathways.

Redesigning secondary structure to invert coenzyme specificity in isopropylmalate dehydrogenase.

Rational engineering of enzymes involves introducing key amino acids guided by a knowledge of protein structure to effect a desirable change in function. To date, all successful attempts to change specificity have been limited to substituting individual amino acids within a protein fold. However, the infant field of protein engineering will only reach maturity when changes in function can be generated by rationally engineering secondary structures. Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to systematically invert the coenzyme specificity of Thermus thermophilus isopropylmalate dehydrogenase from a 100-fold preference for NAD to a 1000-fold preference for NADP. The engineered mutant, which is twice as active as wild type, contains four amino acid substitutions and an alpha-helix and loop that replaces the original beta-turn. These results demonstrate that rational engineering of secondary structures to produce enzymes with novel properties is feasible.

Mutations in copper-zinc superoxide dismutase that cause amyotrophic lateral sclerosis alter the zinc binding site and the redox behavior of the protein.

A series of mutant human and yeast copper-zinc superoxide dismutases has been prepared, with mutations corresponding to those found in familial amyotrophic lateral sclerosis (ALS; also known as Lou Gehrig's disease). These proteins have been characterized with respect to their metal-binding characteristics and their redox reactivities. Replacement of Zn2+ ion in the zinc sites of several of these proteins with either Cu2+ or Co2+ gave metal-substituted derivatives with spectroscopic properties different from those of the analogous derivative of the wild-type proteins, indicating that the geometries of binding of these metal ions to the zinc site were affected by the mutations. Several of the ALS-associated mutant copper-zinc superoxide dismutases were also found to be reduced by ascorbate at significantly greater rate than the wild-type proteins. We conclude that similar alterations in the properties of the zinc binding site can be caused by mutations scattered throughout the protein structure. This finding may help to explain what is perhaps the most perplexing question in copper-zinc superoxide dismutase-associated familial ALS-i.e., how such a diverse set of mutations can result in the same gain of function that causes the disease.

Characterizations of diverse residue clusters in protein three-dimensional structures.

We present new methods for identifying and analyzing statistically significant residue clusters that occur in three-dimensional (3D) protein structures. Residue clusters of different kinds occur in many contexts. They often feature the active site (e.g., in substrate binding), the interface between polypeptide units of protein complexes, regions of protein-protein and protein-nucleic acid interactions, or regions of metal ion coordination. The methods are illustrated with 3D clusters centering on four themes. (i) Acidic or histidine-acidic clusters associated with metal ions. (ii) Cysteine clusters including coordination of metals such as zinc or iron-sulfur structures, cysteine knots prominent in growth factors, multiple sets of buried disulfide pairings that putatively nucleate the hydrophobic core, or cysteine clusters of mostly exposed disulfide bridges. (iii) Iron-sulfur proteins and charge clusters. (iv) 3D environments of multiple histidine residues. Study of diverse 3D residue clusters offers a new perspective on protein structure and function. The algorithms can aid in rapid identification of distinctive sites, suggest correlations among protein structures, and serve as a tool in the analysis of new structures.

Clusters of charged residues in protein three-dimensional structures.

Statistically significant charge clusters (basic, acidic, or of mixed charge) in tertiary protein structures are identified by new methods from a large representative collection of protein structures. About 10% of protein structures show at least one charge cluster, mostly of mixed type involving about equally anionic and cationic residues. Positive charge clusters are very rare. Negative (or histidine-acidic) charge clusters often coordinate calcium, or magnesium or zinc ions [e.g., thermolysin (PDB code: 3tln), mannose-binding protein (2msb), aminopeptidase (1amp)]. Mixed-charge clusters are prominent at interchain contacts where they stabilize quaternary protein formation [e.g., glutathione S-transferase (2gst), catalase (8act), and fructose-1,6-bisphosphate aldolase (1fba)]. They are also involved in protein-protein interaction and in substrate binding. For example, the mixed-charge cluster of aspartate carbamoyl-transferase (8atc) envelops the aspartate carbonyl substrate in a flexible manner (alternating tense and relaxed states) where charge associations can vary from weak to strong. Other proteins with charge clusters include the P450 cytochrome family (BM-3, Terp, Cam), several flavocytochromes, neuraminidase, hemagglutinin, the photosynthetic reaction center, and annexin. In each case in Table 2 we discuss the possible role of the charge clusters with respect to protein structure and function.

Structure determination of murine mitochondrial carbonic anhydrase V at 2.45-A resolution: implications for catalytic proton transfer and inhibitor design.

The three-dimensional structure of murine mitochondrial carbonic anhydrase V has been determined and refined at 2.45-A resolution (crystallographic R factor = 0.187). Significant structural differences unique to the active site of carbonic anhydrase V are responsible for differences in the mechanism of catalytic proton transfer as compared with other carbonic anhydrase isozymes. In the prototypical isozyme, carbonic anhydrase II, catalytic proton transfer occurs via the shuttle group His-64; carbonic anhydrase V has Tyr-64, which is not an efficient proton shuttle due in part to the bulky adjacent side chain of Phe-65. Based on analysis of the structure of carbonic anhydrase V, we speculate that Tyr-131 may participate in proton transfer due to its proximity to zinc-bound solvent, its solvent accessibility, and its electrostatic environment in the protein structure. Finally, the design of isozyme-specific inhibitors is discussed in view of the complex between carbonic anhydrase V and acetazolamide, a transition-state analogue. Such inhibitors may be physiologically important in the regulation of blood glucose levels.

Nuclear magnetic dipole interactions in field-oriented proteins: information for structure determination in solution.

The measurement of dipolar contributions to the splitting of 15N resonances of 1H-15N amide pairs in multidimensional high-field NMR spectra of field-oriented cyanometmyoglobin is reported. The splittings appear as small field-dependent perturbations of normal scalar couplings. Assignment of more than 90 resonances to specific sequential sites in the protein allows correlation of the dipolar contributions with predictions based on the known susceptibility and known structure of the protein. Implications as an additional source of information for protein structure determination in solution are discussed.

The Structure of theCyprinid herpesvirus 3ORF112-Zα·Z-DNA Complex Reveals a Mechanism of Nucleic Acids Recognition Conserved with E3L, a Poxvirus Inhibitor of Interferon Response

In vertebrate species, the innate immune system down-regulates protein translation in response to viral infection through the action of the double-stranded RNA (dsRNA)-activated protein kinase (PKR). In some teleost species another protein kinase, Z-DNA-dependent protein kinase (PKZ), plays a similar role but instead of dsRNA binding domains, PKZ has Zα domains. These domains recognize the left-handed conformer of dsDNA and dsRNA known as Z-DNA/Z-RNA. Cyprinid herpesvirus 3 infects common and koi carp, which have PKZ, and encodes the ORF112 protein that itself bears a Zα domain, a putative competitive inhibitor of PKZ. Here we present the crystal structure of ORF112-Zα in complex with an 18-bp CpG DNA repeat, at 1.5 Å. We demonstrate that the bound DNA is in the left-handed conformation and identify key interactions for the specificity of ORF112. Localization of ORF112 protein in stress granules induced in Cyprinid herpesvirus 3-infected fish cells suggests a functional behavior similar to that of Zα domains of the interferon-regulated, nucleic acid surveillance proteins ADAR1 and DAI.

SVMtm: Support vector machines to predict transmembrane segments

A new method has been developed for prediction of transmembrane helices using support vector machines. Different coding schemes of protein sequences were explored, and their performances were assessed by crossvalidation tests. The best performance method can predict the transmembrane helices with sensitivity of 93.4% and precision of 92.0%. For each predicted transmembrane segment, a score is given to show the strength of transmembrane signal and the prediction reliability. In particular, this method can distinguish transmembrane proteins from soluble proteins with an accuracy of similar to99%. This method can be used to complement current transmembrane helix prediction methods and can be Used for consensus analysis of entire proteomes . The predictor is located at http://genet.imb.uq.edu.au/predictors/ SVMtm. (C) 2004 Wiley Periodicals, Inc.

Structure and folding of potato type II proteinase inhibitors: Circular permutation and intramolecular domain swapping

Potato type II serine proteinase inhibitors are proteins that consist of multiple sequence repeats, and exhibit a multidomain structure. The structural domains are circular permutations of the repeat sequence.. as a result or intramolecular domain swapping. Structural studies give indications for the origins of this folding behaviour, and the evolution of the inhibitor family.