Explore the vast range of titles available.

Background & objectives: The present study proposed a series of computational techniques such as homology modelling, molecular simulation, and molecular docking to be performed to explore the structural features and binding mechanism of Cytochrome c oxidase subunit I (COX1) protein with known inhibitors. Methods: Elucidation of the three-dimensional structure of COX1 protein was carried out by using MODELLER software. The modelled protein was validated using GROMACS, structural qualitative tools and web servers. Finally the model was docked with carbon monoxide (CO) and nitric oxide (NO) using Auto Dock Tools. Results: The three-dimensional structure of mitochondrial transmembrane protein COX1 was built using homology modelling based on high-resolution crystal structures of Bos taurus. Followed by inserting the lipid bilayer, molecular dynamics simulation was performed on the modelled protein structure. The modelled protein was validated using qualitative structural indices. Known inhibitors such as carbon monoxide (CO) and nitric oxide (NO) inhibit their active binding sites of mitochondrial COX1 and the inhibitors were docked into the active site of attained model. A structure-based virtual screening was performed on the basis of the active site inhibition with best scoring hits. The COX1 model was submitted and can be accessible from the Model Archive site through the following link https://www.modelarchive.org/doi/10.5452/ma-at44v. Interpretation & conclusion: Structural characterization and active site identification can be further used as target for the planning of potent mosquitocidal compounds, thereby assisting the information in the field of research.

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.

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 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.

Cholesterol is an integral component of eukaryotic cell membranes and a key molecule in controlling membrane fluidity, organization, and other physicochemical parameters. It also plays a regulatory function in antibiotic drug resistance and the immune response of cells against viruses, by stabilizing the membrane against structural damage. While it iswell understood that, structurally, cholesterol exhibits a densification effect on fluid lipid membranes, its effects on membrane bending rigidity are assumed to be nonuniversal; i.e., cholesterol stiffens saturated lipid membranes, but has no stiffening effect on membranes populated by unsaturated lipids, such as 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). This observation presents a clear challenge to structure–property relationships and to our understanding of cholesterol-mediated biological functions. Here, using a comprehensive approach—combining neutron spin-echo (NSE) spectroscopy, solid-state deuterium NMR (²H NMR) spectroscopy, and molecular dynamics (MD) simulations—we report that cholesterol locally increases the bending rigidity of DOPC membranes, similar to saturated membranes, by increasing the bilayer’s packing density. All three techniques, inherently sensitive to mesoscale bending fluctuations, show up to a threefold increase in effective bending rigidity with increasing cholesterol content approaching a mole fraction of 50%. Our observations are in good agreement with the known effects of cholesterol on the area-compressibility modulus and membrane structure, reaffirming membrane structure–property relationships. The current findings point to a scale-dependent manifestation of membrane properties, highlighting the need to reassess cholesterol’s role in controlling membrane bending rigidity over mesoscopic length and time scales of important biological functions, such as viral budding and lipid–protein interactions.

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.

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

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.

Despite the biological importance of protein–protein complexes, determining their structures and association mechanisms remains an outstanding challenge. Here, we report the results of atomic-level simulations in which we observed five protein–protein pairs repeatedly associate to, and dissociate from, their experimentally determined native complexes using a molecular dynamics (MD)–based sampling approach that does not make use of any prior structural information about the complexes. To study association mechanisms, we performed additional, conventional MD simulations, in which we observed numerous spontaneous association events. A shared feature of native association for these five structurally and functionally diverse protein systems was that if the proteins made contact far from the native interface, the native state was reached by dissociation and eventual reassociation near the native interface, rather than by extensive interfacial exploration while the proteins remained in contact. At the transition state (the conformational ensemble from which association to the native complex and dissociation are equally likely), the protein–protein interfaces were still highly hydrated, and no more than 20% of native contacts had formed.