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

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

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

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

In the present research, molecular simulation and quantum chemistry calculations were combined to investigate the thermal stability of three kinds of insulin aggregations and the effect of Zn (II) ion coordination on these aggregations. The results of molecular simulation indicated that the three insulin dimers in the same sphere closed hexamer had synergistic stability. It is the synergistic stability that enhances the structural and thermal stability of insulin, preserves its bioactivity during production, storage, and delivery of insulin formulations, and prolongs its halflife in human bodies. According to the results of quantum chemistry calculations, each Zn (II)-N (Im-insulin) bond energy can reach 73.610 kJ/mol for insulin hexamer and 79.907 kJ/mol for insulin tetramer. However, the results of Gibbs free energy changes still indicats that the coordination of zinc (II) ions is unfavorable for the formation of insulin hexamer, because the standard Gibbs free energy change of the coordinate reaction of zinc (II) ions associated with the formatting insulin hexamer is positive and increased.
Insulin ; chemistry ; metabolism ; Molecular Dynamics Simulation ; Protein Stability ; Zinc ; chemistry

Glutamate decarboxylase, a unique pyridoxal 5'-phosphate-dependent enzyme, catalyzes α-decarboxylation of L-glutamate to γ-aminobutyrate. However, glutamate decarboxylase from different sources has the common problem of poor thermostability that affects its application in industry. In this study, proline was introduced at 13 different positions in glutamate decarboxylase by using the design strategy of homologous sequence alignment between Thermococcus kodakarensis and Lactobacillus brevis CGMCC No.1306. A mutant enzyme G364P with higher thermostability was obtained. Compared to the wild type, thermostability of the mutant G364P was significantly improved, the half-life time (t1/2) at 55 °C and the semi-inactivation temperature (T₅₀ ¹⁵) of the mutant G364P increased 19.4 min and 5.3 °C, respectively, while kcat/Km of the mutant enzyme remained nearly unchanged. Further analysis of their thermostability by molecular dynamics simulations were performed. The root mean square deviation of G364P and root mean square fluctuation in the loop region including G364 were lower than the wild type at 313 K for 10 ns, and G364P increased one hydrophobic interaction in the loop region. It proves that mutation of flexible 364-Gly to rigid proline endows glutamate decarboxylase with enhanced thermostability.
Glutamate Decarboxylase ; Glutamic Acid ; Lactobacillus brevis ; Molecular Dynamics Simulation ; Proline

More effective production of human insulin is important, because insulin is the main medication that is used to treat multiple types of diabetes and because many people are suffering from diabetes. The current system of insulin production is based on recombinant DNA technology, and the expression vector is composed of a preproinsulin sequence that is a fused form of an artificial leader peptide and the native proinsulin. It has been reported that the sequence of the leader peptide affects the production of insulin. To analyze how the leader peptide affects the maturation of insulin structurally, we adapted several in silico simulations using 13 artificial proinsulin sequences. Three-dimensional structures of models were predicted and compared. Although their sequences had few differences, the predicted structures were somewhat different. The structures were refined by molecular dynamics simulation, and the energy of each model was estimated. Then, protein-protein docking between the models and trypsin was carried out to compare how efficiently the protease could access the cleavage sites of the proinsulin models. The results showed some concordance with experimental results that have been reported; so, we expect our analysis will be used to predict the optimized sequence of artificial proinsulin for more effective production.
Computer Simulation ; DNA, Recombinant ; Humans* ; Insulin ; Molecular Dynamics Simulation* ; Proinsulin ; Protein Sorting Signals ; Trypsin

Today, the understanding of the sequence and structure of biologically relevant targets is growing rapidly and researchers from many disciplines, physics and computational science in particular, are making significant contributions to modern biology and drug discovery. However, it remains challenging to rationally design small molecular ligands with desired biological characteristics based on the structural information of the drug targets, which demands more accurate calculation of ligand binding free-energy. With the rapid advances in computer power and extensive efforts in algorithm development, physics-based computational chemistry approaches have played more important roles in structure-based drug design. Here we reviewed the newly developed computational chemistry methods in structure-based drug design as well as the elegant applications, including binding-site druggability assessment, large scale virtual screening of chemical database, and lead compound optimization. Importantly, here we address the current bottlenecks and propose practical solutions.
Computational Biology ; Drug Design ; Drug Discovery ; High-Throughput Screening Assays ; Molecular Docking Simulation ; Molecular Dynamics Simulation ; Quantitative Structure-Activity Relationship

Fatty acid synthase (FASN, EC 2.3.1.85), is a multi-enzyme dimer complex that plays a critical role in lipogenesis. This lipogenic enzyme has gained importance beyond its physiological role due to its implications in several clinical conditions-cancers, obesity, and diabetes. This has made FASN an attractive pharmacological target. Here, we have attempted to predict the theoretical models for the human enoyl reductase (ER) and beta-ketoacyl reductase (KR) domains based on the porcine FASN crystal structure, which was the structurally closest template available at the time of this study. Comparative modeling methods were used for studying the structure-function relationships. Different validation studies revealed the predicted structures to be highly plausible. The respective substrates of ER and KR domains-namely, trans-butenoyl and beta-ketobutyryl-were computationally docked into active sites using Glide in order to understand the probable binding mode. The molecular dynamics simulations of the apo and holo states of ER and KR showed stable backbone root mean square deviation trajectories with minimal deviation. Ramachandran plot analysis showed 96.0% of residues in the most favorable region for ER and 90.3% for the KR domain, respectively. Thus, the predicted models yielded significant insights into the substrate binding modes of the ER and KR catalytic domains and will aid in identifying novel chemical inhibitors of human FASN that target these domains.
Catalytic Domain* ; Humans ; Lipogenesis ; Models, Theoretical ; Molecular Dynamics Simulation* ; Obesity ; Oxidoreductases*