Browsing by Subject "Molecular Dynamics Simulation"

Now showing 1 - 9 of 9
  • Thumbnail Image
    ItemOpen Access
    Anti-HIV Potential of Beesioside I Derivatives as Maturation Inhibitors: Synthesis, 3D-QSAR, Molecular Docking and Molecular Dynamics Simulations.
    (International journal of molecular sciences, 2023-01) Zhao, Zixuan; Ma, Yinghong; Li, Xiangyuan; Morris-Natschke, Susan L; Sun, Zhaocui; Sun, Zhonghao; Ma, Guoxu; Dong, Zhengqi; Zhao, Xiaohong; Yang, Meihua; Xu, Xudong; Lee, Kuohsiung; Wu, Haifeng; Chen, Chinho
    HIV-1 maturation is the final step in the retroviral lifecycle that is regulated by the proteolytic cleavage of the Gag precursor protein. As a first-in-class HIV-1 maturation inhibitor (MI), bevirimat blocks virion maturation by disrupting capsid-spacer peptide 1 (CA-SP1) cleavage, which acts as the target of MIs. Previous alterations of beesioside I (1) produced (20S,24S)-15ꞵ,16ꞵ-diacetoxy-18,24; 20,24-diepoxy-9,19-cyclolanostane-3ꞵ,25-diol 3-O-3',3'-dimethylsuccinate (3, DSC), showing similar anti-HIV potency compared to bevirimat. To ascertain the binding modes of this derivative, further modification of compound 1 was conducted. Three-dimensional quantitative structure-activity relationship (3D-QSAR) analysis combined with docking simulations and molecular dynamics (MD) were conducted. Five new derivatives were synthesized, among which compound 3b showed significant activity against HIV-1NL4-3 with an EC50 value of 0.28 µM. The developed 3D-QSAR model resulted in great predictive ability with training set (r2 = 0.99, q2 = 0.55). Molecular docking studies were complementary to the 3D-QSAR analysis, showing that DSC was differently bound to CA-SP1 with higher affinity than that of bevirimat. MD studies revealed that the complex of the ligand and the protein was stable, with root mean square deviation (RMSD) values <2.5 Å. The above results provided valuable insights into the potential of DSC as a prototype to develop new antiviral agents.
  • Thumbnail Image
    ItemOpen Access
    Characterizing protein crystal contacts and their role in crystallization: rubredoxin as a case study.
    (Soft Matter, 2014-01-14) Fusco, Diana; Headd, Jeffrey J; De Simone, Alfonso; Wang, Jun; Charbonneau, Patrick
    The fields of structural biology and soft matter have independently sought out fundamental principles to rationalize protein crystallization. Yet the conceptual differences and the limited overlap between the two disciplines have thus far prevented a comprehensive understanding of the phenomenon to emerge. We conduct a computational study of proteins from the rubredoxin family that bridges the two fields. Using atomistic simulations, we characterize the protein crystal contacts, and accordingly parameterize patchy particle models. Comparing the phase diagrams of these schematic models with experimental results enables us to critically examine the assumptions behind the two approaches. The study also reveals features of protein–protein interactions that can be leveraged to crystallize proteins more generally.
  • Thumbnail Image
    ItemOpen Access
    Dimensional study of the caging order parameter at the glass transition.
    (Proc Natl Acad Sci U S A, 2012-08-28) Charbonneau, Patrick; Ikeda, Atsushi; Parisi, Giorgio; Zamponi, Francesco
    The glass problem is notoriously hard and controversial. Even at the mean-field level, little is agreed upon regarding why a fluid becomes sluggish while exhibiting but unremarkable structural changes. It is clear, however, that the process involves self-caging, which provides an order parameter for the transition. It is also broadly assumed that this cage should have a gaussian shape in the mean-field limit. Here we show that this ansatz does not hold. By performing simulations as a function of spatial dimension d, we find the cage to keep a nontrivial form. Quantitative mean-field descriptions of the glass transition, such as mode-coupling theory, density functional theory, and replica theory, all miss this crucial element. Although the mean-field random first-order transition scenario of the glass transition is qualitatively supported here and non-mean-field corrections are found to remain small on decreasing d, reconsideration of its implementation is needed for it to result in a coherent description of experimental observations.
  • Thumbnail Image
    ItemOpen Access
    Drug design from the cryptic inhibitor envelope.
    (Nat Commun, 2016-02-25) Lee, Chul-Jin; Liang, Xiaofei; Wu, Qinglin; Najeeb, Javaria; Zhao, Jinshi; Gopalaswamy, Ramesh; Titecat, Marie; Sebbane, Florent; Lemaitre, Nadine; Toone, Eric J; Zhou, Pei
    Conformational dynamics plays an important role in enzyme catalysis, allosteric regulation of protein functions and assembly of macromolecular complexes. Despite these well-established roles, such information has yet to be exploited for drug design. Here we show by nuclear magnetic resonance spectroscopy that inhibitors of LpxC--an essential enzyme of the lipid A biosynthetic pathway in Gram-negative bacteria and a validated novel antibiotic target--access alternative, minor population states in solution in addition to the ligand conformation observed in crystal structures. These conformations collectively delineate an inhibitor envelope that is invisible to crystallography, but is dynamically accessible by small molecules in solution. Drug design exploiting such a hidden inhibitor envelope has led to the development of potent antibiotics with inhibition constants in the single-digit picomolar range. The principle of the cryptic inhibitor envelope approach may be broadly applicable to other lead optimization campaigns to yield improved therapeutics.
  • Thumbnail Image
    ItemOpen Access
    Evidence for an electrostatic mechanism of force generation by the bacteriophage T4 DNA packaging motor.
    (Nat Commun, 2014-06-17) Migliori, Amy D; Keller, Nicholas; Alam, Tanfis I; Mahalingam, Marthandan; Rao, Venigalla B; Arya, Gaurav; Smith, Douglas E
    How viral packaging motors generate enormous forces to translocate DNA into viral capsids remains unknown. Recent structural studies of the bacteriophage T4 packaging motor have led to a proposed mechanism wherein the gp17 motor protein translocates DNA by transitioning between extended and compact states, orchestrated by electrostatic interactions between complimentarily charged residues across the interface between the N- and C-terminal subdomains. Here we show that site-directed alterations in these residues cause force dependent impairments of motor function including lower translocation velocity, lower stall force and higher frequency of pauses and slips. We further show that the measured impairments correlate with computed changes in free-energy differences between the two states. These findings support the proposed structural mechanism and further suggest an energy landscape model of motor activity that couples the free-energy profile of motor conformational states with that of the ATP hydrolysis cycle.
  • Thumbnail Image
    ItemOpen Access
    Harnessing calcineurin-FK506-FKBP12 crystal structures from invasive fungal pathogens to develop antifungal agents.
    (Nature communications, 2019-09) Juvvadi, Praveen R; Fox, David; Bobay, Benjamin G; Hoy, Michael J; Gobeil, Sophie MC; Venters, Ronald A; Chang, Zanetta; Lin, Jackie J; Averette, Anna Floyd; Cole, D Christopher; Barrington, Blake C; Wheaton, Joshua D; Ciofani, Maria; Trzoss, Michael; Li, Xiaoming; Lee, Soo Chan; Chen, Ying-Lien; Mutz, Mitchell; Spicer, Leonard D; Schumacher, Maria A; Heitman, Joseph; Steinbach, William J
    Calcineurin is important for fungal virulence and a potential antifungal target, but compounds targeting calcineurin, such as FK506, are immunosuppressive. Here we report the crystal structures of calcineurin catalytic (CnA) and regulatory (CnB) subunits complexed with FK506 and the FK506-binding protein (FKBP12) from human fungal pathogens (Aspergillus fumigatus, Candida albicans, Cryptococcus neoformans and Coccidioides immitis). Fungal calcineurin complexes are similar to the mammalian complex, but comparison of fungal and human FKBP12 (hFKBP12) reveals conformational differences in the 40s and 80s loops. NMR analysis, molecular dynamic simulations, and mutations of the A. fumigatus CnA/CnB-FK506-FKBP12-complex identify a Phe88 residue, not conserved in hFKBP12, as critical for binding and inhibition of fungal calcineurin. These differences enable us to develop a less immunosuppressive FK506 analog, APX879, with an acetohydrazine substitution of the C22-carbonyl of FK506. APX879 exhibits reduced immunosuppressive activity and retains broad-spectrum antifungal activity and efficacy in a murine model of invasive fungal infection.
  • Thumbnail Image
    ItemOpen Access
    Infinite swapping replica exchange molecular dynamics leads to a simple simulation patch using mixture potentials.
    (J Chem Phys, 2013-02-28) Lu, Jianfeng; Vanden-Eijnden, Eric
    Replica exchange molecular dynamics (REMD) becomes more efficient as the frequency of swap between the temperatures is increased. Recently Plattner et al. [J. Chem. Phys. 135, 134111 (2011)] proposed a method to implement infinite swapping REMD in practice. Here we introduce a natural modification of this method that involves molecular dynamics simulations over a mixture potential. This modification is both simple to implement in practice and provides a better, energy based understanding of how to choose the temperatures in REMD to optimize efficiency. It also has implications for generalizations of REMD in which the swaps involve other parameters than the temperature.
  • Thumbnail Image
    ItemOpen Access
    Learning about Biomolecular Solvation from Water in Protein Crystals.
    (The journal of physical chemistry. B, 2018-03) Altan, Irem; Fusco, Diana; Afonine, Pavel V; Charbonneau, Patrick
    Water occupies typically 50% of a protein crystal and thus significantly contributes to the diffraction signal in crystallography experiments. Separating its contribution from that of the protein is, however, challenging because most water molecules are not localized and are thus difficult to assign to specific density peaks. The intricateness of the protein-water interface compounds this difficulty. This information has, therefore, not often been used to study biomolecular solvation. Here, we develop a methodology to surmount in part this difficulty. More specifically, we compare the solvent structure obtained from diffraction data for which experimental phasing is available to that obtained from constrained molecular dynamics (MD) simulations. The resulting spatial density maps show that commonly used MD water models are only partially successful at reproducing the structural features of biomolecular solvation. The radial distribution of water is captured with only slightly higher accuracy than its angular distribution, and only a fraction of the water molecules assigned with high reliability to the crystal structure is recovered. These differences are likely due to shortcomings of both the water models and the protein force fields. Despite these limitations, we manage to infer protonation states of some of the side chains utilizing MD-derived densities.
  • Thumbnail Image
    ItemOpen Access
    Solution NMR studies of the plant peptide hormone CEP inform function.
    (FEBS letters, 2013-12) Bobay, Benjamin G; DiGennaro, Peter; Scholl, Elizabeth; Imin, Nijat; Djordjevic, Michael A; Mck Bird, David
    The C-terminally Encoded Peptide (CEP) family of regulatory peptides controls root development in vascular plants. Here, we present the first NMR structures of CEP. We show that root-knot nematode (RKN: Meloidogyne spp.) also encodes CEP, presumably to mimic plant CEP as part of their stereotypic, parasitic interaction with vascular plants. Molecular dynamics simulations of plant- and nematode-encoded CEP displaying known posttranslational modifications (PTM) provided insight into the structural effects of PTM and the conformational plasticity and rigidity of CEP. Potential mechanisms of action are discussed with respect to the structure and sampling of conformational space.