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.

Investigation of the Nucleation and Growth of Colloidal Indium Phosphide: from Molecular Precursors to Semiconductor Nanocrystals through In37P20(O2CR)51 as a Magic-Size Intermediate

Thesis (Ph.D.)--University of Washington, 2016-06

Approximate ab initio calculations and the method of molecular fragments

A two stage approach to performing ab initio calculations on medium and large sized molecules is described. The first step is to perform SCF calculations on small molecules or molecular fragments using the OPIT Program. This employs a small basis set of spherical and p-type Gaussian functions. The Gaussian functions can be identified very closely with atomic cores, bond pairs, lone pairs, etc. The position and exponent of any of the Gaussian functions can be varied by OPIT to produce a small but fully optimised basis set. The second stage is the molecular fragments method. As an example of this, Gaussian exponents and distances are taken from an OPIT calculation on ethylene and used unchanged in a single SCF calculation on benzene. Approximate ab initio calculations of this type give much useful information and are often preferable to semi-empirical approaches, since the nature of the approximations involved is much better defined.

Computing the structural and vibrational properties of polymorphic organic molecular crystals through van der Waals corrected density functional theory and the electronic properties of organic thin films through microelectrostatic calculations.

The present Thesis reports on the various research projects to which I have contributed during my PhD period, working with several research groups, and whose results have been communicated in a number of scientific publications. The main focus of my research activity was to learn, test, exploit and extend the recently developed vdW-DFT (van der Waals corrected Density Functional Theory) methods for computing the structural, vibrational and electronic properties of ordered molecular crystals from first principles. A secondary, and more recent, research activity has been the analysis with microelectrostatic methods of Molecular Dynamics (MD) simulations of disordered molecular systems. While only very unreliable methods based on empirical models were practically usable until a few years ago, accurate calculations of the crystal energy are now possible, thanks to very fast modern computers and to the excellent performance of the best vdW-DFT methods. Accurate energies are particularly important for describing organic molecular solids, since they often exhibit several alternative crystal structures (polymorphs), with very different packing arrangements but very small energy differences. Standard DFT methods do not describe the long-range electron correlations which give rise to the vdW interactions. Although weak, these interactions are extremely sensitive to the packing arrangement, and neglecting them used to be a problem. The calculations of reliable crystal structures and vibrational frequencies has been made possible only recently, thanks to development of some good representations of the vdW contribution to the energy (known as “vdW corrections”).

Core-valence Correlation Effects On Ir Calculations: The Bf3 And Bcl3 Cases.

The first theoretical results of core-valence correlation effects are presented for the infrared wavenumbers and intensities of the BF3 and BCl3 molecules, using (double- and triple-zeta) Dunning core-valence basis sets at the CCSD(T) level. The results are compared with those calculated in the frozen core approximation with standard Dunning basis sets at the same correlation level and with the experimental values. The general conclusion is that the effect of core-valence correlation is, for infrared wavenumbers and intensities, smaller than the effect of adding augmented diffuse functions to the basis set, e.g., cc-pVTZ to aug-cc-pVTZ. Moreover, the trends observed in the data are mainly related to the augmented functions rather than the core-valence functions added to the basis set. The results obtained here confirm previous studies pointing out the large descrepancy between the theoretical and experimental intensities of the stretching mode for BCl3.

Low Molecular Weight Protein Tyrosine Phosphatase (lmwptp) Upregulation Mediates Malignant Potential In Colorectal Cancer.

Phosphatases have long been regarded as tumor suppressors, however there is emerging evidence for a tumor initiating role for some phosphatases in several forms of cancer. Low Molecular Weight Protein Tyrosine Phosphatase (LMWPTP; acid phosphatase 1 [ACP1]) is an 18 kDa enzyme that influences the phosphorylation of signaling pathway mediators involved in cancer and is thus postulated to be a tumor-promoting enzyme, but neither unequivocal clinical evidence nor convincing mechanistic actions for a role of LMWPTP have been identified. In the present study, we show that LMWPTP expression is not only significantly increased in colorectal cancer (CRC), but also follows a step-wise increase in different levels of dysplasia. Chemical inhibition of LMWPTP significantly reduces CRC growth. Furthermore, downregulation of LMWPTP in CRC leads to a reduced migration ability in both 2D- and 3D-migration assays, and sensitizes tumor cells to the chemotherapeutic agent 5-FU. In conclusion, this study shows that LMWPTP is not only overexpressed in colorectal cancer, but it is correlated with the malignant potential of this cancer, suggesting that this phosphatase may act as a predictive biomaker of CRC stage and represents a rational novel target in the treatment of this disease.

Long Terminal Repeat Sequence Analysis of HTLV-2 Molecular Variants Identified in Southern Brazil

In Brazil, human T-lymphotropic virus type 2 (HTLV-2) is endemic in Amerindians and epidemic in intravenous drug users (IDUs). The long terminal repeat (LTR) is the most divergent genomic region of HTLV-2, therefore useful to characterize subtypes. Nucleotide sequence and restriction fragment length polymorphism (RFLP) analysis of LTR genomic segments of fourteen HTLV-2 strains isolated from HIV-infected patients of Londrina, Southern Brazil, were carried out. Molecular analysis disclosed that all HTLV-2 strains belonged to 2a subtype, and RFLP detected the presence of the a4, a5, and a6 subgroups according to Switzer's nomenclature. RFLP correlated with nucleotide sequence, and phylogenetic analysis clustered HTLV-2 sequences of IDUs into subgroups a5 and a6. HTLV-2 sequences from individuals of sexual risk factor clustered into the a4 subgroup. These results extend the knowledge of the genetic diversity of HTLV-2 circulating in Brazil and provide insights into HTLV-2 transmission and virus movement in this geographic area.

Study of the antimicrobial peptide indolicidin and a mutant in micelle medium by molecular dynamics simulation

The antimicrobial peptide indolicidin (IND) and the mutant CP10A in hydrated micelles were studied using molecular dynamics simulations in order to observe whether the molecular dynamics and experimental data could be sufficiently correlated and a detailed description of the interaction of the antimicrobial peptides with a model of the membrane provided by a hydrated micelle system could be obtained. In agreement with the experiments, the simulations showed that the peptides are located near the surface of the micelles. Peptide insertions agree with available experimental data, showing deeper insertion of the mutant compared with the peptide IND. Major insertion into the hydrophobic core of the micelle by all tryptophan and mutated residues of CP10A in relation to IND was observed. The charged residues of the terminus regions of both peptides present similar behavior, indicating that the major differences in the interactions with the micelles of the peptides IND and CP10A occur in the case of the hydrophobic residues.

249 TP53 mutation has high prevalence and is correlated with larger and poorly differentiated HCC in Brazilian patients

Background: Ser-249 TP53 mutation (249(Ser)) is a molecular evidence for aflatoxin-related carcinogenesis in Hepatocellular Carcinoma (HCC) and it is frequent in some African and Asian regions, but it is unusual in Western countries. HBV has been claimed to add a synergic effect on genesis of this particular mutation with aflatoxin. The aim of this study was to investigate the frequency of 249Ser mutation in HCC from patients in Brazil. Methods: We studied 74 HCC formalin fixed paraffin blocks samples of patients whom underwent surgical resection in Brazil. 249Ser mutation was analyzed by RFLP and DNA sequencing. HBV DNA presence was determined by Real-Time PCR. Results: 249Ser mutation was found in 21/74 (28%) samples while HBV DNA was detected in 13/74 (16%). 249Ser mutation was detected in 21/74 samples by RFLP assay, of which 14 were confirmed by 249Ser mutant-specific PCR, and 12 by nucleic acid sequencing. All HCC cases with p53-249ser mutation displayed also wild-type p53 sequences. Poorly differentiated HCC was more likely to have 249Ser mutation (OR = 2.415, 95% CI = 1.001 - 5.824, p = 0.05). The mean size of 249Ser HCC tumor was 9.4 cm versus 5.5 cm on wild type HCC (p = 0.012). HBV DNA detection was not related to 249Ser mutation. Conclusion: Our results indicate that 249Ser mutation is a HCC important factor of carcinogenesis in Brazil and it is associated to large and poorly differentiated tumors.

Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.