density functional tight binding
Dotaz
Zobrazit nápovědu
A new database of nucleic acid base trimers has been developed that includes 141 geometries and stabilization energies obtained at the RI-MP2 level of theory with the TZVPP basis set. Compared to previously compiled biologically oriented databases, this new construct includes considerably more complicated structures; the various intermolecular interactions in the trimers are quite heterogeneous and in particular include simultaneous hydrogen bonding and stacking interactions, which is similar to the situation in actual biopolymers. Validation against these benchmark data is therefore a more demanding task for approximate models, since correct descriptions of all energy terms are unlikely to be accomplished by fortuitous cancellations of systematic errors. The density functionals TPSS (both with and without an empirical dispersion term), PWB6K, M05-2X, and BH&H, and the self-consistent charge density functional tight binding method augmented with an empirical dispersion term (SCC-DFTB-D) were assessed for their abilities accurately to compute structures and energies. The best reproduction of the BSSE corrected RI-MP2 stabilization energies was achieved by the TPSS functional (TZVPP basis set) combined with empirical dispersion; removal of the dispersion correction leads to significantly degraded performance. The M05-2X and PWB6K functionals performed very well in reproducing the RI-MP2 geometries, but showed a systematic moderate underestimation of the magnitude of base stacking interactions. The SCC-DFTB-D method predicts geometries in fair agreement with RI-MP2; given its computational efficiency it represents a good option for initial scanning of analogous biopolymeric potential energy surfaces. BH&H gives geometries of comparable quality to the other functionals but significantly overestimates interaction energies other than stacking.
A detailed quantum chemical study on five peptides (WG, WGG, FGG, GGF and GFA) containing the residues phenylalanyl (F), glycyl (G), tryptophyl (W) and alanyl (A) -- where F and W are of aromatic character -- is presented. When investigating isolated small peptides, the dispersion interaction is the dominant attractive force in the peptide backbone-aromatic side chain intramolecular interaction. Consequently, an accurate theoretical study of these systems requires the use of a methodology covering properly the London dispersion forces. For this reason we have assessed the performance of the MP2, SCS-MP2, MP3, TPSS-D, PBE-D, M06-2X, BH&H, TPSS, B3LYP, tight-binding DFT-D methods and ff99 empirical force field compared to CCSD(T)/complete basis set (CBS) limit benchmark data. All the DFT techniques with a '-D' symbol have been augmented by empirical dispersion energy while the M06-2X functional was parameterized to cover the London dispersion energy. For the systems here studied we have concluded that the use of the ff99 force field is not recommended mainly due to problems concerning the assignment of reliable atomic charges. Tight-binding DFT-D is efficient as a screening tool providing reliable geometries. Among the DFT functionals, the M06-2X and TPSS-D show the best performance what is explained by the fact that both procedures cover the dispersion energy. The B3LYP and TPSS functionals-not covering this energy-fail systematically. Both, electronic energies and geometries obtained by means of the wave-function theory methods compare satisfactorily with the CCSD(T)/CBS benchmark data.
The unwinding free energy of 128 DNA octamers was correlated with the sum of interaction energies among DNA bases and their solvation energies. The former energies were determined by using the recently developed density functional theory procedure augmented by London dispersion energy (RI-DFT-D) that provides accurate hydrogen-bonding and stacking energies highly comparable with CCSD(T)/complete basis set limit benchmark data. Efficient tight-binding DFT covering dispersion energy was also used and yielded satisfactory results. The latter method can be used for extended systems. The solvation energy was determined by using a C-PCM continuum solvent at HF level calculations. Various models were adopted to correlate theoretical energies with experimental unwinding free energies. Unless all energy components (hydrogen-bonding, intra- and interstrand-stacking, and solvation energies) were included and weighted individually, no satisfactory correlation resulted. The most advanced model yielded very close correlation (RMSE=0.32 kcal mol(-1)) fully comparable with the entirely empirical correlation introduced in the original paper. Analysis of the theoretical results shows the importance of inter- and intramolecular stacking energies, and especially the latter term plays a key role in determining DNA-duplex stabilization.
In this work, we investigate the mode of binding of several steroid hormones, namely aldosterone, deoxycorticosterone, and progesterone to the wild-type and S810L mutated mineralocorticoid (MR) receptor using the newly formulated density functional theory with an empirical dispersion term (DFT-D) molecular electronic structure method. It is found that the MR agonists, aldosterone and deoxycorticosterone, form tight hydrogen bonds with residues Thr945 and Asn770, which leads to the formation of hydrogen bond networks near the steroid D-ring, allowing for activation of this transcription factor. Progesterone, an MR antagonist, fails to form the necessary hydrogen bonds near the steroid D-ring. Progesterone is known to be an agonist of the mutated S810L MR receptor. Our studies indicate that this is possible because of a strong hydrogen bond between progesterone and Thr945 and a relatively strong hydrophobic interaction between progesterone and Asn770.
- MeSH
- financování organizované MeSH
- hormony chemie metabolismus MeSH
- leucin genetika metabolismus MeSH
- ligandy MeSH
- molekulární modely MeSH
- mutace genetika MeSH
- receptory mineralokortikoidů genetika chemie metabolismus MeSH
- serin genetika metabolismus MeSH
- steroidy chemie MeSH
- terciární struktura proteinů MeSH
- vodíková vazba MeSH
The geometries of a 13 mer of a DNA double helix (5'-GCGTACACATGCG-3') were determined by molecular dynamics simulations using a Cornell et al. empirical force field. The bases in the central base pair (shown in bold) were replaced (one or both) by a series of hydrophobic base analogues (phenyl, biphenyl, phenylnaphathalene, phenylanthracene and phenylphenanthrene). Due to the large fluctuations of the systems, an average geometry could not be determined. The interaction energies of the Model A, which consisted of three central steps of a duplex without a sugar phosphate backbone, taken from molecular dynamics simulations (geometry sampled every 1 ps), were calculated by the self-consistent charge density functional based tight-binding (SCC-DFTB-D) method and were subsequently averaged. The higher the stability of the systems the higher the aromaticity of the base analogues. To estimate the desolvation energy of the duplex, the COSMO continuum solvent model was used and the calculations were provided on a larger model, Model B (the three central steps of the duplex with a sugar phosphate backbone neutralised by H atoms), taken from molecular dynamics simulations (geometry sampled every 200 ps) and subsequently averaged. The selectivity of the base analogue pairs was ascertained (Model B) by including the desolvation energy and the interaction energy of both strands, as determined by the SCC-DFTB-D method. The highest selectivity was found for a phenylphenanthrene. Replacing the nucleic acid bases with a base analogue leads to structural changes of the central pair. Only with the smallest base analogues (phenyl) does the central base pair stay planar. When passing to larger base analogues the central base pair is usually stacked.
... Adhere to One Another 38 -- Molecular Complementarity Mediated via Noncovalent Interactions Permits Tight ... ... , Highly Specific Binding of Biomolecules 39 -- ŁQ[ Chemical Building Blocks of Cells 40 -- Amino Acids ... ... | Protein Function 78 -- Specific Binding of Ligands Underlies the -- Functions of Most Proteins 78 - ... ... - Enzymes Are Highly Efficient and Specific Catalysts 79 -- An Enzyme\'s Active Site Binds Substrates ... ... 118 -- I Transcription of Protein-Coding Genes and Formation of Functional mRNA 120 -- A Template DNA ...
6th ed. xxxvii, 1150 s. : il., tab. ; 29 cm
- MeSH
- biologie buňky MeSH
- molekulární biologie MeSH
- Publikační typ
- monografie MeSH
- Konspekt
- Biochemie. Molekulární biologie. Biofyzika
- NLK Obory
- biologie
- cytologie, klinická cytologie
