Search

The calibration of torsional interaction terms by fitting relative gas-phase conformational energies against their quantum-mechanical values is a common procedure in force-field development. However, much less attention has been paid to the optimization of third-neighbor nonbonded interaction parameters, despite their strong coupling with the torsions. This article introduces an algorithm termed LLS-SC, aimed at simultaneously parametrizing torsional and third-neighbor interaction terms based on relative conformational energies. It relies on a self-consistent (SC) procedure where each iteration involves a linear least-squares (LLS) regression followed by a geometry optimization of the reference structures. As a proof-of-principle, this method is applied to obtain torsional and third-neighbor interaction parameters for aliphatic chains in the context of the GROMOS 53A6 united-atom force field. The optimized parameter set is compared to the original one, which has been fitted manually against thermodynamic properties for small linear alkanes. The LLS-SC implementation is freely available under http://github.com/mssm-labmmol/profiler. ISSN:0192-8651 ISSN:1096-987X

Hypoxia is an effective preconditioning stimulus and many cellular responses to hypoxia are mediated through a transcription control complex termed the hypoxia-inducible factor (HIF). The stability and activation of HIF are governed by HIF prolyl-4-hydroxylases 2 (PHD2). Hence, the development of a small molecule inhibitor for prolyl hydroxylase has been suggested as a potentially useful therapeutic strategy for the treatment of oxidative/ischemic stress conditions. Thus, to unveil a novel human PHD2 inhibitor, a custom-based virtual screening was carried out to identify the potential inhibitors against PHD2 based on; (1) the per-residue energy decomposition (PRED)-based pharmacophore model, (2) molecular docking, and (3) MD approaches. The PRED analysis was performed to identify the common interaction pattern of HIF fragment (5L9B) and crystallized ligand (4JZR) to develop a relevant accurate allosteric pharmacophore model. The custom pharmacophore model (AAARR) was developed and further used to screen multiple databases. The docking was performed as a secondary strategy for screening the pharmacophore hits. Furthermore, the docked complexes were screened by molecular dynamics (MD) simulation and molecular mechanics/generalized Born surface area (MM-GBSA) based binding free energy calculations to determine the binding energy of the inhibitors and to identify crucial interaction energy fingerprint. One hit has demonstrated good binding free energy and a better binding affinity for PHD2 compared to the other four selected ligands. Thus, the results obtained from pharmacophore, docking, and MD simulations depicted that linker length and metal binding in the scaffold could be effectively used as a potent inhibitor toward human PHD2 in AD therapeutics.

IOData is a free and open-source Python library for parsing, storing, and converting various file formats commonly used by quantum chemistry, molecular dynamics, and plane-wave density-functional-theory software programs. In addition, IOData supports a flexible framework for generating input files for various software packages. While designed and released for stand-alone use, its original purpose was to facilitate the interoperability of various modules in the HORTON and ChemTools software packages with external (third-party) molecular quantum chemistry and solid-state density-functional-theory packages. IOData is designed to be easy to use, maintain, and extend; this is why we wrote IOData in Python and adopted many principles of modern software development, including comprehensive documentation, extensive testing, continuous integration/delivery protocols, and package management. This article is the official release note of the IOData library.

The recently described flavin‐dependent halogenase BrvH is able to catalyse both the bromination and chlorination of indole, but shows significantly higher bromination activity. BrvH was annotated as a tryptophan halogenase, but does not accept tryptophan as a substrate. Its native substrate remains unknown. A predictive model with the data available for BrvH was analysed. A training set of compounds tested in vitro was docked into the active site of a complete protein model based on the X‐ray structure of BrvH. The atoms not resolved experimentally were modelled by using molecular mechanics force fields to obtain this protein model. Furthermore, docking poses for the substrates and known non‐substrates have been calculated. Parameters like distance, partial charge and hybridization state were analysed to derive rules for predicting activity. With this model for activity of the BrvH, a virtual screening suggested several structures for potential substrates. Some of the compounds preselected in this way were tested in vitro, and several could be verified as convertible substrates. Based on information on halogenated natural products, a new dataset was created to specifically search for natural products as substrates/products, and virtual screening in this database yielded further hits., To bind or not to bind? That is the question. The generation of a pharmacophore model adapted from a virtual screening methodology frequently employed in medicinal chemistry helped to suggest and identify new substrates for an orphan flavin‐dependent halogenase for which the native substrate is not yet known.

Prenyl pyrophosphate methyltransferases enhance the structural diversity of terpenoids. However, the molecular basis of their catalytic mechanisms is poorly understood. In this study, using multiple strategies, we characterized a geranyl pyrophosphate (GPP) C6-methyltransferase, BezA. Biochemical analysis revealed that BezA requires Mg2+ and solely methylates GPP. The crystal structures of BezA and its complex with S-adenosyl homocysteine were solved at 2.10 and 2.56 A, respectively. Further analyses using site-directed mutagenesis, molecular docking, molecular dynamics simulations, and quantum mechanics/molecular mechanics calculations revealed the molecular basis of the methylation reaction. Importantly, the function of E170 as a catalytic base to complete the methylation reaction was established. We also succeeded in switching the substrate specificity by introducing a W210A substitution, resulting in an unprecedented farnesyl pyrophosphate C6-methyltransferase.

Parkinson's disease (PD) is one of the most targeted neurodegenerative diseases in clinical research. Awareness of research is due to its increasing number of affected people worldwide. The pathology of PD has been linked to several key proteins upregulation such as the catechol O-Methyltransferase (COMT). Hence, the synthesis of compounds possessing inhibitory capacity has been the frontline of research in recent years. Several compounds have been synthesized among which is the nitrocatechol. However, major limitations associated with the nitrocatechol scaffold include the inability to possess adequate CNS penetration properties and hepatic toxicity associated with the compounds. However, a series of bicyclic hydroxypyridones compounds were synthesized to evaluate their inhibitory potentials on COMT protein with compound 38 (c38) 2-[(2,4-dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one shown to have a 40 fold increase level coverage in its IC50 over brain exposure when compared to the other synthesized compound. The molecular dynamics method was employed to understand the nature of interaction exhibited by c38. Molecular mechanics of c38 revealed a disruptive effect on the secondary structure of COMT protein. Per residue decomposition analysis revealed similar crucial residues involved in the favorable binding of c38 and tolcapone implicated its increased inhibitory capacity on COMT in preventing PD. Free binding energy (ΔGbind ) of c38 further revealed the inhibitory capacity towards COMT protein in comparison to the FDA approved tolcapone. Ligand mobility analysis of both compounds showed a timewise different mobility pattern across the simulation time frame at the active site pocket of the protein connoting the different inhibitory potency exhibited by c38 and tolcapone. Findings from this study revealed optimization of c38 could facilitate the discovery of new compounds with enhanced inhibitory properties towards COMT in treating PD.

Noscapine, an opium alkaloid, was discovered to bind tubulin, arrest dividing cells at mitosis, and selectively induce apoptosis to cancer cells. N-3-Br-Benzyl-Noscapine (Br-Bn-Nos), one of the derivatives of noscapine, was demonstrated to have improved anticancer potential compared with noscapine. We approached to evaluate the single and combined effect of Br-Bn-Nos and docetaxel (DOX) based on molecular modeling and cellular study. The individual predicted free energy of binding (∆Gbind,pred ) for Br-Bn-Nos and DOX with tubulin was found to be -28.89 and -36.07 kcal/mol based on molecular mechanics generalized Born solvation area (MM-GBSA) as well as -26.21 and -34.65 kcal/mol based on molecular mechanics Poisson Boltzmann solvation area (MM-PBSA), respectively. However, the ∆Gbind,pred of Br-Bn-Nos was significantly reduced (-33.02 and -30.24 kcal/mol using MM-GBSA and MM-PBSA) in the presence of DOX on its binding pocket. Parenthetically, the ∆Gbind,pred of DOX was significantly reduced (-37.17 and -32.80 kcal/mol using MM-GBSA and MM-PBSA) in the presence of Br-Bn-Nos on its binding pocket. The reduced ∆Gbind,pred in the presence of Br-Bn-Nos and DOX together indicated a combination effect of both the ligands. The combined interaction of both the agents onto tubulin dimmer was also determined experimentally using purified tubulin, in which a combination regimen of Br-Bn-Nos and DOX reduced the fluorescence intensity of tubulin to a higher value (68%) compared with the single regimen. Further, isobologram analysis revealed the synergistic effect of Br-Bn-Nos and DOX in antiproliferative activity using MCF-7 cell line at 48 hr (sum FIC = 0.49) and at 72 hr (sum FIC = 0.62). The combination dose regimen of Br-Bn-Nos and DOX blocks the cell cycle progression at the G2/M phase and induces apoptosis to cancer cells more effectively compared with the single regimen. Taken together, our study provides compelling evidence that the anticancer potential of noscapine derivatives may be substantially improved when it is used in a combined application with DOX for breast cancer.

Specifically targeting the With-No-Lysine (WNK1) kinase, which is implicated in hypertension, renders a significant challenge in discovering competitive inhibitors due to the highly conserved ATP binding pocket. However, an allosteric inhibitor may impart high specificity against the WNK kinase isoforms since it targets the less conserved site and can provide greater efficacy even under high physiological ATP concentration. In the current study, we have investigated the structural and energetic basis of the specificity of the allosteric inhibitor WNK476 against WNK kinase isoforms by combining molecular dynamics simulations and free energy calculations using molecular mechanics Poisson-Boltzmann surface area (MMPBSA). Our study reveals that the conformational stabilization of αC-helix near the allosteric binding site, including conformational changes in activation and glycine-rich loop regions, favors the specificity of WNK476 towards WNK1. The MM/PBSA calculations suggest that the non-polar contribution from hydrophobic residues and polar solvation energy influences WNK/WNK476 complexation. Despite more favorable electrostatic and van der Waals interactions in WNK2/WNK476, WNK476 is more potent against WNK1 due to the lower contribution of disfavouring components- polar solvation and entropy. Further, we have identified that the hydrophobic residues of DLG, αC-helix, β4 , and β5 regions, and H-bond network near the β4 strand play a critical role in the specificity of WNK476 against WNK1. Finally, our study reveals that residues Leu272 , Val281 , Phe283 , and Leu369 of WNK1 actively contribute to the overall hydrophobic interactions for WNK1/WNK476. Overall, our study might help in the rational design of novel allosteric inhibitors against hypertension.

Biomedicinally important histone lysine methyltransferases (KMTs) catalyze the transfer of a methyl group from S‐adenosylmethionine (AdoMet) cosubstrate to lysine residues in histones and other proteins. Herein, experimental and computational investigations on human KMT‐catalyzed ethylation of histone peptides by using S‐adenosylethionine (AdoEth) and Se‐adenosylselenoethionine (AdoSeEth) cosubstrates are reported. MALDI‐TOF MS experiments reveal that, unlike monomethyltransferases SETD7 and SETD8, methyltransferases G9a and G9a‐like protein (GLP) do have the capacity to ethylate lysine residues in histone peptides, and that cosubstrates follow the efficiency trend AdoMet>AdoSeEth>AdoEth. G9a and GLP can also catalyze AdoSeEth‐mediated ethylation of ornithine and produce histone peptides bearing lysine residues with different alkyl groups, such as H3K9meet and H3K9me2et. Molecular dynamics and free energy simulations based on quantum mechanics/molecular mechanics potential supported the experimental findings by providing an insight into the geometry and energetics of the enzymatic methyl/ethyl transfer process., Transfer window: Histone lysine methyltransferases, which catalyze the methylation of histones with S‐adenosylmethionine cosubstrate, have a limited capacity to catalyze ethylation of histones with S‐adenosylethionine (AdoEth) and Se‐adenosylselenoethionine (AdoSeEth) cosubstrates. Insights into the geometry and energetics of the enzymatic methyl/ethyl transfer process are reported.

Co-chaperon p23 has been well established as molecular chaperon for the heat shock protein 90 (Hsp90) that further leads to immorality in cancer cells by providing defense against Hsp90 inhibitors, and as stimulating agent for generating overexpressed antiapoptotic proteins, that is, Hsp70 and Hsp27. The natural compounds such as catechins from Camellia sinensis (green tea) are also well known for inhibition activity against various cancer. However, molecular interaction profile and potential lead bioactive compounds against co-chaperon p23 from green tea are not yet reported. To this context, we study the various secondary metabolites of green tea against co-chaperon p23 using structure-based virtual screening from Traditional Chinese Medicine (TCM) database. Following 26 compounds were obtained from TCM database and further studied for extra precision molecular docking that showed binding score between -10.221 and -2.276 kcal/mol with co-chaperon p23. However, relative docking score to known inhibitors, that is, ailanthone (-4.54 kcal/mol) and gedunin ( 3.60 kcal/mol) along with ADME profile analysis concluded epicatechin (-7.013 kcal/mol) and cis-theaspirone (-4.495 kcal/mol) as potential lead inhibitors from green tea against co-chaperone p23. Furthermore, molecular dynamics simulation and molecular mechanics generalized born surface area calculations validated that epicatechin and cis-theaspirone have significantly occupied the active region of co-chaperone p23 by hydrogen and hydrophobic interactions with various residues including most substantial amino acids, that is, Thr90, Ala94, and Lys95. Hence, these results supported the fact that green tea contained potential compounds with an ability to inhibit the cancer by disrupting the co-chaperon p23 activity.

This study explores the effects of different force fields and solvent water on the free binding energy calculation of NA‐OS and its mutants using the Gromacs software and g_mmpbsa tool. The recommended calculation models selected for working with small molecule antiviral drug systems (as exampled by NA‐OS) are amber03‐tip4p, amber99sb‐tip3p, and amber99sb‐ildn‐tip3p.

Therapeutic targeting of the adenosine triphosphate (ATP) machinery of Mycobacterium tuberculosis (Mtb) has recently presented a potent and alternative measure to halt the pathogenesis of tuberculosis. This has been potentiated by the development of bedaquiline (BDQ), a novel small molecule inhibitor that selectively inhibits mycobacterial F1 Fo -ATP synthase by targeting its rotor c-ring, resulting in the disruption of ATP synthesis and consequential cell death. Although the structural resolution of the mycobacterial C9 ring in co`mplex with BDQ provided the first-hand detail of BDQ interaction at the c-ring region of the ATP synthase, there still remains a need to obtain essential and dynamic insights into the mechanistic activity of this drug molecule towards crucial survival machinery of Mtb. As such, for the first time, we report an atomistic model to describe the structural dynamics that explicate the experimentally reported antagonistic features of BDQ in halting ion shuttling by the mycobacterial c-ring, using molecular dynamics simulation and the Molecular Mechanics/Poisson-Boltzmann Surface Area methods. Results showed that BDQ exhibited a considerably high ΔG while it specifically maintained high-affinity interactions with Glu65B and Asp32B , blocking their crucial roles in proton binding and shuttling, which is required for ATP synthesis. Moreover, the bulky nature of BDQ induced a rigid and compact conformation of the rotor c-ring, which impedes the essential rotatory motion that drives ion exchange and shuttling. In addition, the binding affinity of a BDQ molecule was considerably increased by the complementary binding of another BDQ molecule, which indicates that an increase in BDQ molecule enhances inhibitory potency against Mtb ATP synthase. Taken together, findings provide atomistic perspectives into the inhibitory mechanisms of BDQ coupled with insights that could enhance the structure-based design of novel ATP synthase inhibitors towards the treatment of tuberculosis.

The mechanism used by the ubiquitin-conjugating enzyme, Ubc13, to catalyze ubiquitination is probed with three computational techniques: Born-Oppenheimer molecular dynamics, single point quantum mechanics/molecular mechanics energies, and classical molecular dynamics. These simulations support a long-held hypothesis and show that Ubc13-catalyzed ubiquitination uses a stepwise, nucleophilic attack mechanism. Furthermore, they show that the first step-the formation of a tetrahedral, zwitterionic intermediate-is rate limiting. However, these simulations contradict another popular hypothesis that supposes that the negative charge on the intermediate is stabilized by a highly conserved asparagine (Asn79 in Ubc13). Instead, calculated reaction profiles of the N79A mutant illustrate how charge stabilization actually increases the barrier to product formation. Finally, an alternate role for Asn79 is suggested by simulations of wild-type, N79A, N79D, and H77A Ubc13: it stabilizes the motion of the electrophile prior to the reaction, positioning it for nucleophilic attack. © 2019 Wiley Periodicals, Inc.

The origin of the spectral shift from a red‐ to a green‐absorbing form in a cyanobacteriochrome, Slr1393g3, was identified by combined quantum mechanics/molecular mechanics simulations. This protein, related to classical phytochromes, carries the open‐chain tetrapyrrole chromophore phycocyanobilin. Our calculations reveal that the effective conjugation length in the chromophore becomes shorter upon conversion from the red to the green form. This is related to the planarity of the entire chromophore. A large distortion was found for the terminal pyrrole rings A and D; however, the D ring contributes more strongly to the photoproduct tuning, despite a larger change in the twist of the A ring. Our findings implicate that the D ring twist can be exploited to regulate the absorption of the photoproduct. Hence, mutations that affect the D ring twist can lead to rational tuning of the photoproduct absorption, allowing the tailoring of cyanobacteriochromes for biotechnological applications such as optogenetics and bioimaging.

The phosphoinositide 3-kinase γ (PI3Kγ) has been verified to be a potential drug target for the treatments of various human physical disorders. Although received lots of attention, the development of PI3Kγ-selective inhibitors is still a challenging subject because of its unique protein structural features. Aiming to uncover the interaction mechanism between the selective inhibitors and PI3Kγ, a series of benzothiazole and thiazolopiperidine PI3Kγ isoform-selective inhibitors were studied with an integrated in silico strategy by combining molecular docking, molecular dynamic simulations, binding free energy calculations, and decomposition analysis. Firstly, three molecular docking models, including rigid receptor docking, induced fit docking (IFD), and quantum mechanical-polarized ligand docking, were respectively, built, and the IFD preliminarily predicted the docking poses of all studied inhibitors and roughly analyzed the binding mechanism. Secondly, four binding complexes with representative inhibitors were selected to perform molecular dynamic simulations and free energy calculations. The predicted binding energies were consistent with the experimental bioactivities and different binding patterns between potent and weak inhibitors were uncovered. Finally, through the Molecular Mechanics/Generalized Born Surface Area binding free energy decomposition, residue-inhibitor interactions spectra were obtained and several key residues contributing to favorable binding were highlighted, which provides valuable information for rational PI3Kγ inhibitor design and modification.

Protein-drug binding mode prediction from the apo-protein structure is challenging because drug binding often induces significant protein conformational changes. Here, the authors report a computational workflow that incorporates a novel pocket generation method. First, the closed protein pocket is expanded by repeatedly filling virtual atoms during molecular dynamics (MD) simulations. Second, after ligand docking toward the prepared pocket structures, binding mode candidates are ranked by MD/Molecular Mechanics Poisson-Boltzmann Surface Area. The authors validated our workflow using CDK2 kinase, which has an especially-closed ATP-binding pocket in the apo-form, and several inhibitors. The crystallographic pose coincided with the top-ranked docking pose for 59% (34/58) of the compounds and was within the top five-ranked ones for 88% (51/58), while those estimated by a conventional prediction protocol were 9% (5/58) and 50% (29/58), respectively. Our study demonstrates that the prediction accuracy is significantly improved by preceding pocket expansion, leading to generation of conformationally-diverse binding mode candidates. © 2018 Wiley Periodicals, Inc.

Although the protein backbone is the most fundamental part of the structure, the fine-tuning of side-chain conformations is important for protein function, for example, in protein-protein and protein-ligand interactions, and also in enzyme catalysis. While several benchmarks testing the performance of protein force fields for side chain properties have already been published, they often considered only a few force fields and were not tested against the same experimental observables; hence, they are not directly comparable. In this work, we explore the ability of twelve force fields, which are different flavors of AMBER, CHARMM, OPLS, or GROMOS, to reproduce average rotamer angles and rotamer populations obtained from extensive NMR studies of the 3 J and residual dipolar coupling constants for two small proteins: ubiquitin and GB3. Based on a total of 196 μs sampling time, our results reveal that all force fields identify the correct side chain angles, while the AMBER and CHARMM force fields clearly outperform the OPLS and GROMOS force fields in estimating rotamer populations. The three best force fields for representing the protein side chain dynamics are AMBER 14SB, AMBER 99SB*-ILDN, and CHARMM36. Furthermore, we observe that the side chain ensembles of buried amino acid residues are generally more accurately represented than those of the surface exposed residues.

Mutation in Cu/Zn superoxide dismutase (SOD1) at position 85 from glycine to arginine was found to be a prominent cause of aggregation characterized by an increased content of β-sheets in familial amyotrophic lateral sclerosis (fALS). Various literatures reported that natural polyphenols could act as a β-sheet breaker and therefore, treated as a potential therapeutics against various aggregated proteins involved in neurodegenerative disorders. Through computational perspective, molecular docking, quantum chemical studies, and discrete molecular dynamics were implemented to study the binding and structural effect of natural polyphenols, kaempferol, and kaempferide on mutant SOD1. Kaempferol exhibited significant binding and greater residual energy contribution with mutant SOD1 than kaempferide. More interestingly, kaempferol was found to reduce the β-sheet content augmenting the mutant conformational stability and flexibility relative to that of kaempferide. Hence, the inhibition of mutant SOD1 aggregation by kaempferol was explored, thereby suggesting kaempferol could act as a drug candidate for the design of the natural therapeutics against fALS. © 2018 BioFactors, 44(5):431-442, 2018.

Static‐lattice atomistic simulations, in vacuum and solvent phase, have been recently employed to quantify the “in vitro—in vivo—in silico” performance‐correlation profile of various drug delivery systems and biomaterial scaffolds. The reactional profile of biopolymers was elucidated by exploring the spatial disposition of the molecular components with respect to the formulation conditions and the final release medium. This manuscript provides a brief overview of recently completed and published studies related to molecular tectonics of: (a) the nanoformation and solvation properties of the surfactant‐emulsified polymeric systems; (b) the formation and chemistry of polyelectrolyte complexes; (c) the effect of a plasticizer and/or drug on the physicomechanical properties of biomedical archetypes; (d) the molecular modeling templates to predict stimuli‐ and environmentally esponsive systems; and (e) the polymer‐mucopeptide complexes and intermacromolecular networks. Furthermore, this report provides a detailed account of the role of molecular mechanics energy relationships toward the interpretation and understanding of the mechanisms that control the formation, fabrication, selection, design, performance, complexation, interaction, stereospecificity, and preference of various biopolymeric systems for biomedical applications.

Leprosy (causative, Mycobacterium leprae) continues to be the persisting public health problem with stable incidence rates, owing to the emergence of dapsone resistance that being the principal drug in the ongoing multidrug therapy. Hence, to overcome the drug resistance, structural modification through medicinal chemistry was used to design newer dapsone derivative(s) (DDs), against folic acid biosynthesis pathway. The approach included theoretical modeling, molecular docking, and molecular dynamic (MD) simulation as well as binding free energy estimation for validation of newly designed seven DDs, before synthesis. Theoretical modeling, docking, and MD simulation studies were used to understand the mode of binding and efficacy of DDs against the wild-type and mutant dihydropteroate synthases (DHPS). Principal component analysis was performed to understand the conformational dynamics of DHPS-DD complexes. Furthermore, the overall stability and negative-binding free energy of DHPS-DD complexes were deciphered using Molecular Mechanics/Poisson-Boltzmann Surface Area technique. Molecular mechanics study revealed that DD3 possesses higher binding free energy than dapsone against mutant DHPS. Energetic contribution analysis portrayed that van der Waals and electrostatic energy contributes profoundly to the overall negative free energy, whereas polar solvation energy opposes the binding. Finally, DD3 was synthesized and characterized using Fourier-transform infrared spectroscopy, UV, liquid chromatography-mass spectrometry, and proton nuclear magnetic resonance techniques. This study suggested that DD3 could be further promoted as newer antileprosy agent. The principles of medicinal chemistry and bioinformatics tools help to locate effective therapeutics to minimize resources and time in current drug development modules.

Aggregatibacter actinomycetemcomitans is a Gram-negative bacteria that has gained wide recognition for its causative role in the development of various immune diseases, which includes localized aggressive periodontitis. Its ability to evade host defense mechanisms is mediated by the secretion of leukotoxin (LtxA), which induces death of white blood cells (leukocytes) by specific binding to their surface-expressed leukocyte function-associated receptor (LFA-1) in its active state. Therapeutic compounds that interfere with this pathogenic process and abrogate A. actinomycetemcomitans virulence have been reported in literature. These include doxycycline, and more recently phytochemical compounds such as hamamelitanin, resveratrol, naringin, and quercetin. However, the question remains how do they work? Therefore, with the aid of computational tools, we explore the molecular mechanisms by which they possibly elicit their therapeutic functions. Molecular mechanics Poisson/Boltzmann surface area analyses revealed that these compounds bind favorably to active LFA-1 with high affinity and considerable stability, indicative of their ability to occupy the LtxA binding site (LBS) and prevent LtxA binding. The conformational transition of open LFA-1 to its closed state further describe the mechanistic activity of these compounds. In addition to notable reductions in structural mobility and flexibility, the burial of surface-exposed interactive side chains at the LBS was observed, an occurrence that could alter the complementary binding of LtxA. It is also important to mention that these occurrences were induced more prominently by the phytochemicals. We believe that these findings will enhance the scope of drug design and discovery for potent LtxA antagonists with improved activities and therapeutic efficacies in the treatment of virulent A. actinomycetemcomitans diseases.

Because of the lack of structural studies on aminoglycoside phosphotransferase (APH) from prevalent volatile human pathogen Bacillus cereus, aminoglycoside resistance therapeutics research remains elusive. Hence, in this computational study, we have performed homology modeling, molecular docking, molecular dynamics (MD), and principal component analysis studies on APH from B. cereus. The structure of APH was predicted by homology modeling using MODELLER 9v12 and validated for its stereochemical qualities. Sequence analysis study of the template (Protein Data Bank ID: 3TDW) and APH from B. cereus sensu lato group showed exact matching of active-site residues. The mechanism of substrate and inhibitor binding to APH was studied using molecular docking, which identified GTP as the more preferred substrate, whereas ZINC71575479 as the most effective inhibitor. The active-site residues, ARG41, TYR90, ASP195, and ASP215 at nucleotide triphosphate-binding cavity of APH were found to be involved in binding with substrate and inhibitor. Molecular dynamics simulation study of APH in apo form and bound form confirmed the stability and effective binding of GTP and ZINC71575479 in a dynamic state. Molecular mechanics Poisson-Boltzmann surface area calculations revealed energetic contributions of active-site residues of APH in binding with GTP and ZINC71575479. The principal component analysis revealed the internal global motion of APH in apo and complex form. Furthermore, experimental studies on APH from B. cereus ATCC 10876 validated the in silico findings for its inhibition. Thus, this study provides more information on structure-function relationships of APH from B. cereus and open avenues for designing effective strategies to overcome antibiotic resistance.

Fucosidosis is a rare lysosomal storage disease which has been classified into two subtypes, depending on the severity of clinical signs and symptoms. Fucosidosis patients’ skin abnormalities include angiokeratoma corporis diffusum, widespread telangiectasia, thick skin, hyperhidrosis and hypohidrosis, acrocyanosis and distal transverse nail bands. It has been described that >50% of fucosidosis patients have angiokeratoma. At molecular level, fucosidosis is caused by lysosomal alpha-L-fucosidase (FUCA1) gene mutations. Obtaining samples for functional studies has been challenging due to the inherent difficulty in finding affected individuals. The effect of FUCA1 dysfunction on gene expression is unknown. The aim of this study was to analyse, in keratinocytes, the transcriptomic effect of FUCA1 knock-down for a better understanding of skin lesions’ pathogenesis affecting fucosidosis patients. FUCA1 knock-down (siRNA) was performed in human HaCaT immortalised keratinocytes. Affymetrix arrays and qPCR were used for analysing gene expression. Bioinformatics was used for functional clustering of modified genes. In total, 387 genes showed differential expression between FUCA1 silenced and non-silenced cells (222 up-regulated and 165 down-regulated). Up-regulated genes belonged to two major groups: keratinocyte differentiation/epidermal development (n = 17) and immune response (n = 61). Several transcription factors were up-regulated in FUCA1-siRNA transfected cells. This effect might partly have been produced by abnormal transcription factor expression, that is FOXN1. We thus propose that fucosidosis-related skin lesions (eg angiokeratoma) and those of other diseases (eg psoriasis) might be caused by dysfunctions in common aetiological overlapping molecular cascades. © 2018 John Wiley and Sons A/S. Published by John Wiley and Sons Ltd

The electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method has been successfully utilized for efficient linear-scaling quantum mechanical (QM) calculation of protein energies. In this work, we applied the EE-GMFCC method for calculation of binding affinity of Endonuclease colicin-immunity protein complex. The binding free energy changes between the wild-type and mutants of the complex calculated by EE-GMFCC are in good agreement with experimental results. The correlation coefficient (R) between the predicted binding energy changes and experimental values is 0.906 at the B3LYP/6-31G*-D level, based on the snapshot whose binding affinity is closest to the average result from the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) calculation. The inclusion of the QM effects is important for accurate prediction of protein-protein binding affinities. Moreover, the self-consistent calculation of PB solvation energy is required for accurate calculations of protein-protein binding free energies. This study demonstrates that the EE-GMFCC method is capable of providing reliable prediction of relative binding affinities for protein-protein complexes. © 2018 Wiley Periodicals, Inc.

Recently, stereoinversions and isomerizations of amino acid residues in the proteins of living beings have been observed. Because isomerized amino acids cause structural changes and denaturation of proteins, isomerizations of amino acid residues are suspected to cause age-related diseases. In this study, AMBER molecular force field parameters were tested by using computationally generated nonapeptides and tripeptides including stereoinverted and/or isomerized amino acid residues. Energy calculations by using density functional theory were also performed for comparison. Although the force field parameters were developed by parameter fitting for l-α-amino acids, the accuracy of the computational results for d-amino acids and β-amino acids was comparable to those for l-α-amino acids. The conformational energies for tripeptides calculated by using density functional theory were reproduced more accurately than those for nonapeptides calculated by using the molecular mechanical force field. The evaluations were performed for the ff99SB, ff03, ff12SB, and the latest ff14SB force field parameters.

Due to the difficulties in monoclonal antibody production specific to mycotoxins, aptameric probes have been considered as suitable alternatives. The low efficiency of the SELEX procedure in screening high affinity aptamers for binding mycotoxins as small molecules can be significantly improved through computational techniques. Previously, we designed five new aptamers to aflatoxin B1 (AFB1) based on a known aptamer sequence (Patent: PCT/ CA2010 / 001 292, Apt1) through a genetic algorithm-based in silico maturation strategy and experimentally measured their affinity to the target toxin. Here, integrated molecular dynamic simulation (MDs) studies with molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis to clarify the binding modes, critical interacting nucleic bases and energy component contributions in the six AFB1-binding aptamers. The aptamer F20, which was selected in the first work, showed the best free binding energy and complex stability compared to other aptamers. The trajectory analysis revealed that AFB1 recognized F20 through the groove binding mode along with precise shape complementarity. The MD simulation results revealed that dynamic water intermediate interactions also play a key role in promoting complex stability. According to the MM-PBSA calculations, Van der Waals contacts were identified as dominant energy components in all complexes. Interestingly, a high consistency is observed between the experimentally obtained binding affinities of the six aptamers with their free energy solvation. The computational findings, confirmed via previous experiments, highlighted the binding modes, the dynamic hydration of complex components and the total free interacting energy as the crucial criteria in discovering high functional aptameric probes. This article is protected by copyright. All rights reserved.

ABCG2 is a substantial member of the ABC transporter superfamily that plays a significant role in multidrug resistance in cancer. Until recently, the 3D structure of ABCG2 has not been resolved, which resulted in the limitation of developing potential ABCG2 inhibitors using structure-based drug discovery. Herein, eMolecules, ChEMBL, and ChEBI databases, containing >25 million compounds, were virtually screened against the ABCG2 transporter in homodimer form. Performance of AutoDock4.2.6 software to predict inhibitor-ABCG2 binding mode and affinity were validated on the basis of available experimental data. The explored databases were filtered based on docking scores. The most potent hits with binding affinities higher than that of experimental bound ligand (MZ29) were then selected and subjected to molecular mechanics minimization, followed by binding energy calculation using molecular mechanics-generalized Born surface area (MM-GBSA) approach. Furthermore, molecular dynamics simulations for 50 ns, followed by MM-GBSA binding energy calculations, were performed for the promising compounds, unveiling eight potential inhibitors with binding affinities

Bovine and camel chymosins are aspartic proteases that are used in dairy food manufacturing. Both enzymes catalyze proteolysis of a milk protein, κ-casein, which helps to initiate milk coagulation. Surprisingly, camel chymosin shows a 70% higher clotting activity than bovine chymosin for bovine milk, while exhibiting only 20% of the unspecific proteolytic activity. By contrast, bovine chymosin is a poor coagulant for camel milk. Although both enzymes are marketed commercially, the disparity in their catalytic activity is not yet well understood at a molecular level, due in part to a lack of atomistic resolution data about the chymosin—κ-casein complexes. Here, we report computational alanine scanning calculations of all four chymosin—κ-casein complexes, allowing us to elucidate the influence that individual residues have on binding thermodynamics. Of the 12 sequence differences in the binding sites of bovine and camel chymosin, eight are shown to be particularly important for understanding differences in the binding thermodynamics (Asp112Glu, Lys221Val, Gln242Arg, Gln278Lys. Glu290Asp, His292Asn, Gln294Glu, and Lys295Leu. Residue in bovine chymosin written first). The relative binding free energies of single-point mutants of chymosin are calculated using the molecular mechanics three dimensional reference interaction site model (MM-3DRISM). Visualization of the solvent density functions calculated by 3DRISM reveals the difference in solvation of the binding sites of chymosin mutants.

BACKGROUND Fusarium graminearum is the main pathogen of Fusarium Head Blight (FHB), a worldwide plant disease and one of the major wheat diseases in China. The control of the FHB is mainly dependent on the application of DMIs fungicides. Fungal sterol 14α-demethylase enzymes (CYP51) are the main target for DMIs fungicides. In order to investigate the binding mechanism between azoles and CYP51B in F.graminearum, the molecular modeling study and biological evaluation were performed. RESULTS Firstly, the homology model based on the crystal structure of Aspergillus fumigatus was built. Secondly, molecular docking and molecular dynamics (MD) simulation were used to identify the optimum binding mode between Propiconazole (PRP), Diniconazole (DIN), Triadimenol (TRL), Tebuconazole (TEC), and Triadimefon (TRN) with FgCYP51B. Furthermore, binding free energy of the five protein-inhibitor complexes was calculated by Molecular Mechanics Generalized Born Surface Area and Poisson–Boltzmann Surface Area (MM-GB/PBSA) methods. The key residues to the selective binding of azoles to FgCYP51B were recognized by per-residue free energy decomposition analysis. The five ligands have similar binding mode in the active pocket. Inhibitors PRP and TEC have more beneficial binding free energy to the enzyme than TRN, TRL and DIN. Furthermore, the Amino acid residues Phe511, Val136, Ile374, Ala308, Ser312 and Try137 of FgCYP51B are the key residues interacting with azoles fungicides. The EC50 value of PRP, TEC, DIN, TRL and TRN is 0.024, 0.047, 0.148, 0.154 and 0.474 mg/L according to the experimental evaluation, respectively. These five molecules exhibited potential inhibitory activity against CYP51B protein from the F. graminearum. CONCLUSION The results suggest that azoles fungicides for FgCYP51B should possess more hydrophobic groups interacting with residues Phe511, Val136, Ile374, Ala308, Ser312 and Tyr137. PRP and TEC should be preferred to control Fusarium Head Blight disease than DIN, TRL and TRN.

Free-energy calculations based on molecular simulations provide access to a wide range of thermodynamic properties such as the solubility and partitioning of a molecule in and between various phases. It is demonstrated how molecular dynamics free-energy simulations may be used to obtain the solubilities of primary alcohols and n-alkanes in water and binding affinities of primary alcohols to α-cyclodextrin. The equivalence of two distinct routes to calculate binding free energies is shown leading to the conclusion that host-guest binding affinities may be used to probe the underlying molecular force field.

Quorum sensing (QS) is a chemical communication between bacteria, with which bacteria sense the population of their own species. Autoinducer-2 (AI-2) is a class of universal quorum-sensing molecules, which is used by both Gram-negative and Gram-positive bacteria. The inhibition of AI-2-mediated QS has various practical applications, including the prevention of the formation of biofilm in dental gums. In this work, we develop a computational protocol for developing AI-2 inhibitors. A challenging aspect of such an endeavor is that the receptor undergoes a large conformational change upon ligand binding. We combine several methods such as molecular docking with multiple conformations, molecular dynamics simulations, and molecular mechanics Poisson–Boltzmann computation, in order to estimate binding affinity of candidate molecules to a quorum-sensing receptor. We apply our method to rank the substances in a chemical library. We indeed find a molecule that has a higher affinity than previously known ligands, thus showing the feasibility of the protocol for the development of quorum-sensing inhibitors.

Prediction of drug metabolism is an important step in the early lead optimization phase and preclinical studies. Its main purpose is to avoid late stage (Phase I–III) or post-marketing withdrawals which usually results in a significant financial loss to a pharmaceutical company. Computational methods for identification of soft-spots in the molecules (site of metabolism; SOM) and CYP450 isoforms responsible for drug metabolism have over the past ten years proved an indispensable tool towards achieving this goal. These methods can be broadly classified into i) ligand based, ii) structure based and iii) hybrid methods. In this review we discuss some of the most successful methods of SOM and isoform specificity prediction, their application domain, advantages, and limitations along with recent developments in the last 5 years in this area. Recent applications of molecular dynamics (MD), quantum mechanics (QM) and quantum mechanics/molecular mechanics (QM/MM) methodology for regioselective and isoform specific drug metabolism predictions are also discussed. Finally, a summary of these methods along with expected developments, future challenges and opportunities are discussed.

To gain more insight into the factors controlling the efficient cysteine arylation by cyclometalated Au(III) complexes, the reaction between selected gold compounds and different peptides was investigated by high‐resolution liquid chromatography electrospray ionization mass spectrometry (HR‐LC‐ESI‐MS). The deducted mechanisms of C–S cross‐coupling, also supported by density functional theory (DFT) and quantum mechanics/molecular mechanics (QM/MM) calculations, evidenced the key role of secondary peptidic gold binding sites in favouring the process of reductive elimination.

An understanding of the fate of organic compounds originating from plants in soil is crucial for determining their persistence and concentrations in the environment. Aristolochic acids are believed to be the causal agents that induce Balkan endemic nephropathy by food contamination through soil adsorption of humic acids, major components of soil. Aristolochic acids are active chemicals in Aristolochia plant species found in endemic villages. In this article, molecular structure interactions between 18 structures of aristolochic acids with an inserted humic acid structure were studied. These structures were optimized in vacuo and by periodic box simulation with water solvate using the computational molecular mechanics MM+ method with HyperChem software. The QSPR models were used for correlation of the relationship between the hydrophobicity values of 18 AA structures coupled with a HA structure by MM+ and QSAR+ properties. Computational hydrophobicity values were considered dependent variables and were related to the structural features obtained by molecular and quantum mechanics calculations by multiple linear regression approaches. The obtained model was validated, and the results indicated differing hydrophobicity between the MM+ and QSAR+ properties.

MexAB-OprM efflux pumps, found in Pseudomonas aeruginosa, play a major role in drug resistance by extruding out drugs and antibiotic molecules from cells. Inhibitors are used to cease the potency of the efflux pumps. In this study, in-silico models are used to investigate the nature of the binding pocket of the MexAB-OprM efflux pump. First, we have performed classical molecular dynamics (MD) simulations to shed light on different aspects of protein-inhibitor interaction in the binding pocket of the pump. Using classical MD simulations, quantum mechanics/molecular mechanics (QM/MM), and various types of analyses, it is found that D13-9001 has a higher binding affinity towards the binding pocket compared to D1 and D2; the results are in sync with the experimental dat. Two stable configurations of D13-9001 are discovered inside the distal pocket which could be one of the primary reasons for the greater efficacy of D13-9001. The free energy barrier upon changing one state to another is calculated by employing umbrella sampling method. Finally, F178 is mutated to have the complete picture as it contributes significantly to the binding energy irrespective of the three inhibitors. Our results may help to design a new generation of inhibitors for such an efflux pump.

Understanding the interactions of drug molecules with biomacromolecules at a micro-scale level is essential to design potent drugs for the treatments of human genome diseases. To unravel the mechanism of binding of aristololactam-β-D-glucoside (ADG) and phenylalanine transfer RNA (tRNAPhe), an integrated computational strategy combining quantum mechanics (QM) calculation, molecular docking and atomistic molecular dynamics (MD) simulation was present in this work. QM calculations were performed to derive the partial charges of ADG, molecular docking was used to determine the binding poses of ADG on the tRNAPhe, and atomistic MD simulations were conducted to examine the thermal stability of five predicted binding poses for the complex of ADG and the tRNAPhe. The binding free energies of the five complexes were then calculated using the molecular mechanics/generalized born surface area approach with the variable internal dielectric constant model. By comparing computed affinities and experimentally-measured values in the binding free energy, we identified a most likely binding structure of ADG and the tRNAPhe. Further analysis of energy of the ADG-tRNA complex revealed that the aristololactam of ADG provides binding specificity to the tRNAPhe, and the D-glucoce contributes to the affinity for binding with the tRNAPhe.

Molecular dynamics simulations and quantum mechanics/molecular mechanics calculations were performed on the in silico Leu597Ala/Ile663Ala double mutant of rabbit ALOX15 (12/15 lipoxygenase). The computational results suggested that subtle steric hindrance by the conserved Leu597 and C-terminal Ile663 residues disturbed H10 abstractions in wildtype ALOX15 (which abstracts H13), but if these two bulky residues were mutated to smaller ones, H10 abstraction was no longer impeded and the regioselectivity of the initial H-abstraction step was changed. However, site-directed mutagenesis with HPLC analysis of the products of the whole oxidation process showed that the regioselectivity of the hydroperoxidation was not altered. This disagreement may be explained by the conformational reorganization of the system needed to rotate the -OO. group from an antarafacial to a suprafacial arrangement prior to back-hydrogen transfer. After H10 abstraction and O2 insertion, the evolution of the peroxy radical at C12 was sterically impeded, whereas peroxyl group rotation at C15 (after H13 abstraction) could easily evolve to a suprafacial arrangement, which thus led to the final product. For this reason, the global regiospecificity was not affected in the mutant. These findings exemplify that the regioselectivity of initial hydrogen abstraction and the regioselectivity of the final product do not necessarily coincide (in fact, they can be opposite) for the hydroperoxidation of arachidonic acid catalyzed by a lipoxygenase.

Histone acetylation is a very important regulatory mechanism in gene expression in the chromatin context. A new protein family-YEATS domains have been found as a novel histone acetylation reader, which could specific recognize the histone lysine acetylation. AF9 is an important one in the YEATS family. Focused on the AF9-H3K9ac (K9 acetylation) complex (ALY) (PDB code: 4TMP) and a serials of mutants, MUT (the acetyllsine of H3K9ac was mutated to lysine), F59A, G77A, and D103A, we applied molecular dynamics simulation and molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculations to examine the role of AF9 protein in recognition interaction. The simulation results and analysis indicate that some residues of the protein have significant influence on recognition and binding to H3K9ac peptides and hydrophobic surface show the hydrophobic interactions play an important role in the binding. Our work can give important information to understand how the protein AF9 recognizes the peptides H3K9ac. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 779-786, 2016.

This paper is devoted to the understanding of the reaction mechanism of mycobacterium tuberculosis glutamine synthetase (mtGS) with atomic detail, using computational quantum mechanics/molecular mechanics (QM/MM) methods at the ONIOM M06-D3/6-311++G(2d,2p):ff99SB//B3LYP/6-31G(d):ff99SB level of theory. The complete reaction undergoes a three-step mechanism: the spontaneous transfer of phosphate from ATP to glutamate upon ammonium binding (ammonium quickly loses a proton to Asp54), the attack of ammonia on phosphorylated glutamate (yielding protonated glutamine), and the deprotonation of glutamine by the leaving phosphate. This exothermic reaction has an activation free energy of 21.5 kcal mol(-1) , which is consistent with that described for Escherichia coli glutamine synthetase (15-17 kcal mol(-1) ). The participating active site residues have been identified and their role and energy contributions clarified. This study provides an insightful atomic description of the biosynthetic reaction that takes place in this enzyme, opening doors for more accurate studies for developing new anti-tuberculosis therapies.

Quantum mechanics/molecular mechanics calculations are employed to assign previously recorded experimental spectroscopic signatures of the intermediates occurring during the photo-induced repair of (6-4) photolesions by photolyases to specific molecular structures. Based on this close comparison of experiment and theory it is demonstrated that the acting repair mechanism involves proton transfer from the protonated His365 to the N3' nitrogen of the lesion, which proceeds simultaneously with intramolecular OH transfer along an oxetane-like transition state.

International audience; More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all-atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. (C) 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Photoswitchable azobenzene cross-linkers can control the folding and unfolding of peptides by photoisomerization and can thus regulate peptide affinities and enzyme activities. Using quantum mechanics/molecular mechanics (QM/MM) methods and classical MM force fields, we report the first molecular dynamics simulations of the photoinduced folding and unfolding processes in the azobenzene cross-linked FK-11 peptide. We find that the interactions between the peptide and the azobenzene cross-linker are crucial for controlling the evolution of the secondary structure of the peptide and responsible for accelerating the folding and unfolding events. They also modify the photoisomerization mechanism of the azobenzene cross-linker compared with the situation in vacuo or in solution.

Nature has developed large groups of enzymatic catalysts with the aim to transfer substrates into useful products, which enables biosystems to perform all their natural functions. As such, all biochemical processes in our body (we drink, we eat, we breath, we sleep, etc.) are governed by enzymes. One of the problems associated with research on biocatalysts is that they react so fast that details of their reaction mechanisms cannot be obtained with experimental work. In recent years, major advances in computational hardware and software have been made and now large (bio)chemical systems can be studied using accurate computational techniques. One such technique is the quantum mechanics/molecular mechanics (QM/MM) technique, which has gained major momentum in recent years. Unfortunately, it is not a black-box method that is easily applied, but requires careful set-up procedures. In this work we give an overview on the technical difficulties and caveats of QM/MM and discuss work-protocols developed in our groups for running successful QM/MM calculations.