Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

On the absorbance changes in the photocycle of the photoactive yellow protein: A quantum-chemical analysis

Spectral changes in the photocycle of the photoactive yellow protein (PYP) are investigated by using ab initio multiconfigurational second-order perturbation theory at the available structures experimentally determined. Using the dark ground-state crystal structure [Genick, U. K., Soltis, S. M., Kuhn, P., Canestrelli, I. L. & Getzoff, E. D. (1998) Nature (London) 392, 206–209], the ππ* transition to the lowest excited state is related to the typical blue-light absorption observed at 446 nm. The different nature of the second excited state (nπ*) is consistent with the alternative route detected at 395-nm excitation. The results suggest the low-temperature photoproduct PYPHL as the most plausible candidate for the assignment of the cryogenically trapped early intermediate (Genick et al.). We cannot establish, however, a successful correspondence between the theoretical spectrum for the nanosecond time-resolved x-ray structure [Perman, B., Šrajer, V., Ren, Z., Teng, T., Pradervand, C., et al. (1998) Science 279, 1946–1950] and any of the spectroscopic photoproducts known up to date. It is fully confirmed that the colorless light-activated intermediate recorded by millisecond time-resolved crystallography [Genick, U. K., Borgstahl, G. E. O., Ng, K., Ren, Z., Pradervand, C., et al. (1997) Science 275, 1471–1475] is protonated, nicely matching the spectroscopic features of the photoproduct PYPM. The overall contribution demonstrates that a combined analysis of high-level theoretical results and experimental data can be of great value to perform assignments of detected intermediates in a photocycle.

An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images

We describe an approach to the high-resolution three-dimensional structural determination of macromolecules that utilizes ultrashort, intense x-ray pulses to record diffraction data in combination with direct phase retrieval by the oversampling technique. It is shown that a simulated molecular diffraction pattern at 2.5-Å resolution accumulated from multiple copies of single rubisco biomolecules, each generated by a femtosecond-level x-ray free electron laser pulse, can be successfully phased and transformed into an accurate electron density map comparable to that obtained by more conventional methods. The phase problem is solved by using an iterative algorithm with a random phase set as an initial input. The convergence speed of the algorithm is reasonably fast, typically around a few hundred iterations. This approach and phasing method do not require any ab initio information about the molecule, do not require an extended ordered lattice array, and can tolerate high noise and some missing intensity data at the center of the diffraction pattern. With the prospects of the x-ray free electron lasers, this approach could provide a major new opportunity for the high-resolution three-dimensional structure determination of single biomolecules.

Coupling the folding of homologous proteins

The empirical observation that homologous proteins fold to similar structures is used to enhance the capabilities of an ab initio algorithm to predict protein conformations. A penalty function that forces homologous proteins to look alike is added to the potential and is employed in the coupled energy optimization of several homologous proteins. Significant improvement in the quality of the computed structures (as compared with the computational folding of a single protein) is demonstrated and discussed.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Proactive and reactive inhibition during overt and covert actions. An electrical neuroimaging study.

Response inhibition is the ability to suppress inadequate but automatically activated, prepotent or ongoing response tendencies. In the framework of motor inhibition, two distinct operating strategies have been described: “proactive” and “reactive” control modes. In the proactive modality, inhibition is recruited in advance by predictive signals, and actively maintained before its enactment. Conversely, in the reactive control mode, inhibition is phasically enacted after the detection of the inhibitory signal. To date, ample evidence points to a core cerebral network for reactive inhibition comprising the right inferior frontal gyrus (rIFG), the presupplementary motor area (pre-SMA) and the basal ganglia (BG). Moreover, fMRI studies showed that cerebral activations during proactive and reactive inhibition largely overlap. These findings suggest that at least part of the neural network for reactive inhibition is recruited in advance, priming cortical regions in preparation for the upcoming inhibition. So far, proactive and reactive inhibitory mechanisms have been investigated during tasks in which the requested response to be stopped or withheld was an “overt” action execution (AE) (i.e., a movement effectively performed). Nevertheless, inhibitory mechanisms are also relevant for motor control during “covert actions” (i.e., potential motor acts not overtly performed), such as motor imagery (MI). MI is the conscious, voluntary mental rehearsal of action representations without any overt movement. Previous studies revealed a substantial overlap of activated motor-related brain networks in premotor, parietal and subcortical regions during overtly executed and imagined movements. Notwithstanding this evidence for a shared set of cerebral regions involved in encoding actions, whether or not those actions are effectively executed, the neural bases of motor inhibition during MI, preventing covert action from being overtly performed, in spite of the activation of the motor system, remain to be fully clarified. Taking into account this background, we performed a high density EEG study evaluating cerebral mechanisms and their related sources elicited during two types of cued Go/NoGo task, requiring the execution or withholding of an overt (Go) or a covert (MI) action, respectively. The EEG analyses were performed in two steps, with different aims: 1) Analysis of the “response phase” of the cued overt and covert Go/NoGo tasks, for the evaluation of reactive inhibitory control of overt and covert actions. 2) Analysis of the “preparatory phase” of the cued overt and covert Go/NoGo EEG datasets, focusing on cerebral activities time-locked to the preparatory signals, for the evaluation of proactive inhibitory mechanisms and their related neural sources. For these purposes, a spatiotemporal analysis of the scalp electric fields was applied on the EEG data recorded during the overt and covert Go/NoGo tasks. The spatiotemporal approach provide an objective definition of time windows for source analysis, relying on the statistical proof that the electric fields are different and thus generated by different neural sources. The analysis of the “response phase” revealed that key nodes of the inhibitory circuit, underpinning inhibition of the overt movement during the NoGo response, were also activated during the MI enactment. In both cases, inhibition relied on the activation of pre-SMA and rIFG, but with different temporal patterns of activation in accord with the intended “covert” or “overt” modality of motor performance. During the NoGo condition, the pre-SMA and rIFG were sequentially activated, pointing to an early decisional role of pre-SMA and to a later role of rIFG in the enactment of inhibitory control of the overt action. Conversely, a concomitant activation of pre-SMA and rIFG emerged during the imagined motor response. This latter finding suggested that an inhibitory mechanism (likely underpinned by the rIFG), could be prewired into a prepared “covert modality” of motor response, as an intrinsic component of the MI enactment. This mechanism would allow the rehearsal of the imagined motor representations, without any overt movement. The analyses of the “preparatory phase”, confirmed in both overt and covert Go/NoGo tasks the priming of cerebral regions pertaining to putative inhibitory network, reactively triggered in the following response phase. Nonetheless, differences in the preparatory strategies between the two tasks emerged, depending on the intended “overt” or “covert” modality of the possible incoming motor response. During the preparation of the overt Go/NoGo task, the cue primed the possible overt response programs in motor and premotor cortex. At the same time, through preactivation of a pre-SMA-related decisional mechanism, it triggered a parallel preparation for the successful response selection and/or inhibition during the subsequent response phase. Conversely, the preparatory strategy for the covert Go/NoGo task was centred on the goal-oriented priming of an inhibitory mechanism related to the rIFG that, being tuned to the instructed covert modality of the motor performance and instantiated during the subsequent MI enactment, allowed the imagined response to remain a potential motor act. Taken together, the results of the present study demonstrate a substantial overlap of cerebral networks activated during proactive recruitment and subsequent reactive enactment of motor inhibition in both overt and covert actions. At the same time, our data show that preparatory cues predisposed ab initio a different organization of the cerebral areas (in particular of the pre-SMA and rIFG) involved with sensorimotor transformations and motor inhibitory control for executed and imagined actions. During the preparatory phases of our cued overt and covert Go/NoGo tasks, the different adopted strategies were tuned to the “how” of the motor performance, reflecting the intended overt and covert modality of the possible incoming action.

Análise Conformacional e estudo das interações eletrônicas de algumas α-fenilseleno-α-dietóxifosforilacetofenonas para-substituídas

A presente Dissertação relata a síntese e o estudo conformacional das α-fenilseleno-α-dietóxifosforilacetofenonas para-substituídas p-X-Φ-C(O)CH[SeΦ][P(O)(OEt2] (X=OMe 1, Me 2, H 3, F 4, Cl 5, Br 6 e NO2 7) através da banda de estiramento da carbonila no infravermelho, em solventes de polaridade crescente apoiado por cálculos ab initio HF/6-3IG**. A comparação entre a freqüência e a intensidade relativa dos componentes do dubleto, para os derivados 6 e 7, e do singleto para os derivados 1-5, no solvente apolar tetracloreto de carbono, e dos componentes do dubleto, nos solventes de polaridade crescente (clorofórmio, diclorometano e acetonitrila), para os derivados 1-7, com os dados do cálculo ab initio de 3 (composto de referência), indicou que ambas as conformações estáveis (g1 e g2) apresentam a ligação C-Se na geometria anti-clinal (gauche) em relação à carbonila (C=O), enquanto que a ligação C-P assume uma geometria sin-periplanar (cis) em relação à carbonila. A análise dos contatos interatômicos de átomos relevante em comparação com a soma de seus raios de van der Waals, indicou que ambas as conformações g1 e g2 são fortemente estabilizadas pelo sinergismo das interações orbitalares e eletrostáticas π*(CO) / nSe e Oδ-[CO].....Pδ+[PO]. Analogamente, as interações mais fracas Oδ-[OR]..... Cδ+[CO], 0-Hδ+[SeΦ]....Oδ-[PO] e 0-Hδ+[ΦC(O)]....Oδ-[CO] estabilizam as conformações g1 e g2, aproximadamente na mesma extensão. No entanto, somente a conformação g1 é estabilizada pela interação eletrostática (ligação de hidrogênio) Hδ+[α-CH].....Oδ-[OR], enquanto que sómente a conformação g2 é desestabilizada pelo Efeito de Campo Repulsivo entre os dipolos Cδ+=.Oδ- e Pδ+-ORδ- Assim sendo, pode-se concluir que no dubleto de VCO no IV, o componente de maior freqüência e de menor intensidade corresponde à conformação menos estável g2 (do cálculo) enquanto que o componente de menor freqüência e mais intenso corresponde à conformação mais estável g1 (do cálculo). Estes dados estão de pleno acordo com os deslocamentos de freqüência mais negativos da carbonila (ΔVCO) do confôrmero mais estável g1 em relação ao menos estável g2.

Mechanical, electrical, and magnetic properties of Ni nanocontacts

The dynamic deformation upon stretching of Ni nanowires as those formed with mechanically controllable break junctions or with a scanning tunneling microscope is studied both experimentally and theoretically. Molecular dynamics simulations of the breaking process are performed. In addition, and in order to compare with experiments, we also compute the transport properties in the last stages before failure using the first-principles implementation of Landauer's formalism included in our transport package ALACANT.

Exponential decay of spin-spin correlation between distant defect states produced by contour hydrogenation of polycyclic aromatic hydrocarbon molecules

The first few low-lying spin states of alternant polycyclic aromatic hydrocarbon (PAH) molecules of several shapes showing defect states induced by contour hydrogenation have been studied both by ab initio methods and by a precise numerical solution of Pariser-Parr-Pople (PPP) interacting model. In accordance with Lieb's theorem, the ground state shows a spin multiplicity equal to one for balanced molecules, and it gets larger values for imbalanced molecules (that is, when the number of π electrons on both subsets is not equal). Furthermore, we find a systematic decrease of the singlet-triplet splitting as a function of the distance between defects, regardless of whether the ground state is singlet or triplet. For example, a splitting smaller than 0.001 eV is obtained for a medium size C46H28 PAH molecule (di-hydrogenated [11]phenacene) showing a singlet ground state. We conclude that π electrons unbound by lattice defects tend to remain localized and unpaired even when long-range Coulomb interaction is taken into account. Therefore they show a biradical character (polyradical character for more than two defects) and should be studied as two or more local doublets. The implications for electron transport are potentially important since these unpaired electrons can trap traveling electrons or simply flip their spin at a very small energy cost.

Theoretical insights on electron donor-acceptor interactions involving carbon dioxide

Electron donor-acceptor (EDA) interactions are widely involved in chemistry and their understanding is essential to design new technological applications in a variety of fields ranging from material sciences and chemical engineering to medicine. In this work, we study EDA complexes of carbon dioxide with ketones using several ab initio and Density Functional Theory methods. Energy contributions to the interaction energy have been analyzed in detail using both variational and perturbational treatments. Dispersion energy has been shown to play a key role in explaining the high stability of a non-conventional structure, which can roughly be described by a cooperative EDA interaction.

Sanctissimi domini nostri Benedicti Papae XIV Bullarium : tomus primus [-tertius]

Contiene: T. 1: In quo continentur constitutiones, epistolae ab initio pontificatus usque ad annum MDCCXLVI -- T. 2: In quo continentur constitutiones, epistolae ab anno MDCCXLVI ad totum annum MDCCXLVIII -- T. 3: In quo continentur constitutiones, epistolae, alique editae ab exitu anni MDCCXLVIII usque ad totum pontificatus annum XII ...

Electronic structure methods applied to gas-carbon reactions

A review is given on the fundamental studies of gas-carbon reactions using electronic structure methods in the last several decades. The three types of electronic structure methods including semi-empirical, ab initio and density functional theory, methods are briefly introduced first, followed by the studies on carbon reactions with hydrogen and oxygen-containing gases (non-catalysed and catalysed). The problems yet to solve and possible promising directions are discussed. (C) 2002 Elsevier Science Ltd. All rights reserved.

Imidoylketene dimerization and rearrangement

FVT of pyrroledione 10 affords the NH-imidoylketene 11, which is characterized by its matrix isolation IR spectrum ( 2117 cm 1). On warming above 170 K, 11 dimerizes to the oxazinone 13, the X-ray crystal structure of which is reported. Imidoylketene 11 also undergoes a (reversible) 1,3-phenyl shift to afford the detectable alpha-oxoketenimine 16 (2062 cm(-1)) which at FVT temperatures above 400degreesC, isomerizes to 2-cyano-2-phenylacetophenone 18 (optimally at 700degreesC). Moreover, imidoylketene 11 can cyclize to azetinone 19, detectable at FVT temperatures up to 570degreesC, which undergoes cycloreversion to diphenylacetylene 20 and isocyanic acid (HNCO) 21. Energy profiles calculated at the B3LYP/6-31G** level for the unsubstituted imidoylketene, the diphenylimidoylketene 11 and the N-tert-butylimidoylketene are also reported.

Rearrangement of 2-quinolyl- and 1-isoquinolylcarbenes to naphthylnitrenes

2-Quinolylcarbene 23 and 1-isoquinolylcarbene 33 are generated by flash vacuum thermolysis (FVT) of the corresponding triazolo[1,5-a]quinoline and triazolo[5,1-a]isoquinoline 19 and 29, as well as 2-(5-tetrazolyl)quinoline and 1-(5-tetrazolyl)isoquinoline 20 and 30, respectively. These carbenes rearrange to 1- and 2-naphthylnitrene 21 and 31, respectively, and the nitrenes are also generated by FVT of 1- and 2-naphthyl azides 18 and 28. The products of FVT of both the nitrene and carbene precursors are the 2- and 3-cyanoindenes 26 and 27 together with the nitrene dimers, viz. azonaphthalenes 25 and 35, and the H-abstraction products, aminonaphthalenes 24 and 34. All the azide, triazole, and tetrazole precursors yield 3-cyanoindene 26 as the principal ring contraction product under conditions of low FVT temperature (340-400 degreesC) and high pressure (1 Torr N-2 as carrier gas for the purpose of collisional deactivation). This ring contraction reaction is strongly subject to chemical activation, which caused extensive isomerization of 3-cyanoindene to 2-cyanoindene under conditions of low pressure (10(-3) Torr). 2-Cyanoindene is calculated to be ca. 1.7 kcal/mol below 3-cyanoindene in energy; accordingly, high-temperature FVT of these cyanoindenes always gives mixtures of the two compounds with the 2-cyano isomer dominating. Photolysis of trizolo[1,5-a]quinoline 19 and triazolo[5,1-a]isoquinoline 29 in Ar matrixes causes partial ring opening to the corresponding 2-diazomethylquinoline 19' and 1-diazomethylisoquinoline 29'. The photolysis of the former gives rise to a small amount of the cyclic ketenimine 22, the intermediate connecting 2-quinolylcarbene and 1-naphthylnitrene.

The rearrangements of naphthyinitrenes: UV/Vis and IR spectra of azirines, cyclic ketenimines, and cyclic nitrile ylides

Ar matrix photolysis of 1- and 2-naphthyl azides 3 and 4 at 313 nm initially affords the singlet naphthyl nitrenes, (1)1 and (1)2. Relaxation to the corresponding lower energy, persistent triplet nitrenes (3)1 and (3)2 competes with cyclization to the azirines 15 and 18, which can also be formed photochemically from the triplet nitrenes. On prolonged irradiation, the azirines can be converted to the seven-membered cyclic ketenimines 10 and 13, respectively, as described earlier by Dunkin and Thomson. However, instead of the o-quinoid ketenimines 16 and 19, which are the expected primary ring-opening products of azirines 15 and 18, respectively, we observed their novel bond-shift isomers 17 and 20, which may be formally regarded as cyclic nitrile ylides. The existence of such ylidic heterocumulenes has been predicted previously, but this work provides the first experimental observation of such species. The factors which are responsible for the special stability of the ylidic species 17 and 20 are discussed.