Toward on-the-fly quantum mechanical/molecular mechanical (QM/MM) docking: development and benchmark of a scoring function.

We address the challenges of treating polarization and covalent interactions in docking by developing a hybrid quantum mechanical/molecular mechanical (QM/MM) scoring function based on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. To benchmark this scoring function within the EADock DSS docking algorithm, we created a publicly available dataset of high-quality X-ray structures of zinc metalloproteins ( http://www.molecular-modelling.ch/resources.php ). For zinc-bound ligands (226 complexes), the QM/MM scoring yielded a substantially improved success rate compared to the classical scoring function (77.0% vs 61.5%), while, for allosteric ligands (55 complexes), the success rate remained constant (49.1%). The QM/MM scoring significantly improved the detection of correct zinc-binding geometries and improved the docking success rate by more than 20% for several important drug targets. The performance of both the classical and the QM/MM scoring functions compare favorably to the performance of AutoDock4, AutoDock4Zn, and AutoDock Vina.

Estudo ab-initio da a-alanina em meio aquoso

Ab initio Hartree-Fock (HF), Density Functional (B3LYP) and electron correlation (MP2) methods have been used to caracterize the aqueous medium intramolecular hydrogen bond in a-alanine. The 6-31G* and 6-31++G** were taken from Gaussian94 library. We were concerned on the structure of three conformers of a-alanine, in their neutral form plus on the structure of the zwitterionic form (Z). The Z structure is a stationary point at the HF/6-31G* level but it is not when diffuse functions and electron correlation are included. This results shows that the Z form does not exist in the gas phase. The inclusion of solvent effects changed significantly the results obtained in gas phase, therefore this inclusion make the Z form a stationary point within all level of theory, and the relative energy depends dramatically on the level of calculation.

Estudo teórico de propriedades ópticas não-lineares de nanotubos de carbono de parede única quimicamente modificados

Structure and first hyperpolarizability for a series of armchair a(5,5) chemically modified carbon nanotubes (CNT) were calculated at semiempirical and density functional levels of theory. The 4,4´-substituted stilbenes were selected as chromophore with substituents at position 4´ set to X=NO2, H, Cl, OH and NH2. The calculated values for static first hyperpolarizability (β) were almost linearly dependent on the electronic effect of the group X, increasing from NO2 to NH2. At DFT level the effect of inserting the chromophore in the CNT surface was to enhance the β value up to 70% relative to the free 4,4´-substituted stilbene.

Etudes ab initio des effets de la température sur le spectre optique des semi-conducteurs

Thèse réalisée en cotutelle avec l'Université Catholique de Louvain (Belgique)

Structure and energetics of mixed 4He-3He drops

Using a finite-range density functional, we have investigated the energetics and structural features of mixed helium clusters. The possibility of doping the cluster with a molecule of sulfur hexafluoride is also considered. It is seen that the repulsion introduced by the impurity strongly modifies the properties of the smallest drops. Although only a qualitative comparison is possible, the gross features displayed by our calculations are in agreement with recent experimental findings.

Effective t-J model Hamiltonian parameters of monolayered cuprate superconductors from ab initio electronic structure calculations

The magnetic coupling constant of selected cuprate superconductor parent compounds has been determined by means of embedded cluster model and periodic calculations carried out at the same level of theory. The agreement between both approaches validates the cluster model. This model is subsequently employed in state-of-the-art configuration interaction calculations aimed to obtain accurate values of the magnetic coupling constant and hopping integral for a series of superconducting cuprates. Likewise, a systematic study of the performance of different ab initio explicitly correlated wave function methods and of several density functional approaches is presented. The accurate determination of the parameters of the t-J Hamiltonian has several consequences. First, it suggests that the appearance of high-Tc superconductivity in existing monolayered cuprates occurs with J/t in the 0.20¿0.35 regime. Second, J/t=0.20 is predicted to be the threshold for the existence of superconductivity and, third, a simple and accurate relationship between the critical temperatures at optimum doping and these parameters is found. However, this quantitative electronic structure versus Tc relationship is only found when both J and t are obtained at the most accurate level of theory.

Radiative efficiencies and global warming potentials using theoretically determined absorption cross-sections for several hydrofluoroethers (HFEs) and hydrofluoropolyethers (HFPEs)

Integrated infrared cross-sections and wavenumber positions for the vibrational modes of a range of hydrofluoroethers (HFEs) and hydrofluoropolyethers (HFPEs) have been calculated. Spectra were determined using a density functional method with an empirically derived correction for the wavenumbers of band positions. Radiative efficiencies (REs) were determined using the Pinnock et al. method and were used with atmospheric lifetimes from the literature to determine global warming potentials (GWPs). For the HFEs and the majority of the molecules in the HG series HFPEs, theoretically determined absorption cross-sections and REs lie within ca. 10% of those determined using measured spectra. For the larger molecules in the HG series and the HG′ series of HFPEs, agreement is less good, with theoretical values for the integrated cross-sections being up to 35% higher than the experimental values; REs are up to 45% higher. Our method gives better results than previous theoretical approaches, because of the level of theory chosen and, for REs, because an empirical wavenumber correction derived for perfluorocarbons is effective in predicting the positions of C–F stretching frequencies at around 1250 cm−1 for the molecules considered here.

Electronic structure of Pd multi-layers on Re(0001): the role of charge transfer

Understanding the origin of the properties of metal-supported metal thin films is important for the rational design of bimetallic catalysts and other applications, but it is generally difficult to separate effects related to strain from those arising from interface interactions. Here we use density functional (DFT) theory to examine the structure and electronic behavior of few-layer palladium films on the rhenium (0001) surface, where there is negligible interfacial strain and therefore other effects can be isolated. Our DFT calculations predict stacking sequences and interlayer separations in excellent agreement with quantitative low-energy electron diffraction experiments. By theoretically simulating the Pd core-level X-ray photoemission spectra (XPS) of the films, we are able to interpret and assign the basic features of both low-resolution and high-resolution XPS measurements. The core levels at the interface shift to more negative energies, rigidly following the shifts in the same direction of the valence d-band center. We demonstrate that the valence band shift at the interface is caused by charge transfer from Re to Pd, which occurs mainly to valence states of hybridized s-p character rather than to the Pd d-band. Since the d-band filling is roughly constant, there is a correlation between the d-band center shift and its bandwidth. The resulting effect of this charge transfer on the valence d-band is thus analogous to the application of a lateral compressive strain on the adlayers. Our analysis suggests that charge transfer should be considered when describing the origin of core and valence band shifts in other metal / metal adlayer systems.

Tightened Lieb-Oxford Bound for Systems of Fixed Particle Number

The Lieb-Oxford bound is a constraint upon approximate exchange-correlation functionals. We explore a nonempirical tightening of that bound in both universal and electron number-dependent form. The test functional is PBE. Regarding both atomization energies (slightly worsened) and bond lengths (slightly improved), we find the PBE functional to be remarkably insensitive to the value of the Lieb-Oxford bound. This both rationalizes the use of the original Lieb-Oxford constant in PBE and suggests that enhancement factors more sensitive to sharpened constraints await discovery.

Azo-hydrazone tautomerism in protonated aminoazobenzenes: Resonance Raman spectroscopy and quantum-chemical calculations

The protonation effect on the vibrational and electronic spectra of 4-aminoazobenzene and 4-(dimethylamino)azobenzene was investigated by resonance Raman spectroscopy, and the results were discussed on the basis of quantum-chemical calculations. Although this class of molecular systems has been investigated in the past concerning the azo-hydrazone tautomerism, the present work is the first to use CASSCF/CASPT2 calculations to unveil the structure of both tautomers as well the nature of the molecular orbitals involved in chromophoric moieties responsible for the resonance Raman enhancement patterns. More specifically both the resonance Raman and theoretical results show clearly that in the neutral species, the charge transfer transition involves mainly the azo moiety, whereas in the protonated forms there is a great difference, depending on the tautomer. In fact, for the azo tautomer the transition is similar to that observed in the corresponding neutral species, whereas in the hydrazone tautomer such a transition is much more delocalized due to the contribution of the quinoid structure. The characterization of protonated species and the understanding of the tautomerization mechanism are crucial for controlling molecular properties depending on the polarity and pH of the medium.

Structural and Electronic Properties of Dipyridamole and Derivatives

Quantum mechanical calculations at the B3LYP theory level, together with the 6-31G* basis set, were employed to obtain the energy, ionization potential, and polarizabilites for dipyridamole and derivatives, which are compared with their biological activity. Density functional calculations of the spin densities were performed for radical formed by electron abstraction of dipyridamole and derivatives. The unpaired electron remains in dipyridamole is localized on the nitrogen atoms in the substituent positions 1, 3, 5, 7, 11, 12, 13, 14, with participation of the 9 and 10 carbons in the pyrimido-pyrimidine ring. The antioxidant activity is related with ionization potential, polarizability and Log P.

Estudo ab-initio da a-alanina em meio aquoso

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelagem em simulação computacional do cristal e superfícies do BaZr'O IND. 3' e SrZrO IND. 3': propriedades eletrônicas e estruturais

Pós-graduação em Ciência e Tecnologia de Materiais - FC

Síntese, caracterização e estudo do comportamento térmico dos oxanatos de lantanídeos(III) e de ítrio(III), no estado sólido

Pós-graduação em Química - IQ

Simulação atomística de sistema em nanoescala

Pós-graduação em Física - IGCE