Drug binding dynamics of the dimeric SARS-CoV-2 main protease, determined by molecular dynamics simulation.

We performed molecular dynamics simulation of the dimeric SARS-CoV-2 (severe acute respiratory syndrome corona virus 2) main protease (Mpro) to examine the binding dynamics of small molecular ligands. Seven HIV inhibitors, darunavir, indinavir, lopinavir, nelfinavir, ritonavir, saquinavir, and tipranavir, were used as the potential lead drugs to investigate access to the drug binding sites in Mpro. The frequently accessed sites on Mpro were classified based on contacts between the ligands and the protein, and the differences in site distributions of the encounter complex were observed among the ligands. All seven ligands showed binding to the active site at least twice in 28 simulations of 200 ns each. We further investigated the variations in the complex structure of the active site with the ligands, using microsecond order simulations. Results revealed a wide variation in the shapes of the binding sites and binding poses of the ligands. Additionally, the C-terminal region of the other chain often interacted with the ligands and the active site. Collectively, these findings indicate the importance of dynamic sampling of protein-ligand complexes and suggest the possibilities of further drug optimisations.

Temperature-Sensitive Substrate and Product Binding Underlie Temperature-Compensated Phosphorylation in the Clock.

Temperature compensation is a striking feature of the circadian clock. Here we investigate biochemical mechanisms underlying temperature-compensated, CKIδ-dependent multi-site phosphorylation in mammals. We identify two mechanisms for temperature-insensitive phosphorylation at higher temperature: lower substrate affinity to CKIδ-ATP complex and higher product affinity to CKIδ-ADP complex. Inhibitor screening of ADP-dependent phosphatase activity of CKIδ identified aurintricarboxylic acid (ATA) as a temperature-sensitive kinase activator. Docking simulation of ATA and mutagenesis experiment revealed K224D/K224E mutations in CKIδ that impaired product binding and temperature-compensated primed phosphorylation. Importantly, K224D mutation shortens behavioral circadian rhythms and changes the temperature dependency of SCN's circadian period. Interestingly, temperature-compensated phosphorylation was evolutionary conserved in yeast. Molecular dynamics simulation and X-ray crystallography demonstrate that an evolutionally conserved CKI-specific domain around K224 can provide a structural basis for temperature-sensitive substrate and product binding. Surprisingly, this domain can confer temperature compensation on a temperature-sensitive TTBK1. These findings suggest the temperature-sensitive substrate- and product-binding mechanisms underlie temperature compensation.

Perturbation analyses of intermolecular interactions.

Conformational fluctuations of a protein molecule are important to its function, and it is known that environmental molecules, such as water molecules, ions, and ligand molecules, significantly affect the function by changing the conformational fluctuations. However, it is difficult to systematically understand the role of environmental molecules because intermolecular interactions related to the conformational fluctuations are complicated. To identify important intermolecular interactions with regard to the conformational fluctuations, we develop herein (i) distance-independent and (ii) distance-dependent perturbation analyses of the intermolecular interactions. We show that these perturbation analyses can be realized by performing (i) a principal component analysis using conditional expectations of truncated and shifted intermolecular potential energy terms and (ii) a functional principal component analysis using products of intermolecular forces and conditional cumulative densities. We refer to these analyses as intermolecular perturbation analysis (IPA) and distance-dependent intermolecular perturbation analysis (DIPA), respectively. For comparison of the IPA and the DIPA, we apply them to the alanine dipeptide isomerization in explicit water. Although the first IPA principal components discriminate two states (the α state and PPII (polyproline II) + ß states) for larger cutoff length, the separation between the PPII state and the ß state is unclear in the second IPA principal components. On the other hand, in the large cutoff value, DIPA eigenvalues converge faster than that for IPA and the top two DIPA principal components clearly identify the three states. By using the DIPA biplot, the contributions of the dipeptide-water interactions to each state are analyzed systematically. Since the DIPA improves the state identification and the convergence rate with retaining distance information, we conclude that the DIPA is a more practical method compared with the IPA. To test the feasibility of the DIPA for larger molecules, we apply the DIPA to the ten-residue chignolin folding in explicit water. The top three principal components identify the four states (native state, two misfolded states, and unfolded state) and their corresponding eigenfunctions identify important chignolin-water interactions to each state. Thus, the DIPA provides the practical method to identify conformational states and their corresponding important intermolecular interactions with distance information.

Perturbational formulation of principal component analysis in molecular dynamics simulation.

Conformational fluctuations of a molecule are important to its function since such intrinsic fluctuations enable the molecule to respond to the external environmental perturbations. For extracting large conformational fluctuations, which predict the primary conformational change by the perturbation, principal component analysis (PCA) has been used in molecular dynamics simulations. However, several versions of PCA, such as Cartesian coordinate PCA and dihedral angle PCA (dPCA), are limited to use with molecules with a single dominant state or proteins where the dihedral angle represents an important internal coordinate. Other PCAs with general applicability, such as the PCA using pairwise atomic distances, do not represent the physical meaning clearly. Therefore, a formulation that provides general applicability and clearly represents the physical meaning is yet to be developed. For developing such a formulation, we consider the conformational distribution change by the perturbation with arbitrary linearly independent perturbation functions. Within the second order approximation of the Kullback-Leibler divergence by the perturbation, the PCA can be naturally interpreted as a method for (1) decomposing a given perturbation into perturbations that independently contribute to the conformational distribution change or (2) successively finding the perturbation that induces the largest conformational distribution change. In this perturbational formulation of PCA, (i) the eigenvalue measures the Kullback-Leibler divergence from the unperturbed to perturbed distributions, (ii) the eigenvector identifies the combination of the perturbation functions, and (iii) the principal component determines the probability change induced by the perturbation. Based on this formulation, we propose a PCA using potential energy terms, and we designate it as potential energy PCA (PEPCA). The PEPCA provides both general applicability and clear physical meaning. For demonstrating its power, we apply the PEPCA to an alanine dipeptide molecule in vacuum as a minimal model of a nonsingle dominant conformational biomolecule. The first and second principal components clearly characterize two stable states and the transition state between them. Positive and negative components with larger absolute values of the first and second eigenvectors identify the electrostatic interactions, which stabilize or destabilize each stable state and the transition state. Our result therefore indicates that PCA can be applied, by carefully selecting the perturbation functions, not only to identify the molecular conformational fluctuation but also to predict the conformational distribution change by the perturbation beyond the limitation of the previous methods.