Predicting crossover generation in DNA shuffling

We introduce a quantitative framework for assessing the generation of crossovers in DNA shuffling experiments. The approach uses free energy calculations and complete sequence information to model the annealing process. Statistics obtained for the annealing events then are combined with a reassembly algorithm to infer crossover allocation in the reassembled sequences. The fraction of reassembled sequences containing zero, one, two, or more crossovers and the probability that a given nucleotide position in a reassembled sequence is the site of a crossover event are estimated. Comparisons of the predictions against experimental data for five example systems demonstrate good agreement despite the fact that no adjustable parameters are used. An in silico case study of a set of 12 subtilases examines the effect of fragmentation length, annealing temperature, sequence identity and number of shuffled sequences on the number, type, and distribution of crossovers. A computational verification of crossover aggregation in regions of near-perfect sequence identity and the presence of synergistic reassembly in family DNA shuffling is obtained.

Site-directed spin labeling and chemical crosslinking demonstrate that helix V is close to helices VII and VIII in the lactose permease of Escherichia coli.

Site-directed chemical cleavage of lactose permease indicates that helix V is in close proximity to helices VII and VIII. To test this conclusion further, permease containing a biotin-acceptor domain and paired Cys residues at positions 148 (helix V) and 228 (helix VII), 148 and 226 (helix VII), or 148 and 275 (helix VIII) was affinity purified and labeled with a sulfhydryl-specific nitroxide spin label. Spin-spin interactions are observed with the 148/228 and 148/275 pairs, indicating close proximity between appropriate faces of helix V and helices VII and VIII. Little or no interaction is evident with the 148/226 pair, in all likelihood because position 226 is on the opposite face of helix VII from position 228. Broadening of the electron paramagnetic resonance spectra in the frozen state was used to estimate distance between the 148/228 and the 148/275 pairs. The nitroxides at positions 148 and 228 or 148 and 275 are within approximately 13-15 A. Finally, Cys residues at positions 148 and 228 are crosslinked by dibromobimane, a bifunctional crosslinker that is approximately 5 A. long, while no crosslinking is detected between Cys residues at positions 148 and 275 or 148 and 226. The results provide strong support for a structure in which helix V is in close proximity to both helices VII and VIII and is oriented in such a fashion that Cys-148 is closer to helix VII.

A method for distance determination in proteins using a designed metal ion binding site and site-directed spin labeling: evaluation with T4 lysozyme.

The use of molecular genetics to introduce both a metal ion binding site and a nitroxide spin label into the same protein opens the use of paramagnetic metalnitroxyl interactions to estimate intramolecular distances in a wide variety of proteins. In this report, a His-Xaa3-His metal ion binding motif was introduced at the N terminus of the long interdomain helix of T4 lysozyme (Lys-65 --> His/Gln-69 --> His) of three mutants, each containing a single nitroxide-labeled cysteine residue at position 71, 76, or 80. The results show that Cu(II)-induced relaxation effects on the nitroxide can be quantitatively analyzed in terms of interspin distance in the range of 10-25 A using Redfield theory, as first suggested by Leigh [Leigh, J.S. (1970) J. Chem. Phys. 52, 2608-2612]. Of particular interest is the observation that distances can be determined both under rigid lattice conditions in frozen solution and in the presence of motion of the spins at room temperature under physiological conditions. The method should be particularly attractive for investigating structure in membrane proteins that are difficult to crystallize. In the accompanying paper, the technique is applied to a polytopic membrane protein, lactose permease.

Distance determination in proteins using designed metal ion binding sites and site-directed spin labeling: application to the lactose permease of Escherichia coli.

As shown in the accompanying paper, the magnetic dipolar interaction between site-directed metal-nitroxide pairs can be exploited to measure distances in T4 lysozyme, a protein of known structure. To evaluate this potentially powerful method for general use, particularly with membrane proteins that are difficult to crystallize, both a paramagnetic metal ion binding site and a nitroxide side chain were introduced at selected positions in the lactose permease of Escherichia coli, a paradigm for polytopic membrane proteins. Thus, three individual cysteine residues were introduced into putative helix IV of a lactose permease mutant devoid of native cysteine residues containing a high-affinity divalent metal ion binding site in the form of six contiguous histidine residues in the periplasmic loop between helices III and IV. In addition, the construct contained a biotin acceptor domain in the middle cytoplasmic loop to facilitate purification. After purification and spin labeling, electron paramagnetic resonance spectra were obtained with the purified proteins in the absence and presence of Cu(II). The results demonstrate that positions 103, 111, and 121 are 8, 14, and > 23 A from the metal binding site. These data are consistent with an alpha-helical conformation of transmembrane domain IV of the permease. Application of the technique to determine helix packing in lactose permease is discussed.

Determination of the distance between two spin labels attached to a macromolecule.

An EPR "spectroscopic ruler" was developed using a series of alpha-helical polypeptides, each modified with two nitroxide spin labels. The EPR line broadening due to electron-electron dipolar interactions in the frozen state was determined using the Fourier deconvolution method. These dipolar spectra were then used to estimate the distances between the two nitroxides separated by 8-25 A. Results agreed well with a simple alpha-helical model. The standard deviation from the model system was 0.9 A in the range of 8-25 A. This technique is applicable to complex systems such as membrane receptors and channels, which are difficult to access with high-resolution NMR or x-ray crystallography, and is expected to be particularly useful for systems for which optical methods are hampered by the presence of light-interfering membranes or chromophores.

Crossover and noncrossover recombination during meiosis: timing and pathway relationships.

During meiosis, crossovers occur at a high level, but the level of noncrossover recombinants is even higher. The biological rationale for the existence of the latter events is not known. It has been suggested that a noncrossover-specific pathway exists specifically to mediate chromosome pairing. Using a physical assay that monitors both crossovers and noncrossovers in cultures of yeast undergoing synchronous meiosis, we find that both types of products appear at essentially the same time, after chromosomes are fully synapsed at pachytene. We have also analyzed a situation in which commitment to meiotic recombination and formation of the synaptonemal complex are coordinately suppressed (mer1 versus mer1 MER2++). We find that suppression is due primarily to restoration of meiosis-specific double-strand breaks, a characteristic of the major meiotic recombination pathway. Taken together, the observations presented suggest that there probably is no noncrossover-specific pathway and that restoration of intermediate events in a single pairing/recombination pathway promotes synaptonemal complex formation. The biological significant of noncrossover recombination remains to be determined, however.

Temperature chaos in 3D Ising spin glasses is driven by rare events

Temperature chaos has often been reported in the literature as a rare-event–driven phenomenon. However, this fact has always been ignored in the data analysis, thus erasing the signal of the chaotic behavior (still rare in the sizes achieved) and leading to an overall picture of a weak and gradual phenomenon. On the contrary, our analysis relies on a largedeviations functional that allows to discuss the size dependences. In addition, we had at our disposal unprecedentedly large configurations equilibrated at low temperatures, thanks to the Janus computer. According to our results, when temperature chaos occurs its effects are strong and can be felt even at short distances.

Correspondence between long-range and short-range spin glasses

We compare the critical behavior of the short-range Ising spin glass with a spin glass with long-range interactions which fall off as a power σ of the distance. We show that there is a value of σ of the long-range model for which the critical behavior is very similar to that of the short range model in four dimensions. We also study a value of σ for which we find the critical behavior to be compatible with that of the three-dimensional model, although we have much less precision than in the four-dimensional case.

Thermodynamic glass transition in a spin glass without time-reversal symmetry

Spin glasses are a longstanding model for the sluggish dynamics that appear at the glass transition. However, spin glasses differ from structural glasses in a crucial feature: they enjoy a time reversal symmetry. This symmetry can be broken by applying an external magnetic field, but embarrassingly little is known about the critical behavior of a spin glass in a field. In this context, the space dimension is crucial. Simulations are easier to interpret in a large number of dimensions, but one must work below the upper critical dimension (i.e., in d < 6) in order for results to have relevance for experiments. Here we show conclusive evidence for the presence of a phase transition in a four-dimensional spin glass in a field. Two ingredients were crucial for this achievement: massive numerical simulations were carried out on the Janus special-purpose computer, and a new and powerful finite-size scaling method.

Sample-to-sample fluctuations of the overlap distributions in the three-dimensional Edwards-Anderson spin glass

We study the sample-to-sample fluctuations of the overlap probability densities from large-scale equilibrium simulations of the three-dimensional Edwards-Anderson spin glass below the critical temperature. Ultrametricity, stochastic stability, and overlap equivalence impose constraints on the moments of the overlap probability densities that can be tested against numerical data. We found small deviations from the Ghirlanda Guerra predictions, which get smaller as system size increases. We also focus on the shape of the overlap distribution, comparing the numerical data to a mean-field-like prediction in which finite-size effects are taken into account by substituting delta functions with broad peaks.

Finite-size scaling analysis of the distributions of pseudo-critical temperatures in spin glasses

Using the results of large scale numerical simulations we study the probability distribution of the pseudo critical temperature for the three dimensional Edwards Anderson Ising spin glass and for the fully connected Sherrington-Kirkpatrick model. We find that the behaviour of our data is nicely described by straightforward finitesize scaling relations.

Static versus dynamic heterogeneities in the D=3 Edwards-Anderson-ising spin glass

We numerically study the aging properties of the dynamical heterogeneities in the Ising spin glass. We find that a phase transition takes place during the aging process. Statics-dynamics correspondence implies that systems of finite size in equilibrium have static heterogeneities that obey finite-size scaling, thus signaling an analogous phase transition in the thermodynamical limit. We compute the critical exponents and the transition point in the equilibrium setting, and use them to show that aging in dynamic heterogeneities can be described by a finite-time scaling ansatz, with potential implications for experimental work.