Regina Monaco

Identification of a glutathione-S-transferase effector domain for inhibition of jun kinase, by molecular dynamics.

We have recently found that the glutathione-S-transferase π-isozyme (GST-π), a cellular detoxification enzyme, potently and selectively inhibits activation of jun protein by its upstream kinase, jun kinase (JNK). This newly identified regulatory activity of GST-π is strongly inhibited by a group of agents that inhibit its enzymatic activity. Since loss of enzymatic activity in general does not correlate with loss of regulatory activity, it is likely that inhibitor binding induces changes in the structure of one or more domains of GST that block its interaction with JNK. To identify regions of GST that change conformation on the binding of inhibitors, we have performed molecular dynamics calculations on GST-π to compute its average structure in the presence and absence of the inhibitor, glutathione sulfonate. Superposition of the two average structures reveals that several regions change local structure depending upon whether the inhibitor is bound or not bound. Two of these regions, residues 36–50 and 194–201, are highly exposed. We have synthesized peptides corresponding to these two segments and find that the 194–201 sequence strongly inhibits the ability of GST-π to block the in vitro phosphorylation of jun by JNK. These results suggest that this region of GST-π is critical to its functioning as a newly discovered regulator of signal transduction.

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.

Inhibition of Oncogenic and Activated Wild-Type ras-p21 Protein-Induced Oocyte Maturation by Peptides from the ras-Binding Domain of the raf-p74 Protein, Identified from Molecular Dynamics Calculations

In the preceding paper we found from molecular dynamics calculations that the structure of the ras-binding domain (RBD) of raf changes predominantly in three regions depending upon whether it binds to ras-p21 protein or to its inhibitor protein, rap-1A. These three regions of the RBD involve residues from the protein–protein interaction interface, e.g., between residues 60 and 72, residues 97–110, and 111–121. Since the rap-1A–RBD complex is inactive, these three regions are implicated in ras-p21-induced activation of raf. We have therefore co-microinjected peptides corresponding to these three regions, 62–76, 97–110, and 111–121, into oocytes with oncogenic p21 and microinjected them into oocytes incubated in in insulin, which activates normal p2l. All three peptides, but not a control peptide, strongly inhibit both oncogenic p21- and insulin-induced oocyte maturation. These findings corroborate our conclusions from the theoretical results that these three regions constitute raf effector domains. Since the 97–110 peptide is the strongest inhibitor of oncogenic p21, while the 111–121 peptide is the strongest inhibitor of insulin-induced oocyte maturation, the possibility exists that oncogenic and activated normal p21 proteins interact differently with the RBD of raf.

Comparison of the Computed Three-Dimensional Structures of Oncogenic Forms (Bound to GDP) of the ras-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.

Conformation of the transmembrane domain of the c-erbB-2 oncogene-encoded protein in its monomeric and dimeric states

The human c-erbB-2 oncogene is homologous to the ratneu oncogene, both encoding transmembrane growth factor receptors. Overexpression and point mutations in the transmembrane domain of the encoded proteins in both cases have been implicated in cell transformation and carcinogenesis. In the case of theneu protein, it has been proposed that these effects are mediated by conformational preferences for anα-helix in the transmembrane domain, which facilitates receptor dimerization, an important step in the signal transduction process. To examine whether this is the case for c-erbB-2 as well, we have used conformational energy analysis to determine the preferred three-dimensional structures for the transmembrane domain of the c-erbB-2 protein from residues 650 to 668 with Val (nontransforming) and Glu (transforming) at position 659. The global minimum energy conformation for the Val-659 peptide from the normal, nontransforming protein was found to contain several bends, whereas the global minimum energy conformation for Glu-659 peptide from the mutant, transforming protein was found to beα-helical. Thus, the difference in conformational preferences for these transmembrane domains may explain the difference in transforming ability of these proteins. The presence of higher-energyα-helical conformations for the transmembrane domain from the normal Val-659 protein may provide an explanation for the presence of a transforming effect from overexpression of c-erbB-2. In addition, docking of the oncogenic sequences in theirα-helical and bend conformations shows that the all-α-helical dimer is clearly favored energetically over the bend dimer.

Structural effects of the binding of GTP to the wild-type and oncogenic forms of the ras-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.

Molecular Dynamics Analysis of the Structures of ras-Guanine Nucleotide Exchange Protein (SOS) Bound to Wild-Type and Oncogenic ras-p21. Identification of Effector Domains of SOS

The X-ray crystal structure of the ras oncogene-encoded p21 protein bound to SOS, the guanine nucleotide exchange-promoting protein, has been determined. We have undertaken to determine if there are differences between the three-dimensional structures of SOS bound to normal and oncogenic (Val 12-p21) proteins. Using molecular dynamics, we have computed the average structures for both complexes and superimposed them. We find four domains of SOS that differ markedly in structure: 631–641, 676–691, 718–729, and 994–1004. Peptides corresponding to these sequences have been synthesized and found to be powerful modulators of oncogenic p21 in cells as described in an accompanying paper. We find that the SOS segment from 809–815 makes contacts with multiple domains of ras-p21 and can facilitate correlated conformational changes in these domains.

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.

Dynamics Simulation of the Interaction Between the Novel Intercalator Diethidium Cation and B-Form DNA

Previous research has described the interaction between the novel molecule diethidium (2,7-diamino 9-[2,7 diamino 10-N-phenanthridium] 10-N-phenanthridium) (Figure 1) and B-form DNA. Our goal is the elucidation of diethidium as the first member of a novel class of drugs which are potential pharmaceutical agents. This class of potential drug molecules differs from previously known intercalators in the following ways: a) Its structure, that of two perpendicular planes, each known to have excellent intercalation properties, is novel b) Unlike known bis-intercalators, the linker region length in diethidium is zero c) The geometry of the drug matches the geometry of the space available in the major groove d) The drug is shown to cause some vectorial disruption of DNA. For this paper, we have performed a series of 200 picosecond dynamics simulations on the complex formed between diethidium in the major groove and a dodecarner of double-stranded B-form DNA, CGCGAATTCGCG, and have shown that this complex has a intricate interaction. The DNA dodecamer is found to be in an intermediate A-B state, but, even in simulations as long as 1 nanosecond, the drug does not back-out or otherwise leave the intercalation site. The drug is found to be mobile within the intercalation site on timescales longer than 1 nanoscale. The mobility of the drug within the intercalation site has been predicted by our previous energy minimization studies.

Computed three-dimensional structures for the ras-binding domain of the raf-p74 protein complexed with ras-p21 and with its suppressor protein, rap-1A.

The three-dimensional structures of theras-p21 protein and its protein inhibitor, rap-1A, have been computed bound to theras-binding domain, RBD (residues 55–131), of theraf-p74 protein, a critical target protein ofras-p21 in theras-induced mitogenic signal transduction pathway. The coordinates of RBD have been reconstructed from the stereoview of an X-ray crystal structure of this domain bound to rap-1A and have been subjected to energy minimization. The energy-minimized structures of bothras- p21 and rap-1A, obtained in previous studies, have been docked against RBD, using the stereo figure of the RBD-rap-1A complex, based on a six-step procedure. The final energy-minimized structure of rap-1A-RBD is identical to the X-ray crystal structure. Comparison of theras-p21- and rap-1A-RBD complexes reveals differences in the structures of effector domains ofras-p21 and rap-1a, including residues 32–47, a domain that directly interacts with RBD, 60–66, 96–110, involved in the interaction ofras-p21 withjun kinase (JNK) andjun protein, and 115–126, involved in the interaction of p21 with JNK. The structure of the RBD remained the same in both complexes with the exception of small deviations in itsβ-2 binding loop (residues 63–71) and residues 89–91, also involved in binding to rap-1A. The results suggest that the binding of these two proteins to RBD may allow them to interact with other cellular target proteins such as JNK andjun.