A Network of Conformational Transitions in the Apo Form of NDM-1 Enzyme Revealed by MD Simulation and a Markov State Model

Molecular insights into the inhibition mechanism of harringtonine against essential proteins associated with SARS-CoV-2 entry

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has recently posed a serious threat to global public health. Harringtonine (HT), as a small-molecule antagonist, has antiviral activity against a variety of viruses. There is evidence that HT can inhibit the SARS-CoV-2 entry into host cells by blocking the Spike protein and transmembrane protease serine 2 (TMPRSS2). However, the molecular mechanism underlying the inhibition effect of HT is largely elusive. Here, docking and all-atom molecular dynamic simulations were used to investigate the mechanism of HT against the receptor binding domain (RBD) of Spike, TMPRSS2, as well as the complex of RBD and angiotensin-converting enzyme 2 complex (RBD-ACE2). The results reveal that HT binds to all proteins primarily through hydrogen bond and hydrophobic interactions. Binding with HT influences the structural stability and dynamic motility processes of each protein. The interactions of HT with residues N33, H34 and K353 of ACE2, and residue K417 and Y453 of RBD contribute to disrupting the binding affinity between RBD and ACE2, which may hinder the virus entry into host cells. Our research provides molecular insights into the inhibition mechanism of HT against SARS-CoV-2 associated proteins, which will help for the novel antiviral drugs development.

Markov State Models Underlying the N-Terminal Premodel of TOPK/PBK

Lymphokine-activated killer T-cell-originated protein kinase (TOPK) is a potential target for cancer therapy. To explore the micromechanism, we proposed the N-terminal premodel (NTPM) of the TOPK monomer via homology modeling and molecular dynamic simulations and analyzed the conformational dynamics by Markov state model analysis. The electronegative insert (ENI) motif of the NTPM can be opened with a small probability under wild type, regulated by the so-called "N-C" interaction zone consisting of the N-terminal head, the coil between β₃-strand and αC-helix, and the ENI motif. Glutamate substitution at threonine residue 9 or tyrosine residue 74 promotes the closed-open transition, revealing the details of phosphorylation. Allosteric effects induce functionally relevant structural changes, such as increased structural flexibility and active sites, which are thought to be necessary for further activation or binding. These findings provide rational structural templates for designing state-dependent inhibitors and give insight into the molecular regulatory mechanisms of TOPK monomers.

Recognition between CD147 and cyclophilin A deciphered by accelerated molecular dynamics simulations

CD147 functions as the receptor of extracellular cyclophilin A (CypA) in various diseases, and CD147-CypA binding ulteriorly underlies the pathological process of various viral infections including HIV-1, SARS, and SARS-CoV-2. Although CyPA has been identified as a key intermediate pro-inflammatory factor, the mechanism by which CD147 cooperates with CypA in the development of the cytokine storm remains largely unknown, and the binding profile of CD147 with CypA remains to be elucidated as well. Here, we prepared three binding models of the CD147-CypA complex, including the active site of CypA severally binding to the groove bound by the Ig1 and Ig2 domains (model-0), P180-G181 (model-1), and P211 (model-2) of CD147, as well as introducing mutations P180A-G181A and P211A individually in each model. All systems were studied using accelerated molecular dynamics simulations and the molecular mechanics generalized Born surface area (MM/GBSA) method. For model-0, CypA bound to the ectodomain of CD147 with the highest binding affinity. Moreover, mutations P180A-G181A of CD147 in model-0 decreased the binding affinity and weakened the dynamic correlation between CD147 and CypA, which resulted in CypA shifting from the initial binding location. Other residue mutations of CD147 did not significantly affect the CD147-CypA binding, as reflected by the energy and structural analyses. Compared with surface plasmon resonance results and nuclear magnetic resonance shift signals, CypA should tend to reciprocally bind to the groove of CD147, and the binding process might be modulated by P180-G181 rather than P211. Besides, residue R201 of CD147 is critical for CD147-CypA binding and needs further experimental verification. These findings further our understanding of the recruitment between CD147 and CypA and its potential role in the development of inflammation and viral infection.

TAB1 binding induced p38α conformation change: an accelerated molecular dynamics simulation study

p38α mitogen-activated protein kinase (MAPK) undergoes autophosphorylation induced by the binding of TGFβ-activated kinase 1 binding protein 1 (TAB1) in myocardial ischemia. Investigation of the conformational transformations in p38α triggered by TAB1 binding is motivated by the need to find selective p38α activation inhibitors to treat myocardial ischemia. Herein, the conformational transformations of p38α were studied via all-atom accelerated molecular dynamics simulations and principal component analysis. With the binding of TAB1, the conformational changes of p38α auto-activation were characterized by the movement of the activation loop (A-loop) away from the αG helix toward the αF, αE helixes and L16-loop. In addition, a diverse intermediate state with an extensional and phosphorylated A-loop different from the transition intermediate state was explored. The conformational changes, including the A-loop alpha-structure breaking and the stronger hydrogen bond network formation, are accompanied by the extension of the A-loop and more intramolecular interactions in p38α. TAB1 correlates with other regions of p38α that are distal from the TAB1-binding site, including the A-loop, αC helix, and L16-loop, which regulates the intramolecular correlation of p38α. And, the phosphorylation further enhances the correlations between the A-loop and the other regions of p38α. The correlation results imply the regulation process of p38α conformational transformations. These findings will improve our understanding of the autophosphorylation of kinase and facilitate the development of selective inhibitors for the treatment of ischemic injury.

Hierarchical Markov State Model Building to Describe Molecular Processes

Markov state models can describe ensembles of pathways via kinetic networks but are difficult to create when large free-energy barriers limit unbiased sampling. Chain-of-states simulations allow sampling over large free-energy barriers but are often constructed using a single pathway that is unlikely to thermodynamically average over orthogonal degrees of freedom in complex systems. Here, we combine the advantages of these two approaches in the form of a Markov state model of Markov state models, which we call a Hierarchical Markov state model. In this approach, independent Markov models are constructed in regions of configuration space that are locally well sampled but are separated by large free-energy barriers from other regions. A string method is used to construct an ensemble of pathways connecting the states of these different local Markov models, and the rate through each pathway is then estimated. These rates are then combined with the rate information from the local Markov models in a master equation to predict global rates, fluxes, and populations. By applying this hierarchical approach to tractable systems, a toy potential and dipeptides, we demonstrate that it is more accurate than the conventional single-pathway description. The advantages of this approach are that it (i) is more realistic than the conventional chain-of-states approach, as an ensemble of pathways rather than a single pathway is used to describe processes in high-dimensional systems, and (ii) it resolves the issue of poor sampling in Markov State model building when large free-energy barriers are present. The divide-and-conquer strategy inherent to this approach should make this procedure straightforward to apply to more complex systems.

Generative network complex (GNC) for drug discovery

It remains a challenging task to generate a vast variety of novel compounds with desirable pharmacological properties. In this work, a generative network complex (GNC) is proposed as a new platform for designing novel compounds, predicting their physical and chemical properties, and selecting potential drug candidates that fulfill various druggable criteria such as binding affinity, solubility, partition coefficient, etc. We combine a SMILES string generator, which consists of an encoder, a drug-property controlled or regulated latent space, and a decoder, with verification deep neural networks, a target-specific three-dimensional (3D) pose generator, and mathematical deep learning networks to generate new compounds, predict their drug properties, construct 3D poses associated with target proteins, and reevaluate druggability, respectively. New compounds were generated in the latent space by either randomized output, controlled output, or optimized output. In our demonstration, 2.08 million and 2.8 million novel compounds are generated respectively for Cathepsin S and BACE targets. These new compounds are very different from the seeds and cover a larger chemical space. For potentially active compounds, their 3D poses are generated using a state-of-the-art method. The resulting 3D complexes are further evaluated for druggability by a championing deep learning algorithm based on algebraic topology, differential geometry, and algebraic graph theories. Performed on supercomputers, the whole process took less than one week. Therefore, our GNC is an efficient new paradigm for discovering new drug candidates.

Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Carbonic Anhydrase II Apo-Enzyme

Human carbonic anhydrase II (HCA II) is an enzyme that catalyzes the reversible hydration of CO2 into bicarbonate (HCO3-) and a proton (H+) as well as other reactions at an extremely high rate. This enzyme plays fundamental roles in human physiology/pathology, such as controlling the pH level in cells and so on. However, the binding mechanism between apo-HCA II and CO2 or other ligands as well as related conformational changes remains poorly understood, and atomic investigation into it could promote our understanding of related internal physiological/pathological mechanisms. In this study, long-time atomic molecular dynamics simulations as well as the clustering and free-energy analysis were performed to reveal the dynamics of apo-HCA II as well as the mechanism upon ligand binding. Our simulations indicate that the crystallographic B-factors considerably underestimate the loop dynamics: multiple conformations can be adopted by loops 1 and 2, especially for loop 1 because loop 1 is one side of the binding pocket, and its left-to-right movement can compress or extend the binding pocket, leading to one inactive (closed) state, three intermediate (semiopen) states, and one active (open) state; CO2 cannot get into the binding pocket of the inactive state but can get into those of intermediate and active states. The coexistence of multiple conformational states proposes a possible conformational selection model for the binding mechanism between apo-HCA II and CO2 or other ligands, revising our previous view of its functional mechanism of conformational change upon ligand binding and offering valuable structural insights into the workings of HCA II.