Clustering Molecular Dynamics trajectories is a common analysis that allows grouping together similar conformations. Several algorithms have been designed and optimized to perform this routine task, and among them, Quality Threshold stands as a very attractive option. This algorithm guarantees that in retrieved clusters no pair of frames will have a similarity value greater than a specified threshold, and hence, a set of strongly correlated frames are obtained for each cluster. In this work, it is shown that various commonly used software implementations are flawed by confusing Quality Threshold with another simplistic well-known clustering algorithm published by Daura et al. (Daura, X.; van Gunsteren, W. F.; Jaun, B.; Mark, A. E.; Gademann, K.; Seebach, D. Peptide Folding: When Simulation Meets Experiment. Angew. Chemie Int. Ed. 1999, 38 (1/2), 236–240). Daura’s algorithm does not impose any quality threshold for the frames contained in retrieved clusters, bringing unrelated structural configurations altogether. The advantages of using Quality Threshold whenever possible to explore Molecular Dynamic trajectories is exemplified. An in-house implementation of the original Quality Threshold algorithm has been developed in order to illustrate our comments, and its code is freely available for further use by the scientific community.

The reaction of C 60 with pregnen-20-carboxaldehyde, a biologically active synthetic steroid, by using a 1,3-dipolar cycloaddition reaction (Prato's protocol) results in the formation of pyrrolidine rings bearing a new stereogenic center on the C2 of the five-membered ring. The formation of the fullerene-steroid hybrids proceeds with preference for the Re face of the 1,3-dipole, with formation of a diastereomeric mixture in 73:15 ratio. The investigation of the chiroptical properties of these conjugates allowed determining the absolute configuration of the new fulleropyrrolidines. In addition, a thorough spectroscopical study permitted to determine the structure of the two mono-cycloadducts. The electrochemical properties of the new hybrids were also evaluated by cyclic voltammetry, both systems exhibit three quasi-reversible reduction waves which are cathodically shifted in regard to the parent C 60 . Theoretical calculations help supporting the experimental data. A conformational study combining semiempirical methods and density functional theory has predicted the most stable diastereomer. On the basis of this agreement, a possible reaction mechanism is presented. Additionally, a molecular docking simulation has been carried out using the HIV-1 protease as receptor, thus paving the way to study the possible application of these stereoisomers in biomedicine.

The growing computational capacity allows the investigation of large biomolecular systems by increasingly extensive molecular dynamics simulations. The resulting huge trajectories demand efficient partition methods to discern relevant structural dissimilarity. Clustering algorithms are available to address this task, but their implementations still need to be improved to gain in computational speed and to reduce the consumption of random access memory. We propose the BitClust code which, based on a combination of Python and C programming languages, performs fast structural clustering of long molecular trajectories. BitClust takes advantage of bitwise operations applied to a bit-encoded pairwise similarity matrix. Our approach allowed us to process a half-million frame trajectory in 6 h using less than 35 GB, a task that is not affordable with any of the similar alternatives.

Computational fragment-based approaches are widely used in drug design and discovery. One of their limitations is the lack of performance of docking methods, mainly the scoring functions. With the emergence of fragment-based approaches for singlestranded RNA ligands, we analyze the performance in docking and screening powers of an MCSS-based approach. The performance is evaluated on a benchmark of proteinnucleotide complexes where the four RNA residues are used as fragments. The screening power can be considered the major limiting factor for the fragment-based modeling or design of sequence-selective oligonucleotides. We show that the MCSS sampling is efficient even for such large and flexible fragments. Hybrid solvent models based on some partial explicit representation improve both the docking and screening powers.