skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Mathematical aspects of molecular replacement. V. Isolating feasible regions in motion spaces
This paper mathematically characterizes the tiny feasible regions within the vast 6D rotation–translation space in a full molecular replacement (MR) search. The capability to a priori isolate such regions is potentially important for enhancing robustness and efficiency in computational phasing in macromolecular crystallography (MX). The previous four papers in this series have concentrated on the properties of the full configuration space of rigid bodies that move relative to each other with crystallographic symmetry constraints. In particular, it was shown that the configuration space of interest in this problem is the right-coset space Γ\ G , where Γ is the space group of the chiral macromolecular crystal and G is the group of rigid-body motions, and that fundamental domains F Γ\ G can be realized in many ways that have interesting algebraic and geometric properties. The cost function in MR methods can be viewed as a function on these fundamental domains. This, the fifth and final paper in this series, articulates the constraints that bodies packed with crystallographic symmetry must obey. It is shown that these constraints define a thin feasible set inside a motion space and that they fall into two categories: (i) the bodies must not interpenetrate, thereby excluding so-called `collision zones' from consideration in MR searches; (ii) the bodies must be in contact with a sufficient number of neighbors so as to form a rigid network leading to a physically realizable crystal. In this paper, these constraints are applied using ellipsoidal proxies for proteins to bound the feasible regions. It is shown that the volume of these feasible regions is small relative to the total volume of the motion space, which justifies the use of ellipsoids as proxies for complex proteins in MR searches, and this is demonstrated with P 1 (the simplest space group) and with P 2 1 2 1 2 1 (the most common space group in MX).  more » « less
Award ID(s):
1640970
PAR ID:
10201757
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Acta Crystallographica Section A Foundations and Advances
Volume:
76
Issue:
2
ISSN:
2053-2733
Page Range / eLocation ID:
145 to 162
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In molecular-replacement (MR) searches, spaces of motions are explored for determining the appropriate placement of rigid-body models of macromolecules in crystallographic asymmetric units. The properties of the space of non-redundant motions in an MR search, called a `motion space', are the subject of this series of papers. This paper, the fourth in the series, builds on the others by showing that when the space group of a macromolecular crystal can be decomposed into a product of two space subgroups that share only the lattice translation group, the decomposition of the group provides different decompositions of the corresponding motion spaces. Then an MR search can be implemented by trading off between regions of the translation and rotation subspaces. The results of this paper constrain the allowable shapes and sizes of these subspaces. Special choices result when the space group is decomposed into a product of a normal Bieberbach subgroup and a symmorphic subgroup (which is a common occurrence in the space groups encountered in protein crystallography). Examples of Sohncke space groups are used to illustrate the general theory in the three-dimensional case (which is the relevant case for MR), but the general theory in this paper applies to any dimension. 
    more » « less
  2. This paper presents a method of computing free motions of a planar assembly of rigid bodies connected by loose joints. Joints are modeled using local distance constraints, which are then linearized with respect to configuration space velocities, yielding a linear programming formulation that allows analysis of systems with thousands of rigid bodies. Potential applications include analysis of collections of modular robots, structural stability perturbation analysis, tolerance analysis for mechanical systems, and formation control of mobile robots. 
    more » « less
  3. This paper presents a method of computing free motions of a planar assembly of rigid bodies connected by loose joints. Joints are modeled using local distance constraints, which are then linearized with respect to configuration space velocities, yielding a linear programming formulation that allows analysis of systems with thousands of rigid bodies. Potential applications include analysis of collections of modular robots, structural stability perturbation analysis, tolerance analysis for mechanical systems, and formation control of mobile robots. 
    more » « less
  4. Structural modelling of octahedral tilts in perovskites is typically carried out using the symmetry constraints of the resulting space group. In most cases, this introduces more degrees of freedom than those strictly necessary to describe only the octahedral tilts. It can therefore be a challenge to disentangle the octahedral tilts from other structural distortions such as cation displacements and octahedral distortions. This paper reports the development of constraints for modelling pure octahedral tilts and implementation of the constraints in diffpy-CMI , a powerful package to analyse pair distribution function (PDF) data. The model in the program allows features in the PDF that come from rigid tilts to be separated from non-rigid relaxations, providing an intuitive picture of the tilting. The model has many fewer refinable variables than the unconstrained space group fits and provides robust and stable refinements of the tilt components. It further demonstrates the use of the model on the canonical tilted perovskite CaTiO 3 which has the known Glazer tilt system α + β − β − . The Glazer model fits comparably to the corresponding space-group model Pnma below r = 14 Å and becomes progressively worse than the space-group model at higher r due to non-rigid distortions in the real material. 
    more » « less
  5. The topological symmetry group of an embedding Γ of an abstract graph γ in S3 is the group of automorphisms of γ that can be realized by homeomorphisms of the pair (S3,Γ). These groups are motivated by questions about the symmetries of molecules in space. The Petersen family of graphs is an important family of graphs for many problems in low-dimensional topology, so it is desirable to understand the possible groups of symmetries of their embeddings in space. In this paper, we find all the groups that can be realized as topological symmetry groups for each of the graphs in the Petersen family. Along the way, we also complete the classification of the realizable topological symmetry groups for K3,3. 
    more » « less