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.

Importance of the embedding environment on the strain within small rings in siliceous materials

The effect of the local environment on the energetic strain within small (SiO)N rings (with N=2,3) in silica materials is investigated via periodic model systems employing density functional calculations. Through comparison of the energies of various nonterminated systems containing small rings in strained and relatively unstrained environments, with alpha quartz, we demonstrate how small ring strain is affected by the nature of the embedding environment. We compare our findings with numerous previously reported calculations, often predicting significantly different small-ring strain energies, leading to a critical assessment of methods of calculating accurate localized ring energies. The results have relevance for estimates of the strain-induced response (e.g., chemical, photo, and radio) of small silica rings, and the propensity for them to form in bulk glasses, thin films, and nanoclusters.

Estudo da adsorção de hidrogênio e sulfeto na superfície de paládio: aspectos experimentais (eletroquímica) e teóricos (ab initio e Teoria do Funcional da Densidade)

The adsorption of H and S2- species on Pd (100) has been studied with ab initio, density-functional calculations and electrochemical methods. A cluster of five Pd atoms with a frozen geometry described the surface. The computational calculations were performed through the GAUSSIAN94 program, and the basis functions adapted to a pseudo-potential obtained by using the Generator Coordinate Method adapted to the this program. Using the cyclic voltammetry technique through a Model 283 Potentiostat/Galvanostat E.G.&G-PAR obtained the electrochemical results. The calculated chemisorption geometry has a Pd-H distance of 1.55Å, and the potential energy surface was calculated using the Becke3P86//(GCM/DFT/SBK) methodology. The adsorption of S2- ions on Pd surface obtained both through comparison between the experimental and theoretical results, at MP2 level, suggest a S2- absorption into the metallic cluster. The produced Pd-(S2-) system was show to be very stable under the employed experimental conditions. The paper has shows the powerful aid of computational methods to interpret adsorption experimental data.

Reduced ceria nanofilms from structure prediction

Experimentally, Ce2O3 films are used to study cerium oxide in its fully or partially reduced state, as present in many applications. We have explored the space of low energy Ce2O3 nanofilms using structure prediction and density functional calculations, yielding more than 30 distinct nanofilm structures. First, our results help to rationalize the roles of thermodynamics and kinetics in the preparation of reduced ceria nanofilms with different bulk crystalline structures (e.g. A-type or bixbyite) depending on the support used. Second, we predict a novel, as yet experimentally unresolved, nanofilm which has a structure that does not correspond to any previously reported bulk A2B3 phase and which has an energetic stability between that of A-type and bixbyite. To assist identification and fabrication of this new Ce2O3 nanofilm we calculate some observable properties and propose supports for its epitaxial growth.

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.

Optical spectroscopy of organic material based on C60

The structure and optical properties of thin films based on C60 materials are studied. Reproducible vacuum method of thin fullerene films production with Cd impurity on Si, glass and mica surfaces is developed. Surface morphology of the films are investigated by AFM and SEM methods. The ab initio quantum - chemical calculations of the geometry, total energy and excited energy states of complex fullerene- cadmium telluride supramolecules are performed. Photoluminescence spectra of composite thin films based on C60 before and after X-ray irradiation were measured. The intensity of additional peaks is defined as the charge composition due to the type of substrate. These results are interpreted as an appearance of the dipole-allowed transitions in the fullerene excited singlet states spectrum cause of an interference with cadmium telluride. X-ray irradiated films were investigated, and additional peaks in photoluminescence spectra were detected. These peaks appear as a result of molecular complexes formation from C60CdTe mixture and dimerization of the films. Density functional B3LYP quantum-chemical calculations for C60CdTe, molecular complexes, (C60)2 and C120O dimers were performed to elucidate some experimental results.

Measurement of repulsive forces between charged phospholipid bilayers

Electrostatic forces between membranes containing charged lipids were assumed to play an important role in influencing interactions between membranes long before quantitative measurements of such forces were available. ~ur measurements were designed to measure electrostatic forces between layers of lecithin charged with lipi~s carrying ionizable head groups. These experiments have shown that the interactions between charged lipid bila.yere are dominated by electrostatic forces only at separations greater than 30 A. At smaller separations the repulsion between charged bilayers is dominated by strong hydration forces. The net repulsive force between egg lecithin bilayers containing various amounts of cherged lipids (phosphatidylglycerol (PG) 5,10 ano 50 mole%, phosphatidyli. nosi tol (PI) 10 mole% and sodium oleate (Na-Ol) 3,5 and 10 mole%, where mole% gives the ratio of the number of moles' of .charged lipid to the total number of moles of all lipids present in the sample) was stuoied with the help ('If the osmotic streas technique described by LeNeveu et aI, (1977). Also, the forces between pure PG were j_nvestigated in the same manner. The results have been plotted showing variation of force as a function of bilay- _ er separation dw• All curVes 90 obtained called force curves, were found to be similar in sha.pe, showing two distinct regions, one when dw<.30 A is a region cf very rapid iiivariation of force with separation ( it is the region dominated by hydre,tion force) and second when dw> 40 A is a region of very slow variation of force with separB.tion ( it is the region dominated by the electrostatic force). Between these two regions there exists a transition area in which, in most systems studied, a phase separation of lipids into fractions containing different amounts of charged groups, was observed. A qualitative analysis showed that our results were v/ell described by the simple electrostatic double -le.yer theory. For quantitative agreement between measured and calculated force curves however, the charge density for the calculations had to be taken as half of that given by the number density of charged lipids present in the lecithin bilayers. It is not clear at the moment what causes such low apparent degree of ionization among the charged head groups, and further study is needed in this area.

Molecular dynamics calculation of mean square displacement in alkali metals and rare gas solids and comparison with lattice dynamics

Molec ul ar dynamics calculations of the mean sq ua re displacement have been carried out for the alkali metals Na, K and Cs and for an fcc nearest neighbour Lennard-Jones model applicable to rare gas solids. The computations for the alkalis were done for several temperatures for temperature vol ume a swell as for the the ze r 0 pressure ze ro zero pressure volume corresponding to each temperature. In the fcc case, results were obtained for a wide range of both the temperature and density. Lattice dynamics calculations of the harmonic and the lowe s t order anharmonic (cubic and quartic) contributions to the mean square displacement were performed for the same potential models as in the molecular dynamics calculations. The Brillouin zone sums arising in the harmonic and the quartic terms were computed for very large numbers of points in q-space, and were extrapolated to obtain results ful converged with respect to the number of points in the Brillouin zone.An excellent agreement between the lattice dynamics results was observed molecular dynamics and in the case of all the alkali metals, e~ept for the zero pressure case of CSt where the difference is about 15 % near the melting temperature. It was concluded that for the alkalis, the lowest order perturbation theory works well even at temperat ures close to the melting temperat ure. For the fcc nearest neighbour model it was found that the number of particles (256) used for the molecular dynamics calculations, produces a result which is somewhere between 10 and 20 % smaller than the value converged with respect to the number of particles. However, the general temperature dependence of the mean square displacement is the same in molecular dynamics and lattice dynamics for all temperatures at the highest densities examined, while at higher volumes and high temperatures the results diverge. This indicates the importance of the higher order (eg. ~* ) perturbation theory contributions in these cases.

Magnetic-field dependence of hole levels in self-assembled InGaAs quantum dots

Recent magnetotransport experiments of holes in InGaAs quantum dots [D. Reuter, P. Kailuweit, A. D. Wieck, U. Zeitler, O. Wibbelhoff, C. Meier, A. Lorke, and J. C. Maan, Phys. Rev. Lett. 94, 026808 (2005)] are interpreted by employing a multiband k¿p Hamiltonian, which considers the interaction between heavy hole and light hole subbands explicitly. No need of invoking an incomplete energy shell filling is required within this model. The crucial role we ascribe to the heavy hole-light hole interaction is further supported by one-band local-spin-density functional calculations, which show that Coulomb interactions do not induce any incomplete hole shell filling and therefore cannot account for the experimental magnetic field dispersion.

Importance of the embedding environment on the strain within small rings in siliceous materials

The effect of the local environment on the energetic strain within small (SiO)N rings (with N=2,3) in silica materials is investigated via periodic model systems employing density functional calculations. Through comparison of the energies of various nonterminated systems containing small rings in strained and relatively unstrained environments, with alpha quartz, we demonstrate how small ring strain is affected by the nature of the embedding environment. We compare our findings with numerous previously reported calculations, often predicting significantly different small-ring strain energies, leading to a critical assessment of methods of calculating accurate localized ring energies. The results have relevance for estimates of the strain-induced response (e.g., chemical, photo, and radio) of small silica rings, and the propensity for them to form in bulk glasses, thin films, and nanoclusters.

Theoretical study of reactivity and dynamics of hybride-bridged diruthenium complexes and silylium

Esta tesis presenta un estudio computacional de los sistemas con hidruros puente. En la primera parte se estudia la química de complejos de dirutenio con cuatro hidruros puente. Esto incluye las siguientes reacciones: el intercambio del hidruro con hidrógeno molecular; la activación del enlace C-H del etileno para formar el complejo de bis(vinilo)-etileno; el acoplamiento C-C entre el etileno coordinado y dos ligandos vinilo para producir el complejo rutenaciclopentadieno. Al final de esta parte, se discuten a detalle los mecanismos de estas reacciones. Además, se demostró la importancia de la flexibilidad de los ligandos hidruro y la cooperación entre los dos centros metálicos. En la segunda parte, se estudió el comportamiento fluxional de dos complejos μ-silileno y de un catión sililio. Con esto, se estableció la ruta más favorable en donde se realiza el intercambio de los ligandos hidruro y de los grupos metilo en los complejos μ-silileno. Finalmente, se encontró que hay dos posibles rutas relativas al cambio en la posición del puente Si-H-Si en el cation sililio poliagóstico, asociadas con la rotación interna de los grupos sililo.

Substituent effects in single helix-double helix interconversions in copper(I) complexes

The copper(I) complex of L, the 1:2 condensate of benzil dihydrazone and 2-formylpyridine, exists as single, helical [CuL](+) and double helical [Cu2L2](2+) in dichloromethane solution but crystallizes only as the double helicate [Cu2L2](ClO4)(2). In contrast, earlier [New J Chem, 27 (2003) 193] it has been found that with L', the 1:2 condensate of benzil dihydrazone and 2-acetylpyridine, only the single helical monomeric species [CuL'](+) is isolable as solid. This contrasting behaviour of the copper(I) complexes of L and L' are scrutinised here by density functional calculations.

Design and theoretical study of nickel catalysts for syndiotactic polyolefins

A nickel catalyst was modeled with ligand L-2, [ NH = CH-CH = CH-O](-), which should have potential use as a syndiotactic polyolefin catalyst, and the reaction mechanism was studied by theoretical calculations using the density functional method at the B3LYP/ LANL2MB level. The mechanism involves the formation of the intermediate [(NiLMe)-Me-2](+), in which the metal occuples a T-shaped geometry. - This intermediate has two possible structures with the methyl group trans either to the oxygen or to the nitrogen atom of L-2. The results show that both structures can lead to the desired product via similar reaction paths, A and B. Thus, the polymerization could be considered as taking place either with the alkyl group occupying the position trans to the Ni-O or trans to the Ni-N bond in the catalyst. The polymerization process thus favors the catalysis of syndiotactic polyolefins. The syndiotactic synthesis effects could also be enhanced by variations in the ligand substituents. From energy considerations, we can conclude that it is more favorable for the methyl group to occupy the trans-O position to form a complex than to occupy the trans-N position. From bond length considerations, it is also more favoured for ethene to occupy the trans-O position than to occupy the trans-N position.

Metallo-organic domino reactions: C-H versus C-C bond breaking

HL and MeL are prepared by condensing benzil dihydrazone with 2-formylpyridine and 2-acetylpyridine, respectively, in 1:2 molar proportions. While in a reaction with [Ru-(C6H6)Cl-2](2), HL yields the cation [Ru(C6H6){5,6-diphenyl-3-(pyridin-2-yl)- 1,2,4-triazine}Cl](+), MeL gives the cation [Ru(C6H6)(MeL)Cl](+). Both the cations are isolated as their hexafluorophosphate salts and characterised by X-ray crystallography. In the case of HL, double domino electrocyclic/elimination reactions are found to occur. The electrocyclic reaction occurs in a C=N-N=C-C=N fragment of HL and the elimination reaction involves breaking of a C-H bond of HL. Density functional calculations on model complexes indicate that the identified electrocyclic reaction is thermochemically as well as kinetically feasible for both HL and MeL in the gas phase. For a double domino reaction, similar to that operative in HL, to occur for MeL, breaking of a C-C bond would be required in the elimination step. Our model calculations show the energy barrier for this elimination step to be much higher (329.1 kJ mol(-1)) for MeL than that for HL (96.3 kJ mol(-1)). Thus, the domino reaction takes place for HL and not for MeL. This accounts for the observed stability of [Ru(C6H6)-(MeL)Cl](+) under the reaction conditions employed.