Predicting and characterizing new molecular species at a high-level theoretical approach: HSeCl, HClSe, and SeCl

New molecular species HSeCl, HClSe, and SeCl were investigated at a high level of theory, CCSD(T), with a series of correlation consistent basis sets with extrapolation to the CBS limit. Account has been taken for valence-only and core-valence correlation effects, and of anharmonic effects on the vibrational frequencies. HSeCl is 43.25 kcal mol (1) more stable than HClSe. A barrier (Delta G(#)) of 47.20 kcal mol (1) separates these species. Internuclear distances are generally overestimated by 0.008 angstrom in the valence-only correlation calculations. Inclusion of anharmonicity leads to much improved vibrational frequencies. For SeCl, we estimate Delta H(f) (0 K) = 23.96 and Delta H(f) (298.15 K) = 24.64 kcal mol (1); for HSeCl, we had 4.20 and 4.97 kcal mol (1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

Exploring new molecular species on the (1)[H, Se, F] singlet potential energy surface: Energetics, structures, IR spectra, and heats of formation

Structural, energetic, and vibrational properties of new molecular species, HSeF and HFSe, the associated transition state, and dissociation fragments are investigated using a state-of-the-art theoretical approach, CCSD(T)/CBS. HSeF is a normal covalently bonded molecule 38.98 kcal mol (1) more stable than the complex HF-Se, which shows an unusual structure with a central fluorine atom and a bond angle of 101.8 degrees.A barrier (Delta G(#)) of 49.01 kcal mol (1) separates the two species. Vibrational frequencies are also quite distinct. Heats of formation are evaluated for the diatomic fragments and HSeF. Final Delta(f)H values depend on the experimental accuracy of those of Se(g) and H(2)Se. (c) 2009 Elsevier B.V. All rights reserved.

Predicting new molecular species of potential interest to atmospheric chemistry: The isomers HSI and HIS

Structural, vibrational, and energetic properties of new molecular species, HSI and HIS are investigated for the first time using a state-of-the-art theoretical approach. These molecules can be easily differentiated by their geometric parameters and vibrational spectra. HSI is much more stable, and a direct unimolecular isomerization is very unlikely. Kinetics estimates predict that only at low temperatures there is a possibility of isolating HIS. For HS-I, we estimate a bond dissociation energy of 46.25 kcal/mol, and a heat of formation at 298.15 K of 12.84 kcal/mol. For the H(2)S + I(2) -> HSI + HI reaction enthalpy, we found 8.40 kcal/ mol. (C) 2008 Elsevier B.V. All rights reserved.

Classic and molecular cytogenetic analyses reveal chromosomal gains and losses correlated with survival in head and neck cancer patients

Purpose: Genetic biomarkers of head and neck tumors could be useful for distinguishing among patients with similar clinical and histopathologic characteristics but having differential probabilities of survival. The purpose of this study was to investigate chromosomal alterations in head and neck carcinomas and to correlate the results with clinical and epidentiologic variables.Experimental Design: Cytogenetic analysis of short-term cultures from 64 primary untreated head and neck squamous cell carcinomas was used to determine the overall pattern of chromosome aberrations. A representative subset of tumors was analyzed in detail by spectral karyotyping and/or confirmatory fluorescence in situ hybridization analysis.Results: Recurrent losses of chromosomes Y (26 cases) and 19 (14 cases), and gains of chromosomes 22 (23 cases), 8 and 20 (11 cases each) were observed. The most frequent structural aberration was del(22)(q13.1) followed by rearrangements involving 6q and 12p. The presence of specific cytogenetic aberrations was found to correlate significantly with an unfavorable outcome. There was a significant association between survival and gains in chromosomes 10 (P = 0.008) and 20 (P = 0.002) and losses of chromosomes 15 (P = 0.005) and 22 (P = 0.021). Univariate analysis indicated that acquisition of monosomy 17 was a significant (P = 0.0012) factor for patients with a previous family history of cancer.Conclusions: the significant associations found in this study emphasize that alterations of distinct regions of the genome may be genetic biomarkers for a poor prognosis. Losses of chromosomes 17 and 22 can be associated with a family history of cancer.

Crystal and molecular structures of 4-substituted 3,4-dihydropyrimidin-2(1H)-ones studied by X-ray and AM1 and B3LYP calculations

(1) C13H13N3O5, Mr = 291.26, P (1) over bar, a = 7.4629(9), b = 7.9203(9), c = 12.126(2) angstrom, alpha = 86.804(5), beta = 78.471(7), gamma = 69.401(8)degrees, V = 657.3(2)angstrom(3), Z = 2, R-1 = 0.0454; (2) C11H12N2O4, Mr=236.23, Pbca, a=7.2713(9), b=14.234(1), c=20.848(3)angstrom, V= 2157.8(4) angstrom(3), Z=8, R-1=0.0504; (3) C13H13N2O3Cl, Mr = 280.70, P2/n, a = 17.344(2), b = 9.237(1), c = 18.398(2) angstrom; beta = 92.61(2)degrees, V = 2944.4(6) angstrom(3), Z = 8, R-1 = 0.0714. The conformational features of three 4-substituted-3-4-dihydropyrimidin-2(1H)-ones were investigated by computational and single crystal X-ray crystallographic studies. The geometries were optimized using semiempirical (AM1) and first principle calculations (B3LYP/6-31G**) methods, the rotational barriers for important functional groups were studied. In all structures the pyrimidinone rings are in a more or less distorted boat conformation. The phenyl and the furane rings are almost perpendicular to the best least-squares plane through the dihydropyrimidinone ring.

Structure characterization of molecular complexes for non-linear optical materials. II. 1 : 1 complexes of 4-methyl-morpholine-N-oxide (1) and 3-picoline-N-oxide (2) with 2,4,6-trinitrophenol, studied by X-ray, AM1 and DFT calculations

(1) C6H2N3O7- center dot C5H12NO2+, Mr = 346.26, P2(1)/c, a = 7.2356(6), b = 10.5765(9), c = 19.593(2) angstrom, 3 beta=95.101(6)degrees, V = 1493.5(2) angstrom(3), Z = 4, R-1 = 0.0414; (2) C6H2N3O7- center dot C6H8NO+, Mr = 38.24, P2(1)/n, a = 7.8713(5), b = 6.1979(7), c = 28.697(3) angstrom, beta = 90.028(7)degrees, V = 1400.0(2) angstrom(3), Z = 4, R-1 = 0.0416. The packing units in both compounds consist of hydrogen bonded cation-anion pairs. The (hyper)polarizabilities have been calculated for the crystallographic and optimized molecules, by AM1 and at the DFT/B3LYP(6-31G**) level.

Fully relativistic calculations on the potential energy surfaces of the lowest 23 states of molecular chlorine

The electronic structure and spectroscopic properties (R(e), omega(e), omega(e)x(e), beta(e), and T(e)) of the ground state and the 22 lowest excited states of chlorine molecule were studied within a four-component relativistic framework using the MOLFDIR program package. The potential energy curves of all possible 23 covalent states were calculated using relativistic complete open shell configuration interaction approach. In addition, four component multireference configuration interaction with single and double excitation calculations were performed in order to infer the effects due to dynamical correlation in vertical excitations. The calculated properties are in good agreement with the available experimental data.

Ab Initio correlated all electron Dirac-Fock calculations for eka-francium fluoride (E119F)

Resultados obtidos por cálculos Dirac-Fock correlacionados de 4 componentes para o fluoreto do elemento E119 (Eka-Frâncio) com base estável e precisa, livre de prolapso variacional, são reportados neste trabalho. No nível CCSD(T), a distância de equilíbrio R_e, frequência harmônica ω_e e energia de dissociação D_e são 2,432 Å, 354,97 cm_-1 e 116,92 kcal mol_-1, respectivamente. Também são reportados base livre de prolapso variacional de 4 componentes para o elemento 119, uma curva analítica de energia potencial precisa e o espectro vibracional a partir dos dados obtidos no nível CCSD(T). Nossos resultados sugerem que a molécula E119F deva ser menos iônica que seus fluoretos alcalinos homólogos mais leves, em contraste com o senso químico comum baseado nas propriedades periódicas - era de se esperar nesta molécula a ligação química mais iônica possível. Também encontramos que a correção do tipo modelo de carga para negligenciar as integrais do tipo SS resulta em erros insignificantes e acelera os cálculos cerca de 3 vezes no nível CCSD(T) e cerca de 4 vezes no nível DFT/B3LYP.

Theoretical study of molecular and electronic structures of 5(1) knot systems: two-layered ONIOM calculations

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

Theoretical study of correlated systems using hybrid functionals

Diese Arbeit unterstreicht das Potential von Hybridfunktionalen (B3LYP) für die Untersuchung einer großen Bandbreite von Systemen. Durch die Einbeziehung der exakten Hartree-Fock Austauschenergie kann B3LYP für molekulare und kristalline Systeme eingesetzt werden. Zum Beispiel können stark korrelierte Systeme mit B3LYP erfolgreich erforscht werden. Die elektronische Struktur von PAHs wurde mit B3LYP Hybriddichtefunktionalen untersucht. Mit der ∆SCF-Methode wurden Elektronenbindungsenergien bestimmt, welche die mit UPS gewonnenen experimentellen Resultate bestätigen und ergänzen. Symmetrieeigenschaften der molekularen Orbitale wurden analysiert, um eine Zuordnung und Einschätzung der zugehörigen Signalstärke zu ermöglichen. Während σ-artige Orbitale nur schwer durch UPS-Messungen an dünnen Filmen detektiert werden können, bieten Rechnungen eine detaillierte Einsicht in die verborgenen Teile der Spektren.rnWeiterhin wurden π−π-Komplexe untersucht, welche von verschiedenen Donor- und Akzeptor-Molekülen gebildet werden. Die Moleküle basieren auf polyzyklischen, aromatischen Kohlenwasserstoffen. Für Ladungstransferkomplexe finden DFT Rechnungen ein Minimum in der Oberfläche der potentiellen Energie. Diese attraktive Wechselwirkung wird durch Coulombanziehung verursacht. Allerdings ist die Coulombanziehung nicht die stärkste Wechselwirkung in Ladungstransferkomplexen. Die Einbeziehung von van der Waals-Korrekturen verbessert den intermolekularen Abstand und die Bindungsenergie.rnEine Verkleinerung der intermolekularen Abstände führt zu einer großen Verschiebung der HOMO- und LUMO-Energie.rnAus der Klasse der kristallinen korrelierten Systeme wurden Rb4O6 und FeSe untersucht. Im Falle von Rb4O6 führen Ladungsordnung und Korrelationen zu einem isolierenden Grundzustand. Das hypothetische druckabhängige Phasendiagramm wurde untersucht. Eine Erhöhung des Drucks führt zu einer vergrößerten Bandlücke. Bei etwa 75 GPa wird die Bandbreite W größer als der Bandabstand U und das System nimmt einen homogen gemischt valenten Zustand mit teilweise besetzten π−π-Orbitalen an. Für Drücke ab 160 GPa wird W sehr viel größer als U und das System wird metallisch.rnIm Fall von FeSe finden wir eine korrelierte und isolierende Phase bei hohen Drücken, während das System bei niedrigen Drücken supraleitendes Verhalten zeigt. Die Berechnungen der Elektronenstruktur mit dem Hybridfunktional B3LYP führt zum korrekten halbleitenden Grundzustand in der NiAs- und MnP-Struktur von FeSe. Die Rolle der Korrelationen, der Stöchiometrie und der Nähe zum Magnetismus wird besprochen. Im Speziellen wird gezeigt, dass die Phase mit NiAs-Struktur starke lokale Korrelationen aufweist, was zu einem halbleitenden Zustand in einem weiten Druckbereich führt.

Excitonic splittings in molecular dimers: why static ab initio calculations cannot match them

After decades of research on molecular excitons, only few molecular dimers are available on which exciton and vibronic coupling theories can be rigorously tested. In centrosymmetric H-bonded dimers consisting of identical (hetero)aromatic chromophores, the monomer electronic transition dipole moment vectors subtract or add, yielding S0 → S1 and S0 → S2 transitions that are symmetry-forbidden or -allowed, respectively. Symmetry breaking by 12C/13C or H/D isotopic substitution renders the forbidden transition weakly allowed. The excitonic coupling (Davydov splitting) can then be measured between the S0 → S1 and S0 → S2 vibrationless bands. We discuss the mass-specific excitonic spectra of five H-bonded dimers that are supersonically cooled to a few K and investigated using two-color resonant two-photon ionization spectroscopy. The excitonic splittings Δcalc predicted by ab initio methods are 5–25 times larger than the experimental excitonic splittings Δexp. The purely electronic ab initio splittings need to be reduced (“quenched”), reflecting the coupling of the electronic transition to the optically active vibrations of the monomers. The so-called quenching factors Γ < 1 can be determined from experiment (Γexp) and/or calculation (Γcalc). The vibronically quenched splittings Γ·Δcalc are found to nicely reproduce the experimental exciton splittings.

MOLECULAR AND GENOMIC BASED INSIGHTS INTO THE EVOLUTION OF ENTEROCOCCUS FAECIUM FROM COMMENSAL TO HOSPITAL-ADAPTED PATHOGEN

The basis for the recent transition of Enterococcus faecium from a primarily commensal organism to one of the leading causes of hospital-acquired infections in the United States is not yet understood. To address this, the first part of my project assessed isolates from early outbreaks in the USA and South America using sequence analysis, colony hybridizations, and minimal inhibitory concentrations (MICs) which showed clinical isolates possess virulence and antibiotic resistance determinants that are less abundant or lacking in community isolates. I also revealed that the level of ampicillin resistance increased over time in clinical strains. By sequencing the pbp5 gene, I demonstrated an ~5% difference in the pbp5 gene between strains with MICs <4ug/ml and those with MICs >4µg/ml, but no specific sequence changes correlated with increases in MICs within the latter group. A 3-10% nucleotide difference was also seen in three other genes analyzed, which suggested the existence of two distinct subpopulations of E. faecium. This led to the second part of my project analyzing concatenated core gene sequences, SNPs, the 16S rRNA, and phylogenetics of 21 E. faecium genomes confirming two distinct clades; a community-associated (CA) clade and hospital-associated (HA) clade. Molecular clock calculations indicate that these two clades likely diverged ~ 300,000 to > 1 million years ago, long before the modern antibiotic era. Genomic analysis also showed that, in addition to core genomic differences, HA E. faecium harbor specific accessory genetic elements that may confer selection advantages over CA E. faecium. The third part of my project discovered 6 E. faecium genes with the newly identified “WxL” domain. My analyses, using RT-PCR, western blots, patient sera, whole-cell ELISA, and immunogold electron microscopy, indicated that E. faecium WxL genes exist in operons, encode bacterial cell surface localized proteins, that WxL proteins are antigenic in humans, and are more exposed on the surface of clinical isolates versus community isolates (even though they are ubiquitous in both clades). ELISAs and BIAcore analyses also showed that proteins encoded by these operons bind several different host extracellular matrix proteins, as well as to each other, suggesting a novel cell-surface complex. In summary, my studies provide new insights into the evolution of E. faecium by showing that there are two distantly related clades; one being more successful in the hospital setting. My studies also identified operons encoding WxL proteins whose characteristics could also contribute to colonization and virulence within this species.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Molecular mechanics calculations of the riboacetal internucleotide linkage in double and triple helices.

Structures of Watson-Crick base paired 15-nucleobase oligomer strands in A-type or B-type conformation in which one strand [a strand of alternating nucleotide and riboacetal thymidine nucleoside (RT) units, RP] is DNA and the other is composed of alternating nucleotides and riboacetal nucleosides have been studied by molecular mechanics. Analogously, oligomer strands of RNA in place of DNA have been modeled. The calculations indicate that the RP strand is more stable when complexed in an A-type duplex relative to a B-type form and that this conformational preference is presumably due to the more uniform nature of the former. Nearly planar ribose rings were more commonly observed in the minimized structures of the B-type DNA.RP duplexes as compared with A-type duplexes, despite the fact that planar ribofuranose rings are known to be energetically unfavorable in oligonucleotides. Computed relative stabilities of all duplexes containing the RP strand suggest that such heteroduplexes are less stable than the corresponding double-stranded DNA and double-stranded RNA species. These findings are in agreement with experimental results which show, when equivalent sequences were compared, that a DNA.RNA control forms a more stable duplex than RP hound to a complementary single-stranded RNA strand. In contrast, molecular mechanics studies of complementary triple-helical (DNA)2.RP, (DNA)2.DNA, and (DNA)2.RNA structures indicate that the binding of RP as a Hoogsteen strand stabilizes the underlying duplex to a greater extent compared with native oligonucleotides. These calculations suggest that puckering of the ribose ring in the riboacetal linkage leads to a more favorable interaction with a complementary nucleic acid target than the proposed planar geometry and that this puckering may account for the enhanced binding of RP to a double-stranded target.