Correlation of structure with transforming activity of the P21 proteins with substitutions of L-amino acids for Gly at position 13

The effects of amino acid substitutions for Gly 13 on the structure of the transforming region (Leu 6-Gly 15) of the P21 proteins have been explored using conformational energy calculations. It has been found that the substitution of Asp for Gly at this position results in a protein capable of transforming cells into malignant ones. Proteins that contain Ser at position 13 (but no other substitutions), however, transform cells with a greatly reduced activity. The transforming peptide with Asp 13 adopts a conformation that is different from the one for the peptide from the normal protein (with Gly 12 and Gly 13) and that may result in expression of a higher energy malignancy-producing form. The Ser-containing peptide adopts as its lowest energy conformation one that is identical to that of the peptide from the normal protein, thus explaining its lack of transforming activity. From analysis of the interactions preventing the Asp 13-containing peptide from adopting the normal conformation, it is predicted that substitutions of amino acids with branched side chains atC , such as Val, Ile, and Thr, should promote cell transformation. This prediction with Val has recently been confirmed in genetic experiments.     

Conformational effects of the substitution of Arg for Gly 13 in the ras oncogene-encoded P21 protein

The effect of the substitution of Arg for Gly 13 on the structure of the transforming region decapeptide (Leu 6-Gly 15) of the ras oncogene encoded P21 protein has been investigated using conformational energy analysis. A human malignancy has been identified that contains a ras gene with a single mutation in the thirteenth codon such that the encoded protein would have Arg substituted for Gly at this position, and transfection of cells in culture with this gene results in malignant transformation. Conformational analysis demonstrates that the Arg 13 decapeptide adopts a conformation identical to that for other peptides with substitutions at position 13 (Asp 13, Val 13) from transforming proteins that is distinctively different from that for peptides (Gly 13, Ser 13) from normal, nontransforming proteins. This is found to be an indirect effect resulting from changes in the conformation of Gly 12 produced by substitutions at position 13. These results are consistent with recent analysis of crystallographic data of proteins on conformational preferences for glycine in tripeptide sequences.     

Correlation of structure with transforming activity of the P21 proteins with substitutions at positions 12 and 16

The structural effects of amino acid substitutions at positions 12 and 16 in the amino-terminal segment (Tyr 4-Ala 18) of the ras-oncogene-encoded P21 proteins have been investigated using conformational energy analysis. The P21 protein with Val at position 12 and Lys at position 16 is known to have high transforming ability, while the P21 protein with Val at position 12 and Asn at position 16 is known to have poor transforming ability, similar to that of the normal protein (with Gly at 12 and Lys at 16.) The current results demonstrate a significant conformational change at position 15 induced by the substitution of Asn for Lys at position 16, which could explain this alteration in transformation potential. These findings are consistent with previous results suggesting the existence of a normal and a malignancy-causing conformation for the P21 proteins and suggest that the critical transforming region may encompass residues 12–15.     

Conformational effects of amino acid substitutions at positions 10, 12, and 13 in the P21 protein

Substitutions of amino acids for Gly 12 or Gly 13 in theras oncogene-encoded P21 proteins have been demonstrated to produce unique structural changes in these proteins that correlate with their ability to produce cell transformation. For example, the P21 proteins with Arg 12 or Val 13 are both known to be actively transforming. Recent site-specific mutagenesis experiments on the transforming Arg 12 protein have found that the substitution of Val for Gly 10 has no effect on transforming activity whereas the substitution of Val for Gly 13 led to a loss of transforming activity. In this study, we examine the structural effects of these substitutions on the amino terminal hydrophobic decapeptide (Leu 6-Gly 15) of P21 using conformational energy analysis. The results show that the transforming proteins with Gly 10 and Arg 12 or Val 10 and Arg 12 can both adopt the putative malignancy-causing conformation, whereas, for the nontransforming protein with Arg 12 and Val 13, this conformation is energetically disallowed. These results further support the theory that due to structural changes the transforming P21 proteins are unable to bind to some regulatory cellular element which may be the recently identified binding protein responsible for the induction of increased GTPase activity in normal P21 compared with transforming mutants.     

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.     

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 X-ray structure of native rubredoxin from Pyrococcus furiosus with the NMR structure of the zinc-substituted protein.

The three-dimensional X-ray structures of the oxidized and reduced forms of rubredoxin from Pyrococcus furiosus, determined at -161 degrees C, and the NMR structure of the zinc-substituted protein, determined in solution at 45 degrees C, are compared. The NMR and X-ray structures, which were determined independently, are very similar and lead to similar conclusions regarding the interactions that confer hyperthermostability.     

Comparison of the low energy conformations of an oncogenic and a non-oncogenic p21 protein, neither of which binds GTP or GDP

Oncogenic p21 protein, encoded by theras-oncogene, that causes malignant transformation of normal cells and many human tumors, is almost identical in sequence to its normal protooncogene-encoded counterpart protein, except for the substitution of arbitrary amino acids for the normally occurring amino acids at critical positions such as Gly 12 and Gin 61. Since p21 is normally activated by the binding of GTP in place of GDP, it has been postulated that oncogenic forms must retain bound GTP for prolonged time periods. However, two multiply substituted p21 proteins have been cloned, neither of which binds GDP or GTP. One of these mutant proteins with Val for Gly 10, Arg for Gly 12, and Thr for Ala 59 causes cell transformation, while the other, similar protein with Gly 10, Arg 12, Val for Gly 13 and Thr 59 does not transform cells. To define the critical conformational changes that occur in the p21 protein that cause it to become oncogenic, we have calculated the low energy conformations of the two multiply substituted mutant p21 proteins using a new adaptation of the electrostatically driven Monte Carlo (EDMC) technique, based on the program ECEPP. We have used this method to explore the conformational space available to both proteins and to compute the average structures for both using statistical mechanical averaging. Comparison of the average structures allows us to detect the major differences in conformation between the two proteins. Starting structures for each protein were calculated using the recently deposited x-ray crystal coordinates for the p21 protein, that was energy-refined using ECEPP, and then perturbed using the EDMC method to compute its average structure. The specific amino acid substitutions for both proteins were then generated into the lowest energy structure generated by this procedure, subjected to energy minimization and then to full EDMC perturbations. We find that both mutant proteins exhibit major differences in conformation in specific regions, viz., residues 35–47, 55–78, 81–93, 96–110, 115–126, and 123–134, compared with the EDMC-refined x-ray structure of the wild-type protein. These regions have been found to be the most flexible in the p21 protein bound to GDP from prior molecular dynamics calculations (Dykeset al., 1993). Comparison of the EDMC-average structure of the transforming mutant with that of the nontransforming mutant reveals major structural differences at residues 10–16, 32–40, and 60–68. These structural differences appear to be the ones that are critical in activation of the p21 protein. Analysis of the correlated motions of the different regions of the two mutant proteins reveals that changes in the conformation of regions in the carboxyl half of the protein are caused by changes in conformation around residues 10–16 and are transmitted by means of residues around Gln 61. The latter region therefore constitutes a molecular switch unit, in agreement with conclusions from prior work.On leave from the Department of Chemistry, University of Gdask, ul. Sobieskiego 18, 80-952 Gdask, Poland.     

Comparison of the NMR solution structure with the X-ray crystal structure of the activation domain from procarboxypeptidase B

Summary The NMR solution structure of the activation domain isolated from porcine procarboxypeptidase B is compared with the X-ray crystal structure of the corresponding segment in the intact proenzyme. For the region of the polypeptide chain that has a well-defined three-dimensional structure in solution, i.e., the backbone atoms of residues 11–76 and 25 amino acid side chains in this segment that form a hydrophobic core in the activation domain, the root-mean-square distance between the two structures is 1.1 Å. There are no significant differences in average atom positions between the two structures, but only the NMR structure shows increased structural disorder in three outlying loops located along the same edge of the activation domain. These regions of increased structural disorder in the free domain coincide only partially with the interface to the enzyme domain in the proenzyme.     

A predicted consensus structure for the N-Terminal fragment of the heat shock protein HSP90 family

A secondary structure has been predicted for the heat shock protein HSP90 family from an aligned set of homologous protein sequences by using a transparent method in both manual and automated implementation that extracts conformational information from patterns of variation and conservation within the family. No statistically significant sequence similarity relates this family to any protein with known crystal structure. However, the secondary structure prediction, together with the assignment of active site positions and possible biochemical properties, suggest that the fold is similar to that seen in N-terminal domain of DNA gyrase B (the ATPase fragment). Proteins 27:450–458, 1997. © 1997 Wiley-Liss, Inc.     

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on Cα atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 Å, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low. © 1995 Wiley-Liss, Inc.     

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein‐ligand crystal structures? This value is important biophysically, for better understanding the protein‐ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh‐resolution crystal structures that can be refined without restraints.     

Optimal ratio between the average conformational entropy and the average energy of interaction between residues for fast protein folding

Based on available experimental data and using a theoretical model of protein folding, we demonstrate that there is an optimal ratio between the average conformational entropy and the average contact energy per residue for fast protein folding. A statistical analysis of the conformational entropy and the number of contacts per residue for 5829 protein domains from four main classes (α, β, α/β, α+β) shows that each class has its own characteristic average number of contacts per residue and average conformational entropy per residue. These class-specific characteristics determine the protein folding rates: α-proteins are the fastest to fold, β-proteins are the second fastest, α+β-proteins are the third, and α/β-proteins are the last to fold.     

Purification and characterization of an activated form of the protein tyrosine kinase Lck from an Escherichia coli expression system

The lymphocyte-specific, nonreceptor protein tyrosine kinase Lck has been purified from an Escherichia coli expression system using a monoclonal antibody column followed by dye-affinity chromatography. Polyacrylamide gel electrophoretic analysis of purified protein revealed a single 56 kDa band, indicating that recombinant Lck was purified to near-homogeneity. The purified enzyme displayed tyrosine kinase activity as measured by both autophosphorylation and phosphorylation of exogenous substrates. Biochemical properties including protein phosphorylation and kinetic characteristics of the enzyme have been assessed. Peptide map analysis revealed that bacterially expressed Lck is phosphorylated predominantly on the autophosphorylation site (tyrosine-394), which is characteristic for activated protein tyrosine kinases. Indeed, we found that the recombinant enzyme is approximately fivefold more active than Lck from resting T cells, which is extensively phosphorylated at the regulatory carboxy-terminal tyrosine residue (tyrosine-505). Thus, we have overproduced recombinant human Lck in E. coli and developed a simple two-step purification procedure which yields highly active enzyme. This will enable the identification and characterization of potential regulators and targets of Lck and thereby greatly facilitate studies which will clarify its role in T cell signal transduction. © 1994 Wiley-Liss, Inc.     

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.     

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.     

Kinetic stability and crystal structure of the viral capsid protein SHP

SHP, the capsid-stabilizing protein of lambdoid phage 21, is highly resistant against denaturant-induced unfolding. We demonstrate that this high functional stability of SHP is due to a high kinetic stability with a half-life for unfolding of 25 days at zero denaturant, while the thermodynamic stability is not unusually high. Unfolding experiments demonstrated that the trimeric state (also observed in crystals and present on the phage capsid) of SHP is kinetically stable in solution, while the monomer intermediate unfolds very rapidly. We also determined the crystal structure of trimeric SHP at 1.5A resolution, which was compared to that of its functional homolog gpD. This explains how a tight network of H-bonds rigidifies crucial interpenetrating residues, leading to the observed extremely slow trimer dissociation or denaturation. Taken as a whole, our results provide molecular-level insights into natural strategies to achieve kinetic stability by taking advantage of protein oligomerization. Kinetic stability may be especially needed in phage capsids to allow survival in harsh environments.     

Backbone dynamics and solution structure refinement of the 15N-labeled human oncogenic protein p13MTCP1: Comparison with X-ray data

Two related oncogenes, TCL1 and MTCP1, are overexpressed in certain T-cell prolymphocytic leukemias as a result of chromosomal rearrangements that involve the translocation of one T-cell receptor gene to either chromosome 14q32 or Xq28, respectively. The human oncoprotein p13 ^MTCP1 is coded by the MTCP1 gene and its primary sequence is highly and only homologous to that of p14 ^TCL1 , the product of TCL1. These two proteins likely represent the first members of a new family of oncogenic proteins. A previous model of the three-dimensional solution structure of p13 ^MTCP1 was determined recently using exclusively homonuclear proton two-dimensional NMR methods and, almost simultaneously, high-resolution crystal structures of p13 ^MTCP1 and p14 ^TCL1 appeared in the literature. In order to gain more insight into the details of the solution structure, we uniformly labeled p13 ^MTCP1 with nitrogen-15. The refined structure benefits from 520 additional NOEs, extracted from either ¹⁵N-edited 3D experiments or homonuclear 2D NOESY recorded at 800 MHz, and from a nearly complete set of angular restraints. Measurements of ¹⁵N spin relaxation times and heteronuclear ¹⁵N{¹H}NOEs at two magnetic field strengths provided additional insights into the dynamics of the protein backbone. On the basis of these new results, a putative binding surface for this particular class of oncogenes is discussed.     

Generation of an antibody with a designed specificity difference for protein C and activated protein C

Rabbit polyclonal antibodies to a synthetic peptide, NH₂-Asp-Thr-Asn-Gln-Val-Asp-Gln-Lys-Asp-Gln-Leu-Asp-Phe-Arg-CONH₂ (APep), have been produced. This sequence is identical to that contained in the tetradecapeptide released from bovine protein C (PC) as a result of its conversion to its activated form (APC), except that Phe₁₃ replaced the normal Pro₁₃, in order to discourage cross-reactivity of antibodies to the carboxylterminal portion of APep with PC. The antibody pool obtained reacted with PC and showed virtually no cross-reactivity toward either APC or several typical plasma proteins. This general approach should serve well as a means of production of antibodies with a designed specificity capable of distinguishing between forms of the same protein that arise by release of peptide material.