Free-energy calculations in multiple dimensions constitute a challenging problem, owing to the significant computational cost incurred to achieve ergodic sampling. The generalized adaptive biasing force (gABF) algorithm calculates n one-dimensional lists of biasing forces to approximate the n-dimensional matrix by ignoring the coupling terms ordinarily taken into account in classical ABF simulations, thereby greatly accelerating sampling in the multidimensional space. This approximation may however occasionally lead to poor, incomplete exploration of the conformational space compared to classical ABF, especially when the selected coarse variables are strongly coupled. It has been found that introducing extended potentials coupled to the coarse variables of interest can virtually eliminate this shortcoming, and, thus, improve the efficiency of gABF simulations. In the present contribution, we propose a new free-energy method, coined extended generalized ABF (egABF), combining gABF with an extended Lagrangian strategy. The results for three illustrative examples indicate that (i) egABF can explore the transition coordinate much more efficiently compared with classical ABF, eABF, and gABF, in both simple and complex cases and (ii) egABF can achieve a higher accuracy than gABF, with a root mean-squared deviation between egABF and eABF free-energy profiles on the order of kT. Furthermore, the new egABF algorithm outruns the previous ABF-based algorithms in high-dimensional free-energy calculations and, hence, represents a powerful importance-sampling alternative for the investigation of complex chemical and biological processes.
In cyclodextrin (CD)-based rotaxanes, the shuttling rate of the macrocycle along the thread is crucial to characterize their function as molecular machines. In general, the composition of the thread and the environment are considered to be important factors affecting the nature of the movement. Yet, the role of ancillary motions on the shuttling rate remains unclear. In the present contribution, two rotaxanes having the same components, yet significantly different shuttling rates between two stable states in an aqueous environment, have been investigated at the atomic level using numerical simulations. These two rotaxanes consist of an axle with two stations linked by a 2-methylpyridinium group and an α-CD sliding on the axle and assuming two different orientations. We found that a number of cyclodextrin glucopyranose units (GLUs) isomerized during shuttling, which we anticipate to affect the shuttling rate. The two-dimensional free-energy landscapes characterizing the isomerization of the GLUs and the shuttling along the thread were mapped and revealed that the energetic barriers hampering spontaneous transition between the two stations significantly differ for the two rotaxanes. Structural analysis shows that this difference mainly arises from steric hindrances caused by the methyl substituent of the pyridinium group, which leads to a different number of the GLUs experiencing conformational change during shuttling. Moreover, the thermodynamic stability of the complex is found to be distinct between the two rotaxanes. This discrepancy may be ascribed to the dipole moment of the complex, which is sensitive to the orientation of CD. It can be concluded that shuttling in the rotaxanes is not only highly coupled with isomerization of GLUs but also affected by thermodynamic stability, resulting in a shuttling rate sensitive to the orientation of the CD. The present results help understand the complex molecular motion in CD-based molecular shuttles, and are expected to serve in the design of molecular filters for selectively screening molecules with a specific orientation.
To enable the fast estimation of atom condensed Fukui functions, machine learning algorithms were trained with databases of DFT pre-calculated values for ca. 23,000 atoms in organic molecules. The problem was approached as the ranking of atom types with the Bradley-Terry (BT) model, and as the regression of the Fukui function. Random Forests (RF) were trained to predict the condensed Fukui function, to rank atoms in a molecule, and to classify atoms as high/low Fukui function. Atomic descriptors were based on counts of atom types in spheres around the kernel atom. The BT coefficients assigned to atom types enabled the identification (93-94 % accuracy) of the atom with the highest Fukui function in pairs of atoms in the same molecule with differences ≥0.1. In whole molecules, the atom with the top Fukui function could be recognized in ca. 50 % of the cases and, on the average, about 3 of the top 4 atoms could be recognized in a shortlist of 4. Regression RF yielded predictions for test sets with R(2) =0.68-0.69, improving the ability of BT coefficients to rank atoms in a molecule. Atom classification (as high/low Fukui function) was obtained with RF with sensitivity of 55-61 % and specificity of 94-95 %.
In order to explore atomic asymmetry and molecular chirality in 2D space, benzenoids composed of 3 to 11 hexagons in 2D space were enumerated in our laboratory. These benzenoids are regarded as planar connected polyhexes and have no internal holes; that is, their internal regions are filled with hexagons. The produced dataset was composed of 357,968 benzenoids, including more than 14 million atoms. Rather than simply labeling the huge number of atoms as being either symmetric or asymmetric, this investigation aims at exploring a quantitative graph theoretical descriptor of atomic asymmetry. Based on the particular characteristics in the 2D plane, we suggested the weighted atomic sum as the descriptor of atomic asymmetry. This descriptor is measured by circulating around the molecule going in opposite directions. The investigation demonstrates that the weighted atomic sums are superior to the previously reported quantitative descriptor, atomic sums. The investigation of quantitative descriptors also reveals that the most asymmetric atom is in a structure with a spiral ring with the convex shape going in clockwise direction and concave shape going in anticlockwise direction from the atom. Based on weighted atomic sums, a weighted F index is introduced to quantitatively represent molecular chirality in the plane, rather than merely regarding benzenoids as being either chiral or achiral. By validating with enumerated benzenoids, the results indicate that the weighted F indexes were in accordance with their chiral classification (achiral or chiral) over the whole benzenoids dataset. Furthermore, weighted F indexes were superior to previously available descriptors. Benzenoids possess a variety of shapes and can be extended to practically represent any shape in 2D space—our proposed descriptor has thus the potential to be a general method to represent 2D molecular chirality based on the difference between clockwise and anticlockwise sums around a molecule.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.