Absorção e Emissão de radiação eletromagnética por moléculas na região do UV-Vis fornece informações sobre os estados eletrônicos excitados, sendo propriedades de grande interesse devido a sua relação com processos biológicos, bem como suas possíveis aplicações em diagnósticos e tecnologia. Essas propriedades são sensíveis ao meio em que as moléculas se encontram, tornando-se assim natural a busca por métodos teóricos que possibilitam descrever essas interações. Neste trabalho, nós usamos a metodologia Sequential-Quantum Mechanics/Molecular Mechanics para estudar o espectro de absorção e de emissão de moléculas de relevância biológica, quando em solução. Simulações clássicas Monte Carlo foram usadas pra construir uma configuração eletrostática média do líquido para posterior cálculo das propriedades, feito através dos métodos multiconfiguracionais CASSCF e CASPT2. Cuidados especiais foram tomados para incluir a polarização eletrônica que o soluto sofre devido à presença do solvente. Nossa contribuição é a adaptação do método do Gradiente de Energia Livre para permitir a obtenção de geometrias de estados excitados do soluto em solução. A técnica foi implementada em um programa e aplicada com sucesso nos sistemas estudados. As perspectivas agora se abrem para a obtenção de intersecções cônicas em meio solvente, permitindo assim o estudo de decaimentos não-radiativos em sistemas solvatados. / Absorption and emission of electromagnetic radiation by molecules in the UV-Vis region yields information about the electronic excited states, being properties of great interest due its relation with biological process, as well as its possible applications in diagnostics and technology. These properties are sensible to the environment in which the molecules are, making natural the search for theoretical methods that describe these interactions. In this work, we have used the Sequential-Quantum Mechanics/Molecular Mechanics methodology to study the spectrum of absorption and emission of molecules of biological relevance, when in solution. Classical Monte Carlo simulations were carried out to construct an average solvent electrostatic configuration and used to represent the liquid in posterior quantum mechanics calculation of the properties, performed by using the multiconfigurational methods CASSCF and CASPT2. Special cares were taken to include the electronic polarization of the solute due the solvent. Our contribution is the adaptation of the Free Energy Gradient method to allow the determination of the solute excited states geometries in solution. The method was implemented in a computer program and successfully applied in the systems studied. The perspectives are now open to the calculation of conical intersection in solvent environment, enabling the study of nonradioactive decays in solvated systems.

La bioluminescence est un processus complexe dans lequel la réaction chimique d'un substrat, catalysée par une protéine, entraîne l'apparition d'une émission de lumière dans le spectre visible. Dans le cas des lucioles, un insecte émettant dans le domaine du jaune-vert, le substrat se nomme luciférine et la protéine luciférase. Cependant, la taille et la complexité de ce système chimique limite sa compréhension, notamment celle du mécanisme réactionnel.L'apport de la chimie théorique dans ce domaine est essentiel et a prouvé son utilité dans de nombreux cas. L'utilisation de la méthode QM/MM, méthode hybride couplant la mécanique quantique et la mécanique moléculaire permet de modéliser et d'étudier ces systèmes biologiques.Cette thèse se focalise sur deux approches différentes de l'étude de la bioluminescence chez les lucioles. La première consiste à étudier l'effet de certaines modifications chimiques sur la couleur de la bioluminescence. On s'intéresse plus particulièrement à un analogue de la luciférine et à certaines luciférases issues d'autres systèmes bioluminescents. Par cette étude on cherche à rationaliser et à prédire l'effet ainsi que l’impact de ces changements sur l’émission. Le deuxième sujet explore deux étapes du mécanisme réactionnel de la bioluminescence: d'une part, la coordination du dioxygène sur un intermédiaire de la réaction, une étape encore non étudiée et d'autre part la réaction de tautomérisation dans l'état excité et au sein de la protéine entre deux formes émissives possibles de la luciférine / The bioluminescence is a complex process that involves the reaction of a substrate, catalysed by an enzyme that sheds light in the visible spectra. In fireflies, the light emitted has a yellow-green tone thanks to the interaction between the substrate luciferin and the protein luciferase. However the size and the complexity of the system prevent its comprehension especially when dealing with the reaction mechanism.The use of computational chemistry is key to understand and improve the comprehension of the bioluminescence. The hybrid QM/MM method that combines quantum mechanics with molecular mechanics is a great tool to model and study bioluminescent systems.This thesis deals with two different approaches of the bioluminescence in fireflies. The first one is related to the study of chemical modifications that tune the emission colour. We will discuss about one analogue of the luciferin and on new luciferases from others bioluminescent species. The goals of this study are to rationalise and predict both the effect and the impact of these modifications on the emission. The second subject deals with two different steps of the bioluminescent mechanism. The first one discusses the binding of the dioxygen to the bioluminescent intermediate, which was so far unstudied and the second one about the tautomerization in the excited state and in the protein of two possible emissive forms of the luciferin

L’obtention, dans des conditions douces, sélectives et de chimie durable de molécules ciblées est aujourd’hui un enjeu majeur. Les métalloenzymesartificielles représentent une voie d’investigation importante, car en jouant, par exemple, sur la seconde sphère de coordination,il est possible de modifier fortement la réactivité de ces systèmes bio-inspirés. Le développement de cette chimie suppose une connaissanceapprofondie des différentes étapes du mécanisme de la réaction envisagée. Pour cela, la chimie théorique est essentielle à la rationalisation dela réactivité chimique mais elle souffre encore de nombreuses insuffisances pour les systèmes que nous nous proposons d’étudier.Dans ce travail, nous avons choisi d’étudier la superoxyde réductase, enzyme détoxifiante du radical superoxyde. Si de nombreuses expérimentalessont disponibles détaillant certains intermédiaires, le mécanisme précis est peu documenté. Le but a été de mettre en place uneméthodologie complète allant du développement de paramètres MM spécifiques à l’étude de la réactivité par métadynamiques QM/MM.Le développement de paramètres MM pour le site actif à fer a permis son étude en dynamique MM donnant des informations sur la conformationsdu squelette peptidique ainsi que l’interaction avec les molécules de solvant. De part la nature du fer, une description QM du site actifà été nécessaire via l’utilisation de DFT hybride. Les métadynamiques QM/MM ont permis quant à elles d’explorer les chemins réactionnelset de caractériser les espèces ainsi formées et les énergies d’activations.Cette méthodologie a permis la compréhension en premier lieu de la réactivité native de la forme sauvage et elle a aussi permis d’explorer lesréactivités nouvelles des mutations de la SOR permettant ainsi de définir le rôle crucial de la seconde sphère de coordination. / Obtaining targeted molecules under gentle, selective and sustainable conditions is still a major challenge. Artificial metalloenzymes are animportant line of enquiry, because by playing, for example, with the second sphere of coordination, it is possible to strongly modify thereactivity of these bio-inspired systems. The development of this chemistry presupposes a thorough knowledge of the different stages of themechanism of the reaction under study. For this reason, theoretical chemistry is essential to rationalize chemical reactivity, but it still suffersfrom many shortcomings for the systems we propose to study.In this work, we study the superoxide reductase, a detoxifying enzyme of the superoxide radical. While many experiments are available detailingsome intermediates, the precise mechanism is not well documented. The aim was to implement a complete methodology ranging from thedevelopment of specific MM parameters to the study of reactivity by QM/MM metadynamics.The development of MM parameters for the iron active site allowed its study by MM dynamics giving informations on the conformation ofthe peptide backbone as well as on the interaction with solvent molecules. Due to the nature of the iron, a QM description of the active sitewas required using hybrid DFT. QM/MM metadynamics have allowed us to explore reaction pathways and to characterize the compoundsformed to obtain the needed activation energies. This methodology made it possible to understand the native reactivity of the wild form ofthe SOR, but also to explore the new reactivity of the mutations of the SOR and thus to define the crucial role of the second coordination sphere.

Den ökande koldioxidhalten i atmosfären har lett till ett ökat intresse för användning av solljus, ett exempel är inom kemi. Att använda ljus i kemi kan leda till nya och intressanta kemiska reaktioner och kan tillämpas i användningsområden som exempelvis fotokatalys och fotovoltaik. Genom att kombinera de kraftfulla egenskaperna av ljus i kemi med den höga selektivitet som karaktäriserar proteiner kan intressanta och kraftfulla maskinerier erhållas, så kallade fotoaktiva proteiner. Att studera de fotoinducerade beteendena i fotoaktiva proteiner är däremot icke-trivialt eftersom många utmaningar uppstår från ett experimentellt och teoretiskt perspektiv. Denna uppsats fokuserar på det teoretiska perspektivet. Simulationer av beteenden hos fotoaktiva proteiner leder till utmaningar relaterade till den stora variationen av tids- och längdskalor involverade, där allt ifrån ultrasnabba och lokaliserade excitationsprocesser till storskaliga strukturförändringar som utspelar sig på större tidsskalor. Detta, tillsammans med det faktum att det inte finns en "svart låda" som vi kan använda för nya teoretiska studier motiverar ett försiktigt tillvägagångssätt för teoretiska studier av fotoaktiva proteiner. I detta projekt görs en teoretisk undersökning av två nyligen framställda fotosensibiliserande proteiner (PSP2 och PSP3). De skapades för att härma huvudprocesserna av fotosyntesen i plantor. Fotosensibiliserande proteinerna skiljer sig med en enda aminosyra (position 203), som är en aspartat i PSP2 och en tyrosin i PSP3. Även fast proteinerna är strukturellt lika så uppvisar de olika fotoinducerade beteenden. Målet med denna studie är att klargöra för de mekanistiska detaljerna bakom dessa olikheter. Som ett första steg i att teoretiskt studera PSP2 och PSP3 upprättar vi ett protokoll för teoretisk förberedelse av proteinerna. Detta protokoll används sedan för att teoretiskt undersöka fotofysiken i dessa två proteiner med hjälp av kvantmekaniska/ molekylär mekaniska simuleringar. Våra resultat ger en första ordningens kartläggning av de exciterade tillstånden i de två proteinerna. Vi erhåller en intressant laddningsöverföring tillstånd i PSP3 som involverar tyrosinen och proteinets kromofor, med en energi lägre än det tillstånd ansvarig för absorption av ljus. Denna laddningsöverföring saknas i PSP2 och kan därför vara anledningen till de olika fotoinducerade beteendena i proteinerna. Detta kräver dock fortsatta studier av de två systemen. / The increase of carbon dioxide levels in the atmosphere have resulted in an increasing interest in using solar light for different purposes, one being chemistry. The introduction of light in chemistry can drive new and exciting chemical reactions, finding applications in many fields such as photocatalysis and photovoltaics. Combining the aspects of light in chemistry together with the high selectivity that characterizes proteins make for interesting and powerful machineries, so-called photoactive proteins. However, studying the behaviors of photoactive proteins is a non-trivial task and many challenges arise from both an experimental and a theoretical point of view. This thesis takes a theoretical perspective. The challenges associated with simulating photoactive protein behavior originate from the wide range of time and length scales involved, ranging from ultrafast and localized excitation processes to large-scale structural changes occurring on longer timescales. This, together with the fact that there is no black box that we can use for novel theoretical studies of photoactive proteins motivates a careful approach for theoretical studies of photoactive proteins. In this project, we conduct a theoretical investigation of two photosensitizer proteins (PSP2 and PSP3), recently engineered to capture the essence of plant photosynthesis. The photosensitizer proteins differ by only one residue (position 203), which is an aspartate in PSP2 while a tyrosine in PSP3. Although structurally similar, the proteins demonstrate different photoinduced behaviors. This study aims to shed light on the mechanistic details underlying these differences. As a first step to study PSP2 and PSP3 computationally, we develop a computational protocol for protein preparation. The protocol is then used for a theoretical investigation of the photophysics of the two proteins using quantum mechanics/molecular mechanics simulations. Our results provide a first mapping of the electronic-state manifold of the two proteins. We find an intriguing charge-transfer state in PSP3, involving the tyrosine and the protein chromophore, located below bright state responsible for light absorption. This state is absent in PSP2 and could therefore be the reason for the different photoinduced behavior of the two proteins. However, this requires further studies of the two systems.

No description available.

The relationship between protein structure, motions, and catalytic activity is an evolving perspective in enzymology. An interactive approach, where experimental and theoretical studies examine the same catalytic mechanism, is instrumental in addressing this issue. We combine various techniques, including steady state and pre-steady state kinetics, temperature dependence of kinetic isotope effects (KIEs), site-directed mutagenesis, X-ray crystallography, and quantum mechanics/molecular mechanics (QM/MM) calculations, to study the catalytic mechanisms of thymidylate synthase (TSase). Since TSase catalyzes the last step of the sole intracellular de novo synthesis of thymidylate (i.e. the DNA base T), it is a common target for antibiotic and anticancer drugs. The proposed catalytic mechanism for TSase comprises a series of bond cleavages and formations including activation of two C-H bonds: a rate-limiting C-H→C hydride transfer and a faster C-H→O proton transfer. This provides an excellent model system to examine the structural and dynamic effects of the enzyme on different C-H cleavage steps in the same catalyzed reaction. Our experiments found that the KIE on the hydride transfer is temperature independent while the KIE on the proton transfer is temperature dependent, implying the protein environment is better organized for H-tunneling in the former. Our QM/MM calculations revealed that the hydride transfer has a transition state (TS) that is invariable with temperature while the proton transfer has multiple subsets of TS structures, which corroborates with our experimental results. The calculations also suggest that collective protein motions rearrange the network of H-bonds to accompany structural changes in the ligands during and between chemical transformations. These computational results not only illustrate functionalities of specific protein residues that reconcile many previous experimental observations, but also provide guidance for future experiments to verify the proposed mechanisms. In addition, we conducted experiments to examine the importance of long-range interactions in TSase-catalyzed reaction, using both kinetic and structural analysis. Those experiments found that a remote mutation affects the hydride transfer by disrupting concerted protein motions, and Mg2+ binds to the surface of TSase and affects the hydride transfer at the interior active site. Both our experiments and computations have exposed interesting features of ecTSase that can potentially provide new targets for antibiotic drugs targeting DNA biosynthesis. The relationship between protein structure, motions, and catalytic activity learned from this project may have general implications to the question of how enzymes work.

Le projet de cette thèse était l'évaluation du processus de réactivation et l'étude du site actif de l'AChE inhibée par un agent neurotoxique par des méthode computationnelles. L'objectif était de guider le design rationnel de nouveau réactivateurs. Une étude initiale avec un modèle QM tronqué a indiqué la nécessité de modéliser l'environnement enzymatique pour compenser la charge du Glu334. Elle a aussi confirmé le rôle du trou oxyanionique dans la stabilisation des états de transition de la réactivation. Des simulations QM/MM de la réactivation par le réactivateur classique 2-PAM, ainsi que par deux réactivateurs au coeur aromatique non chargé ont été effectuées. Il a été démontré que le Glu202, un résidu à proximité de la triade catalytique de l'AChE, doit être protoné pour que la réactivation ait lieu. Ces simulations ont aussi montré que le réactivateur peut être déprotoné dans le site actif de l'AChE par His447. Les réactivateurs au coeur aromatique non chargé sont plus nucléophiles que la 2-PAM et l'un d'entre eux est plus aisément déprotoné dans le site actif. Nos résultats indiquent que la capacité d'un réactivateur à être facilement déprotoné est plus importante que sa nucléophilie. Enfin, un mécanisme de migration de protons a été identifié par des calculs QM/MM et EVB. Il implique deux glutamates derrière le site actif, Glu450 et Glu452. La possibilité que ces deux protons soient temporairement protonés et donc impliqués dans une migration de protons a été confirmé par des calculs CpHMD. La migration de proton passe par la N-protonation d'une liaison amide, ce qui constitue un nouveau mécanisme. / The project of this PhD was to investigate the reactivation process and the active site of nerve agent inhibited AChE by computational methodologies to gain insight about the rational design of new reactivators. An initial truncated QM model study provided some insight in the necessary compensation of Glu334 by the enzyme. It also confirmed the role of the oxyanionic hole in the stabilization of the transition state of the reactivation. QM/MM simulations of the reactivation with classical reactivator 2-PAM, as well as two non-pyridinium reactivators, were performed. It was shown that Glu202, a residue near the catalytic triad of AChE, needs to be protonated for the reactivation to occur. Those simulations also showed that the reactivator can be deprotonated in the active site of AChE by His447. Non-pyridinium reactivator were found to have a greater nucleophilicity than 2-PAM and, for one of them, to be easily deprotonated in the active site. Our results indicate that the capacity of a reactivator to be deprotonated in the active site of the enzyme is more important than its nucleophilicity. Finally, a proton relay mechanism was identified through QM/MM and EVB simulations. It involves two glutamate residues, Glu450 and Glu452, positioned behind the active site. The potential for these two residues to be transiently protonated and thus involved in a proton relay was confirmed by CpHMD simulations. This proton relay mechanism relies on the N-protonation of an amide which is a novel mechanism.

Protein pK(a)s and metalloprotein reduction potentials are studied with computational methodologies based on an ab initio quantum mechanics (QM) description of the protein and a linearized Poisson-Boltzmann Equation (LPBE) description of the solvent. The practical applicability of the QM/LPBE method is extended to proteins by using a QM description of the ionizable residue and a molecular mechanics (MM) description of the rest of the protein. This QM/MM/LPBE method is used to predict the pKa of Lys55 in the serine protease inhibitor turkey ovomucoid third domain (OMTKY3) and the prediction of 11.0 is in good agreement with the experimental value of 11.1. This is the first time a protein pKa value has been predicted with QM/MM methods. The QM/LPBE method is used to predict and interpret the pKa values of the five carboxyl residues (Asp7, Glu10, Glu19, Asp27, and Glu43) in OMTKY3. All the predicted pKa values are within 0.5 pH units of experiment, with a root mean square deviation of 0.31 pH units. We find that the decreased pKa values observed for some of the residues are primarily due to hydrogen bonds to the carboxyl oxygens. Hydrophobic effects are also shown to be important in raising the pKa. Interactions with charged residues are shown to have relatively little effect on the carboxyl pKa values in this protein, in general agreement with experiment. The relative Cu2+/Cu+ reduction potentials of six type-1 copper sites (cucumber stellacyanin, P. aeruginosa azurin, poplar plastocyanin, C. cinereus laccase, T. ferrooxidans rusticyanin and human ceruloplasmin), which lie in a reduction potential range from 260 mV to over 1000 mV, have been studied with the QM/LPBE method. For the first time, the range and relative orderings of the reduction potentials are reproduced well compared to experimental values. The study suggests that the main interactions determing the relative reduction potentials of blue copper sites are located within 6 Å of the Cu atoms. Further analysis suggests that the reduction potential differences of type-1 copper sites are caused by axial ligand interactions, hydrogen bonding to the S(Cys), and protein constraints on the inner sphere ligand orientations.

Ces dernières années, les méthodes hybrides QM/MM combinant la mécanique quantique (QM) et la mécanique moléculaire (MM) se sont revélées des méthodes de choix pour l'étude de systèmes chimiques et biochimiques contenant plus d'une centaine d'atomes. Nous avons mis en évidence les apports et difficultés liés à leur utilisation à travers des systèmes variés: modélisation de ligands phosphines, réactivité d'une protéine de type acide phosphatase pourpre, modélisation de l'interaction protéine-ligand.

Phytochromes constitute one of the six well-characterized families of photosensory proteins in Nature. From the viewpoint of computational modeling, however, phytochromes have been the subject of much fewer studies than most other families of photosensory proteins, which is likely a consequence of relevant high-resolution structural data becoming available only in recent years. In this work, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are used to calculate UV-vis absorption spectra of Deinococcus radiodurans bacteriophytochrome. We investigate how the choice of QM/MM methodology affects the resulting spectra and demonstrate that QM/MM methods can reproduce the experimental absorption maxima of both the Q and Soret bands with an accuracy of about 0.15 eV. Furthermore, we assess how the protein environment influences the intrinsic absorption of the bilin chromophore, with particular focus on the Q band underlying the primary photochemistry of phytochromes.