Combining inference from evolution and geometric probability in protein structure evaluation

Starting from the hypothesis that evolutionarily important residues form a spatially limited cluster in a protein's native fold, we discuss the possibility of detecting a non-native structure based on the absence of such clustering. The relevant residues are determined using the Evolutionary Trace method. We propose a quantity to measure clustering of the selected residues on the structure and show that the exact values for its average and variance over several ensembles of interest can be found. This enables us to study the behavior of the associated z-scores. Since our approach rests on an analytic result, it proves to be general, customizable, and computationally fast. We find that clustering is indeed detectable in a large representative protein set. Furthermore, we show that non-native structures tend to achieve lower residue-clustering z-scores than those attained by the native folds. The most important conclusion that we draw from this work is that consistency between structural and evolutionary information, manifested in clustering of key residues, imposes powerful constraints on the conformational space of a protein.     

Potentiating effect of distant sites in non-phosphorylated cyclic peptide antagonists of the Grb2-SH2 domain

Without the presence of a phosphotyrosyl group, a phage library derived non-phosphorylated cyclic peptide ligand of Grb2-SH2 domain attributed its high affinity and specificity to well-defined and highly favored interactions of its structural elements with the binding pocket of the protein. We have disclosed a significant compensatory role of the Glu(2-) sidechain for the absence of the phosphate functionality on Tyr(0) in the peptide ligand, cyclo(CH(2)CO-Glu(2-)-Leu-Tyr(0)-Glu-Asn-Val-Gly-Met(5+)-Tyr-Cys)-amide (termed G1TE). In this study, we report the importance of hydrophobic residue at the Tyr+5 site in G1TE. Both acidic and basic amino acid substitutes are disfavored at this position, and replacement of Met with beta-tert-butyl-Ala was found to improve the antagonist properties. Besides, the polarity of the cyclization linkage was implicated as important in stabilizing the favored binding conformation. Oxidation of the thioether linkage into sulfoxide facilitated the binding to Grb2-SH2 markedly. Simultaneous modification of the three distant sites within G1TE provided the best agent with an IC(50) of 220 nM, which is among the most potent non-phosphorous- and non-phosphotyrosine-mimic containing Grb2-SH2 domain inhibitors yet reported. This potent peptidomimetic provides a novel template for the development of chemotherapeutic agents for the treatment of erbB2-related cancer. Biological assays on G1TE(Gla(2-)) in which the original residue of Glu(2-) was substituted by gamma-carboxyglutamic acid (Gla) indicated that it could inhibit the interaction between activated GF receptor and Grb2 protein in cell homogenates of MDA-MB-453 breast cancer cells at the 2 microM level. More significantly, both G1TE(Gla(2-)) alone and the conjugate of G1TE(Gla(2-)) with a peptide carrier can effectively inhibit intracellular association of erbB2 and Grb2 in the same cell lines with IC(50) of 50 and 2 microM, respectively.     

Reversible inhibition of the human xenobiotic-metabolizing enzyme arylamine N-acetyltransferase 1 by S-nitrosothiols

Human arylamine N-acetyltransferase 1 (NAT1) is a polymorphic phase II xenobiotic-metabolizing enzyme which catalyzes the biotransformation of primary aromatic amines, hydrazine drugs, and carcinogens. Structural and functional studies have shown that the NAT1 and factor XIII transglutaminase catalytic pockets are structurally related with the existence of a conserved catalytic triad (Cys-His-Asp). In addition, it has been reported that factor XIII transglutaminase activity could be regulated by nitric oxide (NO), in particular S-nitrosothiols (RSNO). We thus tested whether NAT1 could be a target of S-nitrosothiols. We show here that human NAT1 is reversibly inactivated by S-nitrosothiols such as SNAP (S-nitroso-N-acetyl-DL-penicillamine). A second-order rate constant for the inactivation of NAT1 by SNAP was determined (k(inact)=270M(-1)min(-1)) and shown to be in the same range of values reported for other enzymes. The inhibition of NAT1 by S-nitrosothiols was reversed by dithiothreitol and reduced glutathione, but not by ascorbate. As reported for some reactive cysteine-containing enzymes, our results suggest that inactivation of NAT1 by S-nitrosothiols is due to direct attack of the highly reactive cysteine residue in the enzyme active site on the sulfur of S-nitrosothiols to form a mixed disulfide between these NO-derived oxidants and NAT1. Finally, our findings suggest that, in addition to the polymorphic-dependent variation of NAT1 activity, NO-derived oxidants, in particular S-nitrosothiols, could also regulate NAT1 activity.     

Amino-acid substitutions at the fully exposed P1 site of bovine pancreatic trypsin inhibitor affect its stability

It is widely accepted that solvent-exposed sites in proteins play only a negligible role in determining protein energetics. In this paper we show that amino acid substitutions at the fully exposed Lys15 in bovine pancreatic trypsin inhibitor (BPTI) influenced the CD- and DSC-monitored stability: The T(den) difference between the least (P1 Trp) and the most stable (P1 His) mutant is 11.2 degrees C at pH 2.0. The DeltaH(den) versus T(den) plot for all the variants at three pH values (2.0, 2.5, 3.0) is linear (DeltaC(p,den) = 0.41 kcal* mole(-1) * K(-1); 1 cal = 4.18 J) leading to a DeltaG(den) difference of 2.1 kcal*mole(-1). Thermal denaturation of the variants monitored by CD signal at pH 2.0 in the presence of 6 M GdmCl again showed differences in their stability, albeit somewhat smaller (DeltaT(den) =7.1 degrees C). Selective reduction of the Cys14-Cys 38 disulfide bond, which is located in the vicinity of the P1 position did not eliminate the stability differences. A correlation analysis of the P1 stability with different properties of amino acids suggests that two mechanisms may be responsible for the observed stability differences: the reverse hydrophobic effect and amino acid propensities to occur in nonoptimal dihedral angles adopted by the P1 position. The former effect operates at the denatured state level and causes a drop in protein stability for hydrophobic side chains, due to their decreased exposure upon denaturation. The latter factor influences the native state energetics and results from intrinsic properties of amino acids in a way similar to those observed for secondary structure propensities. In conclusion, our results suggest that the protein-stability-derived secondary structure propensity scales should be taken with more caution.     

Modeling alpha-helical coiled-coil interactions: the axial and azimuthal alignment of 1B segments from vimentin intermediate filaments

Attempts at predicting the relative axial alignments of fibrous protein molecules in filamentous structures have relied upon representing the (multichain) molecular structure by a one-dimensional sequence of amino acids. Potential intermolecular ionic and apolar interactions were counted and determined as a function of the relative axial stagger between the molecules. No attempts were made to consider the azimuthal aspect of the interacting molecules and neither were apolar or ionic energy terms used. Surprisingly, this simple approach proved remarkably informative and yielded accurate predictions of the axial periods present. However, a more comprehensive analysis involving the energetics of aggregation taking due regard for the relative azimuths of the molecules as well as their separation should decrease the noise level in the calculations and reveal other pertinent information. Toward that end, we have modeled the interaction between two alpha-helical coiled-coil segments in intermediate filament molecules (1B segments from human vimentin). The relative axial alignment and polarity of the molecules is already known from detailed crosslinking studies and this provides a criterion against which the success (or otherwise) of the modeling can be judged. The results confirm that an antiparallel alignment of two 1B segments is preferred over any of the parallel options (as observed experimentally). The calculated axial alignment, however, is not identical to that observed from detailed crosslinking studies indicating that other parts of the molecule (probably the head and tail domains as well as other coiled-coil segments) have a crucial role in determining the precise mode of axial aggregation. The results also show that the apolar interactions seem to be significantly less important in the alignment process than the ionic ones. This is consistent with the observation of a well-defined period in the linear disposition of the charged (but not apolar) residues along the length of the outer surface of the vimentin molecule.     

Changes in lignin content of leaf litters during mulching

Alkaline nitrobenzene oxidation, ozonation and methoxyl content determinations were applied to decomposing leaf litter of Ginkgo biloba L., Cinnamomum camphora sieb., Zelkova serrata Makino and Firmiana simplex W. F. Wight, respectively, during mulching to investigate the properties and estimate changes in lignin composition and content. Since the Klason lignin residue originated from components highly resistant to degradation by acid, the methoxyl content of Klason residue was used to estimate the lignin content of leaf litter. Quantitative analysis of presumed lignin-derived fragments, by use of alkaline nitrobenzene oxidation and ozonation methods, suggested that the estimated lignin content approximates that of the real lignin content of leaves, which is greatly overestimated by the Klason procedure. The estimated lignin contents ranged from 3.9 to 10.0% while the Klason lignan residue varied from 37.1 to 46.7% in un-mulched leaf litter. The absolute amounts of the measured lignin somewhat decreased during mulching, while the structure of lignin remaining in leaf litters after mulching was considered not to be very different from its original structure.     

ADP-glucose pyrophosphorylase from potato tuber: site-directed mutagenesis of homologous aspartic acid residues in the small and large subunits

Asp142 in the homotetrameric ADP-glucose pyrophosphorylase (ADP-Glc PPase) enzyme from Escherichia coli was demonstrated to be involved in catalysis of this enzyme [Frueauf, J.B., Ballicora, M.A. and Preiss J. (2001) J. Biol. Chem., 276, 46319-46325]. The residue is highly conserved throughout the family of ADP-Glc PPases, as well as throughout the super-family of sugar-nucleotide pyrophosphorylases. In the heterotetrameric ADP-Glc PPase from potato (Solanum tuberosum L.) tuber, the homologous residue is present in both the small (Asp145) and the large (Asp160) subunits. It has been proposed that the small subunit of plant ADP-Glc PPases is catalytic, while the large subunit is modulatory; however, no catalytic residues have been identified. To investigate the function of these conserved Asp residues in the ADP-Glc PPase from potato tuber, we used site-directed mutagenesis to introduce either an Asn or a Glu. Kinetic analysis in the direction of synthesis or pyrophosphorolysis of ADP-Glc showed a significant decrease (more than four orders of magnitude) in the specific activity of the SD145NLwt, SD145NLD160N, and SD145NLD160E mutants, while the effect was smaller (approximately two orders of magnitude) with the SD145ELwt, SD145ELD160N, and SD145ELD160E mutants. By contrast, mutation of the large subunit alone did not affect the specific activity but did alter the apparent affinity for the activator 3-phosphoglycerate, showing two types of apparent roles for this residue in the different subunits. These results show that mutation of Asp160 of the large subunit does not affect catalysis, thus the large subunit is not catalytic, and that the negative charge of Asp145 in the small subunit is necessary for enzyme catalysis.     

Apoflavodoxin (un)folding followed at the residue level by NMR

The denaturant-induced (un)folding of apoflavodoxin from Azotobacter vinelandii has been followed at the residue level by NMR spectroscopy. NH groups of 21 residues of the protein could be followed in a series of 1H-15N heteronuclear single-quantum coherence spectra recorded at increasing concentrations of guanidinium hydrochloride despite the formation of protein aggregate. These NH groups are distributed throughout the whole apoflavodoxin structure. The midpoints of unfolding determined by NMR coincide with the one obtained by fluorescence emission spectroscopy. Both techniques give rise to unfolding curves with transition zones at significantly lower denaturant concentrations than the one obtained by circular dichroism spectroscopy. The NMR (un)folding data support a mechanism for apoflavodoxin folding in which a relatively stable intermediate is involved. Native apoflavodoxin is shown to cooperatively unfold to a molten globule-like state with extremely broadened NMR resonances. This initial unfolding step is slow on the NMR chemical shift timescale. The subsequent unfolding of the molten globule is faster on the NMR chemical shift timescale and the limited appearance of 1H-15N HSQC cross peaks of unfolded apoflavodoxin in the denaturant range studied indicates that it is noncooperative.