Cross-validated maximum likelihood enhances crystallographic simulated annealing refinement

Recently, the target function for crystallographic refinement has been improved through a maximum likelihood analysis, which makes proper allowance for the effects of data quality, model errors, and incompleteness. The maximum likelihood target reduces the significance of false local minima during the refinement process, but it does not completely eliminate them, necessitating the use of stochastic optimization methods such as simulated annealing for poor initial models. It is shown that the combination of maximum likelihood with cross-validation, which reduces overfitting, and simulated annealing by torsion angle molecular dynamics, which simplifies the conformational search problem, results in a major improvement of the radius of convergence of refinement and the accuracy of the refined structure. Torsion angle molecular dynamics and the maximum likelihood target function interact synergistically, the combination of both methods being significantly more powerful than each method individually. This is demonstrated in realistic test cases at two typical minimum Bragg spacings (dmin = 2.0 and 2.8 Å, respectively), illustrating the broad applicability of the combined method. In an application to the refinement of a new crystal structure, the combined method automatically corrected a mistraced loop in a poor initial model, moving the backbone by 4 Å.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

DNA annealing by Rad52 Protein is stimulated by specific interaction with the complex of replication protein A and single-stranded DNA

Homologous recombination in Saccharomyces cerevisiae depends critically on RAD52 function. In vitro, Rad52 protein preferentially binds single-stranded DNA (ssDNA), mediates annealing of complementary ssDNA, and stimulates Rad51 protein-mediated DNA strand exchange. Replication protein A (RPA) is a ssDNA-binding protein that is also crucial to the recombination process. Herein we report that Rad52 protein effects the annealing of RPA–ssDNA complexes, complexes that are otherwise unable to anneal. The ability of Rad52 protein to promote annealing depends on both the type of ssDNA substrate and ssDNA binding protein. RPA allows, but slows, Rad52 protein-mediated annealing of oligonucleotides. In contrast, RPA is almost essential for annealing of longer plasmid-sized DNA but has little effect on the annealing of poly(dT) and poly(dA), which are relatively long DNA molecules free of secondary structure. These results suggest that one role of RPA in Rad52 protein-mediated annealing is the elimination of DNA secondary structure. However, neither Escherichia coli ssDNA binding protein nor human RPA can substitute in this reaction, indicating that RPA has a second role in this process, a role that requires specific RPA–Rad52 protein interactions. This idea is confirmed by the finding that RPA, which is complexed with nonhomologous ssDNA, inhibits annealing but the human RPA–ssDNA complex does not. Finally, we present a model for the early steps of the repair of double-strand DNA breaks in yeast.

DNA-binding and strand-annealing activities of human Mre11: implications for its roles in DNA double-strand break repair pathways

DNA double-strand breaks (DSBs) in eukaryotic cells can be repaired by non-homologous end-joining or homologous recombination. The complex containing the Mre11, Rad50 and Nbs1 proteins has been implicated in both DSB repair pathways, even though they are mechanistically different. To get a better understanding of the properties of the human Mre11 (hMre11) protein, we investigated some of its biochemical activities. We found that hMre11 binds both double- and single-stranded (ss)DNA, with a preference for ssDNA. hMre11 does not require DNA ends for efficient binding. Interestingly, hMre11 mediates the annealing of complementary ssDNA molecules. In contrast to the annealing activity of the homologous recombination protein hRad52, the activity of hMre11 is abrogated by the ssDNA binding protein hRPA. We discuss the possible implications of the results for the role(s) of hMre11 in both DSB repair pathways.

DNA strand annealing is promoted by the yeast Rad52 protein.

The Saccharomyces cerevisiae RAD52 gene plays a pivotal role in genetic recombination. Here we demonstrate that yeast Rad52 is a DNA binding protein. To show that the interaction between Rad52 and DNA is direct and not mediated by other yeast proteins and to facilitate protein purification, a recombinant expression system was developed. The recombinant protein can bind both single- and double-stranded DNA and the addition of either Mg2+ or ATP does not enhance the binding of single-stranded DNA. Furthermore, a DNA binding domain was found in the evolutionary conserved N terminus of the protein. More importantly, we show that the protein stimulates DNA annealing even in the presence of a large excess of nonhomologous DNA. Rad52-promoted annealing follows second-order kinetics and the rate is 3500-fold faster than that of the spontaneous reaction. How this annealing activity relates to the genetic phenotype associated with rad52 mutant cells is discussed.

Chaperonin-facilitated protein folding: optimization of rate and yield by an iterative annealing mechanism.

We develop a heuristic model for chaperonin-facilitated protein folding, the iterative annealing mechanism, based on theoretical descriptions of "rugged" conformational free energy landscapes for protein folding, and on experimental evidence that (i) folding proceeds by a nucleation mechanism whereby correct and incorrect nucleation lead to fast and slow folding kinetics, respectively, and (ii) chaperonins optimize the rate and yield of protein folding by an active ATP-dependent process. The chaperonins GroEL and GroES catalyze the folding of ribulose bisphosphate carboxylase at a rate proportional to the GroEL concentration. Kinetically trapped folding-incompetent conformers of ribulose bisphosphate carboxylase are converted to the native state in a reaction involving multiple rounds of quantized ATP hydrolysis by GroEL. We propose that chaperonins optimize protein folding by an iterative annealing mechanism; they repeatedly bind kinetically trapped conformers, randomly disrupt their structure, and release them in less folded states, allowing substrate proteins multiple opportunities to find pathways leading to the most thermodynamically stable state. By this mechanism, chaperonins greatly expand the range of environmental conditions in which folding to the native state is possible. We suggest that the development of this device for optimizing protein folding was an early and significant evolutionary event.

Toward a mechanism for GroEL.GroES chaperone activity: an ATPase-gated and -pulsed folding and annealing cage.

Free GroEL binds denatured proteins very tightly: it retards the folding of barnase 400-fold and catalyzes unfolding fluctuations in native barnase and its folding intermediate. GroEL undergoes an allosteric transition from its tight-binding T-state to a weaker binding R-state on the cooperative binding of nucleotides (ATP/ADP) and GroES. The preformed GroEL.GroES.nucleotide complex retards the folding of barnase by only a factor of 4, and the folding rate is much higher than the ATPase activity that releases GroES from the complex. Binding of GroES and nucleotides to a preformed GroEL.denatured-barnase complex forms an intermediately fast-folding complex. We propose the following mechanism for the molecular chaperone. Denatured proteins bind to the resting GroEL.GroES.nucleotide complex. Fast-folding proteins are ejected as native structures before ATP hydrolysis. Slow-folding proteins enter chaperoning cycles of annealing and folding after the initial ATP hydrolysis. This step causes transient release of GroES and formation of the GroEL.denatured-protein complexes with higher annealing potential. The intermediately fast-folding complex is formed on subsequent rebinding of GroES. The ATPase activity of GroEL.GroES is thus the gatekeeper that selects for initial entry of slow-folding proteins to the chaperone action and then pumps successive transitions from the faster-folding R-states to the tighter-binding/stronger annealing T-states. The molecular chaperone acts as a combination of folding cage and an annealing machine.