Mutually Beneficial Combination of Molecular Dynamics Computer Simulations and Scattering Experiments.
Ontology highlight
ABSTRACT:
SUBMITTER: Zec N
PROVIDER: S-EPMC8304056 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The combination of molecular dynamics (MD) simulations and small-angle X-ray scattering (SAXS), called the MD-SAXS method, is efficient for investigating protein dynamics. To overcome the time-scale limitation of all-atom MD simulations, coarse-grained (CG) representations are often utilized for biomolecular simulations. In this study, we propose a method to combine CG MD simulations with SAXS, termed the CG-MD-SAXS method. In the CG-MD-SAXS method, the scattering factors of CG particles for proteins and nucleic acids are evaluated using high-resolution structural data in the Protein Data Bank, and the excluded volume and the hydration shell are modeled using two adjustable parameters to incorporate solvent effects. To avoid overfitting, only the two parameters are adjusted for an entire structure ensemble. To verify the developed method, theoretical SAXS profiles for various proteins, DNA/RNA, and a protein-RNA complex are compared with both experimental profiles and theoretical profiles obtained by the all-atom representation. In the present study, we applied the CG-MD-SAXS method to the Swi5-Sfr1 complex and three types of nucleosomes to obtain reliable ensemble models consistent with the experimental SAXS data.
| S-EPMC6976007 | biostudies-literature
Project description:Wide-angle x-ray scattering (WAXS) is emerging as a powerful tool for increasing the resolution of solution structure measurements of biomolecules. Compared to its better known complement, small angle x-ray scattering (SAXS), WAXS targets higher scattering angles and can enhance structural studies of molecules by accessing finer details of solution structures. Although the extension from SAXS to WAXS is easy to implement experimentally, the computational tools required to fully harness the power of WAXS are still under development. Currently, WAXS is employed to study structural changes and ligand binding in proteins; however, the methods are not as fully developed for nucleic acids. Here, we show how WAXS can qualitatively characterize nucleic acid structures as well as the small but significant structural changes driven by the addition of multivalent ions. We show the potential of WAXS to test all-atom molecular dynamics (MD) simulations and to provide insight into understanding how the trivalent ion cobalt(III) hexammine (CoHex) affects the structure of RNA and DNA helices. We find that MD simulations capture the RNA structural change that occurs due to addition of CoHex.
| S-EPMC4884185 | biostudies-literature
Project description:The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.
| S-EPMC4183268 | biostudies-literature
Project description:Time-resolved X-ray solution scattering is an increasingly popular method to measure conformational changes in proteins. Extracting structural information from the resulting difference X-ray scattering data is a daunting task. We present a method in which the limited but precious information encoded in such scattering curves is combined with the chemical knowledge of molecular force fields. The molecule of interest is then refined toward experimental data using molecular dynamics simulation. Therefore, the energy landscape is biased toward conformations that agree with experimental data. We describe and verify the method, and we provide an implementation in GROMACS.
| S-EPMC4325560 | biostudies-literature
Project description:Gramicidin A (gA) channels provide an ideal system to test molecular dynamics (MD) simulations of membrane proteins. The peptide backbone lines a cation-selective pore, and due to the small channel size, the average structure and extent of fluctuations of all atoms in the peptide will influence ion permeation. This raises the question of how well molecular mechanical force fields used in MD simulations and potential of mean force (PMF) calculations can predict structure and dynamics as well as ion permeation. To address this question, we undertook a comparative study of nuclear magnetic resonance (NMR) observables predicted by fully atomistic MD simulations on a gA dimer embedded in a sodium dodecyl sulfate (SDS) micelle with measurements of the gA dimer backbone and tryptophan side chain dynamics using solution-state (15)N NMR on gA dimers in SDS micelles (Vostrikov, V. V.; Gu, H.; Ingólfsson, H. I.; Hinton, J. F.; Andersen, O. S.; Roux, B.; Koeppe, R. E., II. J. Phys. Chem. B2011, DOI 10.1021/jp200906y , accompanying article). This comparison enables us to examine the robustness of the MD simulations done using different force fields as well as their ability to predict important features of the gA channel. We find that MD is able to predict NMR observables, including the generalized order parameters (S(2)), the (15)N spin-lattice (T(1)) and spin-spin (T(2)) relaxation times, and the (1)H-(15)N nuclear Overhauser effect (NOE), with remarkable accuracy. To examine further how differences in the force fields can affect the channel conductance, we calculated the PMF for K(+) and Na(+) permeation through a gA channel in a dimyristoylphosphatidylcholine (DMPC) bilayer. In this case, we find that MD is less successful in quantitatively predicting the single-channel conductance.
| S-EPMC3107394 | biostudies-literature
Project description:Neutron spin-echo spectroscopy is a unique experimental method for the investigation of polymer dynamics. The combination of neutron spin-echo spectroscopy with grazing-incidence geometry (GINSES) opens the possibility to probe the dynamics of soft-matter materials in the vicinity of the solid substrate in the time range up to 100 ns. However, the usage of the GINSES technique has some peculiarities and, due to the novelty of the method and complexity of the scattering geometry, difficulties in further data analysis occur. The current work discusses how virtual experiments within the distorted-wave Born approximation using the BornAgain software can improve GINSES data treatment and aid the understanding of polymer dynamics in the vicinity of the solid surface. With two examples, poly(N-isopropyl acrylamide) brushes and poly(ethylene glycol) microgels on Si surfaces, the simulation as well as the application of the simulation to the GINSES data analysis are presented. The approach allowed a deeper insight to be gained of the background effect and scattering contribution of different layers.
| S-EPMC7941298 | biostudies-literature
Project description:The removal of the heme group from myoglobin (Mb) results in a destabilization of the protein structure. The dynamic basis of the destabilization was followed by comparative measurements on holo- (holo-Mb) and apomyoglobin (apo-Mb). Mean-squared displacements (MSD) and protein resilience on the picosecond-to-nanosecond timescale were measured by elastic incoherent neutron scattering. Differences in thermodynamic parameters, MSD, and resilience were observed for both proteins. The resilience of holo-Mb was significantly lower than that of apo-Mb, indicating entropic stabilization by a higher degree of conformational sampling in the heme-bound folded protein. Molecular dynamics simulations provided site-specific information. Averaged over the whole structure, the molecular dynamics simulations yielded similar MSD and resilience values for the two proteins. The mobility of residues around the heme group in holo-Mb showed a smaller MSD and higher resilience compared to the same residue group in apo-Mb. It is of interest that in holo-Mb, higher MSD values are observed for the residues outside the heme pocket, indicating an entropic contribution to protein stabilization by heme binding, which is in agreement with experimental results.
| S-EPMC3260689 | biostudies-literature
Project description:ORF8a protein is 39 residues long and contains a single transmembrane domain. The protein is synthesized using solid phase peptide synthesis and reconstituted into artificial lipid bilayers that forms cation-selective ion channels with a main conductance level of 8.9±0.8pS at elevated temperature (38.5°C). Computational modeling studies including multi nanosecond molecular dynamics simulations in a hydrated POPC lipid bilayer are done with a 22 amino acid transmembrane helix to predict a putative homooligomeric helical bundle model. A structural model of a pentameric bundle is proposed with cysteines, serines and threonines facing the pore.
| S-EPMC7094593 | biostudies-literature
Project description:Cholesterol and ether lipids are ubiquitous in mammalian cell membranes, and their interactions are crucial in ether lipid mediated cholesterol trafficking. We report on cholesterol's molecular interactions with ether lipids as determined using a combination of small-angle neutron and X-ray scattering, and all-atom molecular dynamics (MD) simulations. A scattering density profile model for an ether lipid bilayer was developed using MD simulations, which was then used to simultaneously fit the different experimental scattering data. From analysis of the data the various bilayer structural parameters were obtained. Surface area constrained MD simulations were also performed to reproduce the experimental data. This iterative analysis approach resulted in good agreement between the experimental and simulated form factors. The molecular interactions taking place between cholesterol and ether lipids were then determined from the validated MD simulations. We found that in ether membranes cholesterol primarily hydrogen bonds with the lipid headgroup phosphate oxygen, while in their ester membrane counterparts cholesterol hydrogen bonds with the backbone ester carbonyls. This different mode of interaction between ether lipids and cholesterol induces cholesterol to reside closer to the bilayer surface, dehydrating the headgroup's phosphate moiety. Moreover, the three-dimensional lipid chain spatial density distribution around cholesterol indicates anisotropic chain packing, causing cholesterol to tilt. These insights lend a better understanding of ether lipid-mediated cholesterol trafficking and the roles that the different lipid species have in determining the structural and dynamical properties of membrane associated biomolecules.
| S-EPMC3539752 | biostudies-literature
Project description:Molecular-dynamics (MD) simulations of Bragg and diffuse X-ray scattering provide a means of obtaining experimentally validated models of protein conformational ensembles. This paper shows that compared with a single periodic unit-cell model, the accuracy of simulating diffuse scattering is increased when the crystal is modeled as a periodic supercell consisting of a 2 × 2 × 2 layout of eight unit cells. The MD simulations capture the general dependence of correlations on the separation of atoms. There is substantial agreement between the simulated Bragg reflections and the crystal structure; there are local deviations, however, indicating both the limitation of using a single structure to model disordered regions of the protein and local deviations of the average structure away from the crystal structure. Although it was anticipated that a simulation of longer duration might be required to achieve maximal agreement of the diffuse scattering calculation with the data using the supercell model, only a microsecond is required, the same as for the unit cell. Rigid protein motions only account for a minority fraction of the variation in atom positions from the simulation. The results indicate that protein crystal dynamics may be dominated by internal motions rather than packing interactions, and that MD simulations can be combined with Bragg and diffuse X-ray scattering to model the protein conformational ensemble.
| S-EPMC5947722 | biostudies-literature