Molecular dynamics-guided optimization of BGM0504 enhances dual-target agonism for combating diabetes and obesity.

The dual activation of glucagon-like peptide-1 receptor (GLP-1R) and glucose-dependent insulinotropic polypeptide receptor (GIPR) has emerged as a promising therapeutic strategy for managing type 2 diabetes and obesity. Tirzepatide, a dual agonist peptide, has exhibited superior clinical efficacy in glycemic and weight control compared to selective GLP-1R agonists. Nevertheless, the structural basis of Tirzepatide's extended half-life, attributed to an acylation side chain on the parent peptide, raises questions regarding its partial agonistic activity. Employing molecular dynamics simulations, we explored the dynamic processes of peptide-receptor interactions. We uncovered a crucial salt bridge between parent peptide and GLP-1R/GIPR at K20, a feature not discernible in cryo-electron microscopy structures. Building upon these insights, we developed an optimization strategy based on the parent peptide which involved repositioning the acylation side chain. The results of both in vitro and in vivo experiments demonstrated that the optimized peptide has twofold to threefold increase in agonistic activity compared to Tirzepatide while maintaining its extended half-life in plasma. This led to the design of BGM0504, which proved to be more effective than its predecessor, Tirzepatide, in both laboratory and animal studies.

Analysis of the initial reaction mechanism of TKX-50 based on Raman intensity.

CONTEXT: Dihydroxylammonium 5,5'-biotetrazolium-1,1'-diolate (TKX-50) has two important properties of typical azole energy-containing ionic salts, including high energy and safety. Therefore, in today's era where more emphasis is placed on explosive performance and explosive detonation control conditions, TKX-50 is a very important object of research, and its reaction process in the initial stage of detonation is gradually receiving more and more attention from researchers in the field of energy-containing materials research. METHODS: In this paper, based on first-principles density-functional theory (DFT), the mechanism of chemical bond breakage of TKX-50 under pressure was determined based on the analysis of the strength and stability of chemical bonds inside the TKX-50 molecules using Raman spectroscopy relative intensity analysis. The results show that TKX-50 is dominated by N-H bond breaking and followed by H-O bond breaking in the initial reaction stage. These reactions lead to the reorganization and structural changes within the molecule, which eventually lead to the decomposition of TKX-50.

Molecular Insights on the Conformational Transitions and Activity Regulation of the c-Met Kinase Induced by Ligand Binding.

The tyrosine-protein kinase Met (c-Met) is an important signaling molecule involved in cellular growth and division. The dysregulation of c-Met may induce many fatal diseases, including non-small cell lung cancer, gastrointestinal cancers, hepatocellular carcinoma, etc. The activation of the c-Met kinase is dominant by the structure and dynamics of many important functional motifs, which are regulated by adenosine triphosphate (ATP) binding. c-Met inhibitors bind to the ATP-binding site or the allosteric pocket to compete with ATP molecules or alter the conformation of the function-related domains. Nevertheless, the mechanisms of ligand binding to c-Met are still unclear, especially the regulation of the functional motifs by different inhibitors. These greatly impede the development of novel drugs to overcome the drug tolerance to the currently marketed c-Met inhibitors. In this study, we used enhanced sampling technology to study the binding and regulation of two specific c-Met inhibitors. The results show that the two ligands adopt different binding processes even though with similar binding affinity. More importantly, our results uncovered different protein conformational features and the correlated motions of functional motifs regulated by the inhibitors, providing the structural basis for the functional suppression of the protein kinases.

Free Energy Calculations Using the Movable Type Method with Molecular Dynamics Driven Protein-Ligand Sampling.

Fast and accurate biomolecular free energy estimation has been a significant interest for decades, and with recent advances in computer hardware, interest in new method development in this field has even grown. Thorough configurational state sampling using molecular dynamics (MD) simulations has long been applied to the estimation of the free energy change corresponding to the receptor-ligand complexing process. However, performing large-scale simulation is still a computational burden for the high-throughput hit screening. Among molecular modeling tools, docking and scoring methods are widely used during the early stages of the drug discovery process in that they can rapidly generate discrete receptor-ligand binding modes and their individual binding affinities. Unfortunately, the lack of thorough conformational sampling in docking and scoring protocols leads to difficulty discovering global minimum binding modes on a complicated energy landscape. The Movable Type (MT) method is a novel absolute binding free energy approach which has demonstrated itself to be robust across a wide range of targets and ligands. Traditionally, the MT method is used with protein-ligand binding modes generated with rigid-receptor or flexible-receptor (induced fit) docking protocols; however, these protocols are by their nature less likely to be effective with more highly flexible targets or with those situations in which binding involves multiple step pathways. In these situations, more thorough samplings are required to better explain the free energy of binding. Therefore, to explore the prediction capability and computational efficiency of the MT method when using more thorough protein-ligand conformational sampling protocols, in the present work, we introduced a series of binding mode modeling protocols ranging from conventional docking routines to single-trajectory conventional molecular dynamics (cMD) and parallel Monte Carlo molecular dynamics (MCMD). Through validation against several structurally and mechanistically diverse protein-ligand test sets, we explore the performance of the MT method as a virtual screening tool to work with the docking protocols and as an MD simulation-based binding free energy tool.

The Conformational Transition Pathways and Hidden Intermediates in DFG-Flip Process of c-Met Kinase Revealed by Metadynamics Simulations.

Protein kinases intrinsically translate their conformations between active and inactive states, which is key to their enzymatic activities. The conformational flipping of the three-residue conservative motif, Asp-Phe-Gly (DFG), is crucial for many kinases' biological functions. Obtaining a detailed demonstration of the DFG flipping process and its corresponding dynamical and thermodynamical features could broaden our understanding of kinases' conformation-activity relationship. In this study, we employed metadynamics simulation, a widely used enhanced sampling technique, to analyze the conformational transition pathways of the DFG flipping for the c-Met kinase. The corresponding free energy landscape suggested two distinct transition pathways between the "DFG-in" and "DFG-out" states of the DFG-flip from c-Met. On the basis of the orientation direction of the F1223 residue, we correspondingly named the two pathways the "DFG-up" path, featuring forming a commonly discovered "DFG-up" transition state, and the "DFG-down" path, a unique transition pathway with F1223 rotating along the opposite direction away from the hydrophobic cavity. The free energies along the two pathways were then calculated using the Path Collective Variable (PCV) metadynamics simulation. The simulation results showed that, though having similar free energy barriers, the free energy cuve for the DFG-down path suggested a two-step conformational transition mechanism, while that for the DFG-up path showed the one-step transition feature. The c-Met DFG flipping mechanism and the new intermediate state discovered in this work could provide a deeper understanding of the conformation-activity relationship for c-Met and, possibly, reveal a new conformational state as the drug target for c-Met and other similar kinases.