Fast, clash-free RNA conformational morphing using molecular junctions

Kinematic Vibrational Entropy Assessment and Analysis of SARS CoV-2 Main Protease

The three-dimensional conformations of a protein influence its function and select for the ligands it can interact with. The total free energy change during protein-ligand complex formation includes enthalphic and entropic components, which together report on the binding affinity and conformational states of the complex. However, determining the entropic contribution is computationally burdensome. Here, we apply kinematic flexibility analysis (KFA) to efficiently estimate vibrational frequencies from static protein and protein-ligand structures. The vibrational frequencies, in turn, determine the vibrational entropies of the structures and their complexes. Our estimates of the vibrational entropy change caused by ligand binding compare favorably to values obtained from a dynamic Normal Mode Analysis (NMA). Higher correlation factors can be achieved by increasing the distance cutoff in the potential energy model. Furthermore, we apply our new method to analyze the entropy changes of the SARS CoV-2 main protease when binding with different ligand inhibitors, which is relevant for the design of potential drugs.

Kinematic Flexibility Analysis: Hydrogen Bonding Patterns Impart a Spatial Hierarchy of Protein Motion

Elastic network models (ENMs) and constraint-based, topological rigidity analysis are two distinct, coarse-grained approaches to study conformational flexibility of macromolecules. In the two decades since their introduction, both have contributed significantly to insights into protein molecular mechanisms and function. However, despite a shared purpose of these approaches, the topological nature of rigidity analysis, and thereby the absence of motion modes, has impeded a direct comparison. Here, we present an alternative, kinematic approach to rigidity analysis, which circumvents these drawbacks. We introduce a novel protein hydrogen bond network spectral decomposition, which provides an orthonormal basis for collective motions modulated by noncovalent interactions, analogous to the eigenspectrum of normal modes. The zero modes decompose proteins into rigid clusters identical to those from topological rigidity, while nonzero modes rank protein motions by their hydrogen bond collective energy penalty. Our kinematic flexibility analysis bridges topological rigidity theory and ENM, enabling a detailed analysis of motion modes obtained from both approaches. Analysis of a large, structurally diverse data set revealed that collectivity of protein motions, reported by the Shannon entropy, is significantly reduced for rigidity theory compared to normal mode approaches. Strikingly, kinematic flexibility analysis suggests that the hydrogen bonding network encodes a protein-fold specific, spatial hierarchy of motions, which goes nearly undetected in ENM. This hierarchy reveals distinct motion regimes that rationalize experimental and simulated protein stiffness variations. Kinematic motion modes highly correlate with reported crystallographic B factors and molecular dynamics simulations of adenylate kinase. A formal expression for changes in free energy derived from the spectral decomposition indicates that motions across nearly 40% of modes obey enthalpy-entropy compensation. Taken together, our results suggest that hydrogen bond networks have evolved to modulate protein structure and dynamics, which can be efficiently probed by kinematic flexibility analysis.