Explore the vast range of titles available.

Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.

The PLUMED consortium unifies developers and contributors to PLUMED, an open-source library for enhanced-sampling, free-energy calculations and the analysis of molecular dynamics simulations. Here, we outline our efforts to promote transparency and reproducibility by disseminating protocols for enhanced-sampling molecular simulations.

Although various higher-order protein structure prediction methods have been developed, almost all of them were developed based on the three-dimensional (3D) structure information of known proteins. Here we predicted the short protein structures by molecular dynamics (MD) simulations in which only Newton’s equations of motion were used and 3D structural information of known proteins was not required. To evaluate the ability of MD simulationto predict protein structures, we calculated seven short test protein (10–46 residues) in the denatured state and compared their predicted and experimental structures. The predicted structure for Trp-cage (20 residues) was close to the experimental structure by 200-ns MD simulation. For proteins shorter or longer than Trp-cage, root-mean square deviation values were larger than those for Trp-cage. However, secondary structures could be reproduced by MD simulations for proteins with 10–34 residues. Simulations by replica exchange MD were performed, but the results were similar to those from normal MD simulations. These results suggest that normal MD simulations can roughly predict short protein structures and 200-ns simulations are frequently sufficient for estimating the secondary structures of protein (approximately 20 residues). Structural prediction method using only fundamental physical laws are useful for investigating non-natural proteins, such as primitive proteins and artificial proteins for peptide-based drug delivery systems.

NMDA-type glutamate receptors (NMDAR) are ligand-gated ion channels that contribute to excitatory neurotransmission in the central nervous system. NMDAR dysfunction has been found to be involved in various neurological disorders. Recent crystallographic and EM studies have shown the static structure of different states of the non-human NMDARs. Here we describe a model of a human NMDA receptor (hNMDAR) and its molecular dynamics (MD) before and after the binding of agonist ligands, glutamate and glycine. It is shown that the binding of ligands promotes a global reduction in molecular flexibility that produces a more tightly packed conformation than the unbound hNMDAR, and a higher cooperative regularity of moving. The ligand-induced synchronization of motion, identified on all structural levels of the modular hNMDA receptor is apparently a fundamental factor in channel gating. Although the time scale of the MD simulations (300 ns) was not sufficient to observe the complete gating event, the obtained data has shown the ligand-induced stabilization of hNMDAR that conforms the \"going to be open state\". We propose a mechanistic dynamic model of the ligand-dependent gating mechanism in the hNMDA receptor. At the binding of the ligands, the differently twisted conformations of the highly flexible receptor are stabilized in unique conformation with a linear molecular axis, which is a condition that is optimal for pore development. By searching the receptor surface, we have identified three new pockets, which are different from the pockets described in the literature as the potential and known positive allosteric modulator binding sites. A successful docking of two NMDAR modulators to their binding sites validates the model of a human NMDA receptor as a biological relevant target.

The belief that we could always stay ahead of the pathogens was forced upon scientists in the whole world by antimicrobial resistance. According to several reports, there are medications that are yet to be made public in the pipeline and there are little motivations to design novel antimicrobials to combat the worldwide drug resistance issues. Presently, the desire to design and develop efficient novel anti-bacterial agents is very high by researchers; thus, this study focuses on identifying the interactions between the studied ligands and Proplasmepsin IV, as well as examining the relationship between the calculated descriptors and binding affinities. This work shows successful prediction of the reacting and inhibiting efficiency of ten (10) cyclic tetra-peptides using insilico method. The optimization of the studied compound revealed the proficiency of methyl (3S,9S,12S)-12-(1,3-dioxoisoindolin-2-yl)-9-(2-(methylthio)ethyl)-5,8,11-trioxo-4,7,10-triaza-1(1,3)-benzenacyclotridecaphane-3-carboxylate ( F5 ) and 2-((3S,9S,12S)-12-(1,3-dioxoisoindolin-2-yl)-3-(methoxycarbonyl)-5,8,11-trioxo-4,7,10-triaza-1(1,3)-benzenacyclotridecaphane-9-yl)acetic acid (F7) to react more than the remaining molecules in term of HOMO and LUMO energies. In comparison, compound F9 demonstrated a higher inhibitory activity than the reference drug, Chloroquine, based on binding affinity. Molecular dynamics simulations over a 100 ns period further explored the binding affinity between F9 and the reference drug. The results showed that the reference drug (− 21.91 ± 1.16 kcal/mol) had a slightly stronger binding affinity than the F9_complex (− 13.85 ± 0.72 kcal/mol). Additionally, pharmacokinetic studies for F9 were compared with those of the reference compound and presented accordingly.

Proper balance between protein-protein and protein-water interactions is vital for atomistic molecular dynamics (MD) simulations of globular proteins as well as intrinsically disordered proteins (IDPs). The overestimation of protein-protein interactions tends to make IDPs more compact than those in experiments. Likewise, multiple proteins in crowded solutions are aggregated with each other too strongly. To optimize the balance, Lennard-Jones (LJ) interactions between protein and water are often increased about 10% (with a scaling parameter, λ = 1.1) from the existing force fields. Here, we explore the optimal scaling parameter of protein-water LJ interactions for CHARMM36m in conjunction with the modified TIP3P water model, by performing enhanced sampling MD simulations of several peptides in dilute solutions and conventional MD simulations of globular proteins in dilute and crowded solutions. In our simulations, 10% increase of protein-water LJ interaction for the CHARMM36m cannot maintain stability of a small helical peptide, (AAQAA)3 in a dilute solution and only a small modification of protein-water LJ interaction up to the 3% increase (λ = 1.03) is allowed. The modified protein-water interactions are applicable to other peptides and globular proteins in dilute solutions without changing thermodynamic properties from the original CHARMM36m. However, it has a great impact on the diffusive properties of proteins in crowded solutions, avoiding the formation of too sticky protein-protein interactions.

The cocrystallization of high-energy explosives has attracted great interests since it can alleviate to a certain extent the power-safety contradiction. 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaaza-isowurtzitane (CL-20), one of the most powerful explosives, has attracted much attention for researchers worldwide. However, the disadvantage of CL-20 has increased sensitivity to mechanical stimuli and cocrystallization of CL-20 with other compounds may provide a way to decrease its sensitivity. The intermolecular interaction of five types of CL-20-based cocrystal (CL-20/TNT, CL-20/HMX, CL-20/FOX-7, CL-20/TKX-50 and CL-20/DNB) by using molecular dynamic simulation was reviewed. The preparation methods and thermal decomposition properties of CL-20-based cocrystal are emphatically analyzed. Special emphasis is focused on the improved mechanical performances of CL-20-based cocrystal, which are compared with those of CL-20. The existing problems and challenges for the future work on CL-20-based cocrystal are discussed.

The prototypical model for NOD-like receptor (NLR) inflammasome assembly includes nucleotide-dependent activation of the NLR downstream of pathogen- or danger-associated molecular pattern (PAMP or DAMP) recognition, followed by nucleation of hetero-oligomeric platforms that lie upstream of inflammatory responses associated with innate immunity. As members of the STAND ATPases, the NLRs are generally thought to share a similar model of ATP-dependent activation and effect. However, recent observations have challenged this paradigm to reveal novel and complex biochemical processes to discern NLRs from other STAND proteins. In this review, we highlight past findings that identify the regulatory importance of conserved ATP-binding and hydrolysis motifs within the nucleotide-binding NACHT domain of NLRs and explore recent breakthroughs that generate connections between NLR protein structure and function. Indeed, newly deposited NLR structures for NLRC4 and NLRP3 have provided unique perspectives on the ATP-dependency of inflammasome activation. Novel molecular dynamic simulations of NLRP3 examined the active site of ADP- and ATP-bound models. The findings support distinctions in nucleotide-binding domain topology with occupancy of ATP or ADP that are in turn disseminated on to the global protein structure. Ultimately, studies continue to reveal how the ATP-binding and hydrolysis properties of NACHT domains in different NLRs integrate with signaling modules and binding partners to control innate immune responses at the molecular level.

A library of modified VEGFR-2 inhibitors was designed as VEGFR-2 inhibitors. Virtual screening was conducted for the hypothetical library using in silico docking, ADMET, and toxicity studies. Four compounds exhibited high in silico affinity against VEGFR-2 and an acceptable range of the drug-likeness. These compounds were synthesised and subjected to in vitro cytotoxicity assay against two cancer cell lines besides VEGFR-2 inhibitory determination. Compound D-1 showed cytotoxic activity against HCT-116 cells almost double that of sorafenib. Compounds A-1, C-6, and D-1 showed good IC 50 values against VEGFR-2. Compound D-1 markedly increased the levels of caspase-8 and BAX expression and decreased the anti-apoptotic Bcl-2 level. Additionally, compound D-1 caused cell cycle arrest at pre-G1 and G2-M phases in HCT-116 cells and induced apoptosis at both early and late apoptotic stages. Compound D-1 decreased the level of TNF-α and IL6 and inhibited TNF-α and IL6. MD simulations studies were performed over 100 ns.

Human α-defensin 5 (HD5) is one of cationic antimicrobial peptides which plays a crucial role in an innate immune system in human body. HD5 shows the killing activity against a broad spectrum of pathogenic bacteria by making a pore in a bacterial membrane and penetrating into a cytosol. Nonetheless, its pore-forming mechanisms remain unclear. Thus, in this work, the constant-velocity steered molecular dynamics (SMD) simulation was used to simulate the permeation of a dimeric HD5 into a gram-negative lipopolysaccharide (LPS) membrane model. Arginine-rich HD5 is found to strongly interact with a LPS surface. Upon arrival, arginines on HD5 interact with lipid A head groups (a top part of LPS) and then drag these charged moieties down into a hydrophobic core resulting in the formation of water-filled pore. Although all arginines are found to interact with a membrane, Arg13 and Arg32 appear to play a dominant role in the HD5 adsorption on a gram-negative membrane. Furthermore, one chain of a dimeric HD5 is required for HD5 adhesion. The interactions of arginine-lipid A head groups play a major role in adhering a cationic HD5 on a membrane surface and retarding a HD5 passage in the meantime. Graphical abstract