Quantitative symmetry and chirality-A fast computational algorithm for large structures: Proteins, macromolecules, nanotubes, and unit cells

Graph-theoretical chirality measure and chirality-property relations for chemical structures with multiscale mirror asymmetries

Chirality is an essential geometric property unifying small molecules, biological macromolecules, inorganic nanomaterials, biological microparticles, and many other chemical structures. Numerous chirality measures have attempted to quantify this geometric property of mirror asymmetry and to correlate these measures with physical and chemical properties. However, their utility has been widely limited because these correlations have been largely notional. Furthermore, chirality measures also require prohibitively demanding computations, especially for chiral structures comprised of thousands of atoms. Acknowledging the fundamental problems with quantification of mirror asymmetry, including the ambiguity of sign-variable pseudoscalar chirality measures, we revisit this subject because of the significance of quantifying chirality for quantitative biomimetics and describing the chirality of nanoscale materials that display chirality continuum and scale-dependent mirror asymmetry. We apply the concept of torsion within the framework of differential geometry to the graph theoretical representation of chiral molecules and nanostructures to address some of the fundamental problems and practical limitations of other chirality measures. Chiral gold clusters and other chiral structures are used as models to elaborate a graph-theoretical chirality (GTC) measure, demonstrating its applicability to chiral materials with different degrees of chirality at different scales. For specific cases, we show that GTC provides an adequate description of both the sign and magnitude of mirror asymmetry. The direct correlations with macroscopic properties, such as chiroptical spectra, are enhanced by using the hybrid chirality measures combining parameters from discrete mathematics and physics. Taking molecular helices as an example, we established a direct relation between GTC and optical activity, indicating that this chirality measure can be applied to chiral metamaterials and complex chiral constructs.

Continuous symmetry and chirality measures: approximate algorithms for large molecular structures

Quantifying imperfect symmetry of molecules can help explore the sources, roles and extent of structural distortion. Based on the established methodology of continuous symmetry and chirality measures, we develop a set of three-dimensional molecular descriptors to estimate distortion of large structures. These three-dimensional geometrical descriptors quantify the gap between the desirable symmetry (or chirality) and the actual one. They are global parameters of the molecular geometry, intuitively defined, and have the ability to detect even minute structural changes of a given molecule across chemistry, including organic, inorganic, and biochemical systems. Application of these methods to large structures is challenging due to countless permutations that are involved in the symmetry operations and have to be accounted for. Our approach focuses on iteratively finding the approximate direction of the symmetry element in the three-dimensional space, and the relevant permutation. Major algorithmic improvements over previous versions are described, showing increased accuracy, reliability and structure preservation. The new algorithms are tested for three sets of molecular structures including pillar[5]arene complexes with Li+, C100 fullerenes, and large unit cells of metal organic frameworks. These developments complement our recent algorithms for calculating continuous symmetry and chirality measures for small molecules as well as protein homomers, and simplify the usage of the full set of measures for various research goals, in molecular modeling, QSAR and cheminformatics.

The SARS-CoV-2 spike protein structure: a symmetry tale on distortion trail

A preliminary step in the SARS-CoV-2 human infection process is a conformational change of the receptor binding domain (RBD) of its spike protein, characterized by a significant loss of symmetry. During this process, the residues which later on bind to the human angiotensin converting enzyme 2 (ACE2) receptor, become exposed at the surface of the protein. Symmetry analysis of a data set of 33 protein structures from experimental measurements and 32 structures from molecular dynamics simulation, show that the initial state carries clear indications on the structure of the final state, with respect to the local distortion along the sequence. This surprising finding implies that this type of analysis predicts the mechanism of change. We further show that the level of local distortion at the initial state increases with variant's transmissibility, for the wild type (WT) along with past and present variants of concern (WT ∼ alpha < beta < delta < Omicron BA.1), in accordance with the trend of their evolutionary path. In other words, the initial structure of the variant which is most infectious is also the most distorted, making its path to the final state shorter. It has been claimed that the RBD migration of the spike protein is allosterically controlled. Our analysis provides a quantitative support to a major theorem in this respect - that information about an allosteric process is encoded in the structure itself, suggesting that the path of local distortion is related to an allosteric information network.

Chirality in the Solid State: Chiral Crystal Structures in Chiral and Achiral Space Groups

Chirality depends on particular symmetries. For crystal structures it describes the absence of mirror planes and inversion centers, and in addition to translations, only rotations are allowed as symmetry elements. However, chiral space groups have additional restrictions on the allowed screw rotations as a symmetry element, because they always appear in enantiomorphous pairs. This study classifies and distinguishes the chiral structures and space groups. Chirality is quantified using Hausdorff distances and continuous chirality measures and selected crystal structures are reported. Chirality is discussed for bulk solids and their surfaces. Moreover, the band structure, and thus, the density of states, is found to be affected by the same crystal parameters as chirality. However, it is independent of handedness. The Berry curvature, as a topological measure of the electronic structure, depends on the handedness but is not proof of chirality because it responds to the inversion of a structure. For molecules, optical circular dichroism is one of the most important measures for chirality. Thus, it is proposed in this study that the circular dichroism in the angular distribution of photoelectrons in high symmetry configurations can be used to distinguish the handedness of chiral solids and their surfaces.

Bottlebrush polymers with flexible enantiomeric side chains display differential biological properties

Chirality and molecular conformation are central components of life: biological systems rely on stereospecific interactions between discrete (macro)molecular conformers, and the impacts of stereochemistry and rigidity on the properties of small molecules and biomacromolecules have been intensively studied. Nevertheless, how these features affect the properties of synthetic macromolecules has received comparably little attention. Here we leverage iterative exponential growth and ring-opening metathesis polymerization to produce water-soluble, chiral bottlebrush polymers (CBPs) from two enantiomeric pairs of macromonomers of differing rigidity. Remarkably, CBPs with conformationally flexible, mirror image side chains show several-fold differences in cytotoxicity, cell uptake, blood pharmacokinetics and liver clearance; CBPs with comparably rigid, mirror image side chains show no differences. These observations are rationalized with a simple model that correlates greater conformational freedom with enhanced chiral recognition. Altogether, this work provides routes to the synthesis of chiral nanostructured polymers and suggests key roles for stereochemistry and conformational rigidity in the design of future biomaterials.

Side chain flexibility and the symmetry of protein homodimers

A comprehensive analysis of crystallographic data of 565 high-resolution protein homodimers comprised of over 250,000 residues suggests that amino acids form two groups that differ in their tendency to distort or symmetrize the structure of protein homodimers. Residues of the first group tend to distort the protein homodimer and generally have long or polar side chains. These include: Lys, Gln, Glu, Arg, Asn, Met, Ser, Thr and Asp. Residues of the second group contribute to protein symmetry and are generally characterized by short or aromatic side chains. These include: Ile, Pro, His, Val, Cys, Leu, Trp, Tyr, Phe, Ala and Gly. The distributions of the continuous symmetry measures of the proteins and the continuous chirality measures of their building blocks highlight the role of side chain geometry and the interplay between entropy and symmetry in dictating the conformational flexibility of proteins.

Chiral Dualism as an Instrument of Hierarchical Structure Formation in Molecular Biology

The origin of chiral asymmetry in biology has attracted the attention of the research community throughout the years. In this paper we discuss the role of chirality and chirality sign alternation (L–D–L–D in proteins and D–L–D–L in DNA) in promoting self-organization in biology, starting at the level of single molecules and continuing to the level of supramolecular assemblies. In addition, we also discuss chiral assemblies in solutions of homochiral organic molecules. Sign-alternating chiral hierarchies created by proteins and nucleic acids are suggested to create the structural basis for the existence of selected mechanical degrees of freedom required for conformational dynamics in enzymes and macromolecular machines.

Improved algorithms for quantifying the near symmetry of proteins: complete side chains analysis

Symmetry of proteins, an important source of their elegant structure and unique functions, is not as perfect as it may seem. In the framework of continuous symmetry, in which symmetry is no longer a binary yes/no property, such imperfections can be quantified and used as a global descriptor of the three-dimensional structure. We present an improved algorithm for calculating the continuous symmetry measure for proteins that takes into account their complete set of atoms including all side chains. Our method takes advantage of the protein sequence and the division into peptides in order to improve the accuracy and efficiency of the calculation over previous methods. The Hungarian algorithm is applied to solve the assignment problem and find the permutation that defines the symmetry operation. Analysis of the symmetry of several sets of protein homomers, with various degrees of rotational symmetry is presented. The new methodology lays the foundations for accurate, efficient and reliable large scale symmetry analysis of protein structure and can be used as a collective variable that describes changes of the protein geometry along various processes, both at the backbone level and for the complete protein structure.

Analytical symmetry detection in protein assemblies. II. Dihedral and cubic symmetries

Protein assemblies are often symmetric, as this organization has many advantages compared to individual proteins. Complex protein structures thus very often possess high-order symmetries. Detection and analysis of these symmetries has been a challenging problem and no efficient algorithms have been developed so far. This paper presents the extension of our cyclic symmetry detection method for higher-order symmetries with multiple symmetry axes. These include dihedral and cubic, i.e., tetrahedral, octahedral, and icosahedral, groups. Our method assesses the quality of a particular symmetry group and also determines all of its symmetry axes with a machine precision. The method comprises discrete and continuous optimization steps and is applicable to assemblies with multiple chains in the asymmetric subunits or to those with pseudo-symmetry. We implemented the method in C++ and exhaustively tested it on all 51,358 symmetric assemblies from the Protein Data Bank (PDB). It allowed us to study structural organization of symmetric assemblies solved by X-ray crystallography, and also to assess the symmetry annotation in the PDB. For example, in 1.6% of the cases we detected a higher symmetry group compared to the PDB annotation, and we also detected several cases with incorrect annotation. The method is available at http://team.inria.fr/nano-d/software/ananas. The graphical user interface of the method built for the SAMSON platform is available at http://samson-connect.net.

Analytical symmetry detection in protein assemblies. I. Cyclic symmetries

Symmetry in protein, and, more generally, in macromolecular assemblies is a key point to understand their structure, stability and function. Many symmetrical assemblies are currently present in the Protein Data Bank (PDB) and some of them are among the largest solved structures, thus an efficient computational method is needed for the exhaustive analysis of these. The cyclic symmetry groups represent the most common assemblies in the PDB. These are also the building blocks for higher-order symmetries. This paper presents a mathematical formulation to find the position and the orientation of the symmetry axis in a cyclic symmetrical protein assembly, and also to assess the quality of this symmetry. Our method can also detect symmetries in partial assemblies. We provide an efficient C++ implementation of the method and demonstrate its efficiency on several examples including partial assemblies and pseudo symmetries. We also compare the method with two other published techniques and show that it is significantly faster on all the tested examples. Our method produces results with a machine precision, its cost function is solely based on 3D Euclidean geometry, and most of the operations are performed analytically. The method is available athttp://team.inria.fr/nano-d/software/ananas. The graphical user interface of the method built for the SAMSON platform is available athttp://samson-connect.net.

The near-symmetry of proteins

The majority of protein oligomers form clusters which are nearly symmetric. Understanding of that imperfection, its origins, and perhaps also its advantages requires the conversion of the currently used vague qualitative descriptive language of the near-symmetry into an accurate quantitative measure that will allow to answer questions such as: "What is the degree of symmetry deviation of the protein?," "how do these deviations compare within a family of proteins?," and so on. We developed quantitative methods to answer this type of questions, which are capable of analyzing the whole protein, its backbone or selected portions of it, down to comparison of symmetry-related specific amino-acids, and which are capable of visualizing the various levels of symmetry deviations in the form of symmetry maps. We have applied these methods on an extensive list of homomers and heteromers and found that apparently all proteins never reach perfect symmetry. Strikingly, even homomeric protein clusters are never ideally symmetric. We also found that the main burden of symmetry distortion is on the amino-acids near the symmetry axis; that it is mainly the more hydrophilic amino-acids that take place in symmetry-distortive interactions; and more. The remarkable ability of heteromers to preserve near-symmetry, despite the different sequences, was also shown and analyzed. The comprehensive literature on the suggested advantages symmetric oligomerizations raises a yet-unsolved key question: If symmetry is so advantageous, why do proteins stop shy of perfect symmetry? Some tentative answers to be tested in further studies are suggested in a concluding outlook.

Continuous symmetry measures for complex symmetry group

Symmetry is a fundamental property of nature, used extensively in physics, chemistry, and biology. The Continuous symmetry measures (CSM) is a method for estimating the deviation of a given system from having a certain perfect symmetry, which enables us to formulate quantitative relation between symmetry and other physical properties. Analytical procedures for calculating the CSM of all simple cyclic point groups are available for several years. Here, we present a methodology for calculating the CSM of any complex point group, including the dihedral, tetrahedral, octahedral, and icosahedral symmetry groups. We present the method and analyze its performances and errors. We also introduce an analytical method for calculating the CSM of the linear symmetry groups. As an example, we apply these methods for examining the symmetry of water, the symmetry maps of AB4 complexes, and the symmetry of several Lennard-Jones clusters.