Browsing by Subject "Molecular Dynamics Simulation"

Now showing 1 - 5 of 5
  • Thumbnail Image
    Item
    Consecutive non-natural PZ nucleobase pairs in DNA impact helical structure as seen in 50 μs molecular dynamics simulations
    (Oxford University Press, 2017-04-20) Molt, Robert W.; Georgiadis, Millie M.; Richards, Nigel G. J.; Department of Biochemistry and Molecular Biology, School of Medicine
    Z: Little is known about the influence of multiple consecutive 'non-standard' ( , 6-amino-5-nitro-2(1H)-pyridone, and , 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) nucleobase pairs on the structural parameters of duplex DNA. nucleobase pairs follow standard rules for Watson-Crick base pairing but have rearranged hydrogen bonding donor and acceptor groups. Using the X-ray crystal structure as a starting point, we have modeled the motions of a DNA duplex built from a self-complementary oligonucleotide (5΄-CTTATPPPZZZATAAG-3΄) in water over a period of 50 μs and calculated DNA local parameters, step parameters, helix parameters, and major/minor groove widths to examine how the presence of multiple, consecutive nucleobase pairs might impact helical structure. In these simulations, the -containing DNA duplex exhibits a significantly wider major groove and greater average values of stagger, slide, rise, twist and h-rise than observed for a 'control' oligonucleotide in which nucleobase pairs are replaced by . The molecular origins of these structural changes are likely associated with at least two differences between and . First, the electrostatic properties of differ from in terms of density distribution and dipole moment. Second, differences are seen in the base stacking of pairs in dinucleotide steps, arising from energetically favorable stacking of the nitro group in with π-electrons of the adjacent base.
  • Thumbnail Image
    Item
    How does Water Pass through a Sugar Transporter?
    (Elsevier B.V., 2014-03-18) Zhu, Fangqiang; Department of Physics, School of Science
  • Thumbnail Image
    Item
    Molecular Recognition in a Diverse Set of Protein-Ligand Interactions Studied with Molecular Dynamics Simulations and End-Point Free Energy Calculations
    (ACS Publications, 2013-10-28) Wang, Bo; Li, Liwei; Hurley, Thomas D.; Meroueh, Samy O.; Department of Biochemistry & Molecular Biology, School of Medicine
    End-point free energy calculations using MM-GBSA and MM-PBSA provide a detailed understanding of molecular recognition in protein-ligand interactions. The binding free energy can be used to rank-order protein-ligand structures in virtual screening for compound or target identification. Here, we carry out free energy calculations for a diverse set of 11 proteins bound to 14 small molecules using extensive explicit-solvent MD simulations. The structure of these complexes was previously solved by crystallography and their binding studied with isothermal titration calorimetry (ITC) data enabling direct comparison to the MM-GBSA and MM-PBSA calculations. Four MM-GBSA and three MM-PBSA calculations reproduced the ITC free energy within 1 kcal•mol−1 highlighting the challenges in reproducing the absolute free energy from end-point free energy calculations. MM-GBSA exhibited better rank-ordering with a Spearman ρ of 0.68 compared to 0.40 for MM-PBSA with dielectric constant (ε = 1). An increase in ε resulted in significantly better rank-ordering for MM-PBSA (ρ = 0.91 for ε = 10). But larger ε significantly reduced the contributions of electrostatics, suggesting that the improvement is due to the non-polar and entropy components, rather than a better representation of the electrostatics. SVRKB scoring function applied to MD snapshots resulted in excellent rank-ordering (ρ = 0.81). Calculations of the configurational entropy using normal mode analysis led to free energies that correlated significantly better to the ITC free energy than the MD-based quasi-harmonic approach, but the computed entropies showed no correlation with the ITC entropy. When the adaptation energy is taken into consideration by running separate simulations for complex, apo and ligand (MM-PBSAADAPT), there is less agreement with the ITC data for the individual free energies, but remarkably good rank-ordering is observed (ρ = 0.89). Interestingly, filtering MD snapshots by pre-scoring protein-ligand complexes with a machine learning-based approach (SVMSP) resulted in a significant improvement in the MM-PBSA results (ε = 1) from ρ = 0.40 to ρ = 0.81. Finally, the non-polar components of MM-GBSA and MM-PBSA, but not the electrostatic components, showed strong correlation to the ITC free energy; the computed entropies did not correlate with the ITC entropy.
  • Thumbnail Image
    Item
    The Role of Semidisorder in Temperature Adaptation of Bacterial FlgM Proteins
    (Elsevier B.V., 2013-12-03) Wang, Jihua; Yang, Yuedong; Cao, Zanxia; Li, Zhixiu; Zhao, Huiying; Zhou, Yaoqi; Department of Biochemistry & Molecular Biology, IU School of Medicine
    Probabilities of disorder for FlgM proteins of 39 species whose optimal growth temperature ranges from 273 K (0°C) to 368 K (95°C) were predicted by a newly developed method called Sequence-based Prediction with Integrated NEural networks for Disorder (SPINE-D). We showed that the temperature-dependent behavior of FlgM proteins could be separated into two subgroups according to their sequence lengths. Only shorter sequences evolved to adapt to high temperatures (>318 K or 45°C). Their ability to adapt to high temperatures was achieved through a transition from a fully disordered state with little secondary structure to a semidisordered state with high predicted helical probability at the N-terminal region. The predicted results are consistent with available experimental data. An analysis of all orthologous protein families in 39 species suggests that such a transition from a fully disordered state to semidisordered and/or ordered states is one of the strategies employed by nature for adaptation to high temperatures.
  • Thumbnail Image
    Item
    α-Tocopherol Is Well Designed to Protect Polyunsaturated Phospholipids: MD Simulations.
    (Elsevier, 2015-10-20) Leng, Xiaoling; Kinnun, Jacob J.; Marquardt, Drew; Ghefli, Mikel; Kučerka, Norbert; Katsaras, John; Atkinson, Jeffrey; Harroun, Thad A.; Feller, Scott E.; Wassall, Stephen R.; Department of Physics, School of Science
    The presumptive function for alpha-tocopherol (αtoc) in membranes is to protect polyunsaturated lipids against oxidation. Although the chemistry of the process is well established, the role played by molecular structure that we address here with atomistic molecular-dynamics simulations remains controversial. The simulations were run in the constant particle NPT ensemble on hydrated lipid bilayers composed of SDPC (1-stearoyl-2-docosahexaenoylphosphatidylcholine, 18:0-22:6PC) and SOPC (1-stearoyl-2-oleoylphosphatidylcholine, 18:0-18:1PC) in the presence of 20 mol % αtoc at 37°C. SDPC with SA (stearic acid) for the sn-1 chain and DHA (docosahexaenoic acid) for the sn-2 chain is representative of polyunsaturated phospholipids, while SOPC with OA (oleic acid) substituted for the sn-2 chain serves as a monounsaturated control. Solid-state 2H nuclear magnetic resonance and neutron diffraction experiments provide validation. The simulations demonstrate that high disorder enhances the probability that DHA chains at the sn-2 position in SDPC rise up to the bilayer surface, whereby they encounter the chromanol group on αtoc molecules. This behavior is reflected in the van der Waals energy of interaction between αtoc and acyl chains, and illustrated by density maps of distribution for acyl chains around αtoc molecules that were constructed. An ability to more easily penetrate deep into the bilayer is another attribute conferred upon the chromanol group in αtoc by the high disorder possessed by DHA. By examining the trajectory of single molecules, we found that αtoc flip-flops across the SDPC bilayer on a submicrosecond timescale that is an order-of-magnitude greater than in SOPC. Our results reveal mechanisms by which the sacrificial hydroxyl group on the chromanol group can trap lipid peroxyl radicals within the interior and near the surface of a polyunsaturated membrane. At the same time, water-soluble reducing agents that regenerate αtoc can access the chromanol group when it locates at the surface.