Molecular dynamics simulation of Escherichia coli dihydrofolate reductase and its protein fragments: relative stabilities in experiment and simulations

We have carried out molecular dynamics simulations of the native dihydrofolate reductase from Escherichia coli and several of its folded protein fragments at standard temperature. The simulations have shown fragments 1--36, 37--88, and 89--159 to be unstable, with a C(alpha)RMSD (C(alpha) root mean squared deviation) >5 A after 3.0 nsec of simulation. The unfolding of fragment 1--36 was immediate, whereas fragments 37--88 and 89--159 gradually unfolded because of the presence of the beta-sheet core structure. In the absence of residues 1--36, the two distinct domains comprising fragment 39--159 associated with each other, resulting in a stable conformation. This conformation retained most of its native structural elements. We have further simulated fragments derived from computational protein cutting. These were also found to be unstable, with the exception of fragment 104--159. In the absence of alpha(4), the loose loop region of residues 120--127 exhibited a beta-strand-like behavior, associating itself with the beta-sheet core of the protein fragment. The current study suggests that the folding of dihydrofolate reductase involves cooperative folding of distinct domains which otherwise would have been unstable as independent folded units in solution. Finally, the critical role of residues 1--36 in allowing the two distinct domains of fragment 104--159 to fold into the final native conformation is discussed.     

Identifying the structural boundaries of independent folding domains in the alpha subunit of tryptophan synthase, a beta/alpha barrel protein.

Two equilibrium intermediates have previously been observed in the urea denaturation of the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli, an eight-stranded beta/alpha barrel protein. In the current study, a series of amino-terminal fragments were characterized to probe the elementary folding units that may be in part responsible for this complex behavior. Stop-codon mutagenesis was used to produce eight fragments ranging in size from 105-214 residues and containing incremental elements of secondary structure. Equilibrium studies by circular dichroism indicate that all of these fragments are capable of adopting secondary structure. All except for the shortest fragment fold cooperatively. The addition of the fourth, sixth, and eighth beta-strands leads to distinct increases in structure, cooperativity, and/or stability, suggesting that folding involves the modular assembly of betaalphabeta supersecondary structural elements. One-dimensional NMR titrations at high concentrations of urea, probing the environment around His92, were also performed to test for the presence of residual structure in the fragments. All fragments that contained the first four betaalpha units of structure exhibited a cooperative unfolding transition at high concentrations of urea with significant but reduced stability relative to the full-length protein. These results suggest that the residual structure in alphaTS requires the participation of hydrophobic residues in multiple beta-strands that span the entire sequence.     

Secondary and tertiary structural changes in gamma delta resolvase: comparison of the wild-type enzyme, the I110R mutant, and the C-terminal DNA binding domain in solution.

gamma delta Resolvase is a site-specific DNA recombinase (M(r) 20.5 kDa) in Escherichia coli that shares homology with a family of bacterial resolvases and invertases. We have characterized the secondary and tertiary structural behavior of the cloned DNA binding domain (DBD) and a dimerization defective mutant in solution. Low-salt conditions were found to destabilize the tertiary structure of the DBD dramatically, with concomitant changes in the secondary structure that were localized near the hinge regions between the helices. The molten tertiary fold appears to contribute significantly to productive DNA interactions and supports a mechanism of DNA-induced folding of the tertiary structure, a process that enables the DBD to adapt in conformation for each of the three imperfect palindromic sites. At high salt concentrations, the monomeric I110R resolvase shows a minimal perturbation to the three helices of the DBD structure and changes in the linker segment in comparison to the cloned DBD containing the linker. Comparative analysis of the NMR spectra suggest that the I110R mutant contains a folded catalytic core of approximately 60 residues and that the segment from residues 100 to 149 are devoid of regular structure in the I110R resolvase. No increase in the helicity of the linker region of I110R resolvase occurs on binding DNA. These results support a subunit rotation model of strand exchange that involves the partial unfolding of the catalytic domains.     

p53 contains large unstructured regions in its native state

The human tumor suppressor protein p53 is understood only to some extent on a structural level. We performed a comprehensive biochemical and biophysical structure-function analysis of p53 full-length protein and p53 fragments. The analysis showed that p53 and the fragments investigated form stable functional units. Full-length p53 and the tetrameric fragment N93p53 (residues 93-393) are, however, destabilized significantly compared to the monomeric core domain (residues 94-312) and the monomeric fragment p53C312 (residues 1-312). At the physiological temperature of 37 degrees C and in the absence of modifications or stabilizing partners, wild-type p53 is more than 50% unfolded correlating with a 75% loss in DNA-binding activity. Furthermore the analysis of CD spectra revealed that full-length p53 contains large unstructured regions in its N and C-terminal parts. Our results indicate that full-length p53 is a modular protein consisting of defined structured and unstructured regions. We propose that p53 belongs to the growing family of loosely folded or partially unstructured native proteins. The lack of a rigid structure combined with the low overall stability may allow the physiological interaction of p53 with a multitude of partner proteins and the regulation of its turnover.     

Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.     

Benchmarking of TASSER in the ab initio limit

A significant number of protein sequences in a given proteome have no obvious evolutionarily related protein in the database of solved protein structures, the PDB. Under these conditions, ab initio or template-free modeling methods are the sole means of predicting protein structure. To assess its expected performance on proteomes, the TASSER structure prediction algorithm is benchmarked in the ab initio limit on a representative set of 1129 nonhomologous sequences ranging from 40 to 200 residues that cover the PDB at 30% sequence identity and which adopt alpha, alpha + beta, and beta secondary structures. For sequences in the 40-100 (100-200) residue range, as assessed by their root mean square deviation from native, RMSD, the best of the top five ranked models of TASSER has a global fold that is significantly close to the native structure for 25% (16%) of the sequences, and with a correct identification of the structure of the protein core for 59% (36%). In the absence of a native structure, the structural similarity among the top five ranked models is a moderately reliable predictor of folding accuracy. If we classify the sequences according to their secondary structure content, then 64% (36%) of alpha, 43% (24%) of alpha + beta, and 20% (12%) of beta sequences in the 40-100 (100-200) residue range have a significant TM-score (TM-score > or = 0.4). TASSER performs best on helical proteins because there are less secondary structural elements to arrange in a helical protein than in a beta protein of equal length, since the average length of a helix is longer than that of a strand. In addition, helical proteins have shorter loops and dangling tails. If we exclude these flexible fragments, then TASSER has similar accuracy for sequences containing the same number of secondary structural elements, irrespective of whether they are helices and/or strands. Thus, it is the effective configurational entropy of the protein that dictates the average likelihood of correctly arranging the secondary structure elements.     

Autonomous folding of the excised coenzyme-binding domain of D-glyceraldehyde 3-phosphate dehydrogenase from Thermotoga maritima.

An important question in protein folding is whether compact substructures or domains are autonomous units of folding and assembly. The protomer of the tetrameric D-glyceraldehyde-3-phosphate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima has a complex coenzyme-binding domain, in which residues 1-146 form a compact substructure with the last 31 residues (313-333). Here it is shown that the gene of a single-chain protein can be expressed in Escherichia coli after deleting the 163 codons corresponding to the interspersed catalytic domain (150-312). The purified gene product is a soluble, monomeric protein that binds both NAD+ and NADH strongly and possesses the same unfolding transition induced by guanidinium chloride as the native tetramer. The autonomous folding of the coenzyme-binding domain has interesting implications for the folding, assembly, function, and evolution of the native enzyme.     

How general is the nucleation–condensation mechanism?

We investigate the structures of the major folding transition states of nine proteins by correlation of published Phi-values with inter-residue contact maps. Combined with previous studies on six proteins, the analysis suggests that at least 10 of the 15 small globular proteins fold via a nucleation-condensation mechanism with a concurrent build-up of secondary and tertiary structure contacts, but a structural consolidation that is clearly nonuniformly distributed over the molecule and most intense in a single structural region suggesting the occurrence of a single folding nucleus. However, on average helix- and sheet-forming residues show somewhat larger Phi-values in the major transition state, suggesting that secondary structure formation is one important driving force in the nucleation-condensation in many proteins and that secondary-structure forming residues tend to be more prominent in folding nuclei. We synthesize the combined information on these 10 of 15 proteins into a unified nucleation-condensation mechanism which also accounts for effects described by the framework, hydrophobic collapse, zipper, and funnel models.     

Testing the role of chain connectivity on the stability and structure of dihydrofolate reductase from E. coli: fragment complementation and circular permutation reveal stable, alternatively folded forms

The effects of chain cleavage and circular permutation on the structure, stability, and activity of dihydrofolate reductase (DHFR) from Escherichia coli were investigated by various spectroscopic and biochemical methods. Cleavage of the backbone after position 86 resulted in two fragments, (1--86) and (87--159) each of which are poorly structured and enzymatically inactive. When combined in a 1 : 1 molar ratio, however, the fragments formed a high-affinity (K(a) = 2.6 x 10(7) M(-1)) complex that displays a weakly cooperative urea-induced unfolding transition at micromolar concentrations. The retention of about 15% of the enzymatic activity of full-length DHFR is surprising, considering that the secondary structure in the complex is substantially reduced from its wild-type counterpart. In contrast, a circularly permuted form with its N-terminus at position 86 has similar overall stability to full-length DHFR, about 50% of its activity, substantial secondary structure, altered side-chain packing in the adenosine binding domain, and unfolds via an equilibrium intermediate not observed in the wild-type protein. After addition of ligand or the tight-binding inhibitor methotrexate, both the fragment complex and the circular permutant adopt more native-like secondary and tertiary structures. These results show that changes in the backbone connectivity can produce alternatively folded forms and highlight the importance of protein-ligand interactions in stabilizing the active site architecture of DHFR.     

Protein structure prediction enhanced with evolutionary diversity: SPEED

For naturally occurring proteins, similar sequence implies similar structure. Consequently, multiple sequence alignments (MSAs) often are used in template‐based modeling of protein structure and have been incorporated into fragment‐based assembly methods. Our previous homology‐free structure prediction study introduced an algorithm that mimics the folding pathway by coupling the formation of secondary and tertiary structure. Moves in the Monte Carlo procedure involve only a change in a single pair of ?,ψ backbone dihedral angles that are obtained from a Protein Data Bank‐based distribution appropriate for each amino acid, conditional on the type and conformation of the flanking residues. We improve this method by using MSAs to enrich the sampling distribution, but in a manner that does not require structural knowledge of any protein sequence (i.e., not homologous fragment insertion). In combination with other tools, including clustering and refinement, the accuracies of the predicted secondary and tertiary structures are substantially improved and a global and position‐resolved measure of confidence is introduced for the accuracy of the predictions. Performance of the method in the Critical Assessment of Structure Prediction (CASP8) is discussed.     

Alpha-hairpin stability and folding of transmembrane segments

Molecular Dynamics (MD) simulations at low dielectric constant have been carried out for peptides matching the double spanning segments of transmembrane proteins. Different folding dynamics have been observed. The peptides folded into the stable helix-turn-helix conformation-alpha-hairpin-with antiparallel-oriented strands or unstable alpha-hairpin conformation that unfolded later into the straight helical structure. The peptide having flexible residues in the TM helices often misfolded into a tangled structure that can be avoided by restricting the flexibility of these residues. General conclusions can be drawn from the observed folding dynamics. The stability and folding of some double spanning transmembrane fragments are self-assembling. The following and/or neighboring peptide chains of the protein may support the stability of the hairpin structure of other fragments. The stability of the TM helices containing flexible residues could be maintained due to contacts with neighboring TM segments.     

Proton NMR and photochemically induced dynamic nuclear polarization studies of peptide fragments obtained by controlled proteolysis of mouse epidermal growth factor

Controlled proteolysis of epidermal growth factor from the mouse leads to fragments of mouse epidermal growth factor containing residues 1-48 and 1-45. The COOH-terminal pentapeptide appears to play a crucial role in determining the hydrophobic interactions between the hormone and the stationary phase during gel chromatography on TSK-125 gel. Proton NMR studies indicate that the overall structure of mouse epidermal growth factor is retained in the protein devoid of the COOH-terminal pentapeptide, while subsequent cleavage of the peptide bond between Arg-45 and Asp-46 starts to perturb the proton resonances most characteristic of the tertiary structure of the hormone, especially those from the aromatic ring protons of Tyr-37. Consequently, photochemically induced dynamic nuclear polarization experiments show an increased exposure of Tyr-37 in the fragment of mouse epidermal growth factor containing residues 1-48. Nuclear Overhauser data suggest that structural changes do occur on fragmentation but seem to be localized in the tiered-beta-sheet domain which contains Tyr-37.     

Autonomous folding of a C-terminal inhibitory fragment of Escherichia coli isoleucine-tRNA synthetase.

We previously reported that C-terminal fragments of Escherichia coli Ile-tRNA synthetase, a monomeric enzyme of 939 amino acids, act as dominant negative inhibitors of the wild-type enzyme in vivo and in vitro. Our experiments suggested that it is possible to block the functional assembly of a monomeric protein by interfering with the folding pathway. We postulated that the inhibitory C-terminal fragments fold autonomously, and in the presence of full-length Ile-tRNA synthetase, trap the N-terminal portion of polypeptide in an unproductive complex. Here, we report the results of experiments aimed at understanding the mechanism of dominant negative inhibition. We have carried out biophysical experiments on fragment 585-939 of Ile-tRNA synthetase, which we previously determined to be the minimal inhibitory unit. Circular dichroism and fluorescence spectroscopy indicate that this fragment forms a compact and stable structure in solution. The secondary structure of this fragment is predominantly alpha-helical, consistent with the crystal structure of Ile-tRNA synthetase from another organism. The C-terminal fragment is capable of forming native-like secondary and tertiary structure after refolding from guanidine HCl. Taken together, the results are consistent with the hypothesis that the inhibitory fragment of Ile-tRNA synthetase forms an independent folding unit.     

Acquisition of native-like interactions in C-terminal fragments of barnase

We have characterised a series of C-terminal fragments of barnase by different biophysical techniques to find out when they acquire secondary and tertiary native-like structure. Fragments B96-110 (which comprises the last 15 residues of the intact protein) up to B37-110 (which involves most of the protein except the two first helices and a loop) were mainly disordered. Only fragment B23-110, which lacks alpha-helix1, showed native-like near and far-UV CD and fluorescence spectra. The intensities of these spectra were lower than those of the full-length protein, which suggests the absence of complete side-chain packing. Urea denaturation followed by fluorescence, far-UV CD and gel-filtration chromatography techniques indicated a co-operative transition only for B23-110. None of the fragments melted co-operatively with temperature. Thus, the formation of secondary and tertiary structure requires most of the polypeptide chain to be present, that is, secondary and tertiary structure are formed in parallel. This agrees with the proposed model for barnase folding, where the residual structure in small fragments is weak and flickering, and it is only consolidated when there are enough tertiary interactions. Thus, the development of structure in the series of C-terminal fragments follows a similar behaviour to that observed in the series of N-terminal fragments of barnase.     

Hierarchic organization of globular proteins: Experimental studies on thermolysin

Previous studies from this laboratory have shown that the thermolysin fragment 121–316, comprising entirely the“all-α” COOH-terminal structural domain 158–316, as well as fragment 206–316 (fragment FII) are able to refold into a native-like, stable structure independently from the rest of the protein molecule. The present report describes conformational properties of fragments 228–316 and 255–316 obtained by chemical and enzymatic cleavage of fragment FII, respectively. These subfragments are able to acquire a stable conformation of native-like characteristics, as judged by quantitative analysis of secondary structure from far-ultra-violet circular dichroism spectra and immunochemical properties using rabbit anti-thermolysin antibodies. Melting curves of the secondary structure of the fragments show cooperativity with a temperature of half-denaturationT _mof 65–66°C. The results of this study provide evidence that it is possible to isolate stable supersecondary structures (folding units) of globular proteins and correlate well with predictions of subdomains of the COOH-terminal structural domain 158–316 of thermolysin.     

Nicked apomyoglobin: a noncovalent complex of two polypeptide fragments comprising the entire protein chain

Limited proteolysis of the 153-residue chain of horse apomyoglobin (apoMb) by thermolysin results in the selective cleavage of the peptide bond Pro88-Leu89. The N-terminal (residues 1-88) and C-terminal (residues 89-153) fragments of apoMb were isolated to homogeneity and their conformational and association properties investigated in detail. Far-UV circular dichroism (CD) measurements revealed that both fragments in isolation acquire a high content of helical secondary structure, while near-UV CD indicated the absence of tertiary structure. A 1:1 mixture of the fragments leads to a tight noncovalent protein complex (1-88/89-153, nicked apoMb), characterized by secondary and tertiary structures similar to those of intact apoMb. The apoMb complex binds heme in a nativelike manner, as given by CD measurements in the Soret region. Second-derivative absorption spectra in the 250-300 nm region provided evidence that the degree of exposure of Tyr residues in the nicked species is similar to that of the intact protein at neutral pH. Also, the microenvironment of Trp residues, located in positions 7 and 14 of the 153-residue chain of the protein, is similar in both protein species, as given by fluorescence emission data. Moreover, in analogy to intact apoMb, the nicked protein binds the hydrophobic dye 1-anilinonaphthalene-8-sulfonate (ANS). Taken together, our results indicate that the two proteolytic fragments 1-88 and 89-153 of apoMb adopt partly folded states characterized by sufficiently nativelike conformational features that promote their specific association and mutual stabilization into a nicked protein species much resembling in its structural features intact apoMb. It is suggested that the formation of a noncovalent complex upon fragment complementation can mimic the protein folding process of the entire protein chain, with the difference that the folding of the complementary fragments is an intermolecular process. In particular, this study emphasizes the importance of interactions between marginally stable elements of secondary structure in promoting the tertiary contacts of a native protein. Considering that apoMb has been extensively used as a paradigm in protein folding studies for the past few decades, the novel fragment complementing system of apoMb here described appears to be very useful for investigating the initial as well as late events in protein folding.     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

Tryptophan replacements in the trp aporepressor from Escherichia coli: probing the equilibrium and kinetic folding models.

Mutants of the dimeric Escherichia coli trp aporepressor are constructed by replacement of the two tryptophan residues in each subunit in order to assess the effects on equilibrium and kinetic fluorescence properties of the folding reaction. The three kinetic phases detected by intrinsic tryptophan fluorescence in refolding of the wild-type aporepressor are also observed in folding of both Trp 19 to Phe and Trp 99 to Phe single mutants, demonstrating that these phases correspond to global rather than local conformational changes. Comparison of equilibrium fluorescence (Royer, C.A., Mann, C.J., & Matthews, C.R., 1993, Protein Sci. 2, 1844-1852) and circular dichroism transition curves induced by urea shows that replacement of either Trp 19 or Trp 99 results in noncoincident behavior. Unlike the wild-type protein (Gittelman, M.S. & Matthews, C.R., 1990, Biochemistry 29, 7011-7020), tertiary and/or quaternary structures are disrupted at lower denaturant concentration than is secondary structure. The equilibrium results can be interpreted in terms of enhancement in the population of a monomeric folding intermediate in which the lone tryptophan residue is highly exposed to solvent, but in which substantial secondary structure is retained. The location of both mutations at the interface between the two subunits (Zhang, R.G., et al., 1987, Nature 327, 591-597) provides a simple explanation for this phenomenon.     

Protein anatomy: structure and function of peptide fragments corresponding to the secondary structure units of barnase

Globular proteins can be decomposed into several modules or secondary structure units. It is useful to investigate the functions of such structural units in order to understand the folding units of proteins. In our previous work, barnase was divided into six peptide fragments corresponding to modules, and some of them were shown to have RNA-binding and RNase activity [Yanagawa, et al. (1993) J. Biol. Chem. 268, 5861-5865]. Barnase mutant proteins obtained by permutation of the structural units also had RNase activity [Tsuji, T. et al. (1999) J. Mol. Biol. 286, 1581-1596]. Here we investigated the structure and function of peptide fragments corresponding to secondary structure units of barnase. The results of circular dichroism spectroscopy indicated that some of the peptide fragments form helical structures in aqueous solutions containing over 30% 2,2,2-trifluoroethanol, and the S6 (94-110) peptide fragment is induced to form a beta-sheet structure in the presence of RNA. The S6 peptide fragment forms aggregate complexes with RNA. Electron microscopic analysis showed that the aggregate complexes were comprised of filaments. These results indicate that not only modules but also secondary structure units dissected from a globular protein have functional and structure-forming capabilities.