Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Modeling molecular kinetics with tICA and the kernel trick.

ABSTRACT: The allure of a molecular dynamics simulation is that, given a sufficiently accurate force field, it can provide an atomic-level view of many interesting phenomena in biology. However, the result of a simulation is a large, high-dimensional time series that is difficult to interpret. Recent work has introduced the time-structure based Independent Components Analysis (tICA) method for analyzing MD, which attempts to find the slowest decorrelating linear functions of the molecular coordinates. This method has been used in conjunction with Markov State Models (MSMs) to provide estimates of the characteristic eigenprocesses contained in a simulation (e.g., protein folding, ligand binding). Here, we extend the tICA method using the kernel trick to arrive at nonlinear solutions. This is a substantial improvement as it allows for kernel-tICA (ktICA) to provide estimates of the characteristic eigenprocesses directly without building an MSM.

SUBMITTER: Schwantes CR

PROVIDER: S-EPMC4610300 | biostudies-literature | 2015 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Modeling molecular kinetics with tICA and the kernel trick.

Schwantes Christian R CR Pande Vijay S VS

Journal of chemical theory and computation 20150201 2

The allure of a molecular dynamics simulation is that, given a sufficiently accurate force field, it can provide an atomic-level view of many interesting phenomena in biology. However, the result of a simulation is a large, high-dimensional time series that is difficult to interpret. Recent work has introduced the time-structure based Independent Components Analysis (tICA) method for analyzing MD, which attempts to find the slowest decorrelating linear functions of the molecular coordinates. Thi ...[more]

PMID: 26528090

Similar Datasets

Nantenine as an acetylcholinesterase inhibitor: SAR, enzyme kinetics and molecular modeling investigations.

Project description:Nantenine, as well as a number of flexible analogs, were evaluated for acetylcholinesterase (AChE) inhibitory activity in microplate spectrophotometric assays based on Ellman's method. It was found that the rigid aporphine core of nantenine is an important structural requirement for its anticholinesterase activity. Nantenine showed mixed inhibition kinetics in enzyme assays. Molecular docking experiments suggest that nantenine binds preferentially to the catalytic site of AChE but is also capable of interacting with the peripheral anionic site (PAS) of the enzyme, thus accounting for its mixed inhibition profile. The aporphine core of nantenine may thus be a useful template for the design of novel PAS or dual-site AChE inhibitors. Inhibiting the PAS is desirable for prevention of aggregation of the amyloid peptide A?, a major causative factor in the progression of Alzheimer's disease (AD).

| S-EPMC4663080 | biostudies-literature

Data on biosynthesis of BPAF glucuronide, enzyme kinetics of BPAF glucuronidation, and molecular modeling.

Project description:Bisphenol AF (BPAF) is in the body mainly metabolized to the corresponding bisphenol AF glucuronide (BPAF-G). While BPAF-G is not commercially available, enzyme-assisted synthesis of BPAF-G using the human recombinant enzyme UGT2A1, purification of BPAF-G by solid phase extraction and semi-preparative HPLC and chemical characterization of BPAF-G by NMR and LC-MS/MS were performed and are described here. Furthermore, BPAF glucuronidation kinetics with the UGT enzymes that showed the highest glucuronidation activity in previous studies (i.e hepatic UGTs 1A3, 2B7, and 2B17, intestinal UGT 1A10 and UGT 2A1 that is present in airways) was performed and data is presented. Hepatic enzymes exhibited high affinities toward BPAF, while extrahepatic UGTs 2A1 and 1A10 showed the high vmax values (3.3 and 3.0?nmol/min/mg, respectively). To understand molecular interactions of BPA, BPAF and BPAF-G with ligand biding sites of several nuclear receptors, molecular modeling was performed and data on the binding modes of BPAF, BPA, and BPAF-G in the ligand-binding sites of nuclear receptors are presented. This article is related to "Endocrine activities and adipogenic effects of bisphenol AF and its main metabolite" (Skledar et al., 2019).

| S-EPMC6356001 | biostudies-literature

Circumventing kinetics in biogeochemical modeling.

Project description:Microbial metabolism drives biogeochemical fluxes in virtually every ecosystem. Modeling these fluxes is challenged by the incredible diversity of microorganisms, whose kinetic parameters are largely unknown. In poorly mixed systems, such as stagnant water columns or sediments, however, long-term bulk microbial metabolism may become limited by physical transport rates of substrates across space. Here we mathematically show that under these conditions, biogeochemical fluxes are largely predictable based on the system's transport properties, chemical boundary conditions, and the stoichiometry of metabolic pathways, regardless of the precise kinetics of the resident microorganisms. We formalize these considerations into a predictive modeling framework and demonstrate its use for the Cariaco Basin subeuphotic zone, one of the largest anoxic marine basins worldwide. Using chemical concentration data solely from the upper boundary (depth 180 m) and lower boundary (depth 900 m), but without a priori knowledge of metabolite fluxes, chemical depth profiles, kinetic parameters, or microbial species composition, we predict the concentrations and vertical fluxes of biologically important substances, including oxygen, nitrate, hydrogen sulfide, and ammonium, across the entire considered depth range (180-900 m). Our predictions largely agree with concentration measurements over a period of 14 years ([Formula: see text] = 0.78-0.92) and become particularly accurate during a period where the system was near biogeochemical steady state (years 2007-2009, [Formula: see text] = 0.86-0.95). Our work enables geobiological predictions for a large class of ecosystems without knowledge of kinetic parameters or geochemical depth profiles. Conceptually, our work provides a possible explanation for the decoupling between microbial species composition and bulk metabolic function, observed in various ecosystems.

| S-EPMC6561284 | biostudies-literature

Modeling inhomogeneous DNA replication kinetics.

Project description:In eukaryotic organisms, DNA replication is initiated at a series of chromosomal locations called origins, where replication forks are assembled proceeding bidirectionally to replicate the genome. The distribution and firing rate of these origins, in conjunction with the velocity at which forks progress, dictate the program of the replication process. Previous attempts at modeling DNA replication in eukaryotes have focused on cases where the firing rate and the velocity of replication forks are homogeneous, or uniform, across the genome. However, it is now known that there are large variations in origin activity along the genome and variations in fork velocities can also take place. Here, we generalize previous approaches to modeling replication, to allow for arbitrary spatial variation of initiation rates and fork velocities. We derive rate equations for left- and right-moving forks and for replication probability over time that can be solved numerically to obtain the mean-field replication program. This method accurately reproduces the results of DNA replication simulation. We also successfully adapted our approach to the inverse problem of fitting measurements of DNA replication performed on single DNA molecules. Since such measurements are performed on specified portion of the genome, the examined DNA molecules may be replicated by forks that originate either within the studied molecule or outside of it. This problem was solved by using an effective flux of incoming replication forks at the model boundaries to represent the origin activity outside the studied region. Using this approach, we show that reliable inferences can be made about the replication of specific portions of the genome even if the amount of data that can be obtained from single-molecule experiments is generally limited.

| S-EPMC3296702 | biostudies-literature

Modeling stochastic kinetics of molecular machines at multiple levels: from molecules to modules.

Project description:A molecular machine is either a single macromolecule or a macromolecular complex. In spite of the striking superficial similarities between these natural nanomachines and their man-made macroscopic counterparts, there are crucial differences. Molecular machines in a living cell operate stochastically in an isothermal environment far from thermodynamic equilibrium. In this mini-review we present a catalog of the molecular machines and an inventory of the essential toolbox for theoretically modeling these machines. The tool kits include 1), nonequilibrium statistical-physics techniques for modeling machines and machine-driven processes; and 2), statistical-inference methods for reverse engineering a functional machine from the empirical data. The cell is often likened to a microfactory in which the machineries are organized in modular fashion; each module consists of strongly coupled multiple machines, but different modules interact weakly with each other. This microfactory has its own automated supply chain and delivery system. Buoyed by the success achieved in modeling individual molecular machines, we advocate integration of these models in the near future to develop models of functional modules. A system-level description of the cell from the perspective of molecular machinery (the mechanome) is likely to emerge from further integrations that we envisage here.

| S-EPMC3672901 | biostudies-other

Omnibus risk assessment via accelerated failure time kernel machine modeling.

Project description:Integrating genomic information with traditional clinical risk factors to improve the prediction of disease outcomes could profoundly change the practice of medicine. However, the large number of potential markers and possible complexity of the relationship between markers and disease make it difficult to construct accurate risk prediction models. Standard approaches for identifying important markers often rely on marginal associations or linearity assumptions and may not capture non-linear or interactive effects. In recent years, much work has been done to group genes into pathways and networks. Integrating such biological knowledge into statistical learning could potentially improve model interpretability and reliability. One effective approach is to employ a kernel machine (KM) framework, which can capture nonlinear effects if nonlinear kernels are used (Scholkopf and Smola, 2002; Liu et al., 2007, 2008). For survival outcomes, KM regression modeling and testing procedures have been derived under a proportional hazards (PH) assumption (Li and Luan, 2003; Cai, Tonini, and Lin, 2011). In this article, we derive testing and prediction methods for KM regression under the accelerated failure time (AFT) model, a useful alternative to the PH model. We approximate the null distribution of our test statistic using resampling procedures. When multiple kernels are of potential interest, it may be unclear in advance which kernel to use for testing and estimation. We propose a robust Omnibus Test that combines information across kernels, and an approach for selecting the best kernel for estimation. The methods are illustrated with an application in breast cancer.

| S-EPMC3869038 | biostudies-literature

Monoamine neurotransmitters as substrates for novel tick sulfotransferases, homology modeling, molecular docking, and enzyme kinetics.

Project description:Blacklegged ticks (Ixodes scapularis) transmit the causative agent of Lyme disease in the Northeastern United States. Current research focuses on elucidating biochemical pathways that may be disrupted to prevent pathogen transmission, thereby preventing disease. Genome screening reported transcripts coding for two putative sulfotransferases in whole tick extracts of the nymphal and larval stages. Sulfotransferases are known to sulfonate phenolic and alcoholic receptor agonists such as 17?-estradiol, thereby inactivating the receptor ligands. We used bioinformatic approaches to predict substrates for Ixosc Sult 1 and Ixosc Sult 2 and tested the predictions with biochemical assays. Homology models of 3D protein structure were prepared, and visualization of the electrostatic surface of the ligand binding cavities showed regions of negative electrostatic charge. Molecular docking identified potential substrates including dopamine, R-octopamine and S-octopamine, which docked into Ixosc Sult 1 with favorable binding affinity and correct conformation for sulfonation. Dopamine, but not R- or S-octopamine, also docked into Ixosc Sult 2 in catalytic binding mode. The predictions were confirmed using cytosolic fractions of whole tick extracts. Dopamine was a good substrate (K(m) = 0.1-0.4 ?M) for the native Ixodes scapularis sulfotransferases from larval and nymphal stages regardless of their fed/unfed status. Octopamine sulfonation was detected only after feeding when gene expression data suggests that Ixosc Sult 1 is present. Because dopamine is known to stimulate salivation in ticks through receptor stimulation, these results imply that the function(s) of Ixosc Sult 1 or 2 may include inactivation of the salivation signal via sulfonation of dopamine and/or octopamine.

| S-EPMC3042050 | biostudies-literature

Inhibition kinetics, molecular docking, and stability studies of the effect of papain-generated peptides from palm kernel cake proteins on angiotensin-converting enzyme (ACE).

Project description:Three novel peptide sequences YGIKVGYAIP, GGIF, and GIFE from papain-generated protein hydrolysate of palm kernel cake proteins were used for stability study against ACE, ACE-inhibition kinetics, and molecular docking studies. Results showed that peptide YGIKVGYAIP was degraded, and its ACE-inhibitory activity decreased after 3 h pre-incubation with ACE, while peptides GGIF and GIFE were resistant. However, although the ACE-inhibitory activity of GIFE increased during this time, the ACE inhibitory activity of GGIF decreased after pre-incubation with ACE, indicating that peptide. YGIKVGYAIP and GGIF are substrate-type, whereas GIFE is a true-inhibitor type. Peptide YGIKVGYAIP showed the lowest Ki (0.054 mM) in the inhibition kinetics study compared to GGIF and GIFE, with Ki of 1.27 m M and 18 mM, respectively. In addition, YGIKVGYAIP revealed the lowest Km and Vmax and higher CE in different peptide concentrations, implying that the enzyme catalysis decreased, and peptides had some binding affinity to the enzyme in lower concentrations, which led to reduced catalytic ability. Furthermore, YGIKVGYAIP showed the lowest docking score of -14.733 and 21 interactions with tACE, while GGIF revealed the higher docking score of -8.006 with 15 interactions with tACE.

| S-EPMC9789325 | biostudies-literature

DNA sequence+shape kernel enables alignment-free modeling of transcription factor binding.

Project description:Transcription factors (TFs) bind to specific DNA sequence motifs. Several lines of evidence suggest that TF-DNA binding is mediated in part by properties of the local DNA shape: the width of the minor groove, the relative orientations of adjacent base pairs, etc. Several methods have been developed to jointly account for DNA sequence and shape properties in predicting TF binding affinity. However, a limitation of these methods is that they typically require a training set of aligned TF binding sites.We describe a sequence?+?shape kernel that leverages DNA sequence and shape information to better understand protein-DNA binding preference and affinity. This kernel extends an existing class of k-mer based sequence kernels, based on the recently described di-mismatch kernel. Using three in vitro benchmark datasets, derived from universal protein binding microarrays (uPBMs), genomic context PBMs (gcPBMs) and SELEX-seq data, we demonstrate that incorporating DNA shape information improves our ability to predict protein-DNA binding affinity. In particular, we observe that (i) the k-spectrum?+?shape model performs better than the classical k-spectrum kernel, particularly for small k values; (ii) the di-mismatch kernel performs better than the k-mer kernel, for larger k; and (iii) the di-mismatch?+?shape kernel performs better than the di-mismatch kernel for intermediate k values.The software is available at https://bitbucket.org/wenxiu/sequence-shape.git.rohs@usc.edu or william-noble@uw.edu.Supplementary data are available at Bioinformatics online.

| S-EPMC5870879 | biostudies-literature

Angiotensin Converting Enzyme (ACE)-Peptide Interactions: Inhibition Kinetics, In Silico Molecular Docking and Stability Study of Three Novel Peptides Generated from Palm Kernel Cake Proteins.

Project description:Three novel peptide sequences identified from palm kernel cake (PKC) generated protein hydrolysate including YLLLK, WAFS and GVQEGAGHYALL were used for stability study against angiotensin-converting enzyme (ACE), ACE-inhibition kinetics and molecular docking studies. Results showed that the peptides were degraded at different cleavage degrees of 94%, 67% and 97% for YLLLK, WAFS and GVQEGAGHYALL, respectively, after 3 h of incubation with ACE. YLLLK was found to be the least stable (decreased ACE-inhibitory activity) compared to WAFS and GVQEGAGHYALL (increased ACE-inhibitory activity). YLLLK showed the lowest Ki (1.51 mM) in inhibition kinetics study when compared to WAFS and GVQEGAGHYALL with Ki of 2 mM and 3.18 mM, respectively. In addition, ACE revealed the lowest Kmapp and Vmaxapp and higher catalytic efficiency (CE) in the presence of YLLLK at different concentrations, implying that the enzyme catalysis decreased and hence the inhibition mode increased. Furthermore, YLLLK showed the lowest docking score of -8.224 and seven interactions with tACE, while peptide GVQEGAGHYALL showed the higher docking score of -7.006 and five interactions with tACE.

| S-EPMC6843816 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data