Machine-learning accelerated structure search for ligand-protected clusters.

Finding low-energy structures of ligand-protected clusters is challenging due to the enormous conformational space and the high computational cost of accurate quantum chemical methods for determining the structures and energies of conformers. Here, we adopted and utilized a kernel rigid regression based machine learning method to accelerate the search for low-energy structures of ligand-protected clusters. We chose the Au25(Cys)18 (Cys: cysteine) cluster as a model system to test and demonstrate our method. We found that the low-energy structures of the cluster are characterized by a specific hydrogen bond type in the cysteine. The different configurations of the ligand layer influence the structural and electronic properties of clusters.

Elucidation of protein-ligand interactions by multiple trajectory analysis methods.

The identification of interaction between protein and ligand including binding positions and strength plays a critical role in drug discovery. Molecular docking and molecular dynamics (MD) techniques have been widely applied to predict binding positions and binding affinity. However, there are few works that describe the systematic exploration of the MD trajectory evolution in this context, potentially leaving out important information. To address the problem, we build a framework, Moira (molecular dynamics trajectory analysis), which enables automating the whole process ranging from docking, MD simulations and various analyses as well as visualizations. We utilized Moira to analyze 400 MD simulations in terms of their geometric features (root mean square deviation and protein-ligand interaction profiler) and energetics (molecular mechanics Poisson-Boltzmann surface area) for these trajectories. Finally, we demonstrate the performance of different analysis techniques in distinguishing native poses among four poses.

Exploring the Conformers of an Organic Molecule on a Metal Cluster with Bayesian Optimization.

Finding low-energy conformers of organic molecules is a complex problem due to the flexibilities of the molecules and the high dimensionality of the search space. When such molecules are on nanoclusters, the search complexity is exacerbated by constraints imposed by the presence of the cluster and other surrounding molecules. To address this challenge, we modified our previously developed active learning molecular conformer search method based on Bayesian optimization and density functional theory. Especially, we have developed and tested strategies to avoid steric clashes between a molecule and a cluster. In this work, we chose a cysteine molecule on a well-studied gold-thiolate cluster as a model system to test and demonstrate our method. We found that cysteine conformers in a cluster inherit the hydrogen bond types from isolated conformers. However, the energy rankings and spacings between the conformers are reordered.

Molecular Conformer Search with Low-Energy Latent Space.

Identifying low-energy conformers with quantum mechanical accuracy for molecules with many degrees of freedom is challenging. In this work, we use the molecular dihedral angles as features and explore the possibility of performing molecular conformer search in a latent space with a generative model named variational auto-encoder (VAE). We bias the VAE towards low-energy molecular configurations to generate more informative data. In this way, we can effectively build a reliable energy model for the low-energy potential energy surface. After the energy model has been built, we extract local-minimum conformations and refine them with structure optimization. We have tested and benchmarked our low-energy latent-space (LOLS) structure search method on organic molecules with 5-9 searching dimensions. Our results agree with previous studies.

Efficient Amino Acid Conformer Search with Bayesian Optimization.

Finding low-energy molecular conformers is challenging due to the high dimensionality of the search space and the computational cost of accurate quantum chemical methods for determining conformer structures and energies. Here, we combine active-learning Bayesian optimization (BO) algorithms with quantum chemistry methods to address this challenge. Using cysteine as an example, we show that our procedure is both efficient and accurate. After only 1000 single-point calculations and approximately 80 structure relaxations, which is less than 10% computational cost of the current fastest method, we have found the low-energy conformers in good agreement with experimental measurements and reference calculations. To test the transferability of our method, we also repeated the conformer search of serine, tryptophan, and aspartic acid. The results agree well with previous conformer search studies.

Lithium adsorption and migration in group IV-VI compounds and GeS/graphene heterostructures: a comparative study.

By means of first-principles calculations, the adsorption and transport properties of lithium (Li) in orthorhombic group IV-VI compounds MX (M = Ge, Sn; X = S, Se) and GeS/graphene heterostructures have been systematically investigated. Strong interactions and distinct charge transfer between Li and compounds MX are observed. The Li diffusion barriers along the zigzag direction are found to be much lower than that along the armchair direction in monolayer and bulk MX, showing distinct anisotropic diffusion features. In particular, monolayer GeS has a lowest barrier of 0.173 eV (zigzag) among them and it will transit from a semiconductor to a metallic state after Li intercalation, indicating fast Li and electron transport properties. As a comparison, the addition of graphene in a GeS/graphene heterostructure could enhance its binding with Li, decrease the Li diffusion barrier and inhibit the volume expansion dramatically, suggesting a potential performance improvement. Our study not only reveals the directional transport properties of Li in MX, but also improves the understanding of the role of graphene in the MX/graphene heterostructure, and shows great potential application in the field of electrode materials.