Balancing Accuracy, Efficiency, and Flexibility in Radiation Calculations for Dynamical Models.
Ontology highlight
ABSTRACT:
SUBMITTER: Pincus R
PROVIDER: S-EPMC6900118 | biostudies-literature | 2019 Oct
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In the companion paper, we presented a set of induced dipole interaction models using four types of screening functions, which include the Applequist (no screening), the Thole linear, the Thole exponential model, and the Thole Tinker-like (another form of exponential screening function) functions. In this work, we evaluate the performance of polarizability models using a large set of amino acid analog pairs in conformations that are frequently observed in protein structures as a benchmark. For each amino acid pair, we calculated quantum mechanical interaction energies at the MP2/aug-cc-pVTZ//MP2/6-311++G(d,p) level with the basis set superposition error (BSSE) correction and compared them with molecular mechanics results. Encouragingly, all polarizable models significantly outperform the additive F94 and F03 models (mimicking AMBER ff94/ff99 and ff03 force fields, respectively) in reproducing the BSSE-corrected quantum mechanical interaction energies. In particular, the root-mean-square errors (RMSEs) for three Thole models in Set A (where the 1-2 and 1-3 interactions are turned off and all 1-4 interactions are included) are 1.456, 1.417, and 1.406 kcal/mol for model AL (Thole Linear), model AE (Thole exponential), and model AT (Thole Tinker-like), respectively. In contrast, the RMSEs are 3.729 and 3.433 kcal/mol for F94 and F03 models, respectively. A similar trend was observed for the average unsigned errors (AUEs), which are 1.057, 1.025, 1.011, 2.219, and 2.070 kcal/mol for AL, AE, AT, F94/ff99, and F03, respectively. Analyses based on the trend line slopes indicate that the two fixed charge models substantially underestimate the relative strengths of noncharge-charge interactions by 24 (F03) and 35% (F94), respectively, whereas the four polarizable models overestimate the relative strengths by 5 (AT), 3 (AL, AE), and 13% (AA), respectively. Agreement was further improved by adjusting the van der Waals parameters. Judging from the notably improved accuracy in comparison with the fixed charge models, the polarizable models are expected to form the foundation for the development of high quality polarizable force fields for protein and nucleic acid simulations.
| S-EPMC3082585 | biostudies-literature
Project description:Understanding the complexity of anisotropic turbulent processes in engineering and environmental fluid flows is a formidable challenge with practical significance because energy often flows intermittently from the smaller scales to impact the largest scales in these flows. Conceptual dynamical models for anisotropic turbulence are introduced and developed here which, despite their simplicity, capture key features of vastly more complicated turbulent systems. These conceptual models involve a large-scale mean flow and turbulent fluctuations on a variety of spatial scales with energy-conserving wave-mean-flow interactions as well as stochastic forcing of the fluctuations. Numerical experiments with a six-dimensional conceptual dynamical model confirm that these models capture key statistical features of vastly more complex anisotropic turbulent systems in a qualitative fashion. These features include chaotic statistical behavior of the mean flow with a sub-Gaussian probability distribution function (pdf) for its fluctuations whereas the turbulent fluctuations have decreasing energy and correlation times at smaller scales, with nearly Gaussian pdfs for the large-scale fluctuations and fat-tailed non-Gaussian pdfs for the smaller-scale fluctuations. This last feature is a manifestation of intermittency of the small-scale fluctuations where turbulent modes with small variance have relatively frequent extreme events which directly impact the mean flow. The dynamical models introduced here potentially provide a useful test bed for algorithms for prediction, uncertainty quantification, and data assimilation for anisotropic turbulent systems.
| S-EPMC4020109 | biostudies-literature
Project description:Spatial systems that experience congestion can be modeled as weighted networks whose weights dynamically change over time with the redistribution of flows. This is particularly true for urban transportation networks. The aim of this work is to find appropriate network measures that are able to detect critical zones for traffic congestion and bottlenecks in a transportation system. We propose for both single and multi-layered networks a path-based measure, called dynamical efficiency, which computes the travel time differences under congested and free-flow conditions. The dynamical efficiency quantifies the reachability of a location embedded in the whole urban traffic condition, in lieu of a myopic description based on the average speed of single road segments. In this way, we are able to detect the formation of congestion seeds and visualize their evolution in time as well-defined clusters. Moreover, the extension to multilayer networks allows us to introduce a novel measure of centrality, which estimates the expected usage of inter-modal junctions between two different transportation means. Finally, we define the so-called dilemma factor in terms of number of alternatives that an interconnected transportation system offers to the travelers in exchange for a small increase in travel time. We find macroscopic relations between the percentage of extra-time, number of alternatives and level of congestion, useful to quantify the richness of trip choices that a city offers. As an illustrative example, we show how our methods work to study the real network of a megacity with probe traffic data.
| S-EPMC8630039 | biostudies-literature
Project description:Accounting for target flexibility and selecting "hot spots" most likely to be able to bind an inhibitor continue to be challenges in the field of structure-based drug design, especially in the case of protein-protein interactions. Computational fragment-based approaches using molecular dynamics (MD) simulations are a promising emerging technology having the potential to address both of these challenges. However, the optimal MD conditions permitting sufficient target flexibility while also avoiding fragment-induced target denaturation remain ambiguous. Using one such technology (Site Identification by Ligand Competitive Saturation, SILCS), conditions were identified to either prevent denaturation or identify and exclude trajectories in which subtle but important denaturation was occurring. The target system used was the well-characterized protein cytokine IL-2, which is involved in a protein-protein interface and, in its unliganded crystallographic form, lacks surface pockets that can serve as small-molecule binding sites. Nonetheless, small-molecule inhibitors have previously been discovered that bind to two "cryptic" binding sites that emerge only in the presence of ligand binding, highlighting the important role of IL-2 flexibility. Using the above conditions, SILCS with hydrophobic fragments was able to identify both sites based on favorable fragment binding while avoiding IL-2 denaturation. An important additional finding was that acetonitrile, a water-miscible fragment, fails to identify either site yet can induce target denaturation, highlighting the importance of fragment choice.
| S-EPMC3438888 | biostudies-literature
Project description:While the spatial topological persistence is naturally constructed from a radius-based filtration, it has hardly been derived from a temporal filtration. Most topological models are designed for the global topology of a given object as a whole. There is no method reported in the literature for the topology of an individual component in an object to the best of our knowledge. For many problems in science and engineering, the topology of an individual component is important for describing its properties. We propose evolutionary homology (EH) constructed via a time evolution-based filtration and topological persistence. Our approach couples a set of dynamical systems or chaotic oscillators by the interactions of a physical system, such as a macromolecule. The interactions are approximated by weighted graph Laplacians. Simplices, simplicial complexes, algebraic groups and topological persistence are defined on the coupled trajectories of the chaotic oscillators. The resulting EH gives rise to time-dependent topological invariants or evolutionary barcodes for an individual component of the physical system, revealing its topology-function relationship. In conjunction with Wasserstein metrics, the proposed EH is applied to protein flexibility analysis, an important problem in computational biophysics. Numerical results for the B-factor prediction of a benchmark set of 364 proteins indicate that the proposed EH outperforms all the other state-of-the-art methods in the field.
| S-EPMC8223814 | biostudies-literature
Project description:For reliable stimulus identification, sensory codes have to be robust by including redundancy to combat noise, but redundancy sacrifices coding efficiency. To address how experience affects the balance between the robustness and efficiency of sensory codes, we probed odor representations in the mouse olfactory bulb during learning over a week, using longitudinal two-photon calcium imaging. When mice learned to discriminate between two dissimilar odorants, responses of mitral cell ensembles to the two odorants gradually became less discrete, increasing the efficiency. In contrast, when mice learned to discriminate between two very similar odorants, the initially overlapping representations of the two odorants became progressively decorrelated, enhancing the robustness. Qualitatively similar changes were observed when the same odorants were experienced passively, a condition that would induce implicit perceptual learning. These results suggest that experience adjusts odor representations to balance the robustness and efficiency depending on the similarity of the experienced odorants.
| S-EPMC5061050 | biostudies-literature
Project description:Dynamics of complex systems is often driven by large and intricate networks of microscopic interactions, whose sheer size obfuscates understanding. With limited experimental data, many parameters of such dynamics are unknown, and thus detailed, mechanistic models risk overfitting and making faulty predictions. At the other extreme, simple ad hoc models often miss defining features of the underlying systems. Here we develop an approach that instead constructs phenomenological, coarse-grained models of network dynamics that automatically adapt their complexity to the available data. Such adaptive models produce accurate predictions even when microscopic details are unknown. The approach is computationally tractable, even for a relatively large number of dynamical variables. Using simulated data, it correctly infers the phase space structure for planetary motion, avoids overfitting in a biological signalling system and produces accurate predictions for yeast glycolysis with tens of data points and over half of the interacting species unobserved.
| S-EPMC4560822 | biostudies-literature
Project description:Proper functioning of biological cells requires that the process of protein expression be carried out with high efficiency and fidelity. Given an amino-acid sequence of a protein, multiple degrees of freedom still remain that may allow evolution to tune efficiency and fidelity for each gene under various conditions and cell types. Particularly, the redundancy of the genetic code allows the choice between alternative codons for the same amino acid, which, although 'synonymous,' may exert dramatic effects on the process of translation. Here we review modern developments in genomics and systems biology that have revolutionized our understanding of the multiple means by which translation is regulated. We suggest new means to model the process of translation in a richer framework that will incorporate information about gene sequences, the tRNA pool of the organism and the thermodynamic stability of the mRNA transcripts. A practical demonstration of a better understanding of the process would be a more accurate prediction of the proteome, given the transcriptome at a diversity of biological conditions.
| S-EPMC3101949 | biostudies-other
Project description:The independent trajectory thermodynamic integration (IT-TI) approach (Lawrenz et. al J. Chem. Theory. Comput. 2009, 5:1106-1116(1)) for free energy calculations with distributed computing is employed to study two distinct cases of protein-ligand binding: first, the influenza surface protein N1 neuraminidase bound to the inhibitor oseltamivir, and second, the M. tuberculosis enzyme RmlC complexed with the molecule CID 77074. For both systems, finite molecular dynamics (MD) sampling and varied molecular flexibility give rise to IT-TI free energy distributions that are remarkably centered on the target experimental values, with a spread directly related to protein, ligand, and solvent dynamics. Using over 2 ?s of total MD simulation, alternative protocols for the practical, general implementation of IT-TI are investigated, including the optimal use of distributed computing, the total number of alchemical intermediates, and the procedure to perturb electrostatics and van der Waals interactions. A protocol that maximizes predictive power and computational efficiency is proposed. IT-TI outperforms traditional TI predictions and allows a straightforward evaluation of the reliability of free energy estimates. Our study has broad implications for the use of distributed computing in free energy calculations of macromolecular systems.
| S-EPMC3146769 | biostudies-literature
Project description:The photochemical dynamics of the pyrrole-pyridine hydrogen-bonded complex has been investigated with computational methods. In this system, a highly polar charge-transfer state of (1)pipi* character drives the proton transfer from pyrrole to pyridine, leading to a conical intersection of S(1) and S(0) energy surfaces. A two-sheeted potential-energy surface including 39 in-plane nuclear degrees of freedom has been constructed on the basis of ab initio multiconfiguration electronic-structure data. The non-Born-Oppenheimer nuclear dynamics has been treated with time-dependent quantum wave-packet methods, including the two or three most relevant nuclear degrees of freedom. The effect of the numerous weakly coupled vibrational modes has been taken into account with reduced-density-matrix methods (multilevel Redfield theory). The results provide insight into the mechanisms of excited-state deactivation of hydrogen-bonded aromatic systems via the electron-driven proton-transfer process. This process is believed to be of relevance for the ultrafast excited-state deactivation of DNA base pairs and may contribute to the photostability of the molecular encoding of the genetic information.
| S-EPMC2529040 | biostudies-literature