Nobel laureates' letter to President Bush.
Explore the source record for details and available documents.
Biomedical subjects
Publications and source records attributed to H A Hauptman.
Explore the source record for details and available documents.
The Shake-and-Bake method, as implemented in the computer program SnB, has been applied to simulated reference-beam data for the small protein crambin at several resolutions in the range 1.5-3.0 A. Sets of triplet invariants were generated having simulated mean triplet-phase errors from 0 to 60 degrees. Provided that these errors were no larger than 40 degrees, it was possible (at all resolutions tested) to find trial sets of individual Bragg phases with mean errors of 40-45 degrees. At 1.5 A, this could be achieved using only a single reference-beam data set. Peak picking provided useful phase constraints even at the lowest resolution tested. These results suggest that direct methods may be useful in conjunction with reference-beam data at resolutions lower than 1.2 A.
For each fixed reciprocal-lattice vector H, the correlation coefficient, rho H, of the pair (magnitude of EK, magnitude of EH-K), where the magnitudes magnitude of EK and magnitude of EH-K are presumed to be known, is itself positively correlated with magnitude of EH. Thus, the correlation coefficients rho H serve to mediate the transfer of information from one region of reciprocal space to another. In particular, the calculated values of the correlation coefficients rho H lead to estimates of the unknown magnitudes magnitude of EH, even when the latter lie outside the sphere of the experimentally observed magnitude of E s (extrapolation). The applications show that the procedure described here to carry out this extrapolation is superior to existing methods.
The alternation of phase refinement with the imposition of real-space constraints is the essence of the Shake-and-Bake procedure. Typically, these constraints prevent trial structures from falling into local minima. Nevertheless, P1 structures appear to migrate to false minima with significant frequency. These false minima are characterized by the presence of a large 'uranium' peak on the corresponding Fourier map. Fortunately, they can be recognized and avoided by considering the values of the minimal function both before and after the application of constraints. However, it appears that finding solutions for large P1 structures is likely also to require parameter-shift conditions different from those that have been found to work well in other space groups. In fact, these conditions often yield an unusually high percentage of solutions.
The multi-trial direct-methods procedure known as Shake-and-Bake has been applied to three small proteins (alpha-1 peptide, vancomycin and lysozyme) that crystallize in space group P1. Phase refinement was accomplished through parameter-shift optimization using both the cosine and exponential forms of the minimal function. By extending error-free data to sufficiently high resolution, 100% convergence of trial structures to solution could be achieved in all three cases by using the exponential minimal function and a shift angle in the range 130-150 degrees. These results suggest optimum parameters for other P1 structures and emphasize the importance of collecting data to the highest possible resolution.
The simple cosine function used in the formulation of the traditional minimal principle and the related Shake-and-Bake algorithm is here replaced by a function of exponential type and its expected value and variance are derived. These lead to the corresponding exponential minimal principle and its associated Exponential Shake-and-Bake algorithm. Recent applications of the exponential function to several protein structures within the Shake-and-Bake framework suggest that this function leads, in general, to significant improvements in the success rate (percentage of trial structures yielding solution) of the Shake-and-Bake procedure. However, only in space group P1 is it presently possible to assign optimal values a priori for the exponential-function parameters.
The positions of eight Se atoms in a selenomethionyl 35 kDa protein were determined at 2.0 and 2.5 A resolution using the direct-methods program SnB. Data at the selenium peak, edge and remote wavelengths were measured and processed independently. Anomalous difference E magnitudes at each wavelength were derived by two different procedures: renormalized diffE values were calculated according to the equation diffE = ¿sigma[(f + f')2 + f"2]||E+|-|E-||¿1/2/2q(sigma f"2)1/2, where q is a least-squares fitted renormalization function of sinstraight theta/lambda such that = 1.0; and difference E magnitudes were calculated from DeltaF2. Locally normalized E magnitudes corresponding to |FA| were also derived from the combination of the data at all three wavelengths through the use of the MADSYS program suite. Each of the independent sets of anomalous difference E magnitudes was capable of producing the correct solution, as did the E data obtained from the FA data. Higher success rates with SnB were observed for the 2.0 A peak and edge diffE data.
Explore the source record for details and available documents.
The structure of toxin II from the scorpion Androctonus australis Hector has been determined ab initio by direct methods using SnB at 0.96 A resolution. For the purpose of this structure redetermination, undertaken as a test of the minimal function and the SnB program, the identity and sequence of the protein was withheld from part of the research team. A single solution obtained from 1 619 random atom trials was clearly revealed by the bimodal distribution of the final value of the minimal function associated with each individual trial. Five peptide fragments were identified from a conservative analysis of the initial E-map, and following several refinement cycles with X-PLOR, a model was built of the complete structure. At the end of the X-PLOR refinement, the sequence was compared with the published sequence and 57 of the 64 residues had been correctly identified. Two errors in sequence resulted from side chains with similar size while the rest of the errors were a result of severe disorder or high thermal motion in the side chains. Given the amino-acid sequence, it is estimated that the initial E-map could have produced a model containing 99% of all main-chain and 81% of side-chain atoms. The structure refinement was completed with PROFFT, including the contributions of protein H atoms, and converged at a residual of 0.158 for 30 609 data with F >or= 2sigma(F) in the resolution range 8.0-0.964 A. The final model consisted of 518 non-H protein atoms (36 disordered), 407 H atoms, and 129 water molecules (43 with occupancies less than unity). This total of 647 non-H atoms represents the largest light-atom structure solved to date.
Shake-and-Bake is a direct-methods procedure in which phase refinement and Fourier refinement are alternated repetitively, unconditionally and automatically. The traditional Shake-and-Bake approach invoked a parameter-shift routine to perform phase refinement in an effort to reduce the value of minimal function. In this paper, parameter shift is replaced with the tangent formula as a means of phase refinement. This study shows that the tangent formula is more efficient than parameter shift for small structures when the number of refinement cycles and number of applications of the tangent formula per Shake-and-Bake cycle are chosen very carefully. For larger structures, including the 400 non-H-atom crambin structure, the two methods generally perform with similar efficiency. However, only parameter shift has successfully produced recognizable solutions for the difficult 317 non-H-atom structure gramicidin A.
Results from probabilistic theory for the single-wavelength anomalous-scattering (SAS) Friedel pair, two-phase structure invariants, psi H = phi H + phi-H, are used to show that the SAS three-phase structure invariants, psi HK = phi H + phi K + phi-H-K, tend to positive values that are easily estimated. Appropriate averages of the estimates provide SAS perturbation corrections in the form of positive origin shifts for the probability distribution of psi HK values and for the tangent formula. The theoretical probabilistic results are verified by empirical statistical analyses of model-calculated phases and experimentally measured structure-factor magnitudes for a small-molecule and a protein crystal structure.
Explore the source record for details and available documents.
The phase problem, when single-wavelength anomalous scattering (SAS) diffraction data are available, is formulated as a problem in global optimization. Although the objective function has a myriad of local maxima, its global maxima, never more than two, are readily accessible and easily identified by virtue of their isolation. The ability to determine the global maxima of the objective function represents the latest and most successful attempt to go directly from the known probabilistic estimates of the three-phase structure invariants to the values of the individual phases. The relationship between the maxima of the objective function and the solutions of the newly formulated system of SAS tangent equations plays a key role in this development.
A procedure is described by which the phase-invariant translation function may be used to provide starting set phases for the minimal function that are significantly better than those generated from random-atom coordinates. Applications to determine the heavy-atom positions for single isomorphous replacement data from macromolecular structures are very encouraging. In the case of chiral space groups, e.g. P4(1)2(1)2 versus P4(3)2(1)2, unlike Patterson functions, these methods provide the correct enantiomorph coordinates for the heavy-atom sites for whichever space group is chosen.
A new strategy for employing three phase triples invariant estimates from Hauptman's single isomorphous replacement (SIR) and anomalous dispersion (SAS) joint probability distribution formulae is outlined which produces a single unique phase-invariant solution in the case where the positions of the heavy-atom scatterers is known. A similar but non-identical result is obtained for the phase invariants of a structure for which a molecular-replacement solution has been obtained. It is important to note that the values of the individual native/derivative phases can be determined directly from the probability distribution formulae without having to utilize the phase-invariant estimates in an active way. Elimination of the multisolution aspect of utilizing phase-invariant estimates should have important repercussions with regard to phasing macromolecular sets of derivatized data. Trial calculations based on experimentally measured 2.5 A data for three derivatives of cytochrome c550 are encouraging. The average of the three SIR maps resolves a number of structural ambiguities seen in the published multiple isomorphous replacement (MIR) map obtained from eight derivatives.
The shake-and-bake procedure, which is based on the minimal function, has been tested and shown to be extremely effective in molecular-fragment recycling applications. Correctly positioned fragments as small as 5% of the scattering power of the structure typically have a 50% chance of producing a solution in a single recycling trial. While starting models for tangent-formula recycling methods normally require an average r.m.s. displacement error of less than approximately 0.25 A from the refined structure to ensure an adequate chance of success, the shake-and-bake method often tolerates r.m.s. model errors well in excess of 0.5 A. Tests indicate that the new method can outperform traditional tangent-formula procedures in difficult structural applications involving multiple copies of pseudosymmetrically related molecules or low-resolution data.
The crystal structure of crambin, a 46-residue protein containing the equivalent of approximately 400 fully occupied non-H-atom positions, was originally solved at 1.5 A by exploiting the anomalous scattering of its six S atoms at a single wavelength far removed from the absorption edge of sulfur. The crambin structure has now been resolved without the use of any anomalous-dispersion measurements. The technique employed was an ab initio 'shake-and-bake' method, consisting of a phase-refinement procedure based on the minimal function alternated with Fourier refinement. This method has successfully yielded solutions for a smaller molecule (28 atoms) using 1.2 A data, and a crambin solution was obtained at 1.1 A.
For diffraction data at two wavelengths from a crystal with anomalous scatterers, there are six types of two-phase structure invariants for Friedel pairs. Two of the six are single-wavelength invariants; the other four are mixed-wavelength invariants. It is shown that the latter can be estimated by a straightforward extension of results from the probabilistic direct-methods theory for the single-wavelength anomalous scattering case described in paper I [Hauptman (1982). Acta Cryst. A38, 632-641]. Statistical tests of the mixed-wavelength estimates are reported for small-molecule and macro-molecular examples.