Recognizing Local and Global Structural Motifs at the Atomic Scale

Beyond Local Structures in Critical Supercooled Water through Unsupervised Learning

The presence of a second critical point in water has been a topic of intense investigation for the last few decades. The molecular origins underlying this phenomenon are typically rationalized in terms of the competition between local high-density (HD) and low-density (LD) structures. Their identification often requires designing parameters that are subject to human intervention. Herein, we use unsupervised learning to discover structures in atomistic simulations of water close to the liquid-liquid critical point (LLCP). Encoding the information on the environment using local descriptors, we do not find evidence for two distinct thermodynamic structures. In contrast, when we deploy nonlocal descriptors that probe instead heterogeneities on the nanometer length scale, this leads to the emergence of LD and HD domains rationalizing the microscopic origins of the density fluctuations close to criticality.

Mapping the electronic transitions of protonation sites in peptides using soft X-ray action spectroscopy

Near-edge X-ray absorption mass spectrometry (NEXAMS) around the nitrogen and oxygen K-edges was employed on gas-phase peptides to probe the electronic transitions related to their protonation sites, namely at basic side chains, the N-terminus and the amide oxygen. The experimental results are supported by replica exchange molecular dynamics and density-functional theory and restricted open-shell configuration with single calculations to attribute the transitions responsible for the experimentally observed resonances. We studied five tailor-made glycine-based pentapeptides, where we identified the signature of the protonation site of N-terminal proline, histidine, lysine and arginine, at 406 eV, corresponding to N 1s → σ*(NHx+) (x = 2 or 3) transitions, depending on the peptides. We compared the spectra of pentaglycine and triglycine to evaluate the sensitivity of NEXAMS to protomers. Separate resonances have been identified to distinguish two protomers in triglycine, the protonation site at the N-terminus at 406 eV and the protonation site at the amide oxygen characterized by a transition at 403.1 eV.

Detecting dynamic domains and local fluctuations in complex molecular systems via timelapse neighbors shuffling

It is known that the behavior of many complex systems is controlled by local dynamic rearrangements or fluctuations occurring within them. Complex molecular systems, composed of many molecules interacting with each other in a Brownian storm, make no exception. Despite the rise of machine learning and of sophisticated structural descriptors, detecting local fluctuations and collective transitions in complex dynamic ensembles remains often difficult. Here, we show a machine learning framework based on a descriptor which we name Local Environments and Neighbors Shuffling (LENS), that allows identifying dynamic domains and detecting local fluctuations in a variety of systems in an abstract and efficient way. By tracking how much the microscopic surrounding of each molecular unit changes over time in terms of neighbor individuals, LENS allows characterizing the global (macroscopic) dynamics of molecular systems in phase transition, phases-coexistence, as well as intrinsically characterized by local fluctuations (e.g., defects). Statistical analysis of the LENS time series data extracted from molecular dynamics trajectories of, for example, liquid-like, solid-like, or dynamically diverse complex molecular systems allows tracking in an efficient way the presence of different dynamic domains and of local fluctuations emerging within them. The approach is found robust, versatile, and applicable independently of the features of the system and simply provided that a trajectory containing information on the relative motion of the interacting units is available. We envisage that "such a LENS" will constitute a precious basis for exploring the dynamic complexity of a variety of systems and, given its abstract definition, not necessarily of molecular ones.

Into the dynamics of rotaxanes at atomistic resolution

Mechanically-interlocked molecules (MIMs) are at the basis of artificial molecular machines and are attracting increasing interest for various applications, from catalysis to drug delivery and nanoelectronics. MIMs are composed of mechanically-interconnected molecular sub-parts that can move with respect to each other, imparting these systems innately dynamical behaviors and interesting stimuli-responsive properties. The rational design of MIMs with desired functionalities requires studying their dynamics at sub-molecular resolution and on relevant timescales, which is challenging experimentally and computationally. Here, we combine molecular dynamics and metadynamics simulations to reconstruct the thermodynamics and kinetics of different types of MIMs at atomistic resolution under different conditions. As representative case studies, we use rotaxanes and molecular shuttles substantially differing in structure, architecture, and dynamical behavior. Our computational approach provides results in agreement with the available experimental evidence and a direct demonstration of the critical effect of the solvent on the dynamics of the MIMs. At the same time, our simulations unveil key factors controlling the dynamics of these systems, providing submolecular-level insights into the mechanisms and kinetics of shuttling. Reconstruction of the free-energy profiles from the simulations reveals details of the conformations of macrocycles on the binding site that are difficult to access via routine experiments and precious for understanding the MIMs' behavior, while their decomposition in enthalpic and entropic contributions unveils the mechanisms and key transitions ruling the intermolecular movements between metastable states within them. The computational framework presented herein is flexible and can be used, in principle, to study a variety of mechanically-interlocked systems.

Integrating Machine Learning in the Coarse-Grained Molecular Simulation of Polymers

Machine learning (ML) is having an increasing impact on the physical sciences, engineering, and technology and its integration into molecular simulation frameworks holds great potential to expand their scope of applicability to complex materials and facilitate fundamental knowledge and reliable property predictions, contributing to the development of efficient materials design routes. The application of ML in materials informatics in general, and polymer informatics in particular, has led to interesting results, however great untapped potential lies in the integration of ML techniques into the multiscale molecular simulation methods for the study of macromolecular systems, specifically in the context of Coarse Grained (CG) simulations. In this Perspective, we aim at presenting the pioneering recent research efforts in this direction and discussing how these new ML-based techniques can contribute to critical aspects of the development of multiscale molecular simulation methods for bulk complex chemical systems, especially polymers. Prerequisites for the implementation of such ML-integrated methods and open challenges that need to be met toward the development of general systematic ML-based coarse graining schemes for polymers are discussed.

Unsupervised Data-Driven Reconstruction of Molecular Motifs in Simple to Complex Dynamic Micelles

The reshuffling mobility of molecular building blocks in self-assembled micelles is a key determinant of many their interesting properties, from emerging morphologies and surface compartmentalization, to dynamic reconfigurability and stimuli-responsiveness. However, the microscopic details of such complex structural dynamics are typically nontrivial to elucidate, especially in multicomponent assemblies. Here we show a machine-learning approach that allows us to reconstruct the structural and dynamic complexity of mono- and bicomponent surfactant micelles from high-dimensional data extracted from equilibrium molecular dynamics simulations. Unsupervised clustering of smooth overlap of atomic position (SOAP) data enables us to identify, in a set of multicomponent surfactant micelles, the dominant local molecular environments that emerge within them and to retrace their dynamics, in terms of exchange probabilities and transition pathways of the constituent building blocks. Tested on a variety of micelles differing in size and in the chemical nature of the constitutive self-assembling units, this approach effectively recognizes the molecular motifs populating them in an exquisitely agnostic and unsupervised way, and allows correlating them to their composition in terms of constitutive surfactant species.

Reconstructing the infrared spectrum of a peptide from representative conformers of the full canonical ensemble

Leucine enkephalin (LeuEnk), a biologically active endogenous opioid pentapeptide, has been under intense investigation because it is small enough to allow efficient use of sophisticated computational methods and large enough to provide insights into low-lying minima of its conformational space. Here, we reproduce and interpret experimental infrared (IR) spectra of this model peptide in gas phase using a combination of replica-exchange molecular dynamics simulations, machine learning, and ab initio calculations. In particular, we evaluate the possibility of averaging representative structural contributions to obtain an accurate computed spectrum that accounts for the corresponding canonical ensemble of the real experimental situation. Representative conformers are identified by partitioning the conformational phase space into subensembles of similar conformers. The IR contribution of each representative conformer is calculated from ab initio and weighted according to the population of each cluster. Convergence of the averaged IR signal is rationalized by merging contributions in a hierarchical clustering and the comparison to IR multiple photon dissociation experiments. The improvements achieved by decomposing clusters containing similar conformations into even smaller subensembles is strong evidence that a thorough assessment of the conformational landscape and the associated hydrogen bonding is a prerequisite for deciphering important fingerprints in experimental spectroscopic data.

Supramolecular Semiconductivity through Emerging Ionic Gates in Ion-Nanoparticle Superlattices

The self-assembly of nanoparticles driven by small molecules or ions may produce colloidal superlattices with features and properties reminiscent of those of metals or semiconductors. However, to what extent the properties of such supramolecular crystals actually resemble those of atomic materials often remains unclear. Here, we present coarse-grained molecular simulations explicitly demonstrating how a behavior evocative of that of semiconductors may emerge in a colloidal superlattice. As a case study, we focus on gold nanoparticles bearing positively charged groups that self-assemble into FCC crystals via mediation by citrate counterions. In silico ohmic experiments show how the dynamically diverse behavior of the ions in different superlattice domains allows the opening of conductive ionic gates above certain levels of applied electric fields. The observed binary conductive/nonconductive behavior is reminiscent of that of conventional semiconductors, while, at a supramolecular level, crossing the "band gap" requires a sufficient electrostatic stimulus to break the intermolecular interactions and make ions diffuse throughout the superlattice's cavities.

Spotting Local Environments in Self-Assembled Monolayer-Protected Gold Nanoparticles

Organic-inorganic (O-I) nanomaterials are versatile platforms for an incredible high number of applications, ranging from heterogeneous catalysis to molecular sensing, cell targeting, imaging, and cancer diagnosis and therapy, just to name a few. Much of their potential stems from the unique control of organic environments around inorganic sites within a single O-I nanomaterial, which allows for new properties that were inaccessible using purely organic or inorganic materials. Structural and mechanistic characterization plays a key role in understanding and rationally designing such hybrid nanoconstructs. Here, we introduce a general methodology to identify and classify local (supra)molecular environments in an archetypal class of O-I nanomaterials, i.e., self-assembled monolayer-protected gold nanoparticles (SAM-AuNPs). By using an atomistic machine-learning guided workflow based on the Smooth Overlap of Atomic Positions (SOAP) descriptor, we analyze a collection of chemically different SAM-AuNPs and detect and compare local environments in a way that is agnostic and automated, i.e., with no need of a priori information and minimal user intervention. In addition, the computational results coupled with experimental electron spin resonance measurements prove that is possible to have more than one local environment inside SAMs, being the thickness of the organic shell and solvation primary factors in the determining number and nature of multiple coexisting environments. These indications are extended to complex mixed hydrophilic-hydrophobic SAMs. This work demonstrates that it is possible to spot and compare local molecular environments in SAM-AuNPs exploiting atomistic machine-learning approaches, establishes ground rules to control them, and holds the potential for the rational design of O-I nanomaterials instructed from data.

Classifying soft self-assembled materials via unsupervised machine learning of defects

Unlike molecular crystals, soft self-assembled fibers, micelles, vesicles, etc., exhibit a certain order in the arrangement of their constitutive monomers but also high structural dynamicity and variability. Defects and disordered local domains that continuously form-and-repair in their structures impart to such materials unique adaptive and dynamical properties, which make them, e.g., capable to communicate with each other. However, objective criteria to compare such complex dynamical features and to classify soft supramolecular materials are non-trivial to attain. Here we show a data-driven workflow allowing us to achieve this goal. Building on unsupervised clustering of Smooth Overlap of Atomic Position (SOAP) data obtained from equilibrium molecular dynamics simulations, we can compare a variety of soft supramolecular assemblies via a robust SOAP metric. This provides us with a data-driven "defectometer" to classify different types of supramolecular materials based on the structural dynamics of the ordered/disordered local molecular environments that statistically emerge within them.

High-Dimensional Fluctuations in Liquid Water: Combining Chemical Intuition with Unsupervised Learning

The microscopic description of the local structure of water remains an open challenge. Here, we adopt an agnostic approach to understanding water's hydrogen bond network using data harvested from molecular dynamics simulations of an empirical water model. A battery of state-of-the-art unsupervised data-science techniques are used to characterize the free-energy landscape of water starting from encoding the water environment using local atomic descriptors, through dimensionality reduction and finally the use of advanced clustering techniques. Analysis of the free energy under ambient conditions was found to be consistent with a rough single basin and independent of the choice of the water model. We find that the fluctuations of the water network occur in a high-dimensional space, which we characterize using a combination of both atomic descriptors and chemical-intuition-based coordinates. We demonstrate that a combination of both types of variables is needed in order to adequately capture the complexity of the fluctuations in the hydrogen bond network at different length scales both at room temperature and also close to the critical point of water. Our results provide a general framework for examining fluctuations in water under different conditions.

Local Kernel Regression and Neural Network Approaches to the Conformational Landscapes of Oligopeptides

The application of machine learning to theoretical chemistry has made it possible to combine the accuracy of quantum chemical energetics with the thorough sampling of finite-temperature fluctuations. To reach this goal, a diverse set of methods has been proposed, ranging from simple linear models to kernel regression and highly nonlinear neural networks. Here we apply two widely different approaches to the same, challenging problem: the sampling of the conformational landscape of polypeptides at finite temperature. We develop a local kernel regression (LKR) coupled with a supervised sparsity method and compare it with a more established approach based on Behler-Parrinello type neural networks. In the context of the LKR, we discuss how the supervised selection of the reference pool of environments is crucial to achieve accurate potential energy surfaces at a competitive computational cost and leverage the locality of the model to infer which chemical environments are poorly described by the DFTB baseline. We then discuss the relative merits of the two frameworks and perform Hamiltonian-reservoir replica-exchange Monte Carlo sampling and metadynamics simulations, respectively, to demonstrate that both frameworks can achieve converged and transferable sampling of the conformational landscape of complex and flexible biomolecules with comparable accuracy and computational cost.

Electrostatic co-assembly of nanoparticles with oppositely charged small molecules into static and dynamic superstructures

Coulombic interactions can be used to assemble charged nanoparticles into higher-order structures, but the process requires oppositely charged partners that are similarly sized. The ability to mediate the assembly of such charged nanoparticles using structurally simple small molecules would greatly facilitate the fabrication of nanostructured materials and harnessing their applications in catalysis, sensing and photonics. Here we show that small molecules with as few as three electric charges can effectively induce attractive interactions between oppositely charged nanoparticles in water. These interactions can guide the assembly of charged nanoparticles into colloidal crystals of a quality previously only thought to result from their co-crystallization with oppositely charged nanoparticles of a similar size. Transient nanoparticle assemblies can be generated using positively charged nanoparticles and multiply charged anions that are enzymatically hydrolysed into mono- and/or dianions. Our findings demonstrate an approach for the facile fabrication, manipulation and further investigation of static and dynamic nanostructured materials in aqueous environments.

Controlling Exchange Pathways in Dynamic Supramolecular Polymers by Controlling Defects

Supramolecular fibers composed of monomers that self-assemble directionally via noncovalent interactions are ubiquitous in nature, and of great interest in chemistry. In these structures, the constitutive monomers continuously exchange in-and-out the assembly according to a well-defined supramolecular equilibrium. However, unraveling the exchange pathways and their molecular determinants constitutes a nontrivial challenge. Here, we combine coarse-grained modeling, enhanced sampling, and machine learning to investigate the key factors controlling the monomer exchange pathways in synthetic supramolecular polymers having an intrinsic dynamic behavior. We demonstrate how the competition of directional vs. nondirectional interactions between the monomers controls the creation/annihilation of defects in the supramolecular polymers, from where monomers exchange proceeds. This competition determines the exchange pathway, dictating whether a fiber statistically swaps monomers from the tips or from all along its length. Finally, thanks to their generality, our models allow the investigation of molecular approaches to control the exchange pathways in these dynamic assemblies.

Physics-Inspired Structural Representations for Molecules and Materials

The first step in the construction of a regression model or a data-driven analysis, aiming to predict or elucidate the relationship between the atomic-scale structure of matter and its properties, involves transforming the Cartesian coordinates of the atoms into a suitable representation. The development of atomic-scale representations has played, and continues to play, a central role in the success of machine-learning methods for chemistry and materials science. This review summarizes the current understanding of the nature and characteristics of the most commonly used structural and chemical descriptions of atomistic structures, highlighting the deep underlying connections between different frameworks and the ideas that lead to computationally efficient and universally applicable models. It emphasizes the link between properties, structures, their physical chemistry, and their mathematical description, provides examples of recent applications to a diverse set of chemical and materials science problems, and outlines the open questions and the most promising research directions in the field.

A Data-Driven Dimensionality Reduction Approach to Compare and Classify Lipid Force Fields

Molecular dynamics simulations of all-atom and coarse-grained lipid bilayer models are increasingly used to obtain useful insights for understanding the structural dynamics of these assemblies. In this context, one crucial point concerns the comparison of the performance and accuracy of classical force fields (FFs), which sometimes remains elusive. To date, the assessments performed on different classical potentials are mostly based on the comparison with experimental observables, which typically regard average properties. However, local differences of the structure and dynamics, which are poorly captured by average measurements, can make a difference, but these are nontrivial to catch. Here, we propose an agnostic way to compare different FFs at different resolutions (atomistic, united-atom, and coarse-grained), by means of a high-dimensional similarity metrics built on the framework of Smooth Overlap of Atomic Position (SOAP). We compare and classify a set of 13 FFs, modeling 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. Our SOAP kernel-based metrics allows us to compare, discriminate, and correlate different FFs at different model resolutions in an unbiased, high-dimensional way. This also captures differences between FFs in modeling nonaverage events (originating from local transitions), for example, the liquid-to-gel phase transition in dipalmitoylphosphatidylcholine (DPPC) bilayers, for which our metrics allows us to identify nucleation centers for the phase transition, highlighting some intrinsic resolution limitations in implicit versus explicit solvent FFs.

Living supramolecular polymerization of fluorinated cyclohexanes

The development of powerful methods for living covalent polymerization has been a key driver of progress in organic materials science. While there have been remarkable reports on living supramolecular polymerization recently, the scope of monomers is still narrow and a simple solution to the problem is elusive. Here we report a minimalistic molecular platform for living supramolecular polymerization that is based on the unique structure of all-cis 1,2,3,4,5,6-hexafluorocyclohexane, the most polar aliphatic compound reported to date. We use this large dipole moment (6.2 Debye) not only to thermodynamically drive the self-assembly of supramolecular polymers, but also to generate kinetically trapped monomeric states. Upon addition of well-defined seeds, we observed that the dormant monomers engage in a kinetically controlled supramolecular polymerization. The obtained nanofibers have an unusual double helical structure and their length can be controlled by the ratio between seeds and monomers. The successful preparation of supramolecular block copolymers demonstrates the versatility of the approach.

Global Free-Energy Landscapes as a Smoothly Joined Collection of Local Maps

Enhanced sampling techniques have become an essential tool in computational chemistry and physics, where they are applied to sample activated processes that occur on a time scale that is inaccessible to conventional simulations. Despite their popularity, it is well known that they have constraints that hinder their application to complex problems. The core issue lies in the need to describe the system using a small number of collective variables (CVs). Any slow degree of freedom that is not properly described by the chosen CVs will hinder sampling efficiency. However, the exploration of configuration space is also hampered by including variables that are not relevant for the activated process under study. This paper presents the Adaptive Topography of Landscape for Accelerated Sampling (ATLAS), a new biasing method capable of working with many CVs. The root idea of ATLAS is to apply a divide-and-conquer strategy, where the high-dimensional CVs space is divided into basins, each of which is described by an automatically determined, low-dimensional set of variables. A well-tempered metadynamics-like bias is constructed as a function of these local variables. Indicator functions associated with the basins switch on and off the local biases so that the sampling is performed on a collection of low-dimensional CV spaces that are smoothly combined to generate an effectively high-dimensional bias. The unbiased Boltzmann distribution is recovered through reweighing, making the evaluation of conformational and thermodynamic properties straightforward. The decomposition of the free-energy landscape in local basins can be updated iteratively as the simulation discovers new (meta)stable states.

Mapping the Structure of Oxygen-Doped Wurtzite Aluminum Nitride Coatings from Ab Initio Random Structure Search and Experiments

Machine learning is changing how we design and interpret experiments in materials science. In this work, we show how unsupervised learning, combined with ab initio random structure searching, improves our understanding of structural metastability in multicomponent alloys. We focus on the case of Al-O-N alloys where the formation of aluminum vacancies in wurtzite AlN upon the incorporation of substitutional oxygen can be seen as a general mechanism of solids where crystal symmetry is reduced to stabilize defects. The ideal AlN wurtzite crystal structure occupation cannot be matched due to the presence of an aliovalent hetero-element into the structure. The traditional interpretation of the c-lattice shrinkage in sputter-deposited Al-O-N films from X-ray diffraction (XRD) experiments suggests the existence of a solubility limit at 8 at % oxygen content. Here, we show that such naive interpretation is misleading. We support XRD data with accurate ab initio modeling and dimensionality reduction on advanced structural descriptors to map structure-property relationships. No signs of a possible solubility limit are found. Instead, the presence of a wide range of non-equilibrium oxygen-rich defective structures emerging at increasing oxygen contents suggests that the formation of grain boundaries is the most plausible mechanism responsible for the lattice shrinkage measured in Al-O-N sputtered films. We further confirm our hypothesis using positron annihilation lifetime spectroscopy.

Big-Data Science in Porous Materials: Materials Genomics and Machine Learning

By combining metal nodes with organic linkers we can potentially synthesize millions of possible metal-organic frameworks (MOFs). The fact that we have so many materials opens many exciting avenues but also create new challenges. We simply have too many materials to be processed using conventional, brute force, methods. In this review, we show that having so many materials allows us to use big-data methods as a powerful technique to study these materials and to discover complex correlations. The first part of the review gives an introduction to the principles of big-data science. We show how to select appropriate training sets, survey approaches that are used to represent these materials in feature space, and review different learning architectures, as well as evaluation and interpretation strategies. In the second part, we review how the different approaches of machine learning have been applied to porous materials. In particular, we discuss applications in the field of gas storage and separation, the stability of these materials, their electronic properties, and their synthesis. Given the increasing interest of the scientific community in machine learning, we expect this list to rapidly expand in the coming years.

Crystal Structure Influences Migration along Li and Mg Surfaces

Dendrite formation on Li metal anodes hinders commercialization of more energy-dense rechargeable batteries. Here, we use the migration energy barrier (MEB) for surface transport as a descriptor for dendrite nucleation and compare Li to Mg. Density functional theory calculations show that the MEB for the hexagonal close-packed structure is 40 and 270 meV lower than that of the body-centered cubic structure for Li and Mg, respectively. This is suggested as a reason why Mg surfaces are less prone to form dendrites than Li. We show that the close-packed facets exhibit lower MEBs because of smaller changes in atomic coordination during migration and thereby less surface distortion.

Identifying and Tracking Defects in Dynamic Supramolecular Polymers

A central paradigm of self-assembly is to create ordered structures starting from molecular monomers that spontaneously recognize and interact with each other via noncovalent interactions. In recent years, great efforts have been directed toward perfecting the design of a variety of supramolecular polymers and materials with different architectures. The resulting structures are often thought of as ideally perfect, defect-free supramolecular fibers, micelles, vesicles, etc., having an intrinsic dynamic character, which are typically studied at the level of statistical ensembles to assess their average properties. However, molecular simulations recently demonstrated that local defects that may be present or may form in these assemblies, and which are poorly captured by conventional approaches, are key to controlling their dynamic behavior and properties. The study of these defects poses considerable challenges, as the flexible/dynamic nature of these soft systems makes it difficult to identify what effectively constitutes a defect and to characterize its stability and evolution. Here, we demonstrate the power of unsupervised machine-learning techniques to systematically identify and compare defects in supramolecular polymer variants in different conditions, using as a benchmark 5 Å resolution coarse-grained molecular simulations of a family of supramolecular polymers. We show that this approach allows a complete data-driven characterization of the internal structure and dynamics of these complex assemblies and of the dynamic pathways for defects formation and resorption. This provides a useful, generally applicable approach to unambiguously identify defects in these dynamic self-assembled materials and to classify them based on their structure, stability, and dynamics.

Deep learning for automated classification and characterization of amorphous materials

It is difficult to quantify structure-property relationships and to identify structural features of complex materials. The characterization of amorphous materials is especially challenging because their lack of long-range order makes it difficult to define structural metrics. In this work, we apply deep learning algorithms to accurately classify amorphous materials and characterize their structural features. Specifically, we show that convolutional neural networks and message passing neural networks can classify two-dimensional liquids and liquid-cooled glasses from molecular dynamics simulations with greater than 0.98 AUC, with no a priori assumptions about local particle relationships, even when the liquids and glasses are prepared at the same inherent structure energy. Furthermore, we demonstrate that message passing neural networks surpass convolutional neural networks in this context in both accuracy and interpretability. We extract a clear interpretation of how message passing neural networks evaluate liquid and glass structures by using a self-attention mechanism. Using this interpretation, we derive three novel structural metrics that accurately characterize glass formation. The methods presented here provide a procedure to identify important structural features in materials that could be missed by standard techniques and give unique insight into how these neural networks process data.

Deep Representation Learning for Complex Free-Energy Landscapes

In this Letter, we analyzed the inductive bias underlying complex free-energy landscapes (FELs) and exploited it to train deep neural networks that yield reduced and clustered representation for the FEL. Our parametric method, called information distilling of metastability (IDM), is end-to-end differentiable and thus scalable to ultralarge data sets. IDM is able to perform clustering in the meantime of reducing the dimensionality. Besides, as an unsupervised learning method, IDM differs from many existing dimensionality reduction and clustering methods in that it requires neither a cherry-picked distance metric nor the ground-true number of clusters defined a priori, and it can be used to unroll and zoom in on the hierarchical FEL with respect to different time scales. Through multiple experiments, we show that IDM can achieve physically meaningful representations that partition the FEL into well-defined metastable states that hence are amenable for downstream tasks such as mechanism analysis and kinetic modeling.

Using Dimensionality Reduction to Analyze Protein Trajectories

In recent years the analysis of molecular dynamics trajectories using dimensionality reduction algorithms has become commonplace. These algorithms seek to find a low-dimensional representation of a trajectory that is, according to a well-defined criterion, optimal. A number of different strategies for generating projections of trajectories have been proposed but little has been done to systematically compare how these various approaches fare when it comes to analysing trajectories for biomolecules in explicit solvent. In the following paper, we have thus analyzed a molecular dynamics trajectory of the C-terminal fragment of the immunoglobulin binding domain B1 of protein G of Streptococcus modeled in explicit solvent using a range of different dimensionality reduction algorithms. We have then tried to systematically compare the projections generated using each of these algorithms by using a clustering algorithm to find the positions and extents of the basins in the high-dimensional energy landscape. We find that no algorithm outshines all the other in terms of the quality of the projection it generates. Instead, all the algorithms do a reasonable job when it comes to building a projection that separates some of the configurations that lie in different basins. Having said that, however, all the algorithms struggle to project the basins because they all have a large intrinsic dimensionality.

Unsupervised machine learning in atomistic simulations, between predictions and understanding

Automated analyses of the outcome of a simulation have been an important part of atomistic modeling since the early days, addressing the need of linking the behavior of individual atoms and the collective properties that are usually the final quantity of interest. Methods such as clustering and dimensionality reduction have been used to provide a simplified, coarse-grained representation of the structure and dynamics of complex systems from proteins to nanoparticles. In recent years, the rise of machine learning has led to an even more widespread use of these algorithms in atomistic modeling and to consider different classification and inference techniques as part of a coherent toolbox of data-driven approaches. This perspective briefly reviews some of the unsupervised machine-learning methods-that are geared toward classification and coarse-graining of molecular simulations-seen in relation to the fundamental mathematical concepts that underlie all machine-learning techniques. It discusses the importance of using concise yet complete representations of atomic structures as the starting point of the analyses and highlights the risk of introducing preconceived biases when using machine learning to rationalize and understand structure-property relations. Supervised machine-learning techniques that explicitly attempt to predict the properties of a material given its structure are less susceptible to such biases. Current developments in the field suggest that using these two classes of approaches side-by-side and in a fully integrated mode, while keeping in mind the relations between the data analysis framework and the fundamental physical principles, will be key to realizing the full potential of machine learning to help understand the behavior of complex molecules and materials.

Atomic Motif Recognition in (Bio)Polymers: Benchmarks From the Protein Data Bank

Rationalizing the structure and structure-property relations for complex materials such as polymers or biomolecules relies heavily on the identification of local atomic motifs, e.g., hydrogen bonds and secondary structure patterns, that are seen as building blocks of more complex supramolecular and mesoscopic structures. Over the past few decades, several automated procedures have been developed to identify these motifs in proteins given the atomic structure. Being based on a very precise understanding of the specific interactions, these heuristic criteria formulate the question in a way that implies the answer, by defining a list of motifs based on those that are known to be naturally occurring. This makes them less likely to identify unexpected phenomena, such as the occurrence of recurrent motifs in disordered segments of proteins, and less suitable to be applied to different polymers whose structure is not driven by hydrogen bonds, or even to polypeptides when appearing in unusual, non-biological conditions. Here we discuss how unsupervised machine learning schemes can be used to recognize patterns based exclusively on the frequency with which different motifs occur, taking high-resolution structures from the Protein Data Bank as benchmarks. We first discuss the application of a density-based motif recognition scheme in combination with traditional representations of protein structure (namely, interatomic distances and backbone dihedrals). Then, we proceed one step further toward an entirely unbiased scheme by using as input a structural representation based on the atomic density and by employing supervised classification to objectively assess the role played by the representation in determining the nature of atomic-scale patterns.

Using Data-Reduction Techniques to Analyze Biomolecular Trajectories

This chapter discusses the way in which dimensionality reduction algorithms such as diffusion maps and sketch-map can be used to analyze molecular dynamics trajectories. The first part discusses how these various algorithms function as well as practical issues such as landmark selection and how these algorithms can be used when the data to be analyzed comes from enhanced sampling trajectories. In the later part a comparison between the results obtained by applying various algorithms to two sets of sample data is performed and discussed. This section is then followed by a summary of how one algorithm in particular, sketch-map, has been applied to a range of problems. The chapter concludes with a discussion on the directions that we believe this field is currently moving.