Comparison of the computed three-dimensional structures of oncogenic forms (bound to GDP) of theras-Gene-Encoded p21 protein with the structure of the normal (non-transforming) wild-type protein

Theras-oncogene-encoded p21 protein becomes oncogenic if amino acid substitutions occur at critical positions in the polypeptide chain. The most commonly found oncogenic forms contain Val in place of Gly 12 or Leu in place of Gln 61. To determine the effects of these substitutions on the three-dimensional structure of the whole p21 protein, we have performed molecular dynamics calculations on each of these three proteins bound to GDP and magnesium ion to compute the average structures of each of the three forms. Comparisons of the computed average structures shows that both oncogenic forms with Val 12 and Leu 61 differ substantially in structure from that of the wild type (containing Gly 12 and Gln 61) in discrete regions: residues 10–16, 32–47, 55–74, 85–89, 100–110, and 119–134. All of these regions occur in exposed loops, and several of them have already been found to be involved in the cellular functioning of the p21 protein. These regions have also previously been identified as the most flexible domains of the wild-type protein and have been bound to be the same ones that differ in conformation between transforming and nontransforming p21 mutant proteins neither of which binds nucleotide. The two oncogenic forms have similar conformations in their carboxyl-terminal domains, but differ in conformation at residues 32–47 and 55–74. The former region is known to be involved in the interaction with at least three downstream effector target proteins. Thus, differences in structure between the two oncogenic proteins may reflect different relative affinities of each oncogenic protein for each of these effector targets. The latter region, 55–74, is known to be a highly mobile segment of the protein. The results strongly suggest that critical oncogenic amino acid substitutions in the p21 protein cause changes in the structures of vital domains of this protein.     

Comparison of the predicted structure for the activated form of the P21 protein with the X-ray crystal structure

The predicted conformation and position of the central transforming region (residues 55–67) of the p21 protein are compared with the conformation and position of this segment in a recently determined X-ray crystal structure of residues 1–166 of this protein in the activated state bound to a nonhydrolyzable GTP derivative. We previously predicted that this segment of the protein would adopt a roughly extended conformation from Ile 55-Thr 58, a reverse turn at Ala 59-Gln 61, followed by an -helix from Glu 62-Met 67. We further predicted that this region of the activated protein occupies a position that is virtually identical to corresponding regions in the homologous purine nucleotide-binding proteins, bacterial elongation factor (EF-tu), and adenylate kinase (ADK). We find that there is a close correspondence between the conformation and position of our predicted structure and those found in the X-ray crystal structure. A mechanism for activation of the protein is proposed and is corroborated by X-ray crystallographic data.     

Structural effects of the binding of GTP to the wild-type and oncogenic forms of theras-gene-encoded p21 proteins

Molecular dynamics calculations have been performed to determine the average structures ofras-gene-encoded p21 proteins bound to GTP, i.e., the normal (wild-type) protein and two oncogenic forms of this protein, the Val 12- and Leu 61-p21 proteins. We find that the average structures for all of these proteins exhibit low coordinate fluctuations (which are highest for the normal protein), indicating convergence to specific structures. From previous dynamics calculations of the average structures of these proteins bound to GDP, major regional differences were found among these proteins (Monacoet al. (1995),J. Protein Chem., in press). We now find that the average structures of the oncogenic proteins are more similar to one another when the proteins are bound to GTP than when they are bound to GDP (Monacoet al. (1995),J. Protein Chem., in press). However, they still differ in structureat specific amino acid residues rather than in whole regions, in contradistinction to the results found for the p21-GDP complexes. Two exceptions are the regions 25–32, in an-helical region, and 97–110. The two oncogenic (Val 12- and Leu 61-) proteins have similar structures which differ significantly in the region of residues 97–110. This region has recently been identified as being critical in the interaction of p21 with kinase target proteins. The differences in structure between the oncogenic proteins suggest the existence of more than one oncogenic form of the p21 protein that can activate different signaling pathways.     

Deletion Mutants of Human Granulocyte-Macrophage Colony-Stimulating Factor

To study the structure–function relationship of the human granulocyte-macrophage colony-stimulating factor (GM-CSF), genes were constructed that encode its three deletion mutants: D1, a mutant with the deletion of six amino acid residues (37–42) some of which are a part of a -structural region; D2, a mutant with the deletion of the unstructured six-aa sequence of a loop (45–50); and D3, a mutant with the deletion of 14 aa residues (37–50) corresponding to the A–B loop and encoded by the second exon of the gmcsf gene. The expression products of these genes in E. coli were accumulated in a fraction of insoluble proteins. The secondary structures of the mutant proteins were similar to that of the full-size GM-CSF, but the biological activity of the deletion mutants was 130 times lower than that of the GM-CSF: they stimulated the proliferation of the TF-1 cell line at 3 ng/ml concentration. The resulting proteins displayed antagonistic properties toward the full-size GM-CSF, with the inhibition degree of its colony-stimulating activity being 27%. A decrease in the mutant activity in the row D2 > D1 > D3 implies the importance of the conserved hydrophobic residues involved in the formation of the -structure for the formation of the GM-CSF functional conformation.     

Prediction of the three-dimensional structure of the rap-1A protein from its homology to theras-gene-encoded p21 protein

rap-1A, an anti-oncogene-encoded protein, is aras-p21-like protein whose sequence is over 80% homologous to p21 and which interacts with the same intracellular target proteins and is activated by the same mechanisms as p21, e.g., by binding GTP in place of GDP. Both interact with effector proteins in the same region, involving residues 32–47. However, activated rap-1A blocks the mitogenic signal transducing effects of p21. Optimal sequence alignment of p21 and rap-1A shows two insertions of rap-1A atras positions 120 and 138. We have constructed the three-dimensional structure of rap-1A bound to GTP by using the energy-minimized three-dimensional structure ofras-p21 as the basis for the modeling using a stepwise procedure in which identical and homologous amino acid residues in rap-1A are assumed to adopt the same conformation as the corresponding residues in p21. Side-chain conformations for homologous and nonhomologous residues are generated in conformations that are as close as possible to those of the corresponding side chains in p21. The entire structure has been subjected to a nested series of energy minimizations. The final predicted structure has an overall backbone deviation of 0.7 å from that ofras-p21. The effector binding domains from residues 32–47 are identical in both proteins (except for different side chains of different residues at position 45). A major difference occurs in the insertion region at residue 120. This region is in the middle of another effector loop of the p21 protein involving residues 115–126. Differences in sequence and structure in this region may contribute to the differences in cellular functions of these two proteins.     

Mechanisms of drug-induced allergic contact dermatitis

Allergic contact dermatitis is induced by a wide variety of drugs that trigger specific immune responses following topical exposure. Identified chemical stuctures involved in such reactions include the mercuric and thiosalicylic acid groups of thimerosal, the diphenylketone group of the anti-inflammatory drug ketoprofen, the amide or ester structure of local anesthetics, and the side-chain and thiazolidine ring of -lactams. The T cell responses to such compounds involve CD4+ and CD8+ + T lymphocytes and also CD4–/CD8– + T cells. Although "T helper 2" cytokine production by drug-specific human T cells from patients with allergic contact dermatitis has been described, T helper 1-like and T cytotoxic 1-like responses clearly play key roles in this cutaneous reaction.     

VPS21 encodes a rab5-like GTP binding protein that is required for the sorting of yeast vacuolar proteins.

Many of the vacuolar protein sorting (vps) mutants of Saccharomyces cerevisiae exhibit severe defects in the sorting of vacuolar proteins but still retain near-normal vacuole morphology. The gene affected in one such mutant, vps21, has been cloned and found to encode a member of the ras-like GTP binding protein family. Sequence comparisons with other known GTP binding proteins indicate that Vps21p is unique but shares striking similarity with mammalian rab5 proteins (> 50% identity and > 70% similarity). Regions with highest similarity are clustered within the putative GTP binding motifs and the proposed effector domains of the Vps21/rab5 proteins. Point mutations constructed within these conserved regions inactivate Vps21p function; the mutant cells missort and secrete the soluble vacuolar hydrolase carboxypeptidase Y (CPY). Cells carrying a complete deletion of the VPS21 coding sequence (i) are viable but exhibit a growth defect at 38 degrees C, (ii) missort multiple vacuolar proteins, (iii) accumulate 40-50 nm vesicles and (iv) contain a large vacuole. VPS21 encodes a 22 kDa protein that binds GTP and fractionates with subcellular membranes. Mutant analysis indicates that the association with a membrane(s) is dependent on geranylgeranylation of the C-terminal cysteine residue(s) of Vps21p. We propose that Vps21p functions in the targeting and/or fusion of transport vesicles that mediate the delivery of proteins to the vacuole.     

Molecular Dynamics on Complexes of ras-p21 and Its Inhibitor Protein, rap-1A, Bound to the ras-Binding Domain of the raf-p74 Protein: Identification of Effector Domains in the raf Protein

We have computed the average structures for the ras-p21 protein and its strongly homologous inhibitor protein, rap-1A, bound to the ras-binding domain (RBD) of the raf protein, using molecular dynamics. Our purpose is to determine the differences in structure between these complexes that would result in no mitogenic activity of rap-1A-RBD but full activity of p21-RBD. We find that despite the similarities of the starting structures for both complexes, the average structures differ considerably, indicating that these two proteins do not interact in the same way with this vital target protein. p21 does not undergo major changes in conformation when bound to the RBD, while rap-1 A undergoes significant changes in structure on binding to the RBD, especially in the critical region around residue 61. The p21 and rap-1A make substantially different contacts with the RBD. For example, the loop region from residues 55–71 of rap-la makes extensive hydrogen-bond contacts with the RBD, while the same residues of p21 do not. Comparison of the structures of the RBD in both complexes reveals that it undergoes considerable changes in structure when its structure bound to p21 is compared with that bound to rap-1A. These changes in structure are due to displacements of regular structure (e.g., -helices and -sheets) rather than to changes in the specific conformations of the segments themselves. Three regions of the RBD have been found to differ significantly from one another in the two complexes: the binding interface between the two proteins at residues 60 and 70, the region around residues 105–106, and 118–120. These regions may constitute effector domains of the RBD whose conformations determine whether or not mitogenic signal transduction will occur.     

The structure of the carboxyl terminus of the p21 protein. Structural relationship to the nucleotide-binding/transforming regions of the protein

The carboxyl-terminal region of theras oncogene-encoded p21 protein is critical to the protein's function, since membrane binding through the C-terminus is necessary for its cellular activity. X-ray crystal structures for truncated p21 proteins are available, but none of these include the C-terminal region of the protein (from residues 172–189). Using conformational energy analysis, we determined the preferred three-dimensional structures for this C-terminal octadecapeptide of the H-ras oncogene p21 protein and generated these structures onto the crystal structure of the remainder of the protein. The results indicate that, like other membrane-associated proteins, the membrane-binding C-terminus of p21 assumes a helical hairpin conformation. In several low-energy orientations, the C-terminal structure is in close proximity to other critical locales of p21. These include the central transforming region (around Gln 61) and the amino terminal transforming region (around Gly 12), indicating that extracellular signals can be transduced through the C-terminal helical hairpin to the effector regions of the protein. This finding is consistent with the results of recent genetic experiments.     

Oligomeric structure of p21 ras proteins as determined by radiation inactivation

Using radiation inactivation we determined that p21 ras proteins exhibit an oligomeric target size when assayed both structurally and functionally. Similar target sizes of p21 in ras-transformed cells and in purified preparations of the protein suggested that its structure is homo-oligomeric. p21 monomers were destroyed by radiation with the same target size as the GTP binding activity, indicating the occurrence of a tight association allowing energy transfer between the monomers. Irradiation in the presence of GTP, dithiothreitol, or EDTA did not change the target size. Normal (Gly12) and transforming (Lys12) forms of the protein exhibited similar target sizes. The homo-oligomeric structure suggests that p21 ras proteins do not conform to the structure of monomeric alpha subunits in classical G proteins (alpha beta gamma heterotrimers) and establishes similarities with other homo-oligomeric proteins (such as Escherichia coli CRP) which acquire the active conformation through subunit reorientation upon nucleotide binding.     

Some aspects of the growth of Gracilaria tenuistipitata in pond culture

The effect of temperature, salinity, nitrogen, culture density and depth on the growth of Gracilaria tenuistipitata were investigated between April 1985 and March 1986 in outdoor ponds in Guangxi Province, South China. The mean annual growth rate was 2.4% per day. Under favourable temperatures of 20–30 °C, daily growth rate may reach as high as 3.3%. Salinity had an obvious effect on growth and photosynthesis and growth peaked at 21, with a broad plateau between 7–27. Growth experiments showed that a total nitrogen (NH₄-N plus NO₃-N) concentration of 4 M was sufficient to enable the plants to maintain a daily growth rate of 2.7%. The best growth of the plant was obtained at a culture density of 0.5–1 kg m^–2 and a culture depth of 30 cm in the pond.     

Post-translational modification of proteins by 15-carbon and 20-carbon isoprenoids in three mammalian cell lines

Summary A number of cellular proteins, including p21^ras, lamin B, and the G-protein subunits, undeDanvillergo post-translational modification by 15-carbon farnesyl or 20-carbon geranylgeranyl isoprenoid moieties derived from pyrophosphate intermediates of the cholesterol biosynthetic pathway. In this study, isoprenylated proteins in three mammalian cell lines (Hela cells, Rat-6 fibroblasts and COS cells) were radiolabeled with an isoprenoid precursor, [³H]mevalonate, and resolved by SDS gel electrophoresis. Groups of proteins with different molecular masses were eluted from the gels and the chain-lengths of the radiolabeled isoprenyl groups, released from the proteins by Raney-nickel-catalyzed desulfurization, were established by gel permeation chromatography. 15-Carbon and 20-carbon isoprenyl groups were found in separate classes of proteins within each cell line. With the exception of p21^ras, which incorporated a 15-carbon group when expressed in COS cells, the proteins in the region of the 21–28 kDa ras-related GTP binding proteins contained mostly 20-carbon isoprenyl chains. In contrast, proteins belonging to the 66–72 kDa nuclear lamin family, as well as unidentified proteins with molecular masses of 41–46 kDa and 53–55 kDa, contained predominantly 15-carbon isoprenyl chains. The chain-lengths of the isoprenoids associated with particular classes of proteins did not vary from one cell line to another, suggesting that the nature of the isoprenoid modification (farnesyl versus geranylgeranyl) is determined by intrinsic structural features of the proteins, rather than the cell type in which the proteins are expressed.Abbreviations MVA Mevalonolactone - SDS Sodium Dodecyl Sulfate - PAGE Polyacrylamide Gel Electrophoresis     

Structural effects of amino acid substitutions on the P21 proteins: Evidence for a malignant conformation

The conformational effects of different amino acid substitutions for Gly at position 12 in theras-oncogene-encoded P21 proteins have been investigated using conformational energy calculations. Mutations that cause amino acid substitutions for Gly 12 result in a protein that produces malignant transformation of cells. It had previously been shown that substitution of Val, Lys, or Ser for Gly at position 12 results in a major conformational change, and that the preferred lowest energy structure for each of the substituted peptides is identical. It is now found that substitution for Gly 12 of other amino acids that have widely disparate helix-nucleating potentials and completely different side chains (Asp, Asn, Cys, Phe, Tle, Leu, and Ala) all produce this identical lowest energy conformation. This finding is consistent with the recent results of site-specific mutagenesis experiments showing that P21 proteins containing these amino acids at position 12 all promote malignant transformation of cells and suggests the existence of a malignancy-causing conformation for the P21 proteins.     

Sequence Pattern for the Occurrence of N-Glycosylation in Proteins

To further understand the occurrence of N-glycosylation, 21 nonhomologous proteins with Asn–x–Ser/Thr sequence were investigated. The results showed that some oligopeptides with Gly residues (G–x–y or y–x–G) are adjacent to the N-glycosylated sequences. These oligopeptides are not only essential for the structure and function of the proteins, but they are also found to be often proteolytic processing sites. These properties suggest that these oligopeptides may be a sequence pattern for the occurrence of N-glycosylation. The implications of the findings for protein structure and function are discussed.     

Mutations in the GTP-binding site of GS alpha alter stimulation of adenylyl cyclase

Mutational replacements of specific residues in the GTP-binding pocket of the 21-kDa ras proteins (p21ras) reduce their GTPase activity. To test the possibility that the cognate regions of G protein alpha chains participate in GTP binding and hydrolysis, we compared signaling functions of normal and mutated alpha chains (termed alpha s) of Gs, the stimulatory regulator of adenylyl cyclase. alpha s chains were expressed in an alpha s-deficient S49 mouse lymphoma cell line, cyc-. alpha s in which leucine replaces glutamine 227 (corresponding to glutamine 61 of p21ras) constitutively activates adenylyl cyclase and reduces the kcat for GTP hydrolysis more than 100-fold. There is a smaller reduction in GTPase activity in another mutant in which valine replaces glycine 49 (corresponding to glycine 12 of p21ras). This mutant alpha s is a poor activator of adenylyl cyclase. Moreover, the glycine 49 protein, unlike normal alpha s, is not protected against tryptic cleavage by hydrolysis resistant GTP analogs; this finding suggests impairment of the mutant protein's ability to attain the active (GTP-bound) conformation. We conclude that alpha s residues near glutamine 227 and glycine 49 participate in binding and hydrolysis of GTP, although the GTP binding regions of alpha s and p21ras are not identical.     

Characterization of NO adducts of the diiron center in protein R2 of<Emphasis Type="Italic"> Escherichia coli</Emphasis> ribonucleotide reductase and site-directed variants; implications for the O<Subscript>2 </Subscript>activation mechanism*

The R2 subunit of Escherichia coli ribonucleotide reductase contains a diiron site that reacts with O₂ to produce a tyrosine radical (Y122·). In wild-type R2 (R2-wt), the first observable reaction intermediate is a high-valent [Fe^III-Fe^IV] state called compound X, but in related diiron proteins such as methane monooxygenase, ⁹-desaturase, and ferritin, peroxodiiron(III) complexes have been characterized. Substitution of iron ligand D84 by E within the active site of R2 allows an intermediate (-1,2-peroxo)diiron species to accumulate. To investigate the possible involvement of a bridging peroxo species within the O₂ activation sequence of R2-wt, we have characterized the iron-nitrosyl species that form at the diiron sites in R2-wt, R2-D84E, and R2-W48F/D84E by using vibrational spectroscopy. Previous work has shown that the diiron center in R2-wt binds one NO per iron to form an antiferromagnetically coupled [{FeNO}⁷]₂ center. In the wt and variant proteins, we also observe that both irons bind one NO to form a {FeNO}⁷ dimer where both Fe–N–O units share a common vibrational signature. In the wt protein, (Fe–NO), (Fe–N–O), and (N–O) bands are observed at 445, 434 and 1742 cm^–1, respectively, while in the variant proteins the (Fe–NO) and (Fe–N–O) bands are observed ~10 cm^–1 higher and the (N–O) ~10 cm^–1 lower at 1735 cm^–1. These results demonstrate that all three proteins accommodate fully symmetric [{FeNO}⁷]₂ species with two identical Fe–N–O units. The formation of equivalent NO adducts in the wt and variant proteins strongly favors the formation of a symmetric bridging peroxo intermediate during the O₂ activation process in R2-wt.     

A statistical analysis of side-chain conformations in proteins: Comparison with ECEPP predictions

A comparison of the statistical distributions of side-chain conformations of 17 amino acids (Gly, Ala, and Pro excluded), observed in 63 nonhomologous globular proteins (covering 10,832 residues), is made with similar distributions calculated from the low-energy conformational states for the same amino acids (blocked with acetyl and N-methylamide groups at the N- and C-termini, respectively) obtained by Vásquezet al. [(1983),Macromolecules 16, 1043–1049 using the ECEPP/2 force field. Those residues (i) with linear side chains (Arg, Lys, Met, Cys, Ser), or those that are unbranched through the-carbon atom (Glu, Gln) show good agreement, whereas (ii) those with side chains that are branched at C or C show poor agreement with ECEPP calculations. A possible explanation for this is shown to be the greater tendency for side-chain atoms in class (ii) to interact with the backbone and/or adjacent side chains. Accordingly, ECEPP/3 calculations, carried out after elongating the backbone chain of the model peptide unit (by adding three Ala residues on each side of the central residue, and then blocking the termini as before), result in distributions that are often closer to the observed side-chain distributions. The implications of these results for the relative importance of short-range versus long-range interactions in determining protein structure are discussed.     

Conformational changes induced by the transforming amino acid substitution in the transmembrane domain of the neu oncogene-encoded p185 protein

Theneu oncogene is frequently found in certain types of human carcinomas and has been shown to be activated in animal models by nitrosourea-induced mutation. The activating mutation in theneu oncogene results in the substitution of a glutamic acid for a valine at position 664 in the transmembrane domain of the encoded protein product of 185 kda (designated p185), which, on the basis of homology studies, is presumed to be a receptor for an as yet unidentified growth factor. It has been proposed that activating amino acid substitutions in this region of p185 lead to a conformational change in the protein which causes signal transduction via an increase in tyrosine kinase activity in the absence of any external signal. Using conformational energy analysis, we have determined the preferred three-dimensional structures for the transmembrane decapeptide (residues 658–667) of the p185 protein with valine and glutamic acid at the critical position 664. The results indicate that the global minimum energy conformation of the decapeptide from the normal protein with Val at position 664 is an -helix with a sharp bend (CD* conformation at residues 664 and 665) in this region, whereas the global minimum conformation for the decapeptide from the mutant transforming protein with Glu at position 664 assumes an all -helical configuration. Furthermore, the second highest energy conformation for the decapeptide from the normal protein is identical to the global minimum energy conformation for the decapeptide from the transforming protein, providing a possible explanation why overexpression of the normal protein also has a transforming effect. These results suggest there may be a normal and a transforming conformation for theneu-encoded p185 proteins which may explain their differences in transforming activity.     

Correlation Between Self-Association Modes and GTPase Activation of Dynamin

The GTPase activity of dynamin is obligatorily coupled, by a mechanism yet unknown, to the internalization of clathrin-coated endocytic vesicles. Dynamin oligomerizes in vitro and in vivo and both its mechanical and enzymatic activities appear to be mediated by this self-assembly. In this study we demonstrate that dynamin is characterized by a tetramer/monomer equilibrium with an equilibrium constant of 1.67 × 10¹⁷ M^–3. Stopped-flow fluorescence experiments show that the association rate constant for 2(3)-O-N-methylanthraniloyl (mant)GTP is 7.0 × 10^–5 M^–1 s^–1 and the dissociation rate constant is 2.1 s^–1, whereas the dissociation rate constant for mantdeoxyGDP is 93 s^–1. We also demonstrate the cooperativity of dynamin binding and GTPase activation on a microtubule lattice. Our results indicate that dynamin self-association is not a sufficient condition for the expression of maximal GTPase activity, which suggests that dynamin molecules must be in the proper conformation or orientation if they are to form an active oligomer.     

Dissection of the GTPase mechanism of Ras protein by MD analysis of Ras mutants

Controlling the hydrolysis rate of GTP bound to the p21ras protein is crucial for the delicate timing of many biological processes. A few mechanisms were suggested for the hydrolysis of GTP. To gain more insight into the individual elementary events of GTP hydrolysis, we carried out molecular dynamic analysis of wild-type p21ras and some of its mutants. It was recently shown that Ras-related proteins and mutants generally follow a linear free energy relationship (LFER) relating the rate of reaction to the pK(a) of the gamma-phosphate group of the bound GTP, indicating that proton transfer from the attacking water to the GTP is the first elementary event in the GTPase mechanism. However, some exceptions were observed. Thus, the Gly12 --> Aspartic p21ras (G12D) mutant had a very low GTPase activity although its pK(a) was very close to that of the wild-type ras. Here we compared the molecular dynamics (MD) of wild-type Ras and G12D, showing that in the mutant the catalytic water molecule is displaced to a position where proton transfer to GTP is unfavorable. These results suggest that the mechanism of GTPase is indeed composed of an initial proton abstraction from water by the GTP, followed by a nucleophilic attack of the hydroxide ion on the gamma-phosphorus of GTP.