Dataset Information

Calculation of absolute molecular entropies and heat capacities made simple.

ABSTRACT: We propose a fully-automated composite scheme for the accurate and numerically stable calculation of molecular entropies by efficiently combining density-functional theory (DFT), semi-empirical methods (SQM), and force-field (FF) approximations. The scheme is systematically expandable and can be integrated seamlessly with continuum-solvation models. Anharmonic effects are included through the modified rigid-rotor-harmonic-oscillator (msRRHO) approximation and the Gibbs-Shannon formula for extensive conformer ensembles (CEs), which are generated by a metadynamics search algorithm and are extrapolated to completeness. For the first time, variations of the ro-vibrational entropy over the CE are consistently accounted-for through a Boltzmann-population average. Extensive tests of the protocol with the two standard DFT approaches B97-3c and B3LYP-D3 reveal an unprecedented accuracy with mean deviations <1 cal mol^-1 K^-1 (about <1-2%) for the total gas phase molecular entropy of medium-sized molecules. Even for the hardship case of extremely flexible linear alkanes (C₁₄H₃₀-C₁₆H₃₄), errors are only about 3 cal mol^-1 K^-1. Comprehensive tests indicate a relatively strong variation of the conformational entropy on the underlying level of theory for typical drug molecules, inferring the complex potential energy surfaces as the main source of error. Furthermore, we show some application examples for the calculation of free energy differences in typical chemical reactions.

SUBMITTER: Pracht P

PROVIDER: S-EPMC8139639 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The hypothetical scanning molecular dynamics (HSMD) method is used here for calculating the absolute free energy of binding, ΔA(0) of the complex of the protein FKBP12 with the ligand SB2 (also denoted L8) - a system that has been studied previously for comparing the performance of different methods. Our preliminary study suggests that considering long-range electrostatics is imperative even for a hydrophobic ligand such as L8. Therefore the system is modeled by the AMBER force field using Particle Mesh Ewald (PME). HSMD consists of three stages applied to both the ligand-solvent and ligand-protein systems. (1) A small set of system configurations (frames) is extracted from an MD trajectory. (2) The entropy of the ligand in each frame is calculated by a reconstruction procedure. (3) The contribution of water and protein to ΔA(0) is calculated for each frame by gradually increasing the ligand-environment interactions from zero to their full value using thermodynamic integration (TI). Unlike the conventional methods, the structure of the ligand is kept fixed during TI, and HSMD is thus free from the end-point problem encountered with the double annihilation method (DAM); therefore, the need for applying restraints is avoided. Furthermore, unlike the conventional methods, the entropy of the ligand and water is obtained directly as a byproduct of the simulation. In this paper, in addition to the difference in the internal entropies of the ligand in the two environments, we calculate for the first time the external entropy of the ligand, which provides a measure for the size of the active site. We obtain ΔA(0) = -10.7 ±1.0 as compared to the experimental values -10.9 and -10.6 kcal/mol. However, a protein/water system treated by periodic boundary conditions grows significantly with increasing protein size and the computation of ΔA(0) would become expensive by all methods. Therefore, we also apply HSMD to FKBP12-L8 described by the GSBP/SSBP model of Roux's group (implemented in the software CHARMM) where only part of the protein and water around the active site are considered and long-range electrostatic effects are taken into account. For comparison this model was also treated by the double decoupling method (DDM). The two methods have led to comparable results for ΔA(0) which are somewhat lower than the experimental value. The ligand was found to be more confined in the active site described by GSBP/SSBP than by PME where its entropy in solvent is larger than in the active site by 1.7 and by 5.5 kcal/mol, respectively.

Project description:BACKGROUND:The 2014/2015 Ebolavirus outbreak resulted in more than 28,000 cases and 11,323 reported deaths, as of March 2016. Domestic transmission of the Guinea strain associated with the outbreak occurred mainly in six African countries, and international transmission was reported in four countries. Outbreak management was limited by the inability to rapidly diagnose infected cases. A further fifteen countries in Africa are predicted to be at risk of Ebolavirus outbreaks in the future as a consequence of climate change and urbanization. Early detection of cases and reduction of transmission rates is critical to prevent and manage future severe outbreaks. We designed a rapid assay for detection of Ebolavirus using recombinase polymerase amplification, a rapid isothermal amplification technology that can be combined with portable lateral flow detection technology. The developed rapid assay operates in 30 min and was comparable with real-time TaqMan™ PCR. METHODS:Designed, screened, selected and optimized oligonucleotides using the NP coding region from Ebola Zaire virus (Guinea strain). We determined the analytical sensitivity of our Ebola rapid molecular test by testing selected primers and probe with tenfold serial dilutions (1.34 × 1010- 1.34 × 101 copies/μL) of cloned NP gene from Mayinga strain of Zaire ebolavirus in pCAGGS vector, and serially diluted cultured Ebolavirus as established by real-time TaqMan™ PCR that was performed using ABI7500 in Fast Mode. We tested extracted and reverse transcribed RNA from cultured Zaire ebolavirus strains - Mayinga, Gueckedou C05, Gueckedou C07, Makona, Kissidougou and Kiwit. We determined the analytical specificity of our assay with related viruses: Marburg, Ebola Reston and Ebola Sudan. We further tested for Dengue virus 1-4, Plasmodium falciparum and West Nile Virus (Kunjin strain). RESULTS:The assay had a detection limit of 134 copies per μL of plasmid containing the NP gene of Ebolavirus Mayinga, and cultured Ebolavirus and was highly specific for the Zaire ebolavirus species, including the Guinea strain responsible for the 2014/2015 outbreak. The assay did not detect related viruses like Marburg, Reston, or Sudan viruses, and other pathogens likely to be isolated from clinical samples. CONCLUSIONS:Our assay could be suitable for implementation in district and primary health laboratories, as only a heating block and centrifuge is required for operation. The technique could provide a pathway for rapid screening of patients and animals for improved management of outbreaks.

Dataset Information

Calculation of absolute molecular entropies and heat capacities made simple.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets