RESUMO

4'-substituted neutral/protonated furfurylidenanilines and trans-styrylfurans are able to exist in two different conformations related to the rotation around the furan ring-bridge double bond. In this work, the equilibrium geometry and the corresponding rotational barrier of the benzene ring for each furan derivative conformation were calculated by DFT methods. The trend and shape of the rotational barrier are rationalized within natural bond orbitals as well as atoms-in-molecules approach. For the corresponding equilibrium geometries, (1)H and (13)C substituent induced shifts (SIS) were calculated and compared with experimental values. Calculated shielding constants are shown to be sensitive to the substituent effect through a linear fit with substituent's Hammett constants. An alternative approach was followed for assessing the effect of substituents over SIS through comparing the differences in isotropic shielding constants with NBO charges as well as with (1)H and (13)C experimental chemical shifts.

RESUMO

Retinitis pigmentosa (RP) is a pathological condition associated with blindness due to progressive retinal degeneration. RP-linked mutations lead to changes at the retinal binding pocket and in the absorption spectra. Here, we evaluate the geometries, electronic effects, and vertical excitation energies in the dark state of mutated human rhodopsins carrying the abnormal substitutions M207R or S186W at the retinal binding pocket. Two models are used, the solvated protein and the protein in a solvated POPC lipid bilayer. We apply homology modeling, classical molecular dynamics simulations, density functional theory (DFT), and quantum mechanical/molecular mechanical (QM/MM) methods. Our results for the wild type bovine and human rhodopsins, used as a reference, are in good agreement with experiment. For the mutants, we find less twisted QM/MM ground-state chromophore geometries around the C(11)-C(12) double bond and substantial blue shifts in the lowest vertical DFT excitation energies. An analysis of the QM energies shows that the chromophore-counterion region is less stable in the mutants compared to the wild type, consistent with recent protein folding studies. The influence of the mutations near the chromophore is discussed in detail to gain more insight into the properties of these mutants. The spectral tuning is mainly associated with counterion effects and structural features of the retinal chain in the case of the hM207R mutant, and with the presence of a neutral chromophore with deprotonated Lys296 in the case of the hS186W mutant.

Assuntos

RESUMO

The behavior of O(2) molecule in models of acid aluminosilicate sites on any kind of material was investigated using reliable QM ab initio calculations. The triplet-singlet energy gap of isolated O(2) was calculated at confident levels of theory with different basis sets as a reference. Models of aluminosilicate active sites interacting with oxygen in their singlet and triplet electronic states were considered for two kinds of O(2) arrangements. Geometry optimizations were performed on both non-corrected and corrected BSSE potential energy surfaces, realizing that good modeling of heavy atom-hydrogen interactions is sensitive to BSSE corrections during these processes. Energies were further evaluated at higher level of theory to test tendencies. Singlet oxygen appears more attractive to active aluminosilicate sites than those calculated with triplet oxygen, indicating a source of oxidative efficiency for designed nanostructures containing such molecular residues. It was clearly seen that aluminosilicate groups, appearing ubiquitously in several materials, could reduce the O(2) triplet-singlets energy gap by at least 10 kJ/mol. Some elegant features of oxygen interactions with such sites were further analyzed by means of the atoms in molecules (AIM) theory.

Assuntos

RESUMO

A set of weak complexes between imidazole (Imi) and nitric oxide (NO) were calculated with UMP2/6-31++G(d,p) and UMP2/6-311++G(2d,2p) levels of theory. Planar and nonplanar geometries were considered. Complexes of NO and protonated imidazole (ImiH(+)) were also studied due to the biological important effect of Imi protonation. It was found that ring protonation increases the stability of planar complexes and does not affect significantly nonplanar complexes. The Z-H...XY (Z = C, N and X, Y = N, O) interactions resulted as hydrogen bonds according to Koch and Popelier criteria based on AIM theory. Charge transfer was also found very important for complex stabilization within our theoretical framework. Planar NO...ImiH(+) complexes are more stable than those obtained with neutral Imi.

RESUMO

Predictive quantitative structure-activity relationship (QSAR) models of anabolic and androgenic activities for the testosterone and dihydrotestosterone steroid analogues were obtained by means of multiple linear regression using quantum and physicochemical molecular descriptors (MD) as well as a genetic algorithm for the selection of the best subset of variables. Quantitative models found for describing the anabolic (androgenic) activity are significant from a statistical point of view: R(2) of 0.84 (0.72 and 0.70). A leave-one-out cross-validation procedure revealed that the regression models had a fairly good predictability [q(2) of 0.80 (0.60 and 0.59)]. In addition, other QSAR models were developed to predict anabolic/androgenic (A/A) ratios and the best regression equation explains 68% of the variance for the experimental values of AA ratio and has a rather adequate q(2) of 0.51. External validation, by using test sets, was also used in each experiment in order to evaluate the predictive power of the obtained models. The result shows that these QSARs have quite good predictive abilities (R(2) of 0.90, 0.72 (0.55), and 0.53) for anabolic activity, androgenic activity, and A/A ratios, respectively. Last, a Williams plot was used in order to define the domain of applicability of the models as a squared area within +/-2 band for residuals and a leverage threshold of h=0.16. No apparent outliers were detected and the models can be used with high accuracy in this applicability domain. MDs included in our QSAR models allow the structural interpretation of the biological process, evidencing the main role of the shape of molecules, hydrophobicity, and electronic properties. Attempts were made to include lipophilicity (octanol-water partition coefficient (logP)) and electronic (hardness (eta)) values of the whole molecules in the multivariate relations. It was found from the study that the logP of molecules has positive contribution to the anabolic and androgenic activities and high values of eta produce unfavorable effects. The found MDs can also be efficiently used in similarity studies based on cluster analysis. Our model for the anabolic/androgenic ratio (expressed by weight of levator ani muscle, LA, and seminal vesicle, SV, in mice) predicts that the 2-aminomethylene-17alpha-methyl-17beta-hydroxy-5alpha-androstan-3-one (43) compound is the most potent anabolic steroid, and the 17alpha-methyl-2beta,17beta-dihydroxy-5alpha-androstane (31) compound is the least potent one of this series. The approach described in this report is an alternative for the discovery and optimization of leading anabolic compounds among steroids and analogues. It also gives an important role to electron exchange terms of molecular interactions to this kind of steroid activity.

Assuntos

RESUMO

A remarkable difference for (2)J(C(2)-H(f)) coupling constant in syn and anti conformers of 5-X-furan-2-carboxaldehydes (X = CH(3), Ph, NO(2), Br) and a rationalization of this difference are reported. On the basis of the current knowledge of the Fermi-contact term transmission, a rather unusual dual-coupling pathway in the syn conformer is presented. The additional coupling pathway resembles somewhat that of the J(H-H) in homoallylic couplings, which are transmitted by hyperconjugative interactions involving the pi(C=C) electronic system. The homoallylic coupling pathway can be labeled as sigma*(C-H) <-- pi(C=C) --> sigma*(C-H). In the present case, this additional coupling pathway, using an analogous notation, can be labeled as sigma*(C(2)-C(C)) <-- LP(1)(O(1))...LP(2)(O(C)) --> sigma*(C(C)-H(f)) (sigma*(C(2)-C(C))) where O(1) and O(C) stand for the ring and carbonyl O atoms, respectively. This additional coupling pathway is not activated in the anti conformers since both oxygen lone pairs do not overlap.

Assuntos

RESUMO

The great cost associated with the development of new anabolic-androgenic steroid (AASs) makes necessary the development of computational methods that shorten the drug discovery pipeline. Toward this end, quantum, and physicochemical molecular descriptors, plus linear discriminant analysis (LDA) were used to analyze the anabolic/androgenic activity of structurally diverse steroids and to discover novel AASs, as well as also to give a structural interpretation of their anabolic-androgenic ratio (AAR). The obtained models are able to correctly classify 91.67% (86.27%) of the AASs in the training (test) sets, respectively. The results of predictions on the 10% full-out cross-validation test also evidence the robustness of the obtained model. Moreover, these classification functions are applied to an "in house" library of chemicals, to find novel AASs. Two new AASs are synthesized and tested for in vivo activity. Although both AASs are less active than some commercially AASs, this result leaves a door open to a virtual variational study of the structure of the two compounds, to improve their biological activity. The LDA-assisted QSAR models presented here, could significantly reduce the number of synthesized and tested AASs, as well as could increase the chance of finding new chemical entities with higher AAR.

Assuntos

RESUMO

Very large molecular systems can be calculated with the so called CNDOL approximate Hamiltonians that have been developed by avoiding oversimplifications and only using a priori parameters and formulas from the simpler NDO methods. A new diagonal monoelectronic term named CNDOL/21 shows great consistency and easier SCF convergence when used together with an appropriate function for charge repulsion energies that is derived from traditional formulas. It is possible to obtain a priori molecular orbitals and electron excitation properties after the configuration interaction of single excited determinants with reliability, maintaining interpretative possibilities even being a simplified Hamiltonian. Tests with some unequivocal gas phase maxima of simple molecules (benzene, furfural, acetaldehyde, hexyl alcohol, methyl amine, 2,5 dimethyl 2,4 hexadiene, and ethyl sulfide) ratify the general quality of this approach in comparison with other methods. The calculation of large systems as porphine in gas phase and a model of the complete retinal binding pocket in rhodopsin with 622 basis functions on 280 atoms at the quantum mechanical level show reliability leading to a resulting first allowed transition in 483 nm, very similar to the known experimental value of 500 nm of "dark state." In this very important case, our model gives a central role in this excitation to a charge transfer from the neighboring Glu(-) counterion to the retinaldehyde polyene chain. Tests with gas phase maxima of some important molecules corroborate the reliability of CNDOL/2 Hamiltonians.

Assuntos

RESUMO

The potential energy surface of the CH(4)...NO van der Waals complexes was explored at the RCCSD(T)/aug-cc-pVTZ level including the full counterpoise correction to the basis set superposition error. The Jahn-Teller distortion of the C(3v) configurations for the CH bonded and CH(3) face complexes was analyzed. From this distortion, two A(') and A(") adiabatic surfaces were considered. The estimated zero point energy of C(s) configurations is above the barrier of the C(3v) ones. Therefore, the CH(3) face complexes are dynamic Jahn-Teller systems. The D(0) (140 cm(-1) for A(") state and 100 cm(-1) for A(')) values obtained are in good agreement with the experimental values (103+/-2 cm(-1)) recently reported.

RESUMO

A systematic study of basis set superposition error (BSSE) behavior in H(3)C-H[ellipsis (horizontal)][NO] complexes for both -H...N- and -H...O- orientations were carried out using MP2 and density-functional theory with Pople's [6-31G(d,p),6-311++G(nd,nd), where n=1,2,3, and 6-311++G(3df,3pd)] and Dunning's augmented correlation consistent basis sets [aug-cc-pVXZ (X=D and T)]. Corrected and uncorrected counterpoise potential-energy surfaces (PESs) were explored and differences obtained between them indicate that reliable optimizations of these molecular interactions must be carried out in a PES free of BSSE, even in the case of large basis sets and popularly used functionals such as B3LYP. Although all basis used could be always considered within a margin of approximation for representing molecular orbitals and show important values of BSSE, 6-311++G(2d,2p) basis set shows the best results in uncorrected PES with respect to the corrected ones. B3LYP functional produces erratic results: complexes appear repulsive and the intermolecular distances are always large, evidencing the lack of a correct dispersive forces treatment in the original parameterization. According to the MP2 results, the -H...N- interactions appear as slightly more stable than those of the -H...O- orientation.