Browse
Submit Data
Databases
API
Help

Dataset Information

22 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Alkane/Water Partition Coefficient Calculation Based on the Modified AM1 Method and Internal Hydrogen Bonding Sampling Using COSMO-RS.

ABSTRACT: We introduce a physics-based model for calculating partition coefficients of solutes between water and alkanes, using a combination of a semi-empirical method for COSMO charge density calculation and statistical sampling of internal hydrogen bonds (IHBs). We validate the model on the experimental partition data (∼3500 molecules) of small organics, drug-like molecules, and statistical assessment of modeling of proteins and ligand drugs. The model combines two novel algorithms: a bond-correction method for improving the calculation of COSMO charge density from AM1 calculations and a sampling method to deal with IHBs. From a comparison of simulated and experimental partition coefficients, we find a root-mean-square deviation of roughly one log 10 unit. From IHB analysis, we know that IHBs can be present in two states: open (in water) and closed (in apolar solvent). The difference can lead to a shift of as much as two log 10 units per IHB; not taking this effect into account can lead to substantial errors. The method takes a few minutes of calculation time on a single core, per molecule. Although this is still much slower than quantitative structure-activity relationship, it is much faster than molecular simulations and can be readily incorporated into any screening method.

SUBMITTER: Petris PC

PROVIDER: S-EPMC8317156 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Effect of Atomic Charges on Octanol-Water Partition Coefficient Using Alchemical Free Energy Calculation.

Project description:The octanol-water partition coefficient (logPow) is an important index for measuring solubility, membrane permeability, and bioavailability in the drug discovery field. In this paper, the logPow values of 58 compounds were predicted by alchemical free energy calculation using molecular dynamics simulation. In free energy calculations, the atomic charges of the compounds are always fixed. However, they must be recalculated for each solvent. Therefore, three different sets of atomic charges were tested using quantum chemical calculations, taking into account vacuum, octanol, and water environments. The calculated atomic charges in the different environments do not necessarily influence the correlation between calculated and experimentally measured ?Gwater values. The largest correlation coefficient values of the solvation free energy in water and octanol were 0.93 and 0.90, respectively. On the other hand, the correlation coefficient of logPow values calculated from free energies, the largest of which was 0.92, was sensitive to the combination of the solvation free energies calculated from the calculated atomic charges. These results reveal that the solvent assumed in the atomic charge calculation is an important factor determining the accuracy of predicted logPow values.

| S-EPMC6017306 | biostudies-literature

Predicting pKa for proteins using COSMO-RS.

Project description:We have used the COSMO-RS implicit solvation method to calculate the equilibrium constants, pKa, for deprotonation of the acidic residues of the ovomucoid inhibitor protein, OMTKY3. The root mean square error for comparison with experimental data is only 0.5 pH units and the maximum error 0.8 pH units. The results show that the accuracy of pKa prediction using COSMO-RS is as good for large biomolecules as it is for smaller inorganic and organic acids and that the method compares very well to previous pKa predictions of the OMTKY3 protein using Quantum Mechanics/Molecular Mechanics. Our approach works well for systems of about 1000 atoms or less, which makes it useful for small proteins as well as for investigating portions of larger proteins such as active sites in enzymes.

| S-EPMC3817581 | biostudies-literature

Screening Deep Eutectic Solvents for CO2 Capture With COSMO-RS.

Project description:In this work, 502 experimental data for CO2 solubilities and 132 for Henry's constants of CO2 in DESs were comprehensively summarized from literatures and used for further verification and development of COSMO-RS. Large systematic deviations of 62. 2, 59.6, 63.0, and 59.1% for the logarithmic CO2 solubilities in the DESs (1:2, 1:3, 1:4, 1:5), respectively, were observed for the prediction with the original COSMO-RS, while the predicted Henry's constants of CO2 in the DESs (1:1.5, 1:2, 1:3, 1:4, 1:5) at temperatures ranging of 293.15-333.15 K are more accurate than the predicted CO2 solubility with the original COSMO-RS. To improve the performance of COSMO-RS, 502 data points of CO2 solubility in the DESs (1:2, 1:3, 1:4, 1:5) were used for correcting COSMO-RS with a temperature-pressure dependent parameter, and the CO2 solubility in the DES (1:6) was predicted to further verify the performance of the corrected model. The results indicate that the corrected COSMO-RS can significantly improve the model performance with the ARDs decreasing down to 6.5, 4.8, 6.5, and 4.5% for the DESs (1:2, 1:3, 1:4, and 1:5), respectively, and the corrected COSMO-RS with the universal parameters can be used to predict the CO2 solubility in DESs with different mole ratios, for example, for the DES (1:6), the corrected COSMO-RS significantly improves the prediction with an ARD of 10.3% that is much lower than 78.2% provided by the original COSMO-RS. Additionally, the result from COSMO-RS shows that the σ-profiles can reflect the strength of molecular interactions between an HBA (or HBD) and CO2, determining the CO2 solubility, and the dominant interactions for CO2 capture in DESs are the H-bond and Van der Waals force, followed by the misfit based on the analysis of the predicted excess enthalpies.

| S-EPMC7031488 | biostudies-literature

A Method of Calculating the Kamlet-Abboud-Taft Solvatochromic Parameters Using COSMO-RS.

Project description:There is demand for safer and bio-based solvents, brought on by legislation and sustainability objectives. The prediction of physical properties is highly desirable to help design new molecules. Here we present an in silico approach to obtain calculated Kamlet-Abboud-Taft solvatochromic parameters using virtual experiments. The tautomerisation equilibrium of methyl acetoacetate and dimedone was calculated in different solvents with COSMO-RS theory and converted into estimates of solvent dipolarity and hydrogen bond accepting ability, respectively. Hydrogen bond donating ability was calculated as a function of the electron deficient surface area on protic solvents. These polarity descriptors correlate with rate constants and equilibria, and so ability of calculated Kamlet-Abboud-Taft solvatochromic parameters to recreate experimental free energy relationships was tested with sixteen case studies taken from the literature. The accuracy of the calculated parameters was also satisfactory for solvent selection, as demonstrated with a 1,4-addition reaction and a multicomponent heterocycle synthesis.

| S-EPMC6630472 | biostudies-literature

Stability of cellulase in ionic liquids: correlations between enzyme activity and COSMO-RS descriptors.

Project description:Ionic liquids (ILs) are effective in pretreating cellulose for enhanced enzymatic saccharification, however ILs can inactivate cellulases. To guide the selection of ILs, the activity of cellulase was correlated with COSMO-RS calculations and descriptors of ILs including hydrogen bond (H-bond) basicity/acidity, polarity and ion size. Trends were deduced using an anion-series and a cation-series of ionic liquids in aqueous solutions. The activity in the cation-series was best correlated with the size of varied cations, whereas the activity in the anion-series showed a pronounced correlation to H-bond basicity and polarity of different anions. COSMO-RS was further used to predict the solubility of cellulose in ILs, which was correlated with cellulase activity on IL-pretreated cellulose. The best correlations were found between the enzyme activity in the anion-series ILs and the logarithmic activity coefficients, the H-bond energy, H-bond basicity and polarizability, underlining that the anion plays a crucial role in cellulose dissolution.

| S-EPMC6877754 | biostudies-literature

Screening ionic liquids for dissolving hemicellulose by COSMO-RS based on the selective model.

Project description:The utilization of biomass resources has attracted more and more attention due to the consumption of non-renewable resources. Compared with cellulose and lignin, hemicellulose has been less studied. Some ionic liquids (ILs) have been proved to be excellent solvents for lignocellulosic pretreatment. However, screening of more efficient ILs is difficult due to numerous possible ILs. Computational chemistry has been proved effective in solvent screening, but a precise model is indispensable. In this work, we focused on building several appropriate models and selected the most suitable one. According to the structure of hemicellulose, six hemicellulose models were constructed and the mid-dimer of the xylan chain hemicellulose (MDXC) model was proved to be the best compared with the reported experimental results. Based on the MDXC model, 1368 ILs were screened to evaluate their ability to dissolve hemicellulose by Conductor-like Screening Model for Real Solvents (COSMO-RS). The activity coefficient (γ), excess enthalpy (HE), and σ-profile indicated that the hydrogen-bond (H-bond) played a vital role in the dissolution of hemicellulose. Anions played a more critical role than cations, where small anions with H-bond acceptor groups could enhance the molecular interactions with hemicellulose. This work provided a thermodynamic understanding of hemicellulose and IL solvent systems. It highlights the importance of building appropriate solute models, which may be necessary to predict of the other thermodynamic properties in the future.

| S-EPMC9169070 | biostudies-literature

MoDoop: An Automated Computational Approach for COSMO-RS Prediction of Biopolymer Solubilities in Ionic Liquids.

Project description:An automated computational framework (MoDoop) was developed to predict the biopolymer solubilities in ionic liquids (ILs) on the basis of conductor-like screening model for real solvents calculations of two thermodynamic properties: logarithmic activity coefficient (ln γ) at infinite dilution and excess enthalpy (H E) of mixture. The calculation was based on the optimized two-dimensional structures of biopolymer models and ILs by searching the lowest-energy conformer and optimizing molecular geometry. Three lignin models together with one IL dataset were used to evaluate the prediction ability of the developed method. The evaluation results show that ln γ is a more reliable property to predict lignin solubilities in ILs and the p-coumaryl alcohol model is considered as the best model to represent lignin molecules. The developed MoDoop approach is efficient for rapid in silico screening of suitable ionic liquids to dissolve biopolymers.

| S-EPMC6648271 | biostudies-literature

A NON-LINEAR STRUCTURE-PROPERTY MODEL FOR OCTANOL-WATER PARTITION COEFFICIENT.

Project description:Octanol-water partition coefficient (K(ow)) is an important thermodynamic property used to characterize the partitioning of solutes between an aqueous and organic phase and has importance in such areas as pharmacology, pharmacokinetics, pharmacodynamics, chemical production and environmental toxicology. We present a non-linear quantitative structure-property relationship model for determining K(ow) values of new molecules in silico. A total of 823 descriptors were generated for 11,308 molecules whose K(ow) values are reported in the PhysProp dataset by Syracuse Research. Optimum network architecture and its associated inputs were identified using a wrapper-based feature selection algorithm that combines differential evolution and artificial neural networks. A network architecture of 50-33-35-1 resulted in the least root-mean squared error (RMSE) in the training set. Further, to improve on single-network predictions, a neural network ensemble was developed by combining five networks that have the same architecture and inputs but differ in layer weights. The ensemble predicted the K(ow) values with RMSE of 0.28 and 0.38 for the training set and internal validation set, respectively. The ensemble performed reasonably well on an external dataset when compared with other popular K(ow) models in the literature.

| S-EPMC3505089 | biostudies-literature

Octanol-water partition coefficient measurements for the SAMPL6 blind prediction challenge.

Project description:Partition coefficients describe the equilibrium partitioning of a single, defined charge state of a solute between two liquid phases in contact, typically a neutral solute. Octanol-water partition coefficients ([Formula: see text]), or their logarithms (log P), are frequently used as a measure of lipophilicity in drug discovery. The partition coefficient is a physicochemical property that captures the thermodynamics of relative solvation between aqueous and nonpolar phases, and therefore provides an excellent test for physics-based computational models that predict properties of pharmaceutical relevance such as protein-ligand binding affinities or hydration/solvation free energies. The SAMPL6 Part II octanol-water partition coefficient prediction challenge used a subset of kinase inhibitor fragment-like compounds from the SAMPL6 [Formula: see text] prediction challenge in a blind experimental benchmark. Following experimental data collection, the partition coefficient dataset was kept blinded until all predictions were collected from participating computational chemistry groups. A total of 91 submissions were received from 27 participating research groups. This paper presents the octanol-water log P dataset for this SAMPL6 Part II partition coefficient challenge, which consisted of 11 compounds (six 4-aminoquinazolines, two benzimidazole, one pyrazolo[3,4-d]pyrimidine, one pyridine, one 2-oxoquinoline substructure containing compounds) with log P values in the range of 1.95-4.09. We describe the potentiometric log P measurement protocol used to collect this dataset using a Sirius T3, discuss the limitations of this experimental approach, and share suggestions for future log P data collection efforts for the evaluation of computational methods.

| S-EPMC7301889 | biostudies-literature

COSMO-RS-Based Descriptors for the Machine Learning-Enabled Screening of Nucleotide Analogue Drugs against SARS-CoV-2.

Project description:Chemical similarity-based approaches employed to repurpose or develop new treatments for emerging diseases, such as COVID-19, correlates molecular structure-based descriptors of drugs with those of a physiological counterpart or clinical phenotype. We propose novel descriptors based on a COSMO-RS (short for conductor-like screening model for real solvents) σ-profiles for enhanced drug screening enabled by machine learning (ML). The descriptors' performance is hereby illustrated for nucleotide analogue drugs that inhibit the ribonucleic acid-dependent ribonucleic acid polymerase, key to viral transcription and genome replication. The COSMO-RS-based descriptors account for both chemical reactivity and structure, and are more effective for ML-based screening than fingerprints based on molecular structure and simple physical/chemical properties. The descriptors are evaluated using principal component analysis, an unsupervised ML technique. Our results correlate with the active monophosphate forms of the leading drug remdesivir and the prospective drug EIDD-2801 with nucleotides, followed by other promising drugs, and are superior to those from molecular structure-based descriptors and molecular docking. The COSMO-RS-based descriptors could help accelerate drug discovery for the treatment of emerging diseases.

| S-EPMC7605243 | 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