Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2020 Jan 17;7(Pt 2):158-165.
doi: 10.1107/S2052252519015975. eCollection 2020 Mar 1.

fragHAR: towards ab initio quantum-crystallographic X-ray structure refinement for polypeptides and proteins

Affiliations

fragHAR: towards ab initio quantum-crystallographic X-ray structure refinement for polypeptides and proteins

Justin Bergmann et al. IUCrJ. .

Abstract

The first ab initio aspherical structure refinement against experimental X-ray structure factors for polypeptides and proteins using a fragmentation approach to break up the protein into residues and solvent, thereby speeding up quantum-crystallographic Hirshfeld atom refinement (HAR) calculations, is described. It it found that the geometric and atomic displacement parameters from the new fragHAR method are essentially unchanged from a HAR on the complete unfragmented system when tested on dipeptides, tripeptides and hexapeptides. The largest changes are for the parameters describing H atoms involved in hydrogen-bond interactions, but it is shown that these discrepancies can be removed by including the interacting fragments as a single larger fragment in the fragmentation scheme. Significant speed-ups are observed for the larger systems. Using this approach, it is possible to perform a highly parallelized HAR in reasonable times for large systems. The method has been implemented in the TONTO software.

Keywords: H atoms; Hirshfeld atom refinement; peptides; quantum crystallography.

PubMed Disclaimer

Figures

Figure 1
Figure 1
The molecular fractionation with conjugate caps (MFCC) procedure for cutting a dipeptide (left) across the peptide bond (shown in black), producing two fragment molecules (right), which are then ‘capped’ with –CH3C=O (red) and –NHCH3 (orange) groups, comprised of the neighbour and next-neighbour non-H atoms.
Figure 2
Figure 2
Crystal structures (100 K) of the three peptide model systems with 50% ADP probability ellipsoids. Hydrogen bonds are shown in green. (a) Gly-Ala (GA), (b) Ala-His-Ala (AHA), (c) cyclo-(Ala)4-(d,l-Pro)2 (A4P2).
Figure 3
Figure 3
Bond lengths between non-H atoms from fragHAR calculations plotted against reference HAR values. Error bars are depicted, but are invisible to the eye on this scale.
Figure 4
Figure 4
X—H bonds (with error bars) in all model compounds for fragHAR calculations versus reference HAR values. Bonds with notable differences are marked with the corresponding H-atom name.
Figure 5
Figure 5
X—H bonds (with error bars) in all model compounds for fragHAR with fragments ‘joined’ across hydrogen bonds versus reference HAR calculations.
Figure 6
Figure 6
Hydrogen ADPs (with error bars) obtained with fragHAR versus those from HAR for the A4P2 system. The hydrogen-bonded H41A atom is labelled.
Figure 7
Figure 7
Timing of the fragHAR (green) and HAR calculations (black) for single-processor serial (square) and parallel (circles) calculations for GA (two processors), AHA (four processors) and A4P2 (four processors).

References

    1. Allen, F. H. & Bruno, I. J. (2010). Acta Cryst. B66, 380–386. - PubMed
    1. Antony, J. & Grimme, S. (2012). J. Comput. Chem. 33, 1730–1739. - PubMed
    1. Becke, A. D. (1988). J. Chem. Phys. 88, 2547–2553.
    1. Borbulevych, O. Y., Plumley, J. A., Martin, R. I., Merz, K. M. & Westerhoff, L. M. (2014). Acta Cryst. D70, 1233–1247. - PMC - PubMed
    1. Brünger, A. T., Kuriyan, J. & Karplus, M. (1987). Science, 235, 458–460. - PubMed