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.

The universal ancestor

A genetic annealing model for the universal ancestor of all extant life is presented; the name of the model derives from its resemblance to physical annealing. The scenario pictured starts when “genetic temperatures” were very high, cellular entities (progenotes) were very simple, and information processing systems were inaccurate. Initially, both mutation rate and lateral gene transfer levels were elevated. The latter was pandemic and pervasive to the extent that it, not vertical inheritance, defined the evolutionary dynamic. As increasingly complex and precise biological structures and processes evolved, both the mutation rate and the scope and level of lateral gene transfer, i.e., evolutionary temperature, dropped, and the evolutionary dynamic gradually became that characteristic of modern cells. The various subsystems of the cell “crystallized,” i.e., became refractory to lateral gene transfer, at different stages of “cooling,” with the translation apparatus probably crystallizing first. Organismal lineages, and so organisms as we know them, did not exist at these early stages. The universal phylogenetic tree, therefore, is not an organismal tree at its base but gradually becomes one as its peripheral branchings emerge. The universal ancestor is not a discrete entity. It is, rather, a diverse community of cells that survives and evolves as a biological unit. This communal ancestor has a physical history but not a genealogical one. Over time, this ancestor refined into a smaller number of increasingly complex cell types with the ancestors of the three primary groupings of organisms arising as a result.

Role of Saccharomyces cerevisiae Msh2 and Msh3 repair proteins in double-strand break-induced recombination

When gene conversion is initiated by a double-strand break (DSB), any nonhomologous DNA that may be present at the ends must be removed before new DNA synthesis can be initiated. In Saccharomyces cerevisiae, removal of nonhomologous ends depends not only on the nucleotide excision repair endonuclease Rad1/Rad10 but also on Msh2 and Msh3, two proteins that are required to correct mismatched bp. These proteins have no effect when DSB ends are homologous to the donor, either in the kinetics of recombination or in the proportion of gene conversions associated with crossing-over. A second DSB repair pathway, single-strand annealing also requires Rad1/Rad10 and Msh2/Msh3, but reveals a difference in their roles. When the flanking homologous regions that anneal are 205 bp, the requirement for Msh2/Msh3 is as great as for Rad1/Rad10; but when the annealing partners are 1,170 bp, Msh2/Msh3 have little effect, while Rad1/Rad10 are still required. Mismatch repair proteins Msh6, Pms1, and Mlh1 are not required. We suggest Msh2 and Msh3 recognize not only heteroduplex loops and mismatched bp, but also branched DNA structures with a free 3′ tail.

Polymer Models of Meiotic and Mitotic Chromosomes

Polymers tied together by constraints exhibit an internal pressure; this idea is used to analyze physical properties of the bottle-brush–like chromosomes of meiotic prophase that consist of polymer-like flexible chromatin loops, attached to a central axis. Using a minimal number of experimental parameters, semiquantitative predictions are made for the bending rigidity, radius, and axial tension of such brushes, and the repulsion acting between brushes whose bristles are forced to overlap. The retraction of lampbrush loops when the nascent transcripts are stripped away, the oval shape of diplotene bivalents between chiasmata, and the rigidity of pachytene chromosomes are all manifestations of chromatin pressure. This two-phase (chromatin plus buffer) picture that suffices for meiotic chromosomes has to be supplemented by a third constituent, a chromatin glue to understand mitotic chromosomes, and explain how condensation can drive the resolution of entanglements. This process resembles a thermal annealing in that a parameter (the affinity of the glue for chromatin and/or the affinity of the chromatin for buffer) has to be tuned to achieve optimal results. Mechanical measurements to characterize this protein–chromatin matrix are proposed. Finally, the propensity for even slightly chemically dissimilar polymers to phase separate (cluster like with like) can explain the apparent segregation of the chromatin into A+T- and G+C-rich regions revealed by chromosome banding.

A green fluorescent protein-reporter mammalian two-hybrid system with extrachromosomal maintenance of a prey expression plasmid: Application to interaction screening

An improved mammalian two-hybrid system designed for interaction trap screening is described in this paper. CV-1/EBNA-1 monkey kidney epithelial cells expressing Epstein–Barr virus nuclear antigen 1 (EBNA-1) were stably transfected with a reporter plasmid for GAL4-dependent expression of the green fluorescent protein (GFP). A resulting clone, GB133, expressed GFP strongly when transfected transiently with transcriptional activators fused to GAL4 DNA-binding domain with minimal background GFP expression. GB133 cells maintained plasmids containing the OriP Epstein–Barr virus replication origin that directs replication of plasmids in mammalian cells in the presence of the EBNA-1 protein. GB133 cells transfected stably with a model bait expressed GFP when further transfected transiently with an expression plasmid for a known positive prey. When the bait-expressing GB133 cells were transfected transiently with an OriP-containing expression plasmid for the positive prey together with excess amounts of empty vector, cells that received the positive prey were readily identified by green fluorescence in cell culture and eventually formed green fluorescent microcolonies, because the prey plasmid was maintained by the EBNA-1/Ori-P system. The green fluorescent microcolonies were harvested directly from the culture dishes under a fluorescence microscope, and total DNA was then prepared. Prey-encoding cDNA was recovered by PCR using primers annealing to the vector sequences flanking the insert-cloning site. This system should be useful in mammalian cells for efficient screening of cDNA libraries by two-hybrid interaction.

A synthetic cryIC gene, encoding a Bacillus thuringiensis δ-endotoxin, confers Spodoptera resistance in alfalfa and tobacco

Spodoptera species, representing widespread polyphagous insect pests, are resistant to Bacillus thuringiensis δ-endotoxins used thus far as insecticides in transgenic plants. Here we describe the chemical synthesis of a cryIC gene by a novel template directed ligation–PCR method. This simple and economical method to construct large synthetic genes can be used when routine resynthesis of genes is required. Chemically phosphorylated adjacent oligonucleotides of the gene to be synthesized are assembled and ligated on a single-stranded, partially homologous template derived from a wild-type gene (cryIC in our case) by a thermostable Pfu DNA ligase using repeated cycles of melting, annealing, and ligation. The resulting synthetic DNA strands are selectively amplified by PCR with short specific flanking primers that are complementary only to the new synthetic DNA. Optimized expression of the synthetic cryIC gene in alfalfa and tobacco results in the production of 0.01–0.2% of total soluble proteins as CryIC toxin and provides protection against the Egyptian cotton leafworm (Spodoptera littoralis) and the beet armyworm (Spodoptera exigua). To facilitate selection and breeding of Spodoptera-resistant plants, the cryIC gene was linked to a pat gene, conferring resistance to the herbicide BASTA.

The nucleocapsid protein specifically anneals tRNALys-3 onto a noncomplementary primer binding site within the HIV-1 RNA genome in vitro

HIV type 1 (HIV-1) specifically uses host cell tRNALys-3 as a primer for reverse transcription. The 3′ 18 nucleotides of this tRNA are complementary to a region on the HIV RNA genome known as the primer binding site (PBS). HIV-1 has a strong preference for maintaining a lysine-specific PBS in vivo, and viral genomes with mutated PBS sequences quickly revert to be complementary to tRNALys-3. To investigate the mechanism for the observed PBS reversion events in vitro, we examined the capability of the nucleocapsid protein (NC) to anneal various tRNA primer sequences onto either complementary or noncomplementary PBSs. We show that NC can anneal different full-length tRNAs onto viral RNA transcripts derived from the HIV-1 MAL or HXB2 isolates, provided that the PBS is complementary to the tRNA used. In contrast, NC promotes specific annealing of only tRNALys-3 onto an RNA template (HXB2) whose PBS sequence has been mutated to be complementary to the 3′ 18 nt of human tRNAPro. Moreover, HIV-1 reverse transcriptase extends this binary complex from the proline-specific PBS. The formation of the noncomplementary binary complex does not occur when a chimeric tRNALys/Pro containing proline-specific D and anticodon domains is used as the primer. Thus, elements outside the acceptor-TΨC domains of tRNALys-3 play an important role in preferential primer use in vitro. Our results support the hypothesis that mutant PBS reversion is a result of tRNALys-3 annealing onto and extension from a PBS that specifies an alternate host cell tRNA.

Solution structure and dynamics of GCN4 cognate DNA: NMR investigations

A 12 bp long GCN4-binding, self-complementary duplex DNA d(CATGACGTCATG)2 has been investigated by NMR spectroscopy to study the structure and dynamics of the molecule in aqueous solution. The NMR structure of the DNA obtained using simulated annealing and iterative relaxation matrix calculations compares quite closely with the X-ray structure of ATF/CREB DNA in complex with GCN4 protein (DNA-binding domain). The DNA is also seen to be curved in the free state and this has a significant bearing on recognition by the protein. The dynamic characteristics of the molecule have been studied by 13C relaxation measurements at natural abundance. A correlation has been observed between sequence-dependent dynamics and recognition by GCN4 protein.

Unwinding of nucleic acids by HCV NS3 helicase is sensitive to the structure of the duplex

Hepatitis C virus (HCV) helicase, non-structural protein 3 (NS3), is proposed to aid in HCV genome replication and is considered a target for inhibition of HCV. In order to investigate the substrate requirements for nucleic acid unwinding by NS3, substrates were prepared by annealing a 30mer oligonucleotide to a 15mer. The resulting 15 bp duplex contained a single-stranded DNA overhang of 15 nt referred to as the bound strand. Other substrates were prepared in which the 15mer DNA was replaced by a strand of peptide nucleic acid (PNA). The PNA–DNA substrate was unwound by NS3, but the observed rate of strand separation was at least 25-fold slower than for the equivalent DNA–DNA substrate. Binding of NS3 to the PNA–DNA substrate was similar to the DNA–DNA substrate, due to the fact that NS3 initially binds to the single-stranded overhang, which was identical in each substrate. A PNA–RNA substrate was not unwound by NS3 under similar conditions. In contrast, morpholino–DNA and phosphorothioate–DNA substrates were utilized as efficiently by NS3 as DNA–DNA substrates. These results indicate that the PNA–DNA and PNA–RNA heteroduplexes adopt structures that are unfavorable for unwinding by NS3, suggesting that the unwinding activity of NS3 is sensitive to the structure of the duplex.