Molecular mechanics (MM4) study of amines

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

On the Relevance of Considering the Intermolecular Interactions on the Prediction of the Vibrational Spectra of Isopropylamine

The effects of implicitly considering the effects of hydrogen bonding on the molecular properties, such as vibrational frequencies, were inferred on the basis of DFT calculations. Several clusters of isopropylamine were assembled and theoretically characterized. The results showed that maximum H-bond cooperativity is achieved when the amine group acts simultaneously as donor and acceptor. The effect of H-bond cooperativity manifests itself in the relative cluster stability and on the structural and vibrational frequency predictions. Referring to the vibrational frequencies it was found that theNH2stretching and torsion vibrational modes are the most affected by the amine involvement in hydrogen bonding. Both stretching modes were found to be significantly redshifted relative to the monomer. TheNH2torsional mode, on the other hand, was found to be blueshifted up to 350 cm-1. Finally, the comparative study between the theory levels performed allows to conclude that the small 6-31G* basis set is able to stabilize weakC–H⋯Ninteractions as long as the new dispersion corrected DFT methods are considered. The impairments observed with conventional DFT methods for describing weak interactions may be overcome with the improvement of basis set, but the associated increase of computational costs may turn the calculations unfeasible.

Spectroscopic properties of neuroleptics: IR and Raman spectra of Risperidone (Risperdal) and of its mono- and di-protonated forms

Structures and IR and Raman spectra of Risperidone in its neutral, mono- and di-protonated forms were calculated in gas phase by DFT-B3LYP/6-31G* level. Mono-protonation occurs at the nitrogen atom of the piperidine ring, while nitrogen atom of the pyrimidine ring is the preferred site for the second protonation. The lowest-energy structure of the mono-protonated Risperidone is characterized by formation of a strong seven-membered O(pyrimidine ring)⋯(+)H-N(piperidine ring) intramolecular hydrogen-bonded cycle. In the high-energy spectral region (3500-2500 cm(-1)), the bands of the N-H(+) stretches and the changes in wavenumbers and IR intensities of the C-H stretches near to the piperidine nitrogen atom (Bohlmann effect) are potentially useful to discriminate conformations and protonation states. Di-protonated structures can be identified by the presence of an isolated absorption peak located in the low-energy IR region (660-690 cm(-1)), attributed to the out-of-plane N-H(+)(pyrimidine ring) bending deformation. The most intense Raman band of neutral Risperidone placed at ca. 1500 cm(-1), assigned to C=C(pyrimidine ring) stretch + C=N(pyrimidine ring) stretch, can be a useful vibrational marker to distinguish the neutral from the protonated forms.

The important role of lone-pairs in force field (MM4) calculations on hydrogen bonding in alcohols

An expanded treatment of hydrogen bonding has been developed for MM4 force field calculations, which is an extension from the traditional van der Waals-electrostatic model. It adds explicit hydrogen-bond angularity by the inclusion of lone-pair directionality. The vectors that account for this directionality are placed along the hydrogen acceptor and its chemically intuitive electron pairs. No physical lone-pairs are used in the calculations. Instead, an H-bond angularity function, and a lone-pair directionality function, are incorporated into the hydrogen-bond term. The inclusion of the lone-pair directionality results in improved accuracy in hydrogen-bonded geometries and interaction energies. In this work is described hydrogen bonding in alcohols, and also in water and hydrogen fluoride dimer. The extension to other compounds such as aldehydes, ketones, amides, and so on is straightforward and will be discussed in future work. The conformational energies of ethylene glycol are discussed.