skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Structure-Factor Extrapolation using the Scalar Approximation: Theory, Applications and Limitations

Abstract

For many experiments in macromolecular crystallography, the overall structure of the protein/nucleic acid is already known and the aim of the experiment is to determine the effect a chemical or physical perturbation/activation has on the structure of the molecule. In a typical experiment, an experimenter will collect a data set from a crystal in the unperturbed state, perform the perturbation (i.e. soaking a ligand into the crystal or activating the sample with light) and finally collect a data set from the perturbed crystal. In many cases the perturbation fails to activate all molecules, so that the crystal contains a mix of molecules in the activated and native states. In these cases, it has become common practice to calculate a data set corresponding to a hypothetical fully activated crystal by linear extrapolation of structure-factor amplitudes. These extrapolated data sets often aid greatly in the interpretation of electron-density maps. However, the extrapolation of structure-factor amplitudes is based on a mathematical shortcut that treats structure factors as scalars, not vectors. Here, a full derivation is provided of the error introduced by this approximation and it is determined how this error scales with key experimental parameters. The perhaps surprising result of this analysis ismore » that for most structural changes encountered in protein crystals, the error introduced by the scalar approximation is very small. As a result, the extrapolation procedure is largely limited by the propagation of experimental uncertainties of individual structure-factor amplitudes. Ultimately, propagation of these uncertainties leads to a reduction in the effective resolution of the extrapolated data set. The program XTRA, which implements SASFE (scalar approximation to structure-factor extrapolation), performs error-propagation calculations and determines the effective resolution of the extrapolated data set, is further introduced.« less

Authors:
Publication Date:
Research Org.:
Brookhaven National Laboratory (BNL) National Synchrotron Light Source
Sponsoring Org.:
Doe - Office Of Science
OSTI Identifier:
959549
Report Number(s):
BNL-82535-2009-JA
TRN: US201016%%693
DOE Contract Number:
DE-AC02-98CH10886
Resource Type:
Journal Article
Resource Relation:
Journal Name: Acta Crystallographica Section D: Biological Crystallography; Journal Volume: 63
Country of Publication:
United States
Language:
English
Subject:
59 BASIC BIOLOGICAL SCIENCES; 99 GENERAL AND MISCELLANEOUS//MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; AMPLITUDES; APPROXIMATIONS; CRYSTALLOGRAPHY; ELECTRON DENSITY; EXTRAPOLATION; PROTEINS; RESOLUTION; SCALARS; STRUCTURE FACTORS; VECTORS; national synchrotron light source

Citation Formats

Genick,U. Structure-Factor Extrapolation using the Scalar Approximation: Theory, Applications and Limitations. United States: N. p., 2007. Web. doi:10.1107/S0907444907038164.
Genick,U. Structure-Factor Extrapolation using the Scalar Approximation: Theory, Applications and Limitations. United States. doi:10.1107/S0907444907038164.
Genick,U. Mon . "Structure-Factor Extrapolation using the Scalar Approximation: Theory, Applications and Limitations". United States. doi:10.1107/S0907444907038164.
@article{osti_959549,
title = {Structure-Factor Extrapolation using the Scalar Approximation: Theory, Applications and Limitations},
author = {Genick,U.},
abstractNote = {For many experiments in macromolecular crystallography, the overall structure of the protein/nucleic acid is already known and the aim of the experiment is to determine the effect a chemical or physical perturbation/activation has on the structure of the molecule. In a typical experiment, an experimenter will collect a data set from a crystal in the unperturbed state, perform the perturbation (i.e. soaking a ligand into the crystal or activating the sample with light) and finally collect a data set from the perturbed crystal. In many cases the perturbation fails to activate all molecules, so that the crystal contains a mix of molecules in the activated and native states. In these cases, it has become common practice to calculate a data set corresponding to a hypothetical fully activated crystal by linear extrapolation of structure-factor amplitudes. These extrapolated data sets often aid greatly in the interpretation of electron-density maps. However, the extrapolation of structure-factor amplitudes is based on a mathematical shortcut that treats structure factors as scalars, not vectors. Here, a full derivation is provided of the error introduced by this approximation and it is determined how this error scales with key experimental parameters. The perhaps surprising result of this analysis is that for most structural changes encountered in protein crystals, the error introduced by the scalar approximation is very small. As a result, the extrapolation procedure is largely limited by the propagation of experimental uncertainties of individual structure-factor amplitudes. Ultimately, propagation of these uncertainties leads to a reduction in the effective resolution of the extrapolated data set. The program XTRA, which implements SASFE (scalar approximation to structure-factor extrapolation), performs error-propagation calculations and determines the effective resolution of the extrapolated data set, is further introduced.},
doi = {10.1107/S0907444907038164},
journal = {Acta Crystallographica Section D: Biological Crystallography},
number = ,
volume = 63,
place = {United States},
year = {Mon Jan 01 00:00:00 EST 2007},
month = {Mon Jan 01 00:00:00 EST 2007}
}
  • A procedure for carrying out the Buenker--Peyerimhoff extrapolation using perturbation theory natural orbitals (PTNO's) is presented. This procedure is then applied to the H/sub 2/O and HeH/sub 2//sup +/ systems within the framework of a configuration space comprising all single and double excitations from a single configuration. The results indicate that significantly smaller secular determinants need to be solved when PTNO's are used; the extrapolated CI energies are within 5 x 10/sup -4/ a.u. of energy obtained from the full calculation within the same configuration space.
  • Based on the random-phase approximation and the transcorrelated (TC) method, we optimize the Jastrow factor together with one-electron orbitals in the Slater determinant in the correlated wave function with a new scheme for periodic systems. The TC method is one of the promising wave function theories for first-principles electronic structure calculation, where the many-body wave function is approximated as a product of a Slater determinant and a Jastrow factor, and the Hamiltonian is similarity-transformed by the Jastrow factor. Using this similarity-transformed Hamiltonian, we can optimize the one-electron orbitals without evaluating 3N-dimensional integrations for the N-electron system. In contrast, optimization ofmore » the Jastrow factor within the framework of the TC method is computationally much more expensive and has not been performed for solid-state calculations before. In this study, we also benefit from the similarity-transformation in optimizing the Jastrow factor. Our optimization scheme is tested in applications to some solids from narrow-gap semiconductors to wide-gap insulators, and it is verified that the band gap of a wide-gap insulator and the lattice constants of some solids are improved by this optimization with reasonable computational cost.« less
  • We present a very general method for calculating the electronic excitation spectra and thermodynamic properties of solids with narrow hybridized f-bands. The theory is developed from the localized, atomic limit using the standard-basis operator (SBO) Green's function formalism and greatly improves previous SBO treatments by going beyond the random-phase approximation (RPA). The present theory includes the scattering and resonance-broadening corrections and thus generalizes Hubbard's s-band alloy approximation by incorporating an arbitrary set of atomic configurations.
  • A recently proposed new perturbation technique is applied to a lambda(phi/sup 2/)/sup 1+//sup delta/ self-interacting scalar field theory in D-dimensional space-time. We calculate, to second order in powers of delta, not only the two-point and the four-point Green's functions but also the (2+2delta)-point Green's function, which can be defined in this approximation. After the renormalization procedures are carried out to this order, we find that the theory is not fully nontrivial for Dgreater than or equal to2.