1
|
Linzer JT, Aminov E, Abdullah AS, Kirkup CE, Diaz Ventura RI, Bijoor VR, Jung J, Huang S, Tse CG, Álvarez Toucet E, Onghai HP, Ghosh AP, Grodzki AC, Haines ER, Iyer AS, Khalil MK, Leong AP, Neuhaus MA, Park J, Shahid A, Xie M, Ziembicki JM, Simmerling C, Nagan MC. Accurately Modeling RNA Stem-Loops in an Implicit Solvent Environment. J Chem Inf Model 2024. [PMID: 39002142 DOI: 10.1021/acs.jcim.4c00756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2024]
Abstract
Ribonucleic acid (RNA) molecules can adopt a variety of secondary and tertiary structures in solution, with stem-loops being one of the more common motifs. Here, we present a systematic analysis of 15 RNA stem-loop sequences simulated with molecular dynamics simulations in an implicit solvent environment. Analysis of RNA cluster ensembles showed that the stem-loop structures can generally adopt the A-form RNA in the stem region. Loop structures are more sensitive, and experimental structures could only be reproduced with modification of CH···O interactions in the force field, combined with an implicit solvent nonpolar correction to better model base stacking interactions. Accurately modeling RNA with current atomistic physics-based models remains challenging, but the RNA systems studied herein may provide a useful benchmark set for testing other RNA modeling methods in the future.
Collapse
Affiliation(s)
- Jason T Linzer
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Ethan Aminov
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Aalim S Abdullah
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Colleen E Kirkup
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Rebeca I Diaz Ventura
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Vinay R Bijoor
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Jiyun Jung
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Sophie Huang
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Chi Gee Tse
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Emily Álvarez Toucet
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Hugo P Onghai
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Arghya P Ghosh
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Alex C Grodzki
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Emilee R Haines
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Aditya S Iyer
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Mark K Khalil
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Alexander P Leong
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Michael A Neuhaus
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Joseph Park
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Asir Shahid
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Matthew Xie
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Jan M Ziembicki
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| | - Carlos Simmerling
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - Maria C Nagan
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
2
|
Stroet M, Setz M, Lee T, Malde AK, van den Bergen G, Sykacek P, Oostenbrink C, Mark AE. On the Validation of Protein Force Fields Based on Structural Criteria. J Phys Chem B 2024; 128:4602-4620. [PMID: 38711373 PMCID: PMC11103706 DOI: 10.1021/acs.jpcb.3c08469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 04/25/2024] [Accepted: 04/29/2024] [Indexed: 05/08/2024]
Abstract
Molecular dynamics simulations depend critically on the quality of the force field used to describe the interatomic interactions and the extent to which it has been validated for use in a specific application. Using a curated test set of 52 high-resolution structures, 39 derived from X-ray diffraction and 13 solved using NMR, we consider the extent to which different parameter sets of the GROMOS protein force field can be distinguished based on comparing a range of structural criteria, including the number of backbone hydrogen bonds, the number of native hydrogen bonds, polar and nonpolar solvent-accessible surface area, radius of gyration, the prevalence of secondary structure elements, J-coupling constants, nuclear Overhauser effect (NOE) intensities, positional root-mean-square deviations (RMSD), and the distribution of backbone ϕ and ψ dihedral angles. It is shown that while statistically significant differences between the average values of individual metrics could be detected, these were in general small. Furthermore, improvements in agreement in one metric were often offset by loss of agreement in another. The work establishes a framework and test set against which protein force fields can be validated. It also highlights the danger of inferring the relative quality of a given force field based on a small range of structural properties or small number of proteins.
Collapse
Affiliation(s)
- Martin Stroet
- The
University of Queensland, St. Lucia, Queensland 4072, Australia
| | - Martina Setz
- Institute
for Molecular Modeling and Simulation, Department of Material Science
and Process Engineering, University of Natural
Resources and Life Sciences, Vienna Muthgasse 18, 1190 Vienna, Austria
| | - Thomas Lee
- The
University of Queensland, St. Lucia, Queensland 4072, Australia
| | - Alpeshkumar K. Malde
- Institute
for Glycomics and School of Environment and Science, Griffith University, Gold Coast, Queensland 4222, Australia
| | | | - Peter Sykacek
- Institute
of Computational Biology, Department of Biotechnology, University of Natural Resources and Life Sciences,
Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Chris Oostenbrink
- Institute
for Molecular Modeling and Simulation, Department of Material Science
and Process Engineering, University of Natural
Resources and Life Sciences, Vienna Muthgasse 18, 1190 Vienna, Austria
- Christian
Doppler Laboratory for Molecular Informatics in the Biosciences, University of Natural Resources and Life Sciences,
Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Alan E. Mark
- The
University of Queensland, St. Lucia, Queensland 4072, Australia
| |
Collapse
|
3
|
Zhang P, Yang W. Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein. J Chem Phys 2023; 159:024118. [PMID: 37431910 PMCID: PMC10481389 DOI: 10.1063/5.0142280] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2023] [Accepted: 06/22/2023] [Indexed: 07/12/2023] Open
Abstract
Molecular dynamics (MD) is an extremely powerful, highly effective, and widely used approach to understanding the nature of chemical processes in atomic details for proteins. The accuracy of results from MD simulations is highly dependent on force fields. Currently, molecular mechanical (MM) force fields are mainly utilized in MD simulations because of their low computational cost. Quantum mechanical (QM) calculation has high accuracy, but it is exceedingly time consuming for protein simulations. Machine learning (ML) provides the capability for generating accurate potential at the QM level without increasing much computational effort for specific systems that can be studied at the QM level. However, the construction of general machine learned force fields, needed for broad applications and large and complex systems, is still challenging. Here, general and transferable neural network (NN) force fields based on CHARMM force fields, named CHARMM-NN, are constructed for proteins by training NN models on 27 fragments partitioned from the residue-based systematic molecular fragmentation (rSMF) method. The NN for each fragment is based on atom types and uses new input features that are similar to MM inputs, including bonds, angles, dihedrals, and non-bonded terms, which enhance the compatibility of CHARMM-NN to MM MD and enable the implementation of CHARMM-NN force fields in different MD programs. While the main part of the energy of the protein is based on rSMF and NN, the nonbonded interactions between the fragments and with water are taken from the CHARMM force field through mechanical embedding. The validations of the method for dipeptides on geometric data, relative potential energies, and structural reorganization energies demonstrate that the CHARMM-NN local minima on the potential energy surface are very accurate approximations to QM, showing the success of CHARMM-NN for bonded interactions. However, the MD simulations on peptides and proteins indicate that more accurate methods to represent protein-water interactions in fragments and non-bonded interactions between fragments should be considered in the future improvement of CHARMM-NN, which can increase the accuracy of approximation beyond the current mechanical embedding QM/MM level.
Collapse
Affiliation(s)
- Pan Zhang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
4
|
Dou B, Zhu Z, Merkurjev E, Ke L, Chen L, Jiang J, Zhu Y, Liu J, Zhang B, Wei GW. Machine Learning Methods for Small Data Challenges in Molecular Science. Chem Rev 2023; 123:8736-8780. [PMID: 37384816 PMCID: PMC10999174 DOI: 10.1021/acs.chemrev.3c00189] [Citation(s) in RCA: 21] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/01/2023]
Abstract
Small data are often used in scientific and engineering research due to the presence of various constraints, such as time, cost, ethics, privacy, security, and technical limitations in data acquisition. However, big data have been the focus for the past decade, small data and their challenges have received little attention, even though they are technically more severe in machine learning (ML) and deep learning (DL) studies. Overall, the small data challenge is often compounded by issues, such as data diversity, imputation, noise, imbalance, and high-dimensionality. Fortunately, the current big data era is characterized by technological breakthroughs in ML, DL, and artificial intelligence (AI), which enable data-driven scientific discovery, and many advanced ML and DL technologies developed for big data have inadvertently provided solutions for small data problems. As a result, significant progress has been made in ML and DL for small data challenges in the past decade. In this review, we summarize and analyze several emerging potential solutions to small data challenges in molecular science, including chemical and biological sciences. We review both basic machine learning algorithms, such as linear regression, logistic regression (LR), k-nearest neighbor (KNN), support vector machine (SVM), kernel learning (KL), random forest (RF), and gradient boosting trees (GBT), and more advanced techniques, including artificial neural network (ANN), convolutional neural network (CNN), U-Net, graph neural network (GNN), Generative Adversarial Network (GAN), long short-term memory (LSTM), autoencoder, transformer, transfer learning, active learning, graph-based semi-supervised learning, combining deep learning with traditional machine learning, and physical model-based data augmentation. We also briefly discuss the latest advances in these methods. Finally, we conclude the survey with a discussion of promising trends in small data challenges in molecular science.
Collapse
Affiliation(s)
- Bozheng Dou
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Zailiang Zhu
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Ekaterina Merkurjev
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Lu Ke
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Long Chen
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Jian Jiang
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Yueying Zhu
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Jie Liu
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Bengong Zhang
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Guo-Wei Wei
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
5
|
Kognole AA, Aytenfisu AH, MacKerell AD. Extension of the CHARMM Classical Drude Polarizable Force Field to N- and O-Linked Glycopeptides and Glycoproteins. J Phys Chem B 2022; 126:6642-6653. [PMID: 36005290 PMCID: PMC9463114 DOI: 10.1021/acs.jpcb.2c04245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Molecular dynamic simulations are an effective tool to study complex molecular systems and are contingent upon the availability of an accurate and reliable molecular mechanics force field. The Drude polarizable force field, which allows for the explicit treatment of electronic polarization in a computationally efficient fashion, has been shown to reproduce experimental properties that were difficult or impossible to reproduce with the CHARMM additive force field, including peptide folding cooperativity, RNA hairpin structures, and DNA base flipping. Glycoproteins are essential components of glycoconjugate vaccines, antibodies, and many pharmaceutically important molecules, and an accurate polarizable force field that includes compatibility between the protein and carbohydrate aspect of the force field is essential to study these types of systems. In this work, we present an extension of the Drude polarizable force field to glycoproteins, including both N- and O-linked species. Parameter optimization focused on the dihedral terms using a reweighting protocol targeting NMR solution J-coupling data for model glycopeptides. Validation of the model include eight model glycopeptides and four glycoproteins with multiple N- and O-linked glycosylations. The new glycoprotein carbohydrate force field can be used in conjunction with the remainder of Drude polarizable force field through a variety of MD simulation programs including GROMACS, OPENMM, NAMD, and CHARMM and may be accessed through the Drude Prepper module in the CHARMM-GUI.
Collapse
Affiliation(s)
| | | | - Alexander D. MacKerell
- Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, United States
| |
Collapse
|
6
|
Molecular Evolution of the Pseudomonas aeruginosa DNA Gyrase gyrA Gene. Microorganisms 2022; 10:microorganisms10081660. [PMID: 36014079 PMCID: PMC9415716 DOI: 10.3390/microorganisms10081660] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Revised: 08/10/2022] [Accepted: 08/13/2022] [Indexed: 11/29/2022] Open
Abstract
DNA gyrase plays important roles in genome replication in various bacteria, including Pseudomonasaeruginosa. The gyrA gene encodes the gyrase subunit A protein (GyrA). Mutations in GyrA are associated with resistance to quinolone-based antibiotics. We performed a detailed molecular evolutionary analyses of the gyrA gene and associated resistance to the quinolone drug, ciprofloxacin, using bioinformatics techniques. We produced an evolutionary phylogenetic tree using the Bayesian Markov Chain Monte Carlo (MCMC) method. This tree indicated that a common ancestor of the gene was present over 760 years ago, and the offspring formed multiple clusters. Quinolone drug-resistance-associated amino-acid substitutions in GyrA, including T83I and D87N, emerged after the drug was used clinically. These substitutions appeared to be positive selection sites. The molecular affinity between ciprofloxacin and the GyrA protein containing T83I and/or D87N decreased significantly compared to that between the drug and GyrA protein, with no substitutions. The rate of evolution of the gene before quinolone drugs were first used in the clinic, in 1962, was significantly lower than that after the drug was used. These results suggest that the gyrA gene evolved to permit the bacterium to overcome quinolone treatment.
Collapse
|
7
|
Nyepetsi M, Mbaiwa F, Oyetunji OA, de Leeuw NH. Understanding the Interactions between Triolein and Cosolvent Binary Mixtures Using Molecular Dynamics Simulations. ACS OMEGA 2022; 7:10212-10224. [PMID: 35382278 PMCID: PMC8973112 DOI: 10.1021/acsomega.1c06762] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Accepted: 03/01/2022] [Indexed: 06/14/2023]
Abstract
Biodiesel is one of the emerging renewable sources of energy to replace fossil-fuel-based resources. It is produced by a transesterification reaction in which a triglyceride reacts with methanol in the presence of a catalyst. The reaction is slow because of the low solubility of methanol in triglycerides, which results in low concentrations of methanol available to react with triglyceride. To speed up the reaction, cosolvents are added to create a single phase which helps to improve the concentration of methanol in the triglyceride phase. In this study, molecular dynamics simulations are used to help understand the role of cosolvents in the solvation of triglyceride (triolein). Six binary mixtures of triolein/cosolvent were used to study the solvation of triolein at 298.15 K. Results of 100 ns simulations at constant temperature and pressure to simulate mixing experiments show that in the first 10 ns all the binary mixtures remain largely unmixed. However, for the cosolvents that are fully miscible with triolein, the partial densities across the simulation boxes show that the systems are fully mixed in the final 10 ns. Some solvents were found to interact strongly with the polar part of triolein, while others interacted with the aliphatic part. The radial distribution functions and clustering of the solvents around triolein were also used as indicators for solvation. The presence of cosolvents also influenced the conformation of triolein molecules. In the presence of solvents that solubilize it, triolein preferred a propeller conformation but took up a trident conformation when there is less or no solubilization. The results show that tetrahydrofuran is the best solvent at solubilizing triolein, followed by cyclopentyl methyl ether, diethyl ether, and hexane. With 1,4-dioxane, the solubility improves with an increase in temperature. The miscibility of a solvent in triolein is aided by its ability to interact with both the polar and nonpolar parts of triolein.
Collapse
Affiliation(s)
- Maipelo Nyepetsi
- Department
of Chemical and Forensic Sciences, Botswana
International University of Science and Technology (BIUST), Palapye, Botswana
| | - Foster Mbaiwa
- Department
of Chemical and Forensic Sciences, Botswana
International University of Science and Technology (BIUST), Palapye, Botswana
| | | | - Nora H. de Leeuw
- School
of Chemistry, Cardiff University, Cardiff CF10 3AT, United Kingdom
- School
of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
8
|
Kognole AA, Lee J, Park SJ, Jo S, Chatterjee P, Lemkul JA, Huang J, MacKerell AD, Im W. CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field. J Comput Chem 2021; 43:359-375. [PMID: 34874077 DOI: 10.1002/jcc.26795] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Revised: 11/10/2021] [Accepted: 11/25/2021] [Indexed: 12/18/2022]
Abstract
Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. In addition, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. The capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.
Collapse
Affiliation(s)
- Abhishek A Kognole
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Jumin Lee
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Sang-Jun Park
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Sunhwan Jo
- Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois, USA
| | - Payal Chatterjee
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, Virginia, USA
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Zhejiang, Hangzhou, China
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, USA
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania, USA
| |
Collapse
|
9
|
Preparing and Analyzing Polarizable Molecular Dynamics Simulations with the Classical Drude Oscillator Model. Methods Mol Biol 2021. [PMID: 34302679 DOI: 10.1007/978-1-0716-1468-6_13] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
Molecular dynamics (MD) simulations performed with force fields that include explicit electronic polarization are becoming more prevalent in the field. The increasing emergence of these simulations is a result of continual refinement against a range of theoretical and empirical target data, optimization of software algorithms for higher performance, and availability of graphical processing unit hardware to further accelerate the simulations. Polarizable MD simulations are likely to be most impactful in biomolecular systems in which heterogeneous environments or unique microenvironments exist that would lead to inaccuracies in simulations performed with fixed-charge, nonpolarizable force fields. The further adoption of polarizable MD simulations will benefit from tutorial material that specifically addresses preparing and analyzing their unique features. In this chapter, we introduce common protocols for preparing routine biomolecular systems containing proteins, including both a globular protein in aqueous solvent and a transmembrane model peptide in a phospholipid bilayer. Details and example input files are provided for preparation of the simulation system using CHARMM, performing the simulations with OpenMM, and analyzing interesting dipole moment properties in CHARMM.
Collapse
|
10
|
Yuan Y, Ma Y, Huo D, Mills MJL, Wei J, Su W, Zhang R. Multipolar Description of Atom-Atom Electrostatic Interaction Energies in Single/Double-Stranded DNAs. J Phys Chem B 2020; 124:10089-10103. [PMID: 33138384 DOI: 10.1021/acs.jpcb.0c06757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular force field simulation is an effective method to explore the properties of DNA molecules in depth. Almost all current popular force fields calculate atom-atom electrostatic interaction energies for DNAs based on the atomic charge and dipole or quadrupole moments, without considering high-rank atomic multipole moments for more accurate electrostatics. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. In this work, a multipole expansion method that allows us to combine polarizability and anisotropy was applied. One single-stranded DNA and one double-stranded DNA were selected as pilot systems. Deoxynucleotides were cut out from pilot systems and capped by mimicking the original DNA environment. Atomic multipole moments were integrated instead of fixed-point charges to calculate atom-atom electrostatic energies to improve the accuracy of force fields for DNA simulations. Also, the applicability of modeling the behavior of both single-stranded and double-stranded DNAs was investigated. The calculation results indicated that the models can be transferred from pilot systems to test systems, which is of great significance for the development of future DNA force fields.
Collapse
Affiliation(s)
- Yongna Yuan
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou 730000, China
| | - Yan Ma
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou 730000, China
| | - Dongxu Huo
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou 730000, China
| | - Matthew J L Mills
- 3M Corporate Research Analytical Laboratory, Saint Paul, Minnesota 55114, United States
| | - Jiaxuan Wei
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou 730000, China
| | - Wei Su
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou 730000, China
| | - Ruisheng Zhang
- School of Information Science & Engineering, Lanzhou University, No. 222 South Tianshui Road, Lanzhou 730000, China
| |
Collapse
|
11
|
Sauceda HE, Gastegger M, Chmiela S, Müller KR, Tkatchenko A. Molecular force fields with gradient-domain machine learning (GDML): Comparison and synergies with classical force fields. J Chem Phys 2020; 153:124109. [DOI: 10.1063/5.0023005] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Huziel E. Sauceda
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- BASLEARN, BASF-TU Joint Lab, Technische Universität Berlin, 10587 Berlin, Germany
| | - Michael Gastegger
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- BASLEARN, BASF-TU Joint Lab, Technische Universität Berlin, 10587 Berlin, Germany
- DFG Cluster of Excellence “Unifying Systems in Catalysis” (UniSysCat), Technische Universität Berlin, 10623 Berlin, Germany
| | - Stefan Chmiela
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
| | - Klaus-Robert Müller
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- Department of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul 136-713, South Korea
- Max Planck Institute for Informatics, Stuhlsatzenhausweg, 66123 Saarbrücken, Germany
- Google Research, Brain Team, Berlin, Germany
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| |
Collapse
|
12
|
Yang S, Liu JDH, Diem M, Wesseling S, Vervoort J, Oostenbrink C, Rietjens IMCM. Molecular Dynamics and In Vitro Quantification of Safrole DNA Adducts Reveal DNA Adduct Persistence Due to Limited DNA Distortion Resulting in Inefficient Repair. Chem Res Toxicol 2020; 33:2298-2309. [PMID: 32786539 DOI: 10.1021/acs.chemrestox.0c00097] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The formation and repair of N2-(trans-isosafrol-3'-yl)-2'-deoxyguanosine (S-3'-N2-dG) DNA adduct derived from the spice and herbal alkenylbenzene constituent safrole were investigated. DNA adduct formation and repair were studied in vitro and using molecular dynamics (MD) simulations. DNA adduct formation was quantified using liquid chromatography-mass spectrometry (LCMS) in wild type and NER (nucleotide excision repair) deficient CHO cells and also in HepaRG cells and primary rat hepatocytes after different periods of repair following exposure to safrole or 1'-hydroxysafrole (1'-OH safrole). The slower repair of the DNA adducts found in NER deficient cells compared to that in CHO wild type cells indicates a role for NER in repair of S-3'-N2-dG DNA adducts. However, DNA repair in liver cell models appeared to be limited, with over 90% of the adducts remaining even after 24 or 48 h recovery. In our further studies, MD simulations indicated that S-3'-N2-dG adduct formation causes only subtle changes in the DNA structure, potentially explaining inefficient activation of NER. Inefficiency of NER mediated repair of S-3'-N2-dG adducts points at persistence and potential bioaccumulation of safrole DNA adducts upon daily dietary exposure.
Collapse
Affiliation(s)
- Shuo Yang
- Division of Toxicology, Wageningen University and Research, Stippeneng 4, 6708 WE Wageningen, The Netherlands
| | - Jakob D H Liu
- Institute of Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria
| | - Matthias Diem
- Institute of Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria
| | - Sebastiaan Wesseling
- Division of Toxicology, Wageningen University and Research, Stippeneng 4, 6708 WE Wageningen, The Netherlands
| | - Jacques Vervoort
- Division of Biochemistry, Wageningen University and Research, Stippeneng 4, 6708 WE Wageningen, The Netherlands
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria
| | - Ivonne M C M Rietjens
- Division of Toxicology, Wageningen University and Research, Stippeneng 4, 6708 WE Wageningen, The Netherlands
| |
Collapse
|
13
|
Petrov D, Tunega D, Gerzabek MH, Oostenbrink C. Molecular modelling of sorption processes of a range of diverse small organic molecules in Leonardite humic acid. EUROPEAN JOURNAL OF SOIL SCIENCE 2020; 71:831-844. [PMID: 33041627 PMCID: PMC7540484 DOI: 10.1111/ejss.12868] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2019] [Revised: 06/27/2019] [Accepted: 07/01/2019] [Indexed: 05/08/2023]
Abstract
Soil organic matter (SOM) is abundant in the environment and plays an important role in several biogeochemical processes, including microbial activity, soil aggregation, plant growth and carbon storage. One of its key functions is the retention and release of various chemical compounds, primarily governed by the sorption process, which strongly affects the environmental fate of nutrients and pollutants. Sorption largely depends on the composition of SOM, as well as its structure, dynamics and the thermodynamic conditions. Although several approaches are available, experimental characterization of sorption mechanisms is not easy. Computational models for predicting sorption coefficients often require a wealth of experimental data for training and are only applicable to compounds and conditions related to the training dataset. Here, we use molecular dynamics (MD) simulations to study the sorption of a range of small organic compounds. As a model SOM system we use the standard Leonardite humic acid (LHA) sample, which physicochemical properties have recently been characterized computationally in detail. This model allowed us to estimate sorption propensities of the systems at two different hydration levels (water activities close to 0 and 1), showing a remarkable correlation with experimental data. Importantly, this molecular modelling approach based on perturbation free-energy calculations is rigorously derived from statistical thermodynamics and requires no experimental sorption data for training. It is therefore in principle applicable to any SOM model or thermodynamic condition. Moreover, the power of MD simulations to provide high-resolution insight into atomistic and molecular interactions was employed to explore how sorbate molecules associate with the LHA matrix and which contacts they form. The heteroatoms of both sorbate and sorbent play an important role and water molecules are identified as further key players in facilitating the sorption process. HIGHLIGHTS Modelling of the sorption processes in soil organic matter at atomistic level.Rigorous, physics-based approach applicable to a range of SOM systems and conditions.Remarkable level of matching with experimental data with additional insight into the molecular mechanism.Interactions between the sorbate and local environment strongly affects the sorption process.
Collapse
Affiliation(s)
- Drazen Petrov
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| | - Daniel Tunega
- Department of Forest and Soil SciencesInstitute of Soil Research, University of Natural Resources and Life Sciences ViennaViennaAustria
- School of Pharmaceutical Science and Technology, Tianjin UniversityTianjinPeople's Republic of China
| | - Martin H. Gerzabek
- Department of Forest and Soil SciencesInstitute of Soil Research, University of Natural Resources and Life Sciences ViennaViennaAustria
| | - Chris Oostenbrink
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| |
Collapse
|
14
|
Ghosh SK, Saha B, Banerjee R. Insight into the sequence-structure relationship of TLR cytoplasm's Toll/Interleukin-1 receptor domain towards understanding the conserved functionality of TLR 2 heterodimer in mammals. J Biomol Struct Dyn 2020; 39:5348-5357. [PMID: 32643540 DOI: 10.1080/07391102.2020.1786457] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The signaling response of TLR2 to ligands has always been as a homodimer or in heterodimerization with TLR1/TLR6. The Toll/Interleukin-1 Receptor (TIR) domain of the TLR cytoplasmic region regulates the dimerization and interactions with adaptor molecules to build an active signaling complex. To understand the conservation of functionality of the TLR2-heterodimers between the distantly related species human(h) and mice(m), the pattern of TIR-TIR interaction in heterodimers has been studied through the sequence-structural point of view. Comparative analysis of primary sequence and structural pattern of TLRs(1/2/6) corroborates higher sequence homology between TLR1 and TLR6. Molecular docking analysis of TLR2-TLR1 and TLR2-TLR6 cytoplasmic dimers in both mouse and human have identified that for interaction the BB loop/near-BB loop residues of TLR2 are involved with the near-DD loop of TLR1 and DD loop residues of TLR6 within the TIR domains, which may cause to differential signaling. Molecular dynamics simulation of dimers for both human and mice species recognize stable interface between near-BB/BB loop region of TLR2 and discrete near-DD and DD loop region of TLR1 and TLR6 respectively. The observed dimerization pattern in both the species is further supported by Alanine scanning mutation study. However, Solvent Accessible Surface Area (SASA) of BB and DD loop regions of the cytoplasmic monomers and the heterodimers suggests that while TLR2 BB loop is actively associated as the dimer interface with its heterodimer partners in both the species, the DD loop acts as the active interfacing region in hTLR1 and mTLR6. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Soumya Kanti Ghosh
- Department of Bioinformatics, Maulana Abul Kalam Azad University of Technology, Kolkata, India
| | | | - Raja Banerjee
- Department of Bioinformatics, Maulana Abul Kalam Azad University of Technology, Kolkata, India.,Department of Biotechnology, Maulana Abul Kalam Azad University of Technology, Kolkata, India
| |
Collapse
|
15
|
Bitencourt-Ferreira G, de Azevedo WF. Molecular Dynamics Simulations with NAMD2. Methods Mol Biol 2020; 2053:109-124. [PMID: 31452102 DOI: 10.1007/978-1-4939-9752-7_8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
X-ray diffraction crystallography is the primary technique to determine the three-dimensional structures of biomolecules. Although a robust method, X-ray crystallography is not able to access the dynamical behavior of macromolecules. To do so, we have to carry out molecular dynamics simulations taking as an initial system the three-dimensional structure obtained from experimental techniques or generated using homology modeling. In this chapter, we describe in detail a tutorial to carry out molecular dynamics simulations using the program NAMD2. We chose as a molecular system to simulate the structure of human cyclin-dependent kinase 2.
Collapse
Affiliation(s)
- Gabriela Bitencourt-Ferreira
- Escola de Ciências da Saúde, Pontifícia Universidade Católica do Rio Grande do Sul-PUCRS, Porto Alegre, RS, Brazil
| | - Walter Filgueira de Azevedo
- Escola de Ciências da Saúde, Pontifícia Universidade Católica do Rio Grande do Sul-PUCRS, Porto Alegre, RS, Brazil.
| |
Collapse
|
16
|
Kognole AA, Aytenfisu AH, MacKerell AD. Balanced polarizable Drude force field parameters for molecular anions: phosphates, sulfates, sulfamates, and oxides. J Mol Model 2020; 26:152. [PMID: 32447472 DOI: 10.1007/s00894-020-04399-0] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2020] [Accepted: 04/28/2020] [Indexed: 12/14/2022]
Abstract
Polarizable force fields are emerging as a more accurate alternative to additive force fields in terms of modeling and simulations of a variety of chemicals including biomolecules. Explicit treatment of induced polarization in charged species such as phosphates and sulfates offers the potential for achieving an improved atomistic understanding of the physical forces driving their interactions with their environments. To help achieve this, in this study we present balanced Drude polarizable force field parameters for molecular ions including phosphates, sulfates, sulfamates, and oxides. Better balance was primarily achieved in the relative values of minimum interaction energies and distances of the anionic model compounds with water at the Drude and quantum mechanical (QM) model chemistries. Parametrization involved reoptimizing available parameters as well as extending the force field to new molecules with the goal of achieving self-consistency with respect to the Lennard-Jones and electrostatic parameters targeting QM and experimental hydration free energies. The resulting force field parameters achieve consistent treatment across the studied anions, facilitating more balanced simulations of biomolecules and small organic molecules in the context of the classical Drude polarizable force field. Graphical abstract.
Collapse
Affiliation(s)
- Abhishek A Kognole
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD, 21201, USA
| | - Asaminew H Aytenfisu
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD, 21201, USA
| | - Alexander D MacKerell
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD, 21201, USA.
| |
Collapse
|
17
|
Aprà E, Bylaska EJ, de Jong WA, Govind N, Kowalski K, Straatsma TP, Valiev M, van Dam HJJ, Alexeev Y, Anchell J, Anisimov V, Aquino FW, Atta-Fynn R, Autschbach J, Bauman NP, Becca JC, Bernholdt DE, Bhaskaran-Nair K, Bogatko S, Borowski P, Boschen J, Brabec J, Bruner A, Cauët E, Chen Y, Chuev GN, Cramer CJ, Daily J, Deegan MJO, Dunning TH, Dupuis M, Dyall KG, Fann GI, Fischer SA, Fonari A, Früchtl H, Gagliardi L, Garza J, Gawande N, Ghosh S, Glaesemann K, Götz AW, Hammond J, Helms V, Hermes ED, Hirao K, Hirata S, Jacquelin M, Jensen L, Johnson BG, Jónsson H, Kendall RA, Klemm M, Kobayashi R, Konkov V, Krishnamoorthy S, Krishnan M, Lin Z, Lins RD, Littlefield RJ, Logsdail AJ, Lopata K, Ma W, Marenich AV, Martin Del Campo J, Mejia-Rodriguez D, Moore JE, Mullin JM, Nakajima T, Nascimento DR, Nichols JA, Nichols PJ, Nieplocha J, Otero-de-la-Roza A, Palmer B, Panyala A, Pirojsirikul T, Peng B, Peverati R, Pittner J, Pollack L, Richard RM, Sadayappan P, Schatz GC, Shelton WA, Silverstein DW, Smith DMA, Soares TA, Song D, Swart M, Taylor HL, Thomas GS, Tipparaju V, Truhlar DG, Tsemekhman K, Van Voorhis T, Vázquez-Mayagoitia Á, Verma P, Villa O, Vishnu A, Vogiatzis KD, Wang D, Weare JH, Williamson MJ, Windus TL, Woliński K, Wong AT, Wu Q, Yang C, Yu Q, Zacharias M, Zhang Z, Zhao Y, Harrison RJ. NWChem: Past, present, and future. J Chem Phys 2020; 152:184102. [PMID: 32414274 DOI: 10.1063/5.0004997] [Citation(s) in RCA: 293] [Impact Index Per Article: 73.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.
Collapse
Affiliation(s)
- E Aprà
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - E J Bylaska
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - W A de Jong
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - N Govind
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - K Kowalski
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - T P Straatsma
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - M Valiev
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - H J J van Dam
- Brookhaven National Laboratory, Upton, New York 11973, USA
| | - Y Alexeev
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - J Anchell
- Intel Corporation, Santa Clara, California 95054, USA
| | - V Anisimov
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - F W Aquino
- QSimulate, Cambridge, Massachusetts 02139, USA
| | - R Atta-Fynn
- Department of Physics, The University of Texas at Arlington, Arlington, Texas 76019, USA
| | - J Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - N P Bauman
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - J C Becca
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - D E Bernholdt
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | | | - S Bogatko
- 4G Clinical, Wellesley, Massachusetts 02481, USA
| | - P Borowski
- Faculty of Chemistry, Maria Curie-Skłodowska University in Lublin, 20-031 Lublin, Poland
| | - J Boschen
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - J Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, 18223 Prague 8, Czech Republic
| | - A Bruner
- Department of Chemistry and Physics, University of Tennessee at Martin, Martin, Tennessee 38238, USA
| | - E Cauët
- Service de Chimie Quantique et Photophysique (CP 160/09), Université libre de Bruxelles, B-1050 Brussels, Belgium
| | - Y Chen
- Facebook, Menlo Park, California 94025, USA
| | - G N Chuev
- Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region 142290, Russia
| | - C J Cramer
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Daily
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M J O Deegan
- SKAO, Jodrell Bank Observatory, Macclesfield SK11 9DL, United Kingdom
| | - T H Dunning
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - M Dupuis
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - K G Dyall
- Dirac Solutions, Portland, Oregon 97229, USA
| | - G I Fann
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - S A Fischer
- Chemistry Division, U. S. Naval Research Laboratory, Washington, DC 20375, USA
| | - A Fonari
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - H Früchtl
- EaStCHEM and School of Chemistry, University of St. Andrews, St. Andrews KY16 9ST, United Kingdom
| | - L Gagliardi
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Garza
- Departamento de Química, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana-Iztapalapa, Col. Vicentina, Iztapalapa, C.P. 09340 Ciudad de México, Mexico
| | - N Gawande
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - S Ghosh
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 5545, USA
| | - K Glaesemann
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A W Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - J Hammond
- Intel Corporation, Santa Clara, California 95054, USA
| | - V Helms
- Center for Bioinformatics, Saarland University, D-66041 Saarbrücken, Germany
| | - E D Hermes
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, USA
| | - K Hirao
- Next-generation Molecular Theory Unit, Advanced Science Institute, RIKEN, Saitama 351-0198, Japan
| | - S Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - M Jacquelin
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - L Jensen
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - B G Johnson
- Acrobatiq, Pittsburgh, Pennsylvania 15206, USA
| | - H Jónsson
- Faculty of Physical Sciences, University of Iceland, Reykjavík, Iceland and Department of Applied Physics, Aalto University, FI-00076 Aalto, Espoo, Finland
| | - R A Kendall
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - M Klemm
- Intel Corporation, Santa Clara, California 95054, USA
| | - R Kobayashi
- ANU Supercomputer Facility, Australian National University, Canberra, Australia
| | - V Konkov
- Chemistry Program, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - S Krishnamoorthy
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M Krishnan
- Facebook, Menlo Park, California 94025, USA
| | - Z Lin
- Department of Physics, University of Science and Technology of China, Hefei, China
| | - R D Lins
- Aggeu Magalhaes Institute, Oswaldo Cruz Foundation, Recife, Brazil
| | | | - A J Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Cardiff, Wales CF10 3AT, United Kingdom
| | - K Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - W Ma
- Institute of Software, Chinese Academy of Sciences, Beijing, China
| | - A V Marenich
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Martin Del Campo
- Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México City, Mexico
| | - D Mejia-Rodriguez
- Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611, USA
| | - J E Moore
- Intel Corporation, Santa Clara, California 95054, USA
| | - J M Mullin
- DCI-Solutions, Aberdeen Proving Ground, Maryland 21005, USA
| | - T Nakajima
- Computational Molecular Science Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - D R Nascimento
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - J A Nichols
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - P J Nichols
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - J Nieplocha
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, 33006 Oviedo, Spain
| | - B Palmer
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A Panyala
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - T Pirojsirikul
- Department of Chemistry, Prince of Songkla University, Hat Yai, Songkhla 90112, Thailand
| | - B Peng
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - R Peverati
- Chemistry Program, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - J Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., 18223 Prague 8, Czech Republic
| | - L Pollack
- StudyPoint, Boston, Massachusetts 02114, USA
| | | | - P Sadayappan
- School of Computing, University of Utah, Salt Lake City, Utah 84112, USA
| | - G C Schatz
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - W A Shelton
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | | | - D M A Smith
- Intel Corporation, Santa Clara, California 95054, USA
| | - T A Soares
- Dept. of Fundamental Chemistry, Universidade Federal de Pernambuco, Recife, Brazil
| | - D Song
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M Swart
- ICREA, 08010 Barcelona, Spain and Universitat Girona, Institut de Química Computacional i Catàlisi, Campus Montilivi, 17003 Girona, Spain
| | - H L Taylor
- CD-adapco/Siemens, Melville, New York 11747, USA
| | - G S Thomas
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - V Tipparaju
- Cray Inc., Bloomington, Minnesota 55425, USA
| | - D G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | | | - T Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Á Vázquez-Mayagoitia
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - P Verma
- 1QBit, Vancouver, British Columbia V6E 4B1, Canada
| | - O Villa
- NVIDIA, Santa Clara, California 95051, USA
| | - A Vishnu
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - K D Vogiatzis
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - D Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong 250014, China
| | - J H Weare
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - M J Williamson
- Department of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - T L Windus
- Department of Chemistry, Iowa State University and Ames Laboratory, Ames, Iowa 50011, USA
| | - K Woliński
- Faculty of Chemistry, Maria Curie-Skłodowska University in Lublin, 20-031 Lublin, Poland
| | - A T Wong
- Qwil, San Francisco, California 94107, USA
| | - Q Wu
- Brookhaven National Laboratory, Upton, New York 11973, USA
| | - C Yang
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Q Yu
- AMD, Santa Clara, California 95054, USA
| | - M Zacharias
- Department of Physics, Technical University of Munich, 85748 Garching, Germany
| | - Z Zhang
- Stanford Research Computing Center, Stanford University, Stanford, California 94305, USA
| | - Y Zhao
- State Key Laboratory of Silicate Materials for Architectures, International School of Materials Science and Engineering, Wuhan University of Technology, Wuhan 430070, China
| | - R J Harrison
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
18
|
Perez SJLP, Claudio GC. Molecular dynamics simulations of two double-helical hexamer fragments of iota-carrageenan in aqueous solution. J Mol Graph Model 2020; 98:107588. [PMID: 32220758 DOI: 10.1016/j.jmgm.2020.107588] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Revised: 03/16/2020] [Accepted: 03/16/2020] [Indexed: 11/16/2022]
Abstract
The gelation of anionic carrageenans is known to occur through a coil-to-helix transition followed by further aggregation or association on which positive counterions play a significant role. In the present work, molecular dynamics (MD) simulations were performed on two double-helical iota-carrageenan hexamer fragments along with their sodium counterions using the Carbohydrate Solution Force Field (CSFF) in an aqueous (TIP4P) solution with the GROMACS molecular dynamics package. Results showed a counterion condensation between the two double helices and that the subsequent forces of interaction between them were predominantly attractive. By varying the distance separating the two double helices, the effect of distance on the counterion distribution and the forces of attraction was also investigated. In the presence of counterions, an increase in the forces of attraction was observed as the distance between the two double helices decreases which can be attributed to the greater counterion density between the two like-charged oligosaccharides.
Collapse
Affiliation(s)
- Ser John Lynon P Perez
- Institute of Chemistry, University of the Philippines, Diliman, Quezon City, Philippines.
| | - Gil C Claudio
- Institute of Chemistry, University of the Philippines, Diliman, Quezon City, Philippines
| |
Collapse
|
19
|
Cellular levels and molecular dynamics simulations of estragole DNA adducts point at inefficient repair resulting from limited distortion of the double-stranded DNA helix. Arch Toxicol 2020; 94:1349-1365. [PMID: 32185416 PMCID: PMC7225201 DOI: 10.1007/s00204-020-02695-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2019] [Accepted: 03/02/2020] [Indexed: 10/25/2022]
Abstract
Estragole, naturally occurring in a variety of herbs and spices, can form DNA adducts after bioactivation. Estragole DNA adduct formation and repair was studied in in vitro liver cell models, and a molecular dynamics simulation was used to investigate the conformation dependent (in)efficiency of N2-(trans-isoestragol-3'-yl)-2'-deoxyguanosine (E-3'-N2-dG) DNA adduct repair. HepG2, HepaRG cells, primary rat hepatocytes and CHO cells (including CHO wild-type and three NER-deficient mutants) were exposed to 50 μM estragole or 1'-hydroxyestragole and DNA adduct formation was quantified by LC-MS immediately following exposure and after a period of repair. Results obtained from CHO cell lines indicated that NER plays a role in repair of E-3'-N2-dG adducts, however, with limited efficiency since in the CHO wt cells 80% DNA adducts remained upon 24 h repair. Inefficiency of DNA repair was also found in HepaRG cells and primary rat hepatocytes. Changes in DNA structure resulting from E-3'-N2-dG adduct formation were investigated by molecular dynamics simulations. Results from molecular dynamics simulations revealed that conformational changes in double-stranded DNA by E-3'-N2-dG adduct formation are small, providing a possible explanation for the restrained repair, which may require larger distortions in the DNA structure. NER-mediated enzymatic repair of E-3'-N2-dG DNA adducts upon exposure to estragole will be limited, providing opportunities for accumulation of damage upon repeated daily exposure. The inability of this enzymatic repair is likely due to a limited distortion of the DNA double-stranded helix resulting in inefficient activation of nucleotide excision repair.
Collapse
|
20
|
Lin FY, MacKerell AD. Improved Modeling of Cation-π and Anion-Ring Interactions Using the Drude Polarizable Empirical Force Field for Proteins. J Comput Chem 2020; 41:439-448. [PMID: 31518010 PMCID: PMC7322827 DOI: 10.1002/jcc.26067] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Revised: 07/15/2019] [Accepted: 08/25/2019] [Indexed: 12/22/2022]
Abstract
Cation-π interactions are noncovalent interactions between a π-electron system and a positively charged ion that are regarded as a strong noncovalent interaction and are ubiquitous in biological systems. Similarly, though less studied, anion-ring interactions are present in proteins along with in-plane interactions of anions with aromatic rings. As these interactions are between a polarizing ion and a polarizable π system, the accuracy of the treatment of these interactions in molecular dynamics (MD) simulations using additive force fields (FFs) may be limited. In the present work, to allow for a better description of ion-π interactions in proteins in the Drude-2013 protein polarizable FF, we systematically optimized the parameters for these interactions targeting model compound quantum mechanical (QM) interaction energies with atom pair-specific Lennard-Jones parameters along with virtual particles as selected ring centroids introduced to target the QM interaction energies and geometries. Subsequently, MD simulations were performed on a series of protein structures where ion-π pairs occur to evaluate the optimized parameters in the context of the Drude-2013 FF. The resulting FF leads to a significant improvement in reproducing the ion-π pair distances observed in experimental protein structures, as well as a smaller root-mean-square differences and fluctuations of the overall protein structures from experimental structures. Accordingly, the optimized Drude-2013 protein polarizable FF is suggested for use in MD simulations of proteins where cation-π and anion-ring interactions are critical. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Fang-Yu Lin
- Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA
| | - Alexander D. MacKerell
- Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA
| |
Collapse
|
21
|
Lemkul JA. Pairwise-additive and polarizable atomistic force fields for molecular dynamics simulations of proteins. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:1-71. [PMID: 32145943 DOI: 10.1016/bs.pmbts.2019.12.009] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Protein force fields have been undergoing continual development since the first complete parameter sets were introduced nearly four decades ago. The functional forms that underlie these models have many common elements for the treatment of bonded and nonbonded forces, which are reviewed here. The most widely used force fields to date use a fixed-charge convention in which electronic polarization effects are treated via a mean-field approximation during partial charge assignment. Despite success in modeling folded proteins over many years, the fixed-charge assumption has limitations that cannot necessarily be overcome within their potential energy equations. To overcome these limitations, several force fields have recently been derived that explicitly treat electronic polarization effects with straightforward extensions of the potential energy functions used by nonpolarizable force fields. Here, we review the history of the most popular nonpolarizable force fields (AMBER, CHARMM, OPLS, and GROMOS) as well as studies that have validated them and applied them to studies of protein folding and misfolding. Building upon these force fields are more recent polarizable interaction potentials, including fluctuating charge models, POSSIM, AMOEBA, and the classical Drude oscillator. These force fields differ in their implementations but all attempt to model electronic polarization in a computationally tractable manner. Despite their recent emergence in the field of protein folding, several studies have already applied these polarizable models to challenging problems in this domain, including the role of polarization in folding free energies and sequence-specific effects on the stability of α-helical structures.
Collapse
Affiliation(s)
- Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, VA, United States.
| |
Collapse
|
22
|
Öhlknecht C, Lier B, Petrov D, Fuchs J, Oostenbrink C. Correcting electrostatic artifacts due to net-charge changes in the calculation of ligand binding free energies. J Comput Chem 2020; 41:986-999. [PMID: 31930547 DOI: 10.1002/jcc.26143] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Revised: 12/16/2019] [Accepted: 12/22/2019] [Indexed: 01/06/2023]
Abstract
Alchemically derived free energies are artifacted when the perturbed moiety has a nonzero net charge. The source of the artifacts lies in the effective treatment of the electrostatic interactions within and between the perturbed atoms and remaining (partial) charges in the simulated system. To treat the electrostatic interactions effectively, lattice-summation (LS) methods or cutoff schemes in combination with a reaction-field contribution are usually employed. Both methods render the charging component of the calculated free energies sensitive to essential parameters of the system like the cutoff radius or the box side lengths. Here, we discuss the results of three previously published studies of ligand binding. These studies presented estimates of binding free energies that were artifacted due to the charged nature of the ligands. We show that the size of the artifacts can be efficiently calculated and raw simulation data can be corrected. We compare the corrected results with experimental estimates and nonartifacted estimates from path-sampling methods. Although the employed correction scheme involves computationally demanding continuum-electrostatics calculations, we show that the correction estimate can be deduced from a small sample of configurations rather than from the entire ensemble. This observation makes the calculations of correction terms feasible for complex biological systems. To show the general applicability of the proposed procedure, we also present results where the correction scheme was used to correct independent free energies obtained from simulations employing a cutoff scheme or LS electrostatics. In this work, we give practical guidelines on how to apply the appropriate corrections easily.
Collapse
Affiliation(s)
- Christoph Öhlknecht
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria.,Austrian Centre of Industrial Biotechnology, Graz, Austria
| | - Bettina Lier
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| | - Drazen Petrov
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| | - Julian Fuchs
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria.,Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| |
Collapse
|
23
|
Diem M, Oostenbrink C. Hamiltonian Reweighing To Refine Protein Backbone Dihedral Angle Parameters in the GROMOS Force Field. J Chem Inf Model 2020; 60:279-288. [PMID: 31873012 DOI: 10.1021/acs.jcim.9b01034] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics simulations of proteins depend critically on the underlying force field, which may be parameterized against experimental data or high-quality quantum calculations. Here, we develop search algorithms based on Monte Carlo and steepest descent calculations to optimize the backbone dihedral angle parameters from a single reference simulation. We apply these tools to improve the agreement between simulations of single, capped amino acids and experimentally determined J values and secondary structure propensities of these molecules. The parameters are further refined based on simulations of a set of seven proteins and finally validated in simulations on a large set of 52 protein structures. Improvements in the dihedral angle distributions are observed, and structural propensities of the proteins are reproduced very well. Overall, the GROMOS 54A8_bb parameter set forms an improvement to previous parameter sets, both for small molecules and for protein simulations.
Collapse
Affiliation(s)
- Matthias Diem
- Institute for Molecular Modeling and Simulation , University of Natural Resources and Life Sciences , Muthgasse 18 , 1190 Vienna , Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation , University of Natural Resources and Life Sciences , Muthgasse 18 , 1190 Vienna , Austria
| |
Collapse
|
24
|
Sauceda HE, Chmiela S, Poltavsky I, Müller KR, Tkatchenko A. Construction of Machine Learned Force Fields with Quantum Chemical Accuracy: Applications and Chemical Insights. MACHINE LEARNING MEETS QUANTUM PHYSICS 2020. [DOI: 10.1007/978-3-030-40245-7_14] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
25
|
Bhatt JS. Solution structure of macromolecules using small angle neutron scattering and molecular simulations. EPJ WEB OF CONFERENCES 2020. [DOI: 10.1051/epjconf/202023603003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An introductory account of using molecular simulations to deduce solution structure of macromolecules using small angle neutron scattering data is presented for biologists. The presence of a liquid solution provides mobility to the molecules, making it difficult to pin down their structure. Here a simple introduction to molecular dynamics and Monte Carlo techniques is followed by a recipe to use the output of the simulations along with the scattering data in order to infer the structure of macromolecules when they are placed in a liquid solution. Some practical issues to be watched for are also highlighted.
Collapse
|
26
|
Sauceda HE, Chmiela S, Poltavsky I, Müller KR, Tkatchenko A. Molecular force fields with gradient-domain machine learning: Construction and application to dynamics of small molecules with coupled cluster forces. J Chem Phys 2019; 150:114102. [DOI: 10.1063/1.5078687] [Citation(s) in RCA: 56] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Huziel E. Sauceda
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
| | - Stefan Chmiela
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
| | - Igor Poltavsky
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Klaus-Robert Müller
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- Department of Brain and Cognitive Engineering, Korea University, Anam-dong, Seongbuk-gu, Seoul 02841, South Korea
- Max Planck Institute for Informatics, Stuhlsatzenhausweg, 66123 Saarbrücken, Germany
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| |
Collapse
|
27
|
Gonçalves YMH, Senac C, Fuchs PFJ, Hünenberger PH, Horta BAC. Influence of the Treatment of Nonbonded Interactions on the Thermodynamic and Transport Properties of Pure Liquids Calculated Using the 2016H66 Force Field. J Chem Theory Comput 2019; 15:1806-1826. [DOI: 10.1021/acs.jctc.8b00425] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Yan M. H. Gonçalves
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Caroline Senac
- Sorbonne
Université,
CNRS, INSERM, Laboratoire d’Imagerie Biomédicale (LIB), F-75006 Paris, France
| | - Patrick F. J. Fuchs
- Sorbonne Université,
École normale supérieure, PSL University, CNRS, Laboratoire
des biomolécules, LBM, 75005 Paris, France
- Université
Paris Diderot, 75013 Paris, France
| | | | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| |
Collapse
|
28
|
Goossens K, De Winter H. Molecular Dynamics Simulations of Membrane Proteins: An Overview. J Chem Inf Model 2018; 58:2193-2202. [PMID: 30336018 DOI: 10.1021/acs.jcim.8b00639] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Simulations of membrane proteins have been rising in popularity in the past decade. Advancements in technology and force fields made it possible to simulate behavior of membrane proteins. Membrane protein simulations can now be used as supporting evidence for experimental findings, for elucidating protein mechanisms, and validating protein crystal structures. Unrelated to experimental data, these simulations can also serve to investigate larger scale processes like protein sorting, protein-membrane interactions, and more. In this review, the history as well as the state-of-the-art methodologies in membrane protein simulations will be summarized. An emphasis will be put on how to set up the system and on the current models for the different components of the simulation system. An overview of the available tools for membrane protein simulation will be given, and current limitations and prospects will also be discussed.
Collapse
Affiliation(s)
- Kenneth Goossens
- Department of Pharmaceutical Sciences, Laboratory of Medicinal Chemistry , University of Antwerp , Universiteitsplein 1 , 2610 Wilrijk , Belgium
| | - Hans De Winter
- Department of Pharmaceutical Sciences, Laboratory of Medicinal Chemistry , University of Antwerp , Universiteitsplein 1 , 2610 Wilrijk , Belgium
| |
Collapse
|
29
|
Orr AA, Gonzalez-Rivera JC, Wilson M, Bhikha PR, Wang D, Contreras LM, Tamamis P. A high-throughput and rapid computational method for screening of RNA post-transcriptional modifications that can be recognized by target proteins. Methods 2018; 143:34-47. [DOI: 10.1016/j.ymeth.2018.01.015] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2017] [Revised: 01/14/2018] [Accepted: 01/26/2018] [Indexed: 12/25/2022] Open
|
30
|
Zhang C, Lu C, Jing Z, Wu C, Piquemal JP, Ponder JW, Ren P. AMOEBA Polarizable Atomic Multipole Force Field for Nucleic Acids. J Chem Theory Comput 2018; 14:2084-2108. [PMID: 29438622 PMCID: PMC5893433 DOI: 10.1021/acs.jctc.7b01169] [Citation(s) in RCA: 152] [Impact Index Per Article: 25.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, 3 J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.
Collapse
Affiliation(s)
- Changsheng Zhang
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Chao Lu
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Zhifeng Jing
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Chuanjie Wu
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC, UMR7616 CNRS, Paris, France
| | - Jay W. Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
31
|
Cacciotto P, Ramaswamy VK, Malloci G, Ruggerone P, Vargiu AV. Molecular Modeling of Multidrug Properties of Resistance Nodulation Division (RND) Transporters. Methods Mol Biol 2018; 1700:179-219. [PMID: 29177832 DOI: 10.1007/978-1-4939-7454-2_11] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Efflux pumps of the resistance nodulation division (RND) superfamily are among the major contributors to intrinsic and acquired multidrug resistance in Gram-negative bacteria. Structural information on AcrAB-TolC and MexAB-OprM, major efflux pumps of Escherichia coli and Pseudomonas aeruginosa respectively, boosted intensive research aimed at understanding the molecular mechanisms ruling the active extrusion processes. In particular, several studies were devoted to the understanding of the determinants behind the extraordinary broad specificity of the RND transporters AcrB and MexB. In this chapter, we discuss the ever-growing role computational methods have been playing in deciphering key structural and dynamical features of these transporters and of their interaction with substrates and inhibitors. We further discuss and illustrate examples from our lab of how molecular docking, homology modeling, all-atom molecular dynamics simulations and in silico free energy estimations can all together give precious insights into the processes of recognition and extrusion of substrates, as well as on the possible inhibition strategies.
Collapse
Affiliation(s)
- Pierpaolo Cacciotto
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Venkata K Ramaswamy
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Giuliano Malloci
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Paolo Ruggerone
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Attilio V Vargiu
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy.
| |
Collapse
|
32
|
Bakarić D, Petrov D, Mouvenchery YK, Heiβler S, Oostenbrink C, Schaumann GE. Ion-induced modification of the sucrose network and its impact on melting of freeze-dried liposomes. DSC and molecular dynamics study. Chem Phys Lipids 2017; 210:38-46. [PMID: 29179944 DOI: 10.1016/j.chemphyslip.2017.11.015] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2017] [Revised: 11/14/2017] [Accepted: 11/23/2017] [Indexed: 02/07/2023]
Abstract
Disaccharides play an important role in survival of anhydrobiotic organisms during extreme environmental conditions. A key protection feature is their capability to form the hydrogen bond (HB) network in a similar fashion as the one made by water. Since various ions also affect the HB network in completely hydrated systems, it is of a great interest to understand how they impact preservation when incorporated in a disaccharide network. To address this, we employ a combination of experimental and modeling techniques to study behavior of multilamellar 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) liposomes freeze-dried with sucrose in presence of NaCl or NaH2PO4·H2O at various concentrations (0.01-1M). Differential scanning calorimetry (DSC) was employed in order to determine the cooperative unit size (CUS), the number of lipid molecules that constitute a domain of cooperative motion in the liposome, and the melting temperature (Tm). In the absence of salt CUS was estimated to be 122±12, whereas in the presence of NaCl CUS increases more (347±34 for c=1M) than for NaH2PO4·H2O (193±26 for 1M). When it comes to Tm, the situation is reversed; NaCl induces increase by about 1K, while NaH2PO4·H2O by about 10K. These findings clearly demonstrate how different interaction forces-hydrogen bonding, charge pairing, and van der Waals interactions between acyl chains-affect CUS and Tm. Their interplay and contribution of particular interaction was further analyzed with molecular dynamics (MD) simulations. This analysis demonstrated that the HB network of DMPC and sucrose is partially disrupted in the presence of NaCl ions, and even to a greater extent in the case of NaH2PO4·H2O ions. Notably, H2PO4- ions outcompete and replace the sucrose molecules at the DMPC surface, which in turn alters the nature of the DMPC-surrounding interactions, from a weaker HB-dominated to a stronger CP-dominated interaction network.
Collapse
Affiliation(s)
- Danijela Bakarić
- University of Koblenz-Landau, Institute for Environmental Sciences, Group of Environmental and Soil Chemistry, Fortstraße 7, D-76829 Landau, Germany.
| | - Dražen Petrov
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences Vienna, Muthgasse 18, A-1190 Vienna, Austria
| | - Yamuna Kunhi Mouvenchery
- University of Koblenz-Landau, Institute for Environmental Sciences, Group of Environmental and Soil Chemistry, Fortstraße 7, D-76829 Landau, Germany
| | - Stefan Heiβler
- Institute for Functional Interfaces, Karlsruhe Institute for Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Chris Oostenbrink
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences Vienna, Muthgasse 18, A-1190 Vienna, Austria
| | - Gabriele E Schaumann
- University of Koblenz-Landau, Institute for Environmental Sciences, Group of Environmental and Soil Chemistry, Fortstraße 7, D-76829 Landau, Germany.
| |
Collapse
|
33
|
Markthaler D, Zeman J, Baz J, Smiatek J, Hansen N. Validation of Trimethylamine-N-oxide (TMAO) Force Fields Based on Thermophysical Properties of Aqueous TMAO Solutions. J Phys Chem B 2017; 121:10674-10688. [DOI: 10.1021/acs.jpcb.7b07774] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Daniel Markthaler
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Johannes Zeman
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Jörg Baz
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Jens Smiatek
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering and ‡Institute for Computational
Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| |
Collapse
|
34
|
Solvation free energy of solvation of biomass model cellobiose molecule: A molecular dynamics analysis. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.06.083] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
|
35
|
Danne R, Poojari C, Martinez-Seara H, Rissanen S, Lolicato F, Róg T, Vattulainen I. doGlycans-Tools for Preparing Carbohydrate Structures for Atomistic Simulations of Glycoproteins, Glycolipids, and Carbohydrate Polymers for GROMACS. J Chem Inf Model 2017; 57:2401-2406. [PMID: 28906114 PMCID: PMC5662928 DOI: 10.1021/acs.jcim.7b00237] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Carbohydrates constitute a structurally
and functionally diverse
group of biological molecules and macromolecules. In cells they are
involved in, e.g., energy storage, signaling, and cell–cell
recognition. All of these phenomena take place in atomistic scales,
thus atomistic simulation would be the method of choice to explore
how carbohydrates function. However, the progress in the field is
limited by the lack of appropriate tools for preparing carbohydrate
structures and related topology files for the simulation models. Here
we present tools that fill this gap. Applications where the tools
discussed in this paper are particularly useful include, among others,
the preparation of structures for glycolipids, nanocellulose, and
glycans linked to glycoproteins. The molecular structures and simulation
files generated by the tools are compatible with GROMACS.
Collapse
Affiliation(s)
- Reinis Danne
- Laboratory of Physics, Tampere University of Technology , P.O. Box 692, FI-33101 Tampere, Finland
| | - Chetan Poojari
- Laboratory of Physics, Tampere University of Technology , P.O. Box 692, FI-33101 Tampere, Finland.,Department of Physics, University of Helsinki , P.O. Box 64, FI-00014 Helsinki, Finland
| | - Hector Martinez-Seara
- Laboratory of Physics, Tampere University of Technology , P.O. Box 692, FI-33101 Tampere, Finland.,Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Sami Rissanen
- Laboratory of Physics, Tampere University of Technology , P.O. Box 692, FI-33101 Tampere, Finland
| | - Fabio Lolicato
- Laboratory of Physics, Tampere University of Technology , P.O. Box 692, FI-33101 Tampere, Finland.,Department of Physics, University of Helsinki , P.O. Box 64, FI-00014 Helsinki, Finland
| | - Tomasz Róg
- Laboratory of Physics, Tampere University of Technology , P.O. Box 692, FI-33101 Tampere, Finland.,Department of Physics, University of Helsinki , P.O. Box 64, FI-00014 Helsinki, Finland
| | - Ilpo Vattulainen
- Laboratory of Physics, Tampere University of Technology , P.O. Box 692, FI-33101 Tampere, Finland.,Department of Physics, University of Helsinki , P.O. Box 64, FI-00014 Helsinki, Finland.,MEMPHYS-Center for Biomembrane Physics, University of Southern Denmark , 5230 Odense, Denmark
| |
Collapse
|
36
|
Medeiros-Silva J, Jekhmane S, Baldus M, Weingarth M. Hydrogen bond strength in membrane proteins probed by time-resolved 1H-detected solid-state NMR and MD simulations. SOLID STATE NUCLEAR MAGNETIC RESONANCE 2017; 87:80-85. [PMID: 28342732 DOI: 10.1016/j.ssnmr.2017.03.003] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/30/2016] [Revised: 03/09/2017] [Accepted: 03/13/2017] [Indexed: 06/06/2023]
Abstract
1H-detected solid-state NMR in combination with 1H/2D exchange steps allows for the direct identification of very strong hydrogen bonds in membrane proteins. On the example of the membrane-embedded potassium channel KcsA, we quantify the longevity of such very strong hydrogen bonds by combining time-resolved 1H-detected solid-state NMR experiments and molecular dynamics simulations. In particular, we show that the carboxyl-side chain of the highly conserved residue Glu51 is involved in ultra-strong hydrogen bonds, which are fully-water-exposed and yet stable for weeks. The astonishing stability of these hydrogen bonds is important for the structural integrity of potassium channels, which we further corroborate by computational studies.
Collapse
Affiliation(s)
- João Medeiros-Silva
- NMR Spectroscopy, Bijvoet Center for Biomolecular Research, Department of Chemistry, Utrecht University, Pandualaan 8, 3584 CH Utrecht, The Netherlands
| | - Shehrazade Jekhmane
- NMR Spectroscopy, Bijvoet Center for Biomolecular Research, Department of Chemistry, Utrecht University, Pandualaan 8, 3584 CH Utrecht, The Netherlands
| | - Marc Baldus
- NMR Spectroscopy, Bijvoet Center for Biomolecular Research, Department of Chemistry, Utrecht University, Pandualaan 8, 3584 CH Utrecht, The Netherlands
| | - Markus Weingarth
- NMR Spectroscopy, Bijvoet Center for Biomolecular Research, Department of Chemistry, Utrecht University, Pandualaan 8, 3584 CH Utrecht, The Netherlands.
| |
Collapse
|
37
|
Santos DES, Pol-Fachin L, Lins RD, Soares TA. Polymyxin Binding to the Bacterial Outer Membrane Reveals Cation Displacement and Increasing Membrane Curvature in Susceptible but Not in Resistant Lipopolysaccharide Chemotypes. J Chem Inf Model 2017; 57:2181-2193. [DOI: 10.1021/acs.jcim.7b00271] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Denys E. S. Santos
- Department
of Fundamental Chemistry, Federal University of Pernambuco, 50740-560 Recife, Brazil
| | - Laércio Pol-Fachin
- Department
of Fundamental Chemistry, Federal University of Pernambuco, 50740-560 Recife, Brazil
- Aggeu Magalhães Institute, Oswaldo Cruz Foundation, 50740-465 Recife, Brazil
| | - Roberto D. Lins
- Aggeu Magalhães Institute, Oswaldo Cruz Foundation, 50740-465 Recife, Brazil
| | - Thereza A. Soares
- Department
of Fundamental Chemistry, Federal University of Pernambuco, 50740-560 Recife, Brazil
- Department
of Chemistry, Umeå Center for Microbial Research, Umeå University, 90.187 Umeå, Sweden
| |
Collapse
|
38
|
Lemkul JA, MacKerell AD. Polarizable Force Field for DNA Based on the Classical Drude Oscillator: II. Microsecond Molecular Dynamics Simulations of Duplex DNA. J Chem Theory Comput 2017; 13:2072-2085. [PMID: 28398748 PMCID: PMC5485260 DOI: 10.1021/acs.jctc.7b00068] [Citation(s) in RCA: 61] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The structure and dynamics of DNA are governed by a sensitive balance between base stacking and pairing, hydration, and interactions with ions. Force-field models that include explicit representations of electronic polarization are capable of more accurately modeling the subtle details of these interactions versus commonly used additive force fields. In this work, we validate our recently refined polarizable force field for DNA based on the classical Drude oscillator model, in which electronic degrees of freedom are represented as negatively charged particles attached to their parent atoms via harmonic springs. The previous version of the force field, called Drude-2013, produced stable A- and B-DNA trajectories on the order of hundreds of nanoseconds, but deficiencies were identified that included weak base stacking ultimately leading to distortion of B-DNA duplexes and unstable Z-DNA. As a result of extensive refinement of base nonbonded terms and bonded parameters in the deoxyribofuranose sugar and phosphodiester backbone, we demonstrate that the new version of the Drude DNA force field is capable of simulating A- and B-forms of DNA on the microsecond time scale and the resulting conformational ensembles agree well with a broad set of experimental properties, including solution X-ray scattering profiles. In addition, simulations of Z-form duplex DNA in its crystal environment are stable on the order of 100 ns. The revised force field is to be called Drude-2017.
Collapse
Affiliation(s)
- Justin A. Lemkul
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| |
Collapse
|
39
|
Lemkul JA, MacKerell AD. Polarizable Force Field for DNA Based on the Classical Drude Oscillator: I. Refinement Using Quantum Mechanical Base Stacking and Conformational Energetics. J Chem Theory Comput 2017; 13:2053-2071. [PMID: 28399366 DOI: 10.1021/acs.jctc.7b00067] [Citation(s) in RCA: 59] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Empirical force fields seek to relate the configuration of a set of atoms to its energy, thus yielding the forces governing its dynamics, using classical physics rather than more expensive quantum mechanical calculations that are computationally intractable for large systems. Most force fields used to simulate biomolecular systems use fixed atomic partial charges, neglecting the influence of electronic polarization, instead making use of a mean-field approximation that may not be transferable across environments. Recent hardware and software developments make polarizable simulations feasible, and to this end, polarizable force fields represent the next generation of molecular dynamics simulation technology. In this work, we describe the refinement of a polarizable force field for DNA based on the classical Drude oscillator model by targeting quantum mechanical interaction energies and conformational energy profiles of model compounds necessary to build a complete DNA force field. The parametrization strategy employed in the present work seeks to correct weak base stacking in A- and B-DNA and the unwinding of Z-DNA observed in the previous version of the force field, called Drude-2013. Refinement of base nonbonded terms and reparametrization of dihedral terms in the glycosidic linkage, deoxyribofuranose rings, and important backbone torsions resulted in improved agreement with quantum mechanical potential energy surfaces. Notably, we expand on previous efforts by explicitly including Z-DNA conformational energetics in the refinement.
Collapse
Affiliation(s)
- Justin A Lemkul
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| |
Collapse
|
40
|
Jedrzejczyk D, Gendaszewska-Darmach E, Pawlowska R, Chworos A. Designing synthetic RNA for delivery by nanoparticles. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:123001. [PMID: 28004640 DOI: 10.1088/1361-648x/aa5561] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The rapid development of synthetic biology and nanobiotechnology has led to the construction of various synthetic RNA nanoparticles of different functionalities and potential applications. As they occur naturally, nucleic acids are an attractive construction material for biocompatible nanoscaffold and nanomachine design. In this review, we provide an overview of the types of RNA and nucleic acid's nanoparticle design, with the focus on relevant nanostructures utilized for gene-expression regulation in cellular models. Structural analysis and modeling is addressed along with the tools available for RNA structural prediction. The functionalization of RNA-based nanoparticles leading to prospective applications of such constructs in potential therapies is shown. The route from the nanoparticle design and modeling through synthesis and functionalization to cellular application is also described. For a better understanding of the fate of targeted RNA after delivery, an overview of RNA processing inside the cell is also provided.
Collapse
Affiliation(s)
- Dominika Jedrzejczyk
- Centre of Molecular and Macromolecular Studies, Polish Academy of Sciences, Sienkiewicza 112, 90-363 Lodz, Poland
| | | | | | | |
Collapse
|
41
|
Pagliai M, Mancini G, Carnimeo I, De Mitri N, Barone V. Electronic absorption spectra of pyridine and nicotine in aqueous solution with a combined molecular dynamics and polarizable QM/MM approach. J Comput Chem 2017; 38:319-335. [PMID: 27910109 PMCID: PMC6680224 DOI: 10.1002/jcc.24683] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2016] [Revised: 11/04/2016] [Accepted: 11/08/2016] [Indexed: 01/06/2023]
Abstract
The electronic absorption spectra of pyridine and nicotine in aqueous solution have been computed using a multistep approach. The computational protocol consists in studying the solute solvation with accurate molecular dynamics simulations, characterizing the hydrogen bond interactions, and calculating electronic transitions for a series of configurations extracted from the molecular dynamics trajectories with a polarizable QM/MM scheme based on the fluctuating charge model. Molecular dynamics simulations and electronic transition calculations have been performed on both pyridine and nicotine. Furthermore, the contributions of solute vibrational effect on electronic absorption spectra have been taken into account in the so called vertical gradient approximation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Marco Pagliai
- Scuola Normale SuperiorePiazza dei Cavalieri 7PisaI‐56126Italy
| | | | - Ivan Carnimeo
- Scuola Normale SuperiorePiazza dei Cavalieri 7PisaI‐56126Italy
- Compunet, Istituto Italiano di Tecnologia (IIT)via Morego 30GenovaI‐16163Italy
| | - Nicola De Mitri
- Scuola Normale SuperiorePiazza dei Cavalieri 7PisaI‐56126Italy
- Present address:
Department of ChemistryUniversity of CambridgeLensfield RoadCambridgeCB2 1EWU.K.
| | - Vincenzo Barone
- Scuola Normale SuperiorePiazza dei Cavalieri 7PisaI‐56126Italy
| |
Collapse
|
42
|
Abstract
Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.
Collapse
Affiliation(s)
| | - Kenneth M. Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, and Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
43
|
Margreitter C, Reif MM, Oostenbrink C. Update on phosphate and charged post-translationally modified amino acid parameters in the GROMOS force field. J Comput Chem 2017; 38:714-720. [DOI: 10.1002/jcc.24733] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2016] [Revised: 12/28/2016] [Accepted: 12/29/2016] [Indexed: 11/10/2022]
Affiliation(s)
- Christian Margreitter
- Department of Material Sciences and Process Engineering; Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences; Vienna Austria
| | - Maria M. Reif
- Physics Department T38; Technical University of Munich; Munich Germany
| | - Chris Oostenbrink
- Department of Material Sciences and Process Engineering; Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences; Vienna Austria
| |
Collapse
|
44
|
Molecular dynamics simulation of the sliding of distamycin anticancer drug along DNA: interactions and sequence selectivity. JOURNAL OF THE IRANIAN CHEMICAL SOCIETY 2016. [DOI: 10.1007/s13738-016-1001-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
45
|
Vermaas JV, Trebesch N, Mayne CG, Thangapandian S, Shekhar M, Mahinthichaichan P, Baylon JL, Jiang T, Wang Y, Muller MP, Shinn E, Zhao Z, Wen PC, Tajkhorshid E. Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation. Methods Enzymol 2016; 578:373-428. [PMID: 27497175 PMCID: PMC6404235 DOI: 10.1016/bs.mie.2016.05.042] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.
Collapse
Affiliation(s)
- J V Vermaas
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - N Trebesch
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - C G Mayne
- University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - S Thangapandian
- University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - M Shekhar
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - P Mahinthichaichan
- University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - J L Baylon
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - T Jiang
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - Y Wang
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - M P Muller
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; College of Medicine, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - E Shinn
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - Z Zhao
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - P-C Wen
- University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
| | - E Tajkhorshid
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; University of Illinois at Urbana-Champaign, Urbana, IL, United States; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States; College of Medicine, University of Illinois at Urbana-Champaign, Urbana, IL, United States.
| |
Collapse
|
46
|
Horta BAC, Merz PT, Fuchs PFJ, Dolenc J, Riniker S, Hünenberger PH. A GROMOS-Compatible Force Field for Small Organic Molecules in the Condensed Phase: The 2016H66 Parameter Set. J Chem Theory Comput 2016; 12:3825-50. [DOI: 10.1021/acs.jctc.6b00187] [Citation(s) in RCA: 88] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bruno A. C. Horta
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Pascal T. Merz
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Patrick F. J. Fuchs
- Institut Jacques Monod, UMR 7592 CNRS, Université Paris-Diderot, Sorbonne Paris Cité, F-75205 Paris, France
| | - Jozica Dolenc
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Chemistry,
Biology and Pharmacy Information Center, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | | |
Collapse
|
47
|
Galindo-Murillo R, Robertson JC, Zgarbová M, Šponer J, Otyepka M, Jurečka P, Cheatham TE. Assessing the Current State of Amber Force Field Modifications for DNA. J Chem Theory Comput 2016; 12:4114-27. [PMID: 27300587 PMCID: PMC4980684 DOI: 10.1021/acs.jctc.6b00186] [Citation(s) in RCA: 308] [Impact Index Per Article: 38.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
![]()
The utility of molecular
dynamics (MD) simulations to model biomolecular
structure, dynamics, and interactions has witnessed enormous advances
in recent years due to the availability of optimized MD software and
access to significant computational power, including GPU multicore
computing engines and other specialized hardware. This has led researchers
to routinely extend conformational sampling times to the microsecond
level and beyond. The extended sampling time has allowed the community
not only to converge conformational ensembles through complete sampling
but also to discover deficiencies and overcome problems with the force
fields. Accuracy of the force fields is a key component, along with
sampling, toward being able to generate accurate and stable structures
of biopolymers. The Amber force field for nucleic acids has been used
extensively since the 1990s, and multiple artifacts have been discovered,
corrected, and reassessed by different research groups. We present
a direct comparison of two of the most recent and state-of-the-art
Amber force field modifications, bsc1 and OL15, that focus on accurate
modeling of double-stranded DNA. After extensive MD simulations with
five test cases and two different water models, we conclude that both
modifications are a remarkable improvement over the previous bsc0
force field. Both force field modifications show better agreement
when compared to experimental structures. To ensure convergence, the
Drew–Dickerson dodecamer (DDD) system was simulated using 100
independent MD simulations, each extended to at least 10 μs,
and the independent MD simulations were concatenated into a single
1 ms long trajectory for each combination of force field and water
model. This is significantly beyond the time scale needed to converge
the conformational ensemble of the internal portions of a DNA helix
absent internal base pair opening. Considering all of the simulations
discussed in the current work, the MD simulations performed to assess
and validate the current force fields and water models aggregate over
14 ms of simulation time. The results suggest that both the bsc1 and
OL15 force fields render average structures that deviate significantly
less than 1 Å from the average experimental structures. This
can be compared to similar but less exhaustive simulations with the
CHARMM 36 force field that aggregate to the ∼90 μs time
scale and also perform well but do not produce structures as close
to the DDD NMR average structures (with root-mean-square deviations
of 1.3 Å) as the newer Amber force fields. On the basis of these
analyses, any future research involving double-stranded DNA simulations
using the Amber force fields should employ the bsc1 or OL15 modification.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, University of Utah , 2000 East 30 South, Skaggs 105, Salt Lake City, Utah 84112, United States
| | - James C Robertson
- Department of Medicinal Chemistry, University of Utah , 2000 East 30 South, Skaggs 105, Salt Lake City, Utah 84112, United States
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17 Listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17 Listopadu 12, 771 46 Olomouc, Czech Republic.,Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17 Listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17 Listopadu 12, 771 46 Olomouc, Czech Republic
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, University of Utah , 2000 East 30 South, Skaggs 105, Salt Lake City, Utah 84112, United States
| |
Collapse
|
48
|
Cunha KC, Rusu VH, Viana IFT, Marques ETA, Dhalia R, Lins RD. Assessing protein conformational sampling and structural stability via de novo design and molecular dynamics simulations. Biopolymers 2016; 103:351-61. [PMID: 25677872 DOI: 10.1002/bip.22626] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2014] [Revised: 01/31/2015] [Accepted: 02/03/2015] [Indexed: 11/10/2022]
Abstract
Molecular dynamics and de novo techniques, associated to quality parameter sets, have excelled at determining the structure of small proteins with high accuracy. To achieve a detailed description of protein conformations, these methods must critically assess the thermodynamic features of the molecular ensembles. Here, a comparison of the conformational ensemble generated by molecular dynamics and de novo techniques were carried out for six Top7-based proteins carrying gp41 HIV-1 epitopes. The native Top7, a highly stable computationally designed protein, was used as benchmark. Structural stability, flexibility, and secondary structure content were assessed. The consistency of the latter was compared to experimental circular dichroism spectra for all proteins. While both methods are capable to identify the stable from unstable chimeric proteins, the sampled conformational space and flexibility differ significantly in both methods. Molecular dynamics simulations seem to better describe secondary structure content and identify regions responsible for conformational instability. The de novo method, as implemented in Rosetta-a prime tool for protein design, overestimates secondary structure content. On the other hand, its empirical energy function is capable to predict the threshold for protein stability.
Collapse
Affiliation(s)
- Keila C Cunha
- Department of Fundamental Chemistry, Federal University of Pernambuco, Recife, PE, 50740-560, Brazil
| | | | | | | | | | | |
Collapse
|
49
|
Lemkul J, Huang J, Roux B, MacKerell AD. An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications. Chem Rev 2016; 116:4983-5013. [PMID: 26815602 PMCID: PMC4865892 DOI: 10.1021/acs.chemrev.5b00505] [Citation(s) in RCA: 371] [Impact Index Per Article: 46.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2015] [Indexed: 11/28/2022]
Abstract
Molecular mechanics force fields that explicitly account for induced polarization represent the next generation of physical models for molecular dynamics simulations. Several methods exist for modeling induced polarization, and here we review the classical Drude oscillator model, in which electronic degrees of freedom are modeled by charged particles attached to the nuclei of their core atoms by harmonic springs. We describe the latest developments in Drude force field parametrization and application, primarily in the last 15 years. Emphasis is placed on the Drude-2013 polarizable force field for proteins, DNA, lipids, and carbohydrates. We discuss its parametrization protocol, development history, and recent simulations of biologically interesting systems, highlighting specific studies in which induced polarization plays a critical role in reproducing experimental observables and understanding physical behavior. As the Drude oscillator model is computationally tractable and available in a wide range of simulation packages, it is anticipated that use of these more complex physical models will lead to new and important discoveries of the physical forces driving a range of chemical and biological phenomena.
Collapse
Affiliation(s)
- Justin
A. Lemkul
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Baltimore, Maryland 21201, United States
| | - Jing Huang
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Baltimore, Maryland 21201, United States
| | - Benoît Roux
- Department
of Biochemistry and Molecular Biology, University
of Chicago, Chicago, Illinois 60637, United
States
| | - Alexander D. MacKerell
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Baltimore, Maryland 21201, United States
| |
Collapse
|
50
|
Validation of polarizable force field parameters for nucleic acids by inter-molecular interactions. Front Chem Sci Eng 2016. [DOI: 10.1007/s11705-016-1572-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|