Relaxation estimation of RMSD in molecular dynamics immunosimulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Schreiner W
PROVIDER: S-EPMC3457668 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:When photoactive molecules interact strongly with confined light modes in optical cavities, new hybrid light-matter states form. They are known as polaritons and correspond to coherent superpositions of excitations of the molecules and of the cavity photon. The polariton energies and thus potential energy surfaces are changed with respect to the bare molecules, such that polariton formation is considered a promising paradigm for controlling photochemical reactions. To effectively manipulate photochemistry with confined light, the molecules need to remain in the polaritonic state long enough for the reaction on the modified potential energy surface to take place. To understand what determines this lifetime, we have performed atomistic molecular dynamics simulations of room-temperature ensembles of rhodamine chromophores strongly coupled to a single confined light mode with a 15 fs lifetime. We investigated three popular experimental scenarios and followed the relaxation after optically pumping (i) the lower polariton, (ii) the upper polariton, or (iii) uncoupled molecular states. The results of the simulations suggest that the lifetimes of the optically accessible lower and upper polaritons are limited by (i) ultrafast photoemission due to the low cavity lifetime and (ii) reversible population transfer into the "dark" state manifold. Dark states are superpositions of molecular excitations but with much smaller contributions from the cavity photon, decreasing their emission rates and hence increasing their lifetimes. We find that population transfer between polaritonic modes and dark states is determined by the overlap between the polaritonic and molecular absorption spectra. Importantly, excitation can also be transferred "upward" from the lower polariton into the dark-state reservoir due to the broad absorption spectra of the chromophores, contrary to the common conception of these processes as a "one-way" relaxation from the dark states down to the lower polariton. Our results thus suggest that polaritonic chemistry relying on modified dynamics taking place within the lower polariton manifold requires cavities with sufficiently long lifetimes and, at the same time, strong light-matter coupling strengths to prevent the back-transfer of excitation into the dark states.
| S-EPMC6914212 | biostudies-literature
Project description:Conformational fluctuations and rotational tumbling of proteins can be experimentally accessed with nuclear spin relaxation experiments. However, interpretation of molecular dynamics from the experimental data is often complicated, especially for molecules with anisotropic shape. Here, we apply classical molecular dynamics simulations to interpret the conformational fluctuations and rotational tumbling of proteins with arbitrarily anisotropic shape. The direct calculation of spin relaxation times from simulation data did not reproduce the experimental data. This was successfully corrected by scaling the overall rotational diffusion coefficients around the protein inertia axes with a constant factor. The achieved good agreement with experiments allowed the interpretation of the internal and overall dynamics of proteins with significantly anisotropic shape. The overall rotational diffusion was found to be Brownian, having only a short subdiffusive region below 0.12 ns. The presented methodology can be applied to interpret rotational dynamics and conformation fluctuations of proteins with arbitrary anisotropic shape. However, a water model with more realistic dynamical properties is probably required for intrinsically disordered proteins.
| S-EPMC6150695 | biostudies-literature
Project description:The estimation of
| S-EPMC9301944 | biostudies-literature
Project description:The kinetic rate constants of binding were estimated for four biochemically relevant molecular systems by a method that uses milestoning theory to combine Brownian dynamics simulations with more detailed molecular dynamics simulations. The rate constants found using this method agreed well with experimentally and theoretically obtained values. We predicted the association rate of a small charged molecule toward both a charged and an uncharged spherical receptor and verified the estimated value with Smoluchowski theory. We also calculated the kon rate constant for superoxide dismutase with its natural substrate, O2-, in a validation of a previous experiment using similar methods but with a number of important improvements. We also calculated the kon for a new system: the N-terminal domain of Troponin C with its natural substrate Ca2+. The kon calculated for the latter two systems closely resemble experimentally obtained values. This novel multiscale approach is computationally cheaper and more parallelizable when compared to other methods of similar accuracy. We anticipate that this methodology will be useful for predicting kinetic rate constants and for understanding the process of binding between a small molecule and a protein receptor.
| S-EPMC4624728 | biostudies-literature
Project description:Analyzing nonlinear conformational relaxation dynamics in elastic networks corresponding to two classical motor proteins, we find that they respond by well defined internal mechanical motions to various initial deformations and that these motions are robust against external perturbations. We show that this behavior is not characteristic for random elastic networks. However, special network architectures with such properties can be designed by evolutionary optimization methods. Using them, an example of an artificial elastic network, operating as a cyclic machine powered by ligand binding, is constructed.
| S-EPMC1868896 | biostudies-literature
Project description:Amyloid-β (Aβ) aggregation is a hallmark of Alzheimer's disease. As an intrinsically disordered protein, Aβ undergoes extensive dynamics on multiple length and time scales. Access to a comprehensive picture of the reorientational dynamics in Aβ requires therefore the combination of complementary techniques. Here, we integrate 15N spin relaxation rates at three magnetic fields with microseconds-long molecular dynamics simulation, ensemble-based hydrodynamic calculations, and previously published nanosecond fluorescence correlation spectroscopy to investigate the reorientational dynamics of Aβ1-40 (Aβ40) at single-residue resolution. The integrative analysis shows that librational and dihedral angle fluctuations occurring at fast and intermediate time scales are not sufficient to decorrelate orientational memory in Aβ40. Instead, slow segmental motions occurring at ∼5 ns are detected throughout the Aβ40 sequence and reach up to ∼10 ns for selected residues. We propose that the modulation of time scales of reorientational dynamics with respect to intra- and intermolecular diffusion plays an important role in disease-related Aβ aggregation.
| S-EPMC6598774 | biostudies-literature
Project description:Interactions between hydrogen protons of water molecules and macromolecules within the myelin sheath surrounding the axons are a major factor influencing the magnetic resonance (MR) contrast in white matter (WM) regions. In past decades, several studies have investigated the underlying effects and reported a wide range of R1 rates for the myelin associated compartments at different field strengths. However, it was also shown that the experimental quantification of the compartment-specific R1 rates is associated with large uncertainties. The current study therefore investigates the longitudinal relaxation rates within the myelin sheath using a molecular dynamic (MD) simulation. For this purpose, a realistic molecular model of the myelin sheath was employed to determine the dipole-dipole induced R1 relaxation rate of the hydrogen protons at clinically relevant field strengths. The results obtained clearly reflect the spatial heterogeneity of R1 with a increased relaxivity of myelin water due to a reduced molecular mobility near the membrane surface. Moreover, the calculated R1 rates for both myelin water and macromolecules are in excellent agreement with experimental findings from the literature at different field strengths.
| S-EPMC6794311 | biostudies-literature
Project description:Molecular recognition by enzymes is a complicated process involving thermodynamic energies governing protein-ligand interactions. In order to aid the estimation of inhibitory activity of compounds targeting an enzyme, several computational methods can be employed to dissect this intermolecular contact. Herein, we report a structural dynamics investigation of an epigenetic enzyme HDAC2 in differentiating its binding to various inhibitors within the sub-sites of its active site. Molecular dynamics (MD) simulation was employed to elucidate the intermolecular interactions as well as the dynamics behavior of ligand binding. MD trajectories of five distinct HDAC2-inhibitor complexes reveal that compounds lacking adequate contacts with the opening rim of the active site possess high fluctuation along the cap portion, thus weakening the overall affinity. Key intermolecular interactions determining the effective binding of inhibitors include hydrogen bonds with Gly154, Asp181, and Tyr308; hydrophobic interactions between Phe155/Phe210 and the linker region; and a pi-stacking with Arg39 at the foot pocket. Decomposition of the binding free energy calculated per-residue by MM/PBSA also indicates that the interactions within the internal foot pocket, especially with residues Met35, Leu144, Gly305, and Gly306, can contribute significantly to the ligand binding. Additionally, configurational entropy of the binding was estimated and compared to the scale of the binding free energy in order to assess its contribution to the binding and to differentiate various ligand partners. It was found that the levels of entropic contribution are comparable among a set of structurally similar carbamide ligands, while it is greatly different for the set of unrelated ligands, ranging from 2.75 to 16.38 kcal/mol for the five inhibitors examined. These findings exemplify the importance of assessing molecular dynamics as well as estimating the entropic contribution in evaluating the ligand binding mechanism.
| S-EPMC9387782 | biostudies-literature
Project description:We study the dynamics of semiflexible hyperbranched macromolecules having only dendritic units and no linear spacers, while the structure of these macromolecules is modeled through T-fractals. We construct a full set of eigenmodes of the dynamical matrix, which couples the set of Langevin equations. Based on the ensuing relaxation spectra, we analyze the mechanical relaxation moduli. The fractal character of the macromolecules reveals itself in the storage and loss moduli in the intermediate region of frequencies through scaling, whereas at higher frequencies, we observe the locally-dendritic structure that is more pronounced for higher stiffness.
| S-EPMC6432473 | biostudies-literature
Project description:PURPOSE:To estimate multiple components within a single voxel in magnetic resonance fingerprinting when the number and types of tissues comprising the voxel are not known a priori. THEORY:Multiple tissue components within a single voxel are potentially separable with magnetic resonance fingerprinting as a result of differences in signal evolutions of each component. The Bayesian framework for inverse problems provides a natural and flexible setting for solving this problem when the tissue composition per voxel is unknown. Assuming that only a few entries from the dictionary contribute to a mixed signal, sparsity-promoting priors can be placed upon the solution. METHODS:An iterative algorithm is applied to compute the maximum a posteriori estimator of the posterior probability density to determine the magnetic resonance fingerprinting dictionary entries that contribute most significantly to mixed or pure voxels. RESULTS:Simulation results show that the algorithm is robust in finding the component tissues of mixed voxels. Preliminary in vivo data confirm this result, and show good agreement in voxels containing pure tissue. CONCLUSIONS:The Bayesian framework and algorithm shown provide accurate solutions for the partial-volume problem in magnetic resonance fingerprinting. The flexibility of the method will allow further study into different priors and hyperpriors that can be applied in the model. Magn Reson Med 80:159-170, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
| S-EPMC5876128 | biostudies-literature