Abstract:Various kinds of flexibility have been observed in metal-organic frameworks, which may originate from the topology of the material or the presence of flexible ligands. The construction of free energy profiles describing the full dynamical behavior along the phase transition path is challenging since it is not trivial to identify collective variables able to identify all metastable states along the reaction path. In this work, a systematic three-step protocol to uniquely identify the dominant order parameters for structural transformations in flexible metal-organic frameworks and subsequently construct accurate free energy profiles is presented. Methodologically, this protocol is rooted in the time-structure based independent component analysis (tICA), a well-established statistical modeling technique embedded in the Markov state model methodology and often employed to study protein folding, that allows for the identification of the slowest order parameters characterizing the structural transformation. To ensure an unbiased and systematic identification of these order parameters, the tICA decomposition is performed based on information from a prior replica exchange (RE) simulation, as this technique enhances the sampling along all degrees of freedom of the system simultaneously. From this simulation, the tICA procedure extracts the order parameters-often structural parameters-that characterize the slowest transformations in the material. Subsequently, these order parameters are adopted in traditional enhanced sampling methods such as umbrella sampling, thermodynamic integration, and variationally enhanced sampling to construct accurate free energy profiles capturing the flexibility in these nanoporous materials. In this work, the applicability of this tICA-RE protocol is demonstrated by determining the slowest order parameters in both MIL-53(Al) and CAU-13, which exhibit a strongly different type of flexibility. The obtained free energy profiles as a function of this extracted order parameter are furthermore compared to the profiles obtained when adopting less-suited collective variables, indicating the importance of systematically selecting the relevant order parameters to construct accurate free energy profiles for flexible metal-organic frameworks, which is in correspondence with experimental findings. The method succeeds in mapping the full free energy surface in terms of appropriate collective variables for MOFs exhibiting linker flexibility. For CAU-13, we show the decreased stability of the closed pore phase by systematically adding adsorbed xylene molecules in the framework.

Fast and Anisotropic Flexibility-Rigidity Index

Generalized flexibility-rigidity index

An automated protocol to construct flexibility parameters for classical forcefields: applications to metal-organic frameworks

Buried and accessible surface area control intrinsic protein flexibility

Protein flexibility upon ligand binding: Docking predictions and statistical analysis

Assessment of Four Theoretical Approaches to Predict Protein Flexibility in the Crystal Phase and Solution

Multiscale differential geometry learning for protein flexibility analysis

Dimension Conversion and Scaling of Disordered Protein Chains

Identifying the minimal sets of distance restraints for FRET-assisted protein structural modeling

Advanced atom-level representations for protein flexibility prediction utilizing graph neural networks

A-RFP: An Adaptive Residue Flexibility Prediction Method Improving Protein-ligand Docking Based on Homologous Proteins

Protein flexibility prediction by an all-atom mean-field statistical theory

Flexibility of $Α$-Helices: Results of a Statistical Analysis of Database Protein Structures

3DFlex: determining structure and motion of flexible proteins from cryo-EM

Protocol for Identifying Accurate Collective Variables in Enhanced Molecular Dynamics Simulations for the Description of Structural Transformations in Flexible Metal–Organic Frameworks

Accurate Prediction of Protein Structural Flexibility by Deep Learning Integrating Intricate Atomic Structures and Cryo-EM Density Information

FLEXR GUI: a graphical user interface for multi‐conformer modeling of proteins

A Finite Element framework for computation of protein normal modes and mechanical response

Accurate calculation of absolute protein-ligand binding free energies

Exploring protein functions from structural flexibility using CABS‐flex modeling

A Systematic Analysis of the Effect of Small-Molecule Binding on Protein Flexibility of the Ligand-Binding Sites