The platinum-olefin binding energy in series of (PH<Subscript>3</Subscript>)<Subscript>2</Subscript>Pt(olefin) complexes - a theoretical study

Theoretical investigation of Pt(0)-olefin organometallic complexes containing tertiary phosphine ligands was focused on the strength of platinum-olefin electronic interaction. DFT theoretical study of electronic effects in a substantial number of ethylene derivatives was evaluated in terms of the Pt-olefin binding energy using MP2 correlation theory. Organometallics bearing coordinated olefins with general formula (R¹R²C = CR³R⁴)Pt(PH₃)₂ [R = various substituents] had been selected, including olefins containing both electron-donor substituents as well as electron-withdrawing groups. The stability of the corresponding complexes increases with a strengthening electron-withdrawal ability of the olefin substituents. Figure Representation of (CH₂ = CHR)Pt(PPh₃)₂ and the stability chart     

Theoretical study of the novel sandwich compound [Au<Subscript>3</Subscript>Cl<Subscript>3</Subscript>Tr<Subscript>2</Subscript>]<Superscript>2+</Superscript>

A theoretical study of a sandwich compound with a metal monolayer sheet between two aromatic ligands is presented. A full geometry optimization of the [Au₃Cl₃Tr₂]²⁺ (1) compound, which is a triangular gold(I) monolayer sheet capped by chlorines and bounded to two cycloheptatrienyl (Tr) ligands was carried out using perturbation theory at the MP2 computational level and DFT. Compound (1) is in agreement with the 18–electron rule, the bonding nature in the complex may be interpreted from the donation interaction coming from the Tr rings to the Au array, and from the back-donation from the latter to the former. NICS calculations show a strong aromatic character in the gold monolayer sheet and Tr ligands; calculations done with HOMA, also report the same aromatic behavior on the cycloheptatrienyl fragments giving us an insight on the stability of (1). The Au –Au bond lengths indicate that an intramolecular aurophilic interaction among the Au(I) cations plays an important role in the bonding of the central metal sheet. Figure (a) Ground state geometry of complex 1; (b) Top view of compound 1 and Wiberg bond orders computed with the MP2/B1 computational method; (c) Lateral view of compound 1 and NICS values calculated with the MP2/B1 method; the values in parenthesis were obtained at the VWN/TZP level     

Structure-based 3D-QSAR studies on thiazoles as 5-HT<Subscript>3</Subscript> receptor antagonists

Structure-based 3D-QSAR studies were performed on 20 thiazoles against their binding affinities to the 5-HT₃ receptor with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The thiazoles were initially docked into the binding pocket of a human 5-HT_3A receptor homology model, constructed on the basis of the crystal structure of the snail acetylcholine binding protein (AChBP), using the GOLD program. The docked conformations were then extracted and used to build the 3D-QSAR models, with cross-validated values 0.785 and 0.744 for CoMFA and CoMSIA, respectively. An additional five molecules were used to validate the models further, giving satisfactory predictive values of 0.582 and 0.804 for CoMFA and CoMSIA, respectively. The results would be helpful for the discovery of new potent and selective 5-HT₃ receptor antagonists.      

Molecular modeling and biophysical analysis of the c-MYC NHE-III<Subscript>1</Subscript> silencer element

G-Quadruplex and i-Motif-forming sequences in the promoter regions of several oncogenes show promise as targets for the regulation of oncogenes. In this study, molecular models were created for the c-MYC NHE-III₁ (nuclease hypersensitivity element III₁) from two 39-base complementary sequences. The NHE modeled here consists of single folded conformers of the polypurine intramolecular G-Quadruplex and the polypyrimidine intramolecular i-Motif structures, flanked by short duplex DNA sequences. The G-Quadruplex was based on published NMR structural data for the c-MYC 1:2:1 loop isomer. The i-Motif structure is theoretical (with five cytosine–cytosine pairs), where the central intercalated cytosine core interactions are based on NMR structural data obtained for a tetramolecular [d(A₂C₄)₄] model i-Motif. The loop structures are in silico predictions of the c-MYC i-motif loops. The porphyrin meso-tetra(N-methyl-4-pyridyl)porphine (TMPyP4), as well as the ortho and meta analogs TMPyP2 and TMPyP3, were docked to six different locations in the complete c-MYC NHE. Comparisons are made for drug binding to the NHE and the isolated G-Quadruplex and i-Motif structures. NHE models both with and without bound cationic porphyrin were simulated for 100 ps using molecular dynamics techniques, and the non-bonded interaction energies between the DNA and porphyrins calculated for all of the docking interactions. Figure Molecular models of the average structure of the final 20 ps of the molecular dynamics simulation of the c-MYC NHE-III₁ (nuclease hypersensitivity element III₁) “silencer” element. The G-Quadruplex structure is at the top-center, and the i-Motif is at the bottom-center of each picture. a “Rotation #1” of the G-Quadruplex, with the T15 loop at the top and rear and the G19/A20 loop at the top and front of the picture. b “Rotation #2” of the G-Quadruplex, with the T15 loop at the top and front of the image, and the G19/A20 loop at the front and adjacent to the G-Quadruplex/i-Motif interface     

G2, G3, and complete basis set calculations on the thermodynamic properties of triazane

As a follow-up study to our study on tetrazane (N₄H₆), we present computed thermodynamic properties of triazane (N₃H₅). Calculated properties include optimized geometries, infrared vibrations, enthalpy of formation, enthalpy of combustion, and proton affinities. We have also mapped the potential energy surface as the molecule is rotated about the N-N bond. We have predicted a specific enthalpy of combustion for triazane of about -20 kJ g⁻¹. Figure Schematic diagram of the dielectric barrier discharge (left) and typical temporal profiles of voltage and current, as obtained from the simulations (right)     

Computational studies of the binding site of α<Subscript>1A</Subscript>-adrenoceptor antagonists

Aimed at achieving a good understanding of the 3-dimensional structures of human α_1A-adrenoceptor (α_1A-AR), we have successfully developed its homology model based on the crystal structure of β₂-AR. Subsequent structural refinements were performed to mimic the receptor’s natural membrane environment by using molecular mechanics (MM) and molecular dynamics (MD) simulations in the GBSW implicit membrane model. Through molecular docking and further simulations, possible binding modes of subtype-selective α_1A-AR antagonists, Silodosin, RWJ-69736 and (+)SNAP-7915, were examined. Results of the modeling and docking studies are qualitatively consistent with available experimental data from mutagenesis studies. The homology model built should be very useful for designing more potent subtype-selective α_1A-AR antagonists and for guiding further mutagenesis studies. Figure The superposition of β₂-AR crystal structure (gold ribbons) and α_1A-AR homology model (blue ribbons)     

Aromatic C20F20 cage and its endohedral complexes X@C20F20 (X = H-, F-, Cl-, Br-, H, He)

The structure and stability of endohedral X@C₂₀F₂₀ complexes (X = H⁻, F⁻, Cl⁻, Br⁻, H, He) have been computed at the B3LYP level of theory. All complexes in I _h symmetry were found to be energy minimum structures. H⁻@C₂₀F₂₀ and F⁻@C₂₀F₂₀ complexes have negative inclusion energies, while other complexes have positive inclusion energies. Similarity between C₂₀F₂₀ and C₂₀H₂₀ has been found for X = H and He. On the basis of the computed nucleus independent chemical shift values at the cage center, both C₂₀F₂₀ and C₂₀F₂₀ are aromatic. Figure Endohedral X@C₂₀F₂₀ complexes     

A molecular-dynamics simulation study of diffusion of a single model carbonic chain on a graphite (001) surface

Molecular-dynamics simulations have been used to study the diffusion of a short single model carbonic chain on the graphite (001) surface. The calculated diffusion coefficient (D) first increases, then decreases with increasing chain length (N). This abnormal behavior is similar to polymer lateral diffusion at the solid–liquid interface. Furthermore, we have studied the relation between the mean-square gyration radius and N. Figure Log–log plot of the self-diffusion coefficient D versus the chain length N. The error bars are the standard deviation measured in three repeated simulations     

Exploring the binding features of rimonabant analogues and acyclic CB<Subscript>1</Subscript> antagonists: docking studies and QSAR analysis

In order to elucidate the structural requirements for human CB₁ receptor antagonism, 78 antagonists belonging to five different chemical classes were selected from the literature and docked into the receptor binding site, built by homology modeling techniques. To further explore the structure-activity relationships within the considered chemical classes, a pharmacophore model and a QSAR analysis were developed. In a first step five alignments, one for each group of compounds were generated. All of them were then submitted to a MOE pharmacophore search in order to obtain a final pharmacophore model representative of the whole dataset which was used to elaborate the following 3D-QSAR analysis, by means of the CoMFA methodology. The results of these investigations are expected to be useful in the process of design and development of new potent CB₁ antagonists. Figure Compounds 1-78 are aligned into the putative CB1 receptor binding site. The three key features shared by all of them are reported in coloured spheres. The hydrophobic/aromatic ones are depicted in purple while the acceptor functions are coloured in blue.     

G2, G3, and complete basis set calculations of the thermodynamic properties of cis- and trans-triazene

Following our recent study on triazane, we present a follow-up study on the thermodynamic properties of triazane’s unsaturated analog, triazene. We predict optimized structural parameters, vibrational frequencies, enthalpies of formation, enthalpies of combustion, specific enthalpies of combustion, and proton affinities. Our results indicate that the cis form of triazene has a specific enthalpy of combustion of −15.2 kJ g⁻¹ and the trans form has a specific enthalpy of combustion of −14.7 kJ g⁻¹. Figure Structures of cis- and trans-triazane, N₃H₃     

Theoretical study of [ XN5 ]− (X=O, S, Se, Te) systems

A series of [XN₅]⁻ (X=O, S, Se, Te) compounds has been examined with ab initio and Density Functional Theory (DFT) methods. The five-membered nitrogen ring series of structures are global minima and may exist or be characterized due to their significant dissociation barriers (29.7–32.7 kcal mol⁻¹). Nucleus-independent chemical shifts (NICS) criteria and the presence of (4n+2) π-electrons confirmed that the five-membered nitrogen ring in their structures exhibits characteristics of aromaticity. Thus, the strong stability of the five-membered nitrogen ring structures may be attributed partially to their aromaticity.      

Ligation of Aza bases to the AgF<Subscript>2</Subscript> molecule: a theoretical study

We have analyzed the electronic structure and chemical bonding for molecular adducts of the Ag(II)F₂ molecule with various aza Lewis bases including ammonia, nitriles, secondary amines, and their derivatives exhibiting various degrees of fluorination. Density functional theory calculations indicate that a progressive shift occurs of the spin density from the Ag center towards the coordinating nitrogen atoms of aza ligands, as the ligation energy increases. Chemistry of Ag(II) might be extended with little effort beyond the known aza connections, to include nitriles, perfluorinated nitriles and perfluorinated amines. Figure Properties of a variety of novel adducts of the AgF₂ molecule with two aza bases (L), possible precursors of the AgF₂L₂ extended solids, were assessed by the DFT calculations Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. This work is dedicated to memory of Wojciech Ochmański, unforgettable person, good-hearted man, whose craftsmanship in work was second-to-none.     

Conformational stability and three-dimensional model of the δ-opioid pharmacophore for the extended antiparallel dimer structure of Met-enkephalin in water

The conformational stability of the extended antiparallel dimer structure of Met-enkephalin in water was analyzed by examining the hydration structure of enkephalin using molecular dynamics simulations. The result shows that, despite of the hydrophicility of the terminal atoms in the pentapeptide, the main contributor for the stability of the dimer in water is the four intermolecular hydrogen bonds between the Gly² and Phe⁴ groups. The three-dimensional model of the δ-opioid pharmacophore for this dimer structure was also established. Such a model was demonstrated to match the δ-opioid pharmacophore query derived from the non-peptides SIOM, TAN-67, and OMI perfectly. This result thus strongly supports the assumption that the dimer structure of Met-enkephalin is a possible δ-receptor binding conformation. Figure Schematic model of the extended antiparallel dimer structure of Met-enkephalin     

Aromaticity and electronic properties of Heterosuperbenzene (Heterohexabenzocoronene)

A global electrophilicity parameter and the aromaticity of some heterocyclic polyaromatic hydrocarbons were evaluated on the basis of DFT calculations. The substitution of carbon atoms by nitrogen atoms dramatically changes the global electrophilicity of the molecules, with the fully substituted molecule being the most electrophilic with a reactivity very close to that of fullerene. Figure Fully substituted heterohexabenzocoronene Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Conformation and tautomerizm of the 2-methyl-4-pyridin-2′-yl-1,5-benzodiazepine molecule. An ab initio study

The protomeric tautomerizm and conformation of the 2-methyl-4-pyridin-2′-yl-1,5-benzodiazepine molecule were investigated, and its three neutral tautomers (B₁,B₂,B₃) and their rotamers (C₁,C₂,C₃) were considered. Full geometry optimizations were carried out at the HF/6-31G* and B3LYP/6-31G* levels in gas phase and in water. The tautomerization processes in water (ɛ = 78.54) were studied by using self-consistent reaction field theory. The calculation showed that the boat conformation is dominant for the seven-membered diazepine ring in all of the structures, even with different double bond positions. The calculated relative free energies (ΔG) showed that the tautomer C₁ was the most stable structure, and its conformer B₁ was the second most stable in the gas phase and in water. Figure 2-Methyl-4-pyridin-2′-yl-1,5-benzodiazepine     

Investigation of the intermolecular proton transfer in the supersystems adenine-methanol/ethanol/i-propanol: MP2 and DFT levels study

Twelve H-bonded supersystems constructed between the adenine tautomers and methanol, ethanol, and i-propanol were studied at the B3LYP and MP2 levels of theory using 6-311G(d,p) and 6-311++G(d,p) basis functions. The thermodynamic parameters of the complex formations were calculated in order to estimate the exact stability of the supersystems. It was proven that the calculated energy barriers of the alcohol-assisted proton transfers are about 60% lower than those of the intramolecular proton transfers in adenine found earlier (Gu and Leszczynski in J Phys Chem A 103:2744–2750, 1999). Figure H-bonded complex between i-propanol and adenine     

Molecular dynamics study of the initial stages of catalyzed single-wall carbon nanotubes growth: force field development

Effective force fields for Ni-C interactions developed by Yamaguchi and Maruyama for the formation of metallofullerenes are modified to simulate the catalyzed growth of single-wall carbon nanotubes on Ni_n clusters with n >20, and the reactive empirical bond order Brenner potential for C-C interactions is also revised to include the effect of the metal atoms on such interactions. Figure Force field parameters for carbon-metal interactions obtained from DFT calculations in small clusters.     

Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements

Several modifications that have been made to the NDDO core-core interaction term and to the method of parameter optimization are described. These changes have resulted in a more complete parameter optimization, called PM6, which has, in turn, allowed 70 elements to be parameterized. The average unsigned error (AUE) between calculated and reference heats of formation for 4,492 species was 8.0 kcal mol⁻¹. For the subset of 1,373 compounds involving only the elements H, C, N, O, F, P, S, Cl, and Br, the PM6 AUE was 4.4 kcal mol⁻¹. The equivalent AUE for other methods were: RM1: 5.0, B3LYP 6–31G*: 5.2, PM5: 5.7, PM3: 6.3, HF 6–31G*: 7.4, and AM1: 10.0 kcal mol⁻¹. Several long-standing faults in AM1 and PM3 have been corrected and significant improvements have been made in the prediction of geometries. Figure Calculated structure of the complex ion [Ta₆Cl₁₂]²⁺ (footnote): Reference value in parenthesis Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Oxidation mechanism in the metabolism of (S)-N-[1-(3-morpholin-4-ylphenyl)ethyl]-3-phenylacrylamide on oxyferryl active site in CYP3A4 Cytochrome: DFT modeling

The metabolism mechanism of (S)-N-[1-(3-morpholin-4ylphenyl)ethyl]-3-phenylacrylamide, mediated by CYP3A4 Cytochrome has been investigated by density functional QM calculations aided with molecular mechanics/molecular dynamics simulations. Two different orientations of phenyl ring for substrate approach toward oxyferryl center, imposing two subsequent rearrangement pathways have been investigated. Starting from σ-complex in perpendicular orientation enzymatic mechanism involves consecutive proton shuttle intermediate, which further leads to the formation of alcohol and ketone. Parallel conformation leads solely to ketone product by 1,2 hydride shift. Although parallel and perpendicular σ-complexes are energetically equivalent both for the gas phase or PCM solvent model, molecular dynamics studies in full CYP3A4 environment show that perpendicular conformation of the σ-complex should be privileged, stabilized by hydrophobic interactions of phenylacrylamide chain. After assessing probability of the two conformations we postulate that the alcohol, accessible with the lowest energy barriers should be the major metabolite for studied substrate and CYP3A4 enzyme. Figure Orientation of phenyl ring towards porphyrin plane selected by substrate interaction with enzymatic cavity channels enzymatic reaction     

Ab initio and DFT study of the aromaticity of some Fulvalenes derived from Methylidenecyclopropabenzene

Methylidencyclopropabenzene (MCPB) 1 and Fulvalenes 2–4 are molecules of special interest due to the relation between structure and aromaticity. The aim of this work was to analyze this relation and to quantify the aromaticity in 1–4 using different methods. Magnetic properties are directly related with aromaticity; here we studied the magnetic susceptibility and the anisotropy of the magnetic susceptibility. Nucleus indepedent chemical shift (NICS) and the anisotropy of the induced current density (ACID) were also employed. Tools of very different nature, geometric indexes HOMA and Bird, were determinated too for 1–4. All of these measures were found to be in agreement. Figure Both spatial NICS and ACID plot allow to show the aromaticity/antiaromaticity of a ring