Objectives: This study aims to assess the drug-likeness and binding of nucleobase-substituted ponatinib analogues towards wild-type and T315I mutant BCR-ABL tyrosine kinases. Methodology: A total of 415 ponatinib analogues, encompassing single and combinatorial modifications on five parts of the drug were generated, profiled in SwissADME, and subjected to molecular docking using AutoDock4. Complexes formed by the top analogues then underwent a 100-ns molecular dynamics simulation with GROMACS. Results: Analogues featuring the replacement of the imidazo[1,2b]pyridazine with adenine and cytosine exhibited promising binding free energies, attributed to the presence of primary amines that facilitate crucial hydrogen bond interactions in the hinge region. RMSD, RMSF, and atomic distance analyses of the MD trajectories revealed that the six top analogues formed stable complexes in their inactive DFG-out conformations. Changes in the MMPBSA and MMGBSA-calculated free energies were mainly driven by changes in hydrogen bonds. Furthermore, drug-likeness predictions supported the formulation of most analogues for oral administration. Conclusion Among the top analogues, VP10004 and VP81014 exhibited the most favorable binding free energies and interactions with the target models, while VP10312 was identified as the most feasible candidate for synthesis.
Hydrogen Bonding ; Molecular Dynamics Simulation ; Molecular Docking Simulation

The binding of talin-F0 domain to ras-related protein 1b (Rap1b) plays an important role in the formation of thrombosis. However, since talin is a force-sensitive protein, it remains unclear whether and how force regulates the talin-F0/Rap1b interaction. To explore the effect of force on the binding affinity and the dynamics mechanisms of talin-F0/Rap1b, molecular dynamics simulation was used to observe and compare the changes in functional and conformational information of the complex under different forces. Our results showed that when the complex was subjected to tensile forces, there were at least two dissociation pathways with significantly different mechanical strengths. The key event determining the mechanical strength difference between the two pathways was whether the β4 sheet of the F0 domain was pulled away from the original β1-β4 parallel structure. As the force increased, the talin-F0/Rap1b interaction first strengthened and then weakened, exhibiting the signature of a transition from catch bonds to slip bonds. The mechanical load of 20 pN increased the interaction index of two residue pairs, ASP ⁵⁴-ARG ⁴¹ and GLN ¹⁸-THR ⁶⁵, which resulted in a significant increase in the affinity of the complex. This study predicts the regulatory mechanism of the talin-F0/Rap1b interaction by forces in the intracellular environment and provides novel ideas for the treatment of related diseases and drug development.
Molecular Dynamics Simulation ; Talin

As one of the traditional computer simulation techniques, molecular simulation can intuitively display and quantify molecular structure and explain experimental phenomena from the microscopic molecular level. When the simulation system increases, the amount of calculation will also increase, which will cause a great burden on the simulation system. Coarse-grained molecular dynamics is a method of mesoscopic molecular simulation, which can simplify the molecular structure and improve computational efficiency, as a result, coarse-grained molecular dynamics is often used when simulating macromolecular systems such as drug carrier materials. In this article, we reviewed the recent research results of using coarse-grained molecular dynamics to simulate drug carriers, in order to provide a reference for future pharmaceutical preparation research and accelerate the entry of drug research into the era of precision drug design.
Molecular Dynamics Simulation ; Drug Carriers

In resting platelets, the 17 ^th domain of filamin a (FLNa17) constitutively binds to the platelet membrane glycoprotein Ibα (GPIbα) at its cytoplasmic tail (GPIbα-CT) and inhibits the downstream signal activation, while the binding of ligand and blood shear force can activate platelets. To imitate the pull force transmitted from the extracellular ligand of GPIbα and the lateral tension from platelet cytoskeleton deformation, two pulling modes were applied on the GPIbα-CT/FLNa17 complex, and the molecular dynamics simulation method was used to explore the mechanical regulation on the affinity and mechanical stability of the complex. In this study, at first, nine pairs of key hydrogen bonds on the interface between GPIbα-CT and FLNa17 were identified, which was the basis for maintaining the complex structural stability. Secondly, it was found that these hydrogen bonding networks would be broken down and lead to the dissociation of FLNa17 from GPIbα-CT only under the axial pull force; but, under the lateral tension, the secondary structures at both terminals of FLNa17 would unfold to protect the interface of the GPIbα-CT/FLNa17 complex from mechanical damage. In the range of 0~40 pN, the increase of pull force promoted outward-rotation of the nitrogen atom of the 563 ^rd phenylalanine (PHE ⁵⁶³-N) at GPIbα-CT and the dissociation of the complex. This study for the first time revealed that the extracellular ligand-transmitted axial force could more effectively relieve the inhibition of FLNa17 on the downstream signal of GPIbα than pure mechanical tension at the atomic level, and would be useful for further understanding the platelet intracellular force-regulated signal pathway.
Filamins/metabolism* ; Platelet Glycoprotein GPIb-IX Complex/metabolism* ; Molecular Dynamics Simulation ; Ligands ; Protein Binding ; Blood Platelets/metabolism* ; von Willebrand Factor/metabolism*

This study aims to gain insight into the DNA-specific recognition mechanism of c-Myb transcription factor during the regulation of cell early differentiation and proliferation. Therefore, we chose the chicken myeloid gene, mitochondrial import protein 1 (mim-1), as a target to study the binding specificity between potential dual-Myb-binding sites. The c-Myb-binding site in mim-1 is a pseudo-palindromic sequence AACGGTT, which contains two AACNG consensuses. Simulation studies in different biological scenarios revealed that c-Myb binding with mim-1 in the forward strand (complex F) ismore stable than that inthereverse strand (complex R). The principal component analysis (PCA) dynamics trajectory analyses suggested an opening motion of the recognition helices of R2 and R3 (R2R3), resulting in the dissociation of DNA from c-Myb in complex R at 330 K, triggered by the reduced electrostatic potential on the surface of R2R3. Furthermore, the DNA confirmation and hydrogen-bond interaction analyses indicated that the major groove width of DNA increased in complex R, which affected on the hydrogen-bond formation ability between R2R3 and DNA, and directly resulted in the dissociation of DNA from R2R3. The steered molecular dynamics (SMD) simulation studies also suggested that the electrostatic potential, major groove width, and hydrogen bonds made major contribution to the DNA‍-specific recognition. In vitro trials confirmed the simulation results that c-Myb specifically bound to mim-1 in the forward strand. This study indicates that the three-dimensional (3D) structure features play an important role in the DNA-specific recognition mechanism by c-Myb besides the AACNG consensuses, which is beneficial to understanding the cell early differentiation and proliferation regulated by c-Myb, as well as the prediction of novel c-Myb-binding motifs in tumorigenesis.
Molecular Dynamics Simulation ; Consensus ; DNA ; Hydrogen

We aim to screen out the active components that may have therapeutic effect on coronavirus disease 2019 (COVID-19) from the severe and critical cases' prescriptions in the "Coronavirus Disease 2019 Diagnosis and Treatment Plan (Trial Ninth Edition)" issued by the National Health Commission of the People's Republic of China and explain its mechanism through the interactions with proteins. The ETCM database and SwissADME database were used to screen the active components contained in 25 traditional Chinese medicines in 3 prescriptions, and the PDB database was used to obtain the crystal structures of 4 proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Molecular docking was performed using Autodock Vina and molecular dynamics simulations were performed using GROMACS. Binding energy results showed that 44 active ingredients including xambioona, gancaonin L, cynaroside, and baicalin showed good binding affinity with multiple targets of SARS-CoV-2, while molecular dynamics simulations analysis showed that xambioona bound more tightly to the nucleocapsid protein of SARS-CoV-2 and exerted a potent inhibitory effect. Modern technical methods are used to study the active components of traditional Chinese medicine and show that xambioona is an effective inhibitor of SARS-CoV-2 nucleocapsid protein, which provides a theoretical basis for the development of new anti-SARS-CoV-2 drugs and their treatment methods.
Humans ; SARS-CoV-2 ; Molecular Docking Simulation ; Medicine, Chinese Traditional ; Molecular Dynamics Simulation ; Nucleocapsid Proteins ; Antiviral Agents/pharmacology* ; COVID-19 Drug Treatment

Antiviral Agents/pharmacology* ; COVID-19 ; Coronavirus 3C Proteases ; Humans ; Lactams ; Leucine ; Molecular Docking Simulation ; Molecular Dynamics Simulation ; Nitriles ; Proline ; Protease Inhibitors/pharmacology* ; SARS-CoV-2

Proteus mirabilis lipase (PML) features tolerance to organic solvents and great potential for biodiesel synthesis. However, the thermal stability of the enzyme needs to be improved before it can be used industrially. Various computational design strategies are emerging methods for the modification of enzyme thermal stability. In this paper, the complementary algorithm-based ABACUS, PROSS, and FoldX were employed for positive selection of PML mutations, and their pairwise intersections were further subjected to negative selection by PSSM and G_REMLIN to narrow the mutation library. Thereby, 18 potential single-point mutants were screened out. According to experimental verification, 7 mutants had melting temperature (T_m) improved, and the ΔT_m of K208G and G206D was the highest, which was 3.75 ℃ and 3.21 ℃, respectively. Five mutants with activity higher than the wild type (WT) were selected for combination by greedy accumulation. Finally, the T_m of the five-point combination mutant M10 increased by 10.63 ℃, and the relative activity was 140% that of the WT. K208G and G206D exhibited certain epistasis during the combination, which made a major contribution to the improvement of the thermal stability of M10. Molecular dynamics simulation indicated that new forces were generated at and around the mutation sites, and the rearrangement of forces near G206D/K208G might stabilize the Ca²⁺ binding site which played a key role in the stabilization of PML. This study provides an efficient and user-friendly computational design scheme for the thermal stability modification of natural enzymes and lays a foundation for the modification of PML and the expansion of its industrial applications.
Enzyme Stability ; Lipase/chemistry* ; Molecular Dynamics Simulation ; Proteus mirabilis/metabolism* ; Solvents/chemistry*

A new coronavirus (SARS-CoV-2) has been identified as the etiologic agent for the COVID-19 outbreak. Currently, effective treatment options remain very limited for this disease; therefore, there is an urgent need to identify new anti-COVID-19 agents. In this study, we screened over 6,000 compounds that included approved drugs, drug candidates in clinical trials, and pharmacologically active compounds to identify leads that target the SARS-CoV-2 papain-like protease (PLpro). Together with main protease (M
Antiviral Agents/therapeutic use* ; Binding Sites ; COVID-19/virology* ; Coronavirus Papain-Like Proteases/metabolism* ; Crystallography, X-Ray ; Drug Evaluation, Preclinical ; Drug Repositioning ; High-Throughput Screening Assays/methods* ; Humans ; Imidazoles/therapeutic use* ; Inhibitory Concentration 50 ; Molecular Dynamics Simulation ; Mutagenesis, Site-Directed ; Naphthoquinones/therapeutic use* ; Protease Inhibitors/therapeutic use* ; Protein Structure, Tertiary ; Recombinant Proteins/isolation & purification* ; SARS-CoV-2/isolation & purification*

Aims: The surplus use of herbicide Dalapon® contains 2,2-dichloropropionic acid (2,2-DCP) poses great danger to human and ecosystem due to its toxicity. Hence, this study focused on the isolation and characterization of a dehalogenase producing bacteria from Sungai Skudai, Johor, capable of utilizing 2,2-DCP as a carbon source and in silico analysis of its putative dehalogenase. Methodology and results: Isolation of the target bacteria was done by using 2,2-DCP-enriched culture as the sole carbon source that allows a bacterium to grow in 20 mM of 2,2-DCP at 30 °C with the corresponding doubling time of 8.89 ± 0.03 h. The isolated bacterium was then designated as Klebsiella pneumoniae strain YZ based on biochemical tests and basic morphological examination. The full genome of K. pneumoniae strain KLPN_25 (accession number: RRE04903) which obtained from NCBI database was screened for the presence of dehalogenase gene, assuming both strains YZ and KLPN_25 were the same organisms. A putative dehalogenase gene was then identified as type II dehalogenase from the genome sequence of strain KLPN_25. The protein structure of the type II dehalogenase of KLPN_25 strain was then pairwise aligned with the crystal structure of L-2-haloacid dehalogenase (L-DEX) Pseudomonas sp. strain YL as the template, revealing the existence of conserved amino acids residues, uniquely known to participate in the dehalogenation mechanism. The finding thus implies that the amino acid residues of type II dehalogenase possibly shares similar catalytic functions with the L-DEX. Conclusion, significance and impact of the study In conclusion, this study confirmed the presence of new dehalogenase from the genus Klebsiella with potential to degrade 2,2-DCP from the river water. The structural information of type II dehalogenase provides insights for future work in designing haloacid dehalogenases.
Klebsiella pneumoniae--isolation & ; purification ; Computer Simulation ; Molecular Dynamics Simulation