Search

Hybrid quantum mechanics/molecular mechanics (QM/MM) is one of the most reliable approaches for accurately modeling and studying the complex pharmaceutical discovery system. Classical mechanics has significantly accelerated the drug discovery process in the past decade. However, the current challenge is the large pool of false positives, which require extensive validation. Hybrid QM/MM is an effective solution for accurately studying ligand binding, structural mechanisms, free energy evaluation, and spectroscopic characterization. This article highlights the methodological details relevant to cost-effective hybrid QM/MM methods. This approach, combined with traditional pharmacoinformatics methods, could be a reliable strategy to balance the cost and accuracy of the calculations., Competing Interests: Declaration of Competing Interest The author have no conflicts of interest to declare., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

The formation of amyloid fibers has been implicated in a number of neurodegenerative diseases. The growth of amyloid fibers is strongly thermodynamically favorable, but kinetic traps exist where the incoming monomer binds in an incompatible conformation that blocks further elongation. Unfortunately, this process is difficult to follow experimentally at the atomic level. It is also too complex to simulate in full detail and to date has been explored either through coarse-grained simulations, which may miss many important interactions, or full atomic simulations, in which the incoming peptide is constrained to be near the ideal fiber geometry. Here we use an alternate approach starting from a docked complex in which the monomer is from an experimental NMR structure of one of the major conformations in the unbound ensemble, a largely unstructured peptide with the central hydrophobic region in a 3 10 helix. A 1000 ns full atomic simulation in explicit solvent shows the formation of a metastable intermediate by sequential, concerted movements of both the fiber and the monomer. A Markov state model shows that the unfolded monomer is trapped at the end of the fiber in a set of interconverting antiparallel β-hairpin conformations. The simulation here may serve as a model for the binding of other non-β-sheet conformations to amyloid fibers.

Multidrug resistance against the existing antibiotics is one of the most challenging threats across the globe. Antimicrobial peptides (AMPs), in this regard, are considered to be one of the effective alternatives that can overcome bacterial resistance. MSI-594, a 24-residue linear alpha-helical cationic AMP, has been shown to function via the carpet mechanism to disrupt bacterial membrane systems. To better understand the role of lipid composition in the function of MSI-594, in the present study, eight different model membrane systems have been studied using accelerated molecular dynamics (aMD) simulations. The simulated results are helpful in discriminating the particular effects of cationic MSI-594 against zwitterionic POPC, anionic POPG and POPS, and neutral POPE lipid moieties. Additionally, the effects of various heterogeneous POPC/POPG (7 : 3), POPC/POPS (7 : 3), and POPG/POPE (1 : 3 and 3 : 1) bilayer systems on the dynamic interaction of MSI-594 have also been investigated. The effect on the lipid bilayer due to the interaction with the peptide is characterized by lipid acyl-chain order, membrane thickness, and acyl-chain dynamics. Our simulation results show that the lipid composition affects the membrane interaction of MSI-594, suggesting that membrane selectivity is crucial to its mechanism of action. The results reported in this study are helpful to obtain accurate atomistic-level information governing MSI-594 and its membrane disruptive antimicrobial mechanism of action, and to design next generation potent antimicrobial peptides.

Antifreeze proteins (AFPs) are the key biomolecules that enable species to survive under subzero temperature conditions. The physiologically relevant activities of AFPs are based on the adsorption to ice crystals, followed by the inhibition of subsequent crystal layer growth of ice, routed with depression in freezing point in a noncolligative manner. The functional attributes governing the mechanism by which AFPs inhibit freezing of body fluids in bacteria, fungi, plants, and fishes are mainly attributed to their adsorption onto the surface of ice within the physiological system. Importantly, AFPs are also known for their application in cryopreservation of biological samples that might be related to membrane interaction. To date, there is a paucity of information detailing the interaction of AFPs with membrane structures. Here, we focus on elucidating the biophysical properties of the interactions between AFPs and micelle models that mimic the membrane system. Micelle model systems of zwitterionic DPC and negatively charged SDS were utilized in this study, against which a significant interaction is experienced by two AFP molecules, namely, Peptide 1m and wfAFP (the popular AFP sourced from winter flounder). Using low- and high-resolution biophysical characterization techniques, such as circular dichroism (CD) and NMR spectroscopy, a strong evidence for the interactions of these AFPs with the membrane models is revealed in detail and is corroborated by in-depth residue-specific information derived from molecular dynamics simulation. Altogether, these results not only strengthen the fact that AFPs interact actively with membrane systems, but also demonstrate that membrane-associated AFPs are dynamic and capable of adopting a number of conformations rendering fluidity to the system.

Antifreeze proteins (AFPs) enable the polar living species to survive subzero temperature conditions through effective lowering of the freezing point of body fluids. At the molecular level, AFPs directly interact with the growing seeds of ice crystals to inhibit their formation. To understand the structural and dynamic aspects of this interaction at the atomistic level, molecular dynamics (MD) simulations were carried out on several type I AFPs at multiple temperatures, including the physiologically relevant temperature of 273 K, a lower temperature of 227 K, and the conventional 300 K. A comparison of the principal component analysis (PCA) and mean squared deviation plots for Winter flounder AFP, HPLC6 (mutant of winter flounder AFP), Sculpin, and peptide 1m AFPs reveals that simulations at 273 and 227 K result in the formation of more conserved metastable states than at 300 K. Other parameters such as root-mean-square deviation (rmsd), solvent accessibility surface area (SASA), H-bonding and residual density function (RDF) also suggest the same. MD simulations with ice crystal, where AFPs are complexed to ice plane with TIP4P/ice water model, help in finding relevance of dynamic behavior, and physiological temperature becomes more pronounced. Additionally, a control study on a nonantifreeze protein (LL37) is included, which aids in exploring significant information. On the basis of this approach, it was found that AFPs at 273 and 227 K display relevant dynamic properties that appear at 300 K for nonantifreeze proteins. The present study hence emphasizes the importance of performing computational simulations for antifreeze proteins at the physiologically relevant temperature (273 K), and even at lower temperatures (like 227 K), rather than at room temperatures (300 K).

Enzyme adenosine kinase is responsible for phosphorylation of adenosine to AMP and is crucial for parasites which are purine auxotrophs. The present study describes development of robust homology model of Leishmania donovani adenosine kinase to forecast interaction phenomenon with inhibitory molecules using structure-based drug designing strategy. Docking calculation using reported organic small molecules and natural products revealed key active site residues such as Arg131 and Asp16 for ligand binding, which is consistent with previous studies. Molecular dynamics simulation of ligand protein complex revealed the importance of hydrogen bonding with active site residues and solvent molecules, which may be crucial for successful development of drug candidates. Precise role of Phe168 residue in the active site was elucidated in this report that provided stability to ligand-protein complex via aromatic- π contacts. Overall, the present study is believed to provide valuable information to design a new compound with improved activity for antileishmanial therapeutics development.

Approximately 27 million people are suffering from cancer that contains either an inactivating missense mutation of TP53 gene or partially abrogated p53 signaling pathway. Concerted action of folded and intrinsically disordered domains accounts for multi-faceted role of p53. The intricacy of dynamic p53 structure is believed to shed light on its cellular activity for developing new cancer therapies. In this review, insights into structural details of p53, diverse single point mutations affecting its core domain, thermodynamic understanding and therapeutic strategies for pharmacological rescue of p53 function has been illustrated. An effort has been made here to bridge the structural and sequential evidence of p53 from experimental to computational studies. First, we focused on the individual domains and the crucial protein–protein or DNA–protein contacts that determine conformation and dynamic behavior of p53. Next, the oncogenic mutations associated with cancer and its contribution to thermodynamic fluctuation has been discussed. Thus the emerging anti-cancer strategies include targeting of destabilized cancer mutants with selective inhibition of its negative regulators. Recent advances in development of small molecule inhibitors and peptides exploiting p53–MDM2 interaction has been included. In a nutshell, this review attempts to describe structural biology of p53 which provide new openings for structure-guided rescue. [ABSTRACT FROM AUTHOR]

Microtubule affinity regulating kinase (MARK4) is considered as a potential drug target for diabetes, cancer, and neurodegenerative diseases. Since the role of MARK4 in the phosphorylation of tau protein and subsequently Alzheimer's disease has been established, therefore, we have investigated the binding affinity and MARK4 inhibitory potential of cholic acid (CHA) using both computational and spectroscopic methods. Molecular docking suggested a strong binding of CHA to the functionally important residues of MARK4. We further performed 500 ns molecular dynamics simulation which suggested the MARK4-CHA system was quite stable throughout the simulation trajectory. CHA potential binds to the MARK4 with a binding constant (K) of 107 M−1 at 288 K. Further, MARK4 activity was inhibited by CHA with an IC 50 = 5.5 μM. Further insights into the thermodynamic parameters suggested that MARK4-CHA complex formation is driven by both electrostatic and van der Waals interactions. Overall study provides a rationale to use CHA in the drug development via MARK4 inhibition, towards possible therapeutic implications in neurodegenerative diseases. [ABSTRACT FROM AUTHOR]