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

Direct inhibition of oncogenic KRAS by Bacillus pumilus ribonuclease (binase)

RAS proteins function as molecular switches that transmit signals from cell surface receptors into specific cellular responses via activation of defined signaling pathways (Fang, 2015). Aberrant constitutive RAS activation occurs with high incidence in different types of cancer (Bos, 1989). Thus, inhibition of RAS-mediated signaling is extremely important for therapeutic approaches against cancer. Here we showed that the ribonuclease (RNase) binase, directly interacts with endogenous KRAS. Further, molecular structure models suggested an inhibitory nature of binase-RAS interaction involving regions of RAS that are important for different aspects of its function. Consistent with these models, phosphorylation analysis of effectors of RAS-mediated signaling revealed that binase inhibits the MAPK/ERK signaling pathway. Interestingly, RAS activation assays using a non-hydrolysable GTP analog (GTPγS) demonstrated that binase interferes with the exchange of GDP by GTP. Furthermore, we showed that binase reduced the interaction of RAS with the guanine nucleotide exchange factor (GEF), SOS1. Our data support a model in which binase-KRAS interaction interferes with the function of GEFs and stabilizes the inactive GDP-bound conformation of RAS thereby inhibiting MAPK/ERK signaling. This model plausibly explains the previously reported, antitumor-effect of binase specific towards RAS-transformed cells and suggests the development of anticancer therapies based on this ribonuclease.

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 the ras-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 Gln 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 (Dykes et 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.(ABSTRACT TRUNCATED AT 400 WORDS)

Activated conformations of the ras-gene-encoded p21 protein. 2. Comparison of the computed and high-resolution x-ray crystallographic structures of Gly-12 p21

The ras-oncogene-encoded p21 protein is a G-protein that has been shown to cause the malignant transformation of normal cells and has been implicated in causing human tumors. p21 is thought to be activated by the binding of GTP in place of GDP to the protein. We have previously constructed the three-dimensional structure of the p21 protein bound to GDP from an available alpha-carbon tracing of this protein using a combination of molecular dynamics and energy minimization (Dykes, et al., J. Biomol. Struct. Dynamics, 9:1025-1044). Until the recent publication of the all-heavy-atom x-ray crystallographic molecular coordinates of p21 residues 1-166 bound to a non-hydrolyzable GTP derivative (GppNp), no all-atom structure of the p21 protein has been available in the Brookhaven National Laboratories Protein Data Bank (PDB). In this communication we compare our computed structure for the p21-GDP complex to this x-ray crystal structure. We find that the two structures agree quite closely with one another, the overall RMS deviation for the backbone being 1.47 A and 2.71 A for all of the atoms. We have identified the regions of the protein that are responsible for the most significant deviations between the two structures, i.e., residues 32-40 and 61-74. We find that the main chain in the 32-40 segment deviates significantly from residue 32 to residue 36 and the side chain phenolic rings of residue 32 differ greatly between the two structures. The 61-74 region is the least-well defined region in the whole protein crystallographically having, by far, the highest temperature factor (B-factor). The backbone and side chain conformations in the 61-74 segment differ markedly, the overall RMS deviation being 3.1 A for the backbone and 5.7 A for all atoms. Both of these regions have been found in x-ray crystallographic studies of p21-GDP and p21-GTP complexes to undergo significant changes in conformation upon the binding of GTP in place of GDP to the protein. We have further compared our computed structure of the p21 protein with the x-ray crystal structure with regard to the conformations of individual segments, in particular, structurally conserved sequences (SCR), i.e., those sequences that have structural and sequence homology to corresponding sequences in the related G-protein, bacterial elongation factor Tu(EF-Tu), and variable loop regions.(ABSTRACT TRUNCATED AT 400 WORDS)

Activated conformations of the ras-gene-encoded p21 protein. 1. An energy-refined structure for the normal p21 protein complexed with GDP

A complete three-dimensional structure for the ras-gene-encoded p21 protein with Gly 12 and Gln 61, bound to GDP, has been constructed in four stages using the available alpha-carbon coordinates as deposited in the Brookhaven National Laboratories Protein Data Bank. No all-atom structure has been made available despite the fact that the first crystallographic structure for the p21 protein was reported almost four years ago. In the p21 protein, if amino acid substitutions are made at any one of a number of different positions in the amino acid sequence, the protein becomes permanently activated and causes malignant transformation of normal cells or, in some cell lines, differentiation and maturation. For example, all amino acids except Gly and Pro at position 12 result in an oncogenic protein; all amino acids except Gln, Glu and Pro at position 61 likewise cause malignant transformation of cells. We have constructed our all-atom structure of the non-oncogenic protein from the x-ray structure in order to determine how oncogenic amino acid substitutions affect the three-dimensional structure of this protein. In Stage 1 we generated a poly-alanine backbone (except at Gly and Pro residues) through the alpha-carbon structure, requiring the individual Ala, Pro or Gly residues to conform to standard amino acid geometry and to form trans-planar peptide bonds. Since no alpha-carbon coordinates for residues 60-65 have been determined, these residues were modeled by generating them in the extended conformation and then subjecting them to molecular dynamics using the computer application DISCOVER and energy minimization using DISCOVER and the ECEPP (Empirical Conformational Energies for Peptides Program). In Stage 2, the positions of residues that are homologous to corresponding residues of bacterial elongation factor Tu (EF-Tu) to which p21 bears an overall 40% sequence homology, were determined from their corresponding positions in a high-resolution structure of EF-Tu. Non-homologous loops were taken from the structure generated in Stage 1 and were placed between the appropriate homologous segments so as to connect them. In Stage 3, all bad contacts that occurred in this resulting structure were removed, and the coordinates of the alpha-carbon atoms were forced to superimpose as closely as possible on the corresponding atoms of the reference (x-ray) structure. Then the side chain positions of residues of the non-homologous loop regions were modeled using a combination of molecular dynamics and energy minimization using DISCOVER and ECEPP respectively.(ABSTRACT TRUNCATED AT 400 WORDS)