During in silico crystal structure prediction of organic molecules, millions of candidate structures are often generated. These candidates must be compared to remove duplicates prior to further analysis (e.g., optimization with electronic structure methods) and ultimately compared with structures determined experimentally. The agreement of predicted and experimental structures forms the basis of evaluating the results from the Cambridge Crystallographic Data Centre (CCDC) blind assessment of crystal structure prediction, which further motivates the importance of rigorous alignments. Evaluating crystal structure packings in a reproducible manner requires not only calculating a coordinate root-mean-square deviation (RSMD) for N molecules (or N asymmetric units), but we argue should also include metrics to describe the shape of the compared molecular clusters to account for alternative approaches used to prioritize selection of molecules. Here we describe a flexible algorithm called Progressive Alignment of Crystals (PAC) to evaluate crystal packing similarity using coordinate RMSD and introduce radius of gyration (Rg) as a metric to quantify the shape of the superimposed clusters. We show that the absence of metrics to describe cluster shape adds ambiguity to the results of the CCDC blind assessments because it is not possible to determine whether the superposition algorithm prioritized tightly packed molecular clusters (i.e., to minimize Rg) or prioritized reduced RMSD (i.e., via possibly elongated clusters with relatively larger Rg). For example, we show that when the PAC algorithm described here uses “single linkage” to prioritize molecules for inclusion in the superimposed clusters, our results are nearly identical to those calculated by the widely used program COMPACK. However, we favor the lower Rg values obtained by use of “average linkage” for molecule prioritization because the resulting RMSDs more equally reflect the importance of packing along each dimension. We conclude by showing that the PAC algorithm is faster than COMPACK when using a single process, demonstrate its utility for biomolecular crystals, and finally present parallel scaling up to 64 processes in the open-source code Force Field X.

Autism spectrum disorders (ASD) involve brain wide abnormalities that contribute to a constellation of symptoms including behavioral inflexibility, cognitive dysfunction, learning impairments, altered social interactions, and perceptive time difficulties. Although a single genetic variation does not cause ASD, genetic variations such as one involving a non-canonical Wnt signaling gene, Prickle2, has been found in individuals with ASD. Previous work looking into phenotypes of Prickle2 knock-out (Prickle2−/−) and heterozygous mice (Prickle2−/+) suggest patterns of behavior similar to individuals with ASD including altered social interaction and behavioral inflexibility. Growing evidence implicates the cerebellum in ASD. As Prickle2 is expressed in the cerebellum, this animal model presents a unique opportunity to investigate the cerebellar contribution to autism-like phenotypes. Here, we explore cerebellar structural and physiological abnormalities in animals with Prickle2 disruption using immunohistochemistry, whole-cell patch clamp electrophysiology, and several cerebellar-associated motor and timing tasks, including interval timing and eyeblink conditioning. Histologically, Prickle2−/− mice have significantly more empty spaces or gaps between Purkinje cells in the posterior lobules and a decreased propensity for Purkinje Cells to fire action potentials. These structural cerebellar abnormalities did not impair cerebellar-associated behaviors as eyeblink conditioning and interval timing remained intact. Therefore, although Prickle−/− mice show classic phenotypes of ASD, they do not recapitulate the involvement of the adult cerebellum and may not represent the pathophysiological heterogeneity of the disorder.