DFT method
Dotaz
Zobrazit nápovědu
In this work we investigate the performance of the DFT method, augmented with an empirical dispersion function (DFT-D), paired with the PCM implicit solvation model, for the computation of noncovalent interaction energies of biologically-relevant, solvated model complexes. It is found that this method describes intermolecular interactions within water and ether (protein-like) environments with roughly the same accuracy as in the gas phase. Another important finding is that, when environmental effects are taken into account, the empirical dispersion term associated with the DFT-D method need be modified very little (or not at all), in order to obtain the optimum, most well balanced, performance.
- MeSH
- algoritmy MeSH
- biologické modely MeSH
- chemické modely MeSH
- DNA chemie MeSH
- hydrofobní a hydrofilní interakce MeSH
- kvantová teorie MeSH
- makromolekulární látky chemie MeSH
- povrchové vlastnosti MeSH
- přenos energie MeSH
- proteiny chemie MeSH
- RNA chemie MeSH
- rozpustnost MeSH
- roztoky MeSH
- termodynamika MeSH
- vodíková vazba MeSH
Oxidative stress can lead to various derivatives of the tyrosine residue in peptides and proteins. A typical product is 3-nitro-L-tyrosine residue (Nit), which can affect protein behavior during neurodegenerative processes, such as those associated with Alzheimer's and Parkinson's diseases. Surface enhanced Raman spectroscopy (SERS) is a technique with potential for detecting peptides and their metabolic products at very low concentrations. To explore the applicability to Nit, we use SERS to monitor tyrosine nitration in Met-Enkephalin, rev-Prion protein, and α-synuclein models. Useful nitration indicators were the intensity ratio of two tyrosine marker bands at 825 and 870 cm-1 and a bending vibration of the nitro group. During the SERS measurement, a conversion of nitrotyrosine to azobenzene containing peptides was observed. The interpretation of the spectra has been based on density functional theory (DFT) simulations. The CAM-B3LYP and ωB97XD functionals were found to be most suitable for modeling the measured data. The secondary structure of the α-synuclein models was monitored by electronic and vibrational circular dichroism (ECD and VCD) spectroscopies and modeled by molecular dynamics (MD) simulations. The results suggest that the nitration in these peptides has a limited effect on the secondary structure, but may trigger their aggregation.
- MeSH
- azosloučeniny chemie MeSH
- cirkulární dichroismus MeSH
- peptidy chemická syntéza chemie MeSH
- Ramanova spektroskopie metody MeSH
- sekundární struktura proteinů MeSH
- simulace molekulární dynamiky MeSH
- teorie funkcionálu hustoty MeSH
- tyrosin analogy a deriváty analýza MeSH
- Publikační typ
- časopisecké články MeSH
We designed 0D, 1D, and 2D supramolecular assemblies made of diaryliodonium salts (functioning as double σ-hole donors) and carboxylates (as σ-hole acceptors). The association was based on two charge-supported halogen bonds (XB), which occurred between IIII sites of the iodonium cations and the carboxylate anions. The sequential introduction of the carboxylic groups in the aryl ring of the benzoic acid added a dimension to the 0D supramolecular organization of the benzoate, which furnished 1D-chained and 2D-layered structures when terephthalate and trimesate anions, correspondingly, were applied as XB acceptors. The structure-directing XB were studied using DFT calculations under periodic boundary conditions and were followed by the one-electron-potential analysis and the Bader atoms-in-molecules topological analysis of electron density. These theoretical methods confirmed the existence of the XB and verified the philicities of the interaction partners in the designed solid-state structures.
To find and calibrate a robust and reliable computational protocol for mapping conformational space of medium-sized molecules, exhaustive conformational sampling has been carried out for a series of seven macrocyclic compounds of varying ring size and one acyclic analogue. While five of them were taken from the MD/LLMOD/force field study by Shelley and co-workers ( Watts , K. S. ; Dalal , P. ; Tebben , A. J. ; Cheney , D. L. ; Shelley , J. C. Macrocycle Conformational Sampling with MacroModel . J. Chem. Inf. MODEL: 2014 , 54 , 2680 - 2696 ), three represent potential macrocyclic inhibitors of human cyclophilin A. The free energy values (GDFT/COSMO-RS) for all of the conformers of each compound were obtained by a composite protocol based on in vacuo quantum mechanics (DFT-D3 method in a large basis set), standard gas-phase thermodynamics, and the COSMO-RS solvation model. The GDFT/COSMO-RS values were used as the reference for evaluating the performance of conformational sampling algorithms: standard and extended MD/LLMOD search (simulated-annealing molecular dynamics with low-lying eigenvector following algorithms, employing the OPLS2005 force field plus GBSA solvation) available in MacroModel and plain molecular dynamics (MD) sampling at high temperature (1000 K) using the semiempirical quantum mechanics (SQM) potential SQM(PM6-D3H4/COSMO) followed by energy minimization of the snapshots. It has been shown that the former protocol (MD/LLMOD) may provide a more complete set of initial structures that ultimately leads to the identification of a greater number of low-energy conformers (as assessed by GDFT/COSMO-RS) than the latter (i.e., plain SQM MD). The CPU time needed to fully evaluate one medium-sized compound (∼100 atoms, typically resulting in several hundred or a few thousand conformers generated and treated quantum-mechanically) is approximately 1,000-100,000 CPU hours on today's computers, which transforms to 1-7 days on a small-sized computer cluster with a few hundred CPUs. Finally, our data sets based on the rigorous quantum-chemical approach allow us to formulate a composite conformational sampling protocol with multiple checkpoints and truncation of redundant structural data that offers superior insights at affordable computational cost.
- MeSH
- algoritmy MeSH
- kalibrace MeSH
- krystalografie MeSH
- kvantová teorie MeSH
- makrocyklické sloučeniny chemie MeSH
- molekulární konformace * MeSH
- rychlé screeningové testy MeSH
- simulace molekulární dynamiky MeSH
- termodynamika MeSH
- vysoká teplota MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- srovnávací studie MeSH
AIMS: The aim of the present research was to synthesize glycoluril derivative 2,4-Bis(4- cyanobenzyl)glycoluril through a convergent scheme. BACKGROUND: For this purpose, Sandmeyer reaction procedure was employed for the synthesis of said compound. The structure of the pure compound was confirmed by using different spectroscopic techniques, such as 1HNMR, 13C-NMR and (HR-MS) Mass spectrometry. OBJECTIVE: Convergent synthesis of 2,4-BIS (4-CYANOBENZYL)GLYCOLURIL USING SANDMEYER REACTION and urease inhibition study. METHODS: The structure of the pure compound was confirmed by using different spectroscopic techniques such as 1H-NMR, 13C-NMR and (HR-MS) Mass spectrometry. The electronic properties of the newly synthesized compound and thiourea were determined by using density functional theory. RESULTS: Furthermore, the compound was evaluated against urease enzyme and was found to be potent inhibitors with an IC50 value of 11.5 ± 1.50 μM when compared with standard inhibitor thiourea (IC50 = 21.0 ± 1.90 μM). The compound may serve as a lead compound to synthesize new cyano-based bambusuril in the future with enhanced biological properties. CONCLUSION: We have synthesized a new glycoluril derivative 2,4-Bis(4-cyanobenzyl)glycoluril by the sandmeyer reaction. It has been obtained in the form of light yellowish powder in good yield (96%). Glycoluril based macrocycles have been used in various fields; starting from the 2,4-Bis(4-nitrobenzyl)glycoluril (already reported compound), which has undergone reduction (CH3OH,Pt/C) , diazotization (NaNO2/HCl), cyanation (CuCl/KCN), respectively in order to synthesize the desired new glycoluril derivative. The obtained product will be used as a building block for the synthesis of the cyano based bambusuril marcocycle in the future. The yield of the obtained product has been monitored by using different amounts of cyanating reagent, but the best results are shown by the use of 4 mmol of CuCl/KCN. KCN with CuCl assisted the conversion of diazo group into the cyano group with enhanced yield when used in excess amount. It acts as a catalyst. The solubility characteristic of 2,4-Bis(4-cyanobenzyl)glycoluril has also been determined in different organic solvents. 1H NMR technique proved to be very helpful for the structure determination of our desired product. Benzylic protons give signals at 7.5 ppm and 7.8 ppm, respectively. The downfield peaks confirm the presence of CN group near the benzylic protons. Methine protons show a signal at 5.2 ppm, which ensures the basic skeleton of glycoluril. Ureidyl protons also confirm the synthesis of the heterocyclic 2,4-Bis(4-cyanobenzyl)glycoluril compound. The negative and positive electrostatic potential sites, molecular descriptors, and charge density distribution of frontier molecular orbitals are revealing that 4a with promising sites for electrophilic and nucleophilic attacks would result to enhance the urease inhibition, which is in good agreement with the experimental data.
A series of triterpenoids of the lupane, taraxastane, friedelane and baccharane type were oxidized using selenium dioxide (SeO2) and benzeneseleninic anhydride (BSA) under various conditions. Depending on the reaction conditions, different reaction pathways were observed, including dehydrogenation, allylic oxidation, and 1,2-diketone formation. In this way, derivatives functionalized in the triterpene core (especially in rings A, D, and E), difficult to obtain by other methods, can be easily prepared. In some cases, rarely observed α-phenylseleno-ketones were isolated. An unexpected reaction involving the cleavage of the carbon-carbon double bond was observed in the presence of stoichiometric amounts of osmium tetroxide. Further transformations of selected intermediates facilitated the synthesis of new, functionally enriched derivatives. The key reaction pathways were investigated using density functional theory (DFT), focusing on bond length variations and transition states, revealing energetically favored pathways and critical transition structures, including covalent and noncovalent interactions. Solvent and isomerization equilibrium effects were proposed to explain the experimentally observed discrepancies. Cytotoxic activity of selected derivatives was investigated. Derivatives 4 and 38 showed strongest cytotoxicity in cancer cells and fibroblasts (IC50 2.6-26.4 μM); some compounds were selective for G-361 or HeLa cells. These results suggest that they may find application in pharmaceuticals.
- MeSH
- lidé MeSH
- molekulární struktura MeSH
- nádorové buněčné linie MeSH
- oxidace-redukce MeSH
- pentacyklické triterpeny MeSH
- proliferace buněk účinky léků MeSH
- protinádorové látky * farmakologie chemie chemická syntéza MeSH
- screeningové testy protinádorových léčiv MeSH
- selen * chemie MeSH
- teorie funkcionálu hustoty MeSH
- triterpeny * chemie farmakologie chemická syntéza MeSH
- vztah mezi dávkou a účinkem léčiva MeSH
- vztahy mezi strukturou a aktivitou MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
The CCSD(T) interaction energies for the H-bonded and stacked structures of the uracil dimer are determined at the aug-cc-pVDZ and aug-cc-pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug-cc-pVDZ and aug-cc-pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug-cc-pVTZ and aug-cc-pVQZ basis sets) and extrapolated DeltaCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H-bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS-MP2 (SCS: spin-component-scaled), SCS(MI)-MP2 (MI: molecular interaction), MP3, dispersion-augmented DFT (DFT-D), M06-2X, and DFT-SAPT (SAPT: symmetry-adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)-MP2 method. Remarkably good binding energies are also obtained with the DFT-SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.
A new and more aggressive strain of coronavirus, known as SARS-CoV-2, which is highly contagious, has rapidly spread across the planet within a short period of time. Due to its high transmission rate and the significant time-space between infection and manifestation of symptoms, the WHO recently declared this a pandemic. Because of the exponentially growing number of new cases of both infections and deaths, development of new therapeutic options to help fight this pandemic is urgently needed. The target molecules of this study were the nitro derivatives of quinoline and quinoline N-oxide. Computational design at the DFT level, docking studies, and molecular dynamics methods as a well-reasoned strategy will aid in elucidating the fundamental physicochemical properties and molecular functions of a diversity of compounds, directly accelerating the process of discovering new drugs. In this study, we discovered isomers based on the nitro derivatives of quinoline and quinoline N-oxide, which are biologically active compounds and may be low-cost alternatives for the treatment of infections induced by SARS-CoV-2.
- MeSH
- chinoliny chemie terapeutické užití MeSH
- COVID-19 MeSH
- farmakoterapie COVID-19 MeSH
- počítačová simulace MeSH
- preklinické hodnocení léčiv MeSH
- SARS-CoV-2 chemie MeSH
- simulace molekulární dynamiky MeSH
- simulace molekulového dockingu MeSH
- teorie funkcionálu hustoty MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Accurate estimation of protein-ligand binding affinity is the cornerstone of computer-aided drug design. We present a universal physics-based scoring function, named SQM2.20, addressing key terms of binding free energy using semiempirical quantum-mechanical computational methods. SQM2.20 incorporates the latest methodological advances while remaining computationally efficient even for systems with thousands of atoms. To validate it rigorously, we have compiled and made available the PL-REX benchmark dataset consisting of high-resolution crystal structures and reliable experimental affinities for ten diverse protein targets. Comparative assessments demonstrate that SQM2.20 outperforms other scoring methods and reaches a level of accuracy similar to much more expensive DFT calculations. In the PL-REX dataset, it achieves excellent correlation with experimental data (average R2 = 0.69) and exhibits consistent performance across all targets. In contrast to DFT, SQM2.20 provides affinity predictions in minutes, making it suitable for practical applications in hit identification or lead optimization.
- MeSH
- ligandy MeSH
- proteiny * metabolismus MeSH
- racionální návrh léčiv * MeSH
- termodynamika MeSH
- vazba proteinů MeSH
- Publikační typ
- časopisecké články MeSH
Resolution of identity standard density functional theory augmented with a damped empirical dispersion term (RI-DFT-D) calculations have been carried out on a set of lowest energy minima of tryptophyl-glycine (Trp-Gly) and tryptophyl-glycyl-glycine (Trp-Gly-Gly) peptides. RI-DFT-D (TPSS/TZVP) results are in excellent agreement with benchmark data based on the CCSD(T) method. Experimental spectra could be assigned according to the calculated IR frequencies. Central processing unit (CPU) time requirements are only slightly higher than those needed for the DFT calculations. Consequently, RI-DFT-D theory seems to be a promising methodology for studying oligopeptides with accuracy comparable to ab initio quantum chemical calculations.
