A theoretical study of the fluorine valence shell in methyl fluoride

Molecular modeling study of dihydrofolate reductase inhibitors. Molecular dynamics simulations, quantum mechanical calculations, and experimental corroboration

A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

Topology of Charge Density of Flucytosine and Related Molecules and Characteristics of their Bond Charge Distributions

The molecular charge distribution of flucytosine (4-amino-5-fluoro-2-pyrimidone), uracil, 5-fluorouracil, and thymine was studied by means of density functional theory calculations (DFT). The resulting distributions were analyzed by means of the atoms in molecules (AIM) theory. Bonds were characterized through vectors formed with the charge density value, its Laplacian, and the bond ellipticity calculated at the bond critical point (BCP). Within each set of C=O, C-H, and N-H bonds, these vectors showed little dispersion. C-C bonds formed three different subsets, one with a significant degree of double bonding, a second corresponding to single bonds with a finite ellipticity produced by hyperconjugation, and a third one formed by a pure single bond. In N-C bonds, a decrease in bond length (an increase in double bond character) was not reflected as an increase in their ellipticity, as in all C-C bonds studied. It was also found that substitution influenced the N-C, C-O, and C-C bond ellipticity much more than density and its Laplacian at the BCP. The Laplacian of charge density pointed to the existence of both bonding and nonbonding maxima in the valence shell charge concentration of N, O, and F, while only bonding ones were found for the C atoms. The nonbonding maxima related to the sites for electrophilic attack and H bonding in O and N, while sites of nucleophilic attack were suggested by the holes in the valence shell of the C atoms of the carbonyl groups.