Ultrafast photooxidation of protein-bound anionic flavin radicals.

The photophysical properties of anionic semireduced flavin radicals are largely unknown despite their importance in numerous biochemical reactions. Here, we studied the photoproducts of these intrinsically unstable species in five different flavoprotein oxidases where they can be stabilized, including the well-characterized glucose oxidase. Using ultrafast absorption and fluorescence spectroscopy, we unexpectedly found that photoexcitation systematically results in the oxidation of protein-bound anionic flavin radicals on a time scale of less than ∼100 fs. The thus generated photoproducts decay back in the remarkably narrow 10- to 20-ps time range. Based on molecular dynamics and quantum mechanics computations, positively charged active-site histidine and arginine residues are proposed to be the electron acceptor candidates. Altogether, we established that, in addition to the commonly known and extensively studied photoreduction of oxidized flavins in flavoproteins, the reverse process (i.e., the photooxidation of anionic flavin radicals) can also occur. We propose that this process may constitute an excited-state deactivation pathway for protein-bound anionic flavin radicals in general. This hitherto undocumented photochemical reaction in flavoproteins further extends the family of flavin photocycles.

Catalytic Mechanism of Fatty Acid Photodecarboxylase: On the Detection and Stability of the Initial Carbonyloxy Radical Intermediate.

In fatty acid photodecarboxylase (FAP), light-induced formation of the primary radical product RCOO⋅ from fatty acid RCOO- occurs in 300 ps, upon which CO2 is released quasi-immediately. Based on the hypothesis that aliphatic RCOO⋅ (spectroscopically uncharacterized because unstable) absorbs in the red similarly to aromatic carbonyloxy radicals such as 2,6-dichlorobenzoyloxy radical (DCB⋅), much longer-lived linear RCOO⋅ has been suggested recently. We performed quantum chemical reaction pathway and spectral calculations. These calculations are in line with the experimental DCB⋅ decarboxylation dynamics and spectral properties and show that in contrast to DCB⋅, aliphatic RCOO⋅ radicals a) decarboxylate with a very low energetic barrier and on the timescale of a few ps and b) exhibit little red absorption. A time-resolved infrared spectroscopy experiment confirms very rapid, ≪300 ps RCOO⋅ decarboxylation in FAP. We argue that this property is required for the observed high quantum yield of hydrocarbons formation by FAP.

Photoconduction and Electroluminescence of Copper (II) Protoporphyrin and Chlorin Cu-C-e6.

Cu (II) protoporphyrin Cu-PP-IX and chlorin Cu-C-e6 were found to have both thin solid film formation and charge carrier transport abilities. In the layers deposited by resistive thermal evaporation, the mobilities of holes and electrons are on the order of 10-5 cm2 V-1 s-1. Organic light-emitting diodes incorporating the dye molecules as emitting dopants demonstrate electroluminescence in the UV and near-IR ranges.

Photochemical and Molecular Dynamics Studies of Halide Binding in Flavoenzyme Glucose Oxidase.

Glucose oxidase (GOX), a characteristic flavoprotein oxidase with widespread industrial applications, binds fluoride (F- ) and chloride (Cl- ). We investigated binding properties of halide inhibitors of GOX through time-resolved spectral characterization of flavin-related photochemical processes and molecular dynamic simulations. Cl- and F- bind differently to the protein active site and have substantial but opposite effects on the population and decay of the flavin excited state. Cl- binds closer to the flavin, whose excited-state decays in <100 fs due to anion-π interactions. Such interactions appear absent in F- binding, which, however, significantly increases the active-site rigidity leading to more homogeneous, picosecond fluorescence decay kinetics. These findings are discussed in relation to the mechanism of halide inhibition of GOX by occupying the accommodation site of catalytic intermediates and increasing the active-site rigidity.

Parametrization of Force Field Bonded Terms under Structural Inconsistency.

Parametrization of the bonded part of molecular mechanics (MM) force fields (FFs) is typically done by fitting reference quantum mechanical (QM) energies or forces of representative structures. FFs for small molecules are constructed in incremental parametrization procedures, where parameters developed previously are retained for novel molecules, followed by optimization of missing, not previously optimized parameters. Equilibrium QM and MM geometries of molecules can deviate due to parameters transferred from existing molecules in the FF. In this work, we demonstrate that conventional parametrization methods based on fitting QM energies and/or forces to derive parameters for bond and angle terms produce largely suboptimal force constants when MM and QM equilibrium structures deviate. We further developed and tested a new method to derive CHARMM FF parameters based on the potential energy surface scans where a structural deviation between QM and MM optimized geometries is explicitly allowed during parametrization. The test of the new method was performed on a diverse set of 32 molecules. The results show that without any need for additional restraints, the new method produces robust and largely transferable parameters for bond and angle terms. The new method also improves the agreement for the normal modes for all molecules in the test set, reducing the average error in the reproduction of QM normal mode frequencies from 9.5% computed with CGenFF parameters to 6.8% computed with the new parameters. The new method will allow parametrization of molecules under structural deviations, common for force fields for small molecules, producing robust and transferable parameters.

The effect of PbS quantum dots on molecular dynamics and conductivity of PTB7:PC71BM bulk heterojunction as revealed by dielectric spectroscopy.

A ternary photovoltaic blend containing the PTB7 donor component, the PC71BM acceptor component, and colloidal quantum dots of lead sulfide (PbS QDs) was investigated using broadband dielectric spectroscopy. In the dielectric loss spectrum of PTB7:PC71BM:PbS QDs, γ- and ß-relaxation processes in PTB7 were recognized and analyzed in terms of Arrhenius-type equations. To elucidate the effect of PbS QDs on molecular dynamics of PTB7, the activation energies of both processes were evaluated and compared with those obtained for the binary PTB7:PC71BM blend. Using the CELIV method, the charge carrier mobility was estimated. The PbS QD incorporation into the binary blend was shown to decrease both electron and hole mobility in the ternary PTB7:PC71BM:PbS QD blend. For evaluating the charge carrier lifetime in the ternary blend, the Cole-Cole diagrams derived from the dc conductivity data were plotted. The charge carrier lifetime was found to be much less than the hole extraction time, thus providing effective accumulation of charge carries at the electrodes in the ternary blend under investigation.

Copper(II) Etioporphyrinate as a Promising Photoluminescent and Electroluminescent Temperature Sensor.

Luminescent temperature sensors are of great interest because they allow remote determination of temperature in transparent media, such as living tissues, as well as on scattering or transparent surfaces of materials. This study analyzes the luminescent properties of copper(II) etioporphyrinate (Cu-EtioP) in a polystyrene film upon variation of temperature from -195 °C to +65 °C in a cryostat. It is shown that the ratio of intensities of phosphorescence transitions in the red spectral region of such a material varies significantly, that is, the material has thermosensory properties. The phosphorescence decay curves of copper(II) etioporphyrinate in a polystyrene film are analyzed. The quantum yield of phosphorescence of copper(II) etioporphyrinate determined by the absolute method was 3.15%. It was also found that the electroluminescence (EL) spectra of copper(II) etioporphyrinate in a poly(9-vinylcarbazole) (PVK) matrix demonstrated a similar change in the spectra in the temperature range -3 °C to +80 °C. That is, copper(II) etioporphyrinate can also be used as a luminescent temperature sensor as part of an active OLED layer.

Cyclodipeptide Synthases of the NYH Subfamily Recognize tRNA Using an α-Helix Enriched with Positive Residues.

Cyclodipeptide synthases (CDPSs) perform nonribosomal protein synthesis using two aminoacyl-tRNA substrates to produce cyclodipeptides. At present, there are no structural details of the CDPS:tRNA interaction available. Using AlbC, a CDPS that produces cyclo(l-Phe-l-Phe), the interaction between AlbC and its Phe-tRNA substrate was investigated. Simulations of models of the AlbC:tRNA complex, proposed by rigid-body docking or homology modeling, demonstrated that interactions with residues of an AlbC α-helix, α4, significantly contribute to the free energy of binding of AlbC to tRNA. Individual residue contributions to the tRNA binding free energy of the discovered binding mode explain well the available biochemical data, and the results of in vivo assay experiments performed in this work and guided by simulations. In molecular dynamics simulations, the phenylalanyl group predominantly occupied the two positions observed in the experimental structure of AlbC in the dipeptide intermediate state, suggesting that tRNAs of the first and second substrates interact with AlbC in a similar manner. Overall, given the high degree of sequence and structural similarity among the members of the CDPS NYH protein subfamily, the mechanism of the protein:tRNA interaction is expected to be pertinent to a wide range of proteins interacting with tRNA.

Characterization of Light-Induced, Short-Lived Interacting Radicals in the Active Site of Flavoprotein Ferredoxin-NADP+ Oxidoreductase.

Radicals of flavin adenine dinucleotide (FAD), as well as tyrosine and tryptophan, are widely involved as key reactive intermediates during electron-transfer (ET) reactions in flavoproteins. Due to the high reactivity of these species and their corresponding short lifetime, characterization of these intermediates in functional processes of flavoproteins is usually challenging but can be achieved by ultrafast spectroscopic studies of light-activatable flavoproteins. In ferredoxin-NADP+ oxidoreductase from Bacillus subtilis (BsFNR), fluorescence of the FAD cofactor that very closely interacts with a neighboring tyrosine residue (Tyr50) is strongly quenched. Here we study short-lived photoproducts of this enzyme and its variants, with Tyr50 replaced by tryptophan or glycine. Using time-resolved fluorescence and absorption spectroscopies, we show that, upon the excitation of WT BsFNR, ultrafast ET from Tyr50 to the excited FAD cofactor occurs in ∼260 fs, an order of magnitude faster than the decay by charge recombination, facilitating the characterization of the reaction intermediates in the charge-separated state with respect to other recently studied systems. These studies are corroborated by experiments on the Y50W mutant protein, which yield photoproducts qualitatively similar to those observed in other tryptophan-bearing flavoproteins. By combining the experimental results with molecular dynamics simulations and quantum mechanics calculations, we investigate in detail the effects of protein environment and relaxations on the spectral properties of those radical intermediates and demonstrate that the spectral features of radical anionic FAD are highly sensitive to its environment, and in particular to the dynamics and nature of the counterions formed in the photoproducts. Altogether, comprehensive characterizations are provided for important radical intermediates that are generally involved in functional processes of flavoproteins.

Photochemical processes in flavo-enzymes as a probe for active site dynamics: TrmFO of Thermus thermophilus.

Quenching of flavin fluorescence by electron transfer from neighboring aromatic residues is ubiquitous in flavoproteins. Apart from constituting a functional process in specific light-active systems, time-resolved spectral characterization of the process can more generally be employed as a probe for the active site configuration and dynamics. In the C51A variant of the bacterial RNA-transforming flavoenzyme TrmFO from the bacterium Thermus thermophilus, fluorescence is very short-lived (~ 1 ps), and close-by Tyr343 is known to act as the main quencher, as confirmed here by the very similar dynamics observed in protein variants with modified other potential quenchers, Trp283 and Trp214. When Tyr343 is modified to redox-inactive phenylalanine, slower and highly multiphasic kinetics are observed on the picosecond-nanosecond timescale, reflecting heterogeneous electron donor-acceptor configurations. We demonstrate that Trp214, which is located on a potentially functional flexible loop, contributes to electron donor quenching in this variant. Contrasting with observations in other nucleic acid-transforming enzymes, these kinetics are strikingly temperature-independent. This indicates (a) near-barrierless electron transfer reactions and (b) no exchange between different configurations on the timescale up to at least 2 ns, despite the presumed flexibility of Trp214. Results of extensive molecular dynamics simulations are presented to explain this unexpected finding in terms of slowly exchanging protein configurations.

Ion-Driven Self-Assembly of Lanthanide Bis-phthalocyaninates into Conductive Quasi-MOF Nanowires: an Approach toward Easily Recyclable Organic Electronics.

Controlled self-assembly and rapid disintegration of supramolecular nanowires is potentially useful for ecology-friendly organic electronics. Herein, a novel method exploiting the binding between crown-substituted double-decker lanthanide phthalocyaninates (ML2, M = Lu, Ce, Tb) and K+ ions is applied for the one-step fabrication of macroscopically long conductive one-dimensional quasi-metal-organic frameworks. Their properties are controlled by the size of the lanthanide ion guiding the assembly through either intra- or intermolecular interactions. A LuL2 linker with a small interdeck distance yields fully conjugated intermolecular-bonded K+-LuL2 nanowires with a thickness of 10-50 nm, a length of up to 50 µm, and a conductivity of up to 11.4 S cm-1, the highest among them being reported for phthalocyanine assemblies. The large size of CeL2 and TbL2 leads to the formation of mixed intra- and intermolecular K+-ML2 phases with poor electric properties. A field-assisted method is developed to deposit aligned conductive K+-LuL2 assemblies on solids. The solid-supported nanowires can be disintegrated into starting components in a good aprotic solvent for further recycling.

New Unsymmetrically Substituted Benzothiadiazole-Based Luminophores: Synthesis, Optical, Electrochemical Studies, Charge Transport, and Electroluminescent Characteristics.

Three new benzothiadiazole (BTD)-containing luminophores with different configurations of aryl linkers have been prepared via Pd-catalyzed cross-coupling Suzuki and Buchwald-Hartwig reactions. Photophysical and electroluminescent properties of the compounds were investigated to estimate their potential for optoelectronic applications. All synthesized structures have sufficiently high quantum yields in film. The BTD with aryl bridged carbazole unit demonstrated the highest electrons and holes mobility in a series. OLED with light-emitting layer (EML) based on this compound exhibited the highest brightness, as well as current and luminous efficiency. The synthesized compounds are not only luminophores with a high photoluminescence quantum yield, but also active transport centers for charge carriers in EML of OLED devices.

Mechanism of Naphthoquinone Selectivity of Thymidylate Synthase ThyX.

Naphthoquinones (NQs) are natural and synthetic compounds with a wide range of biological activities commonly attributed to their redox activity and/or chemical reactivity. However, genetic and biochemical experiments have recently demonstrated that 2-hydroxy-NQs (2-OH-NQs) act as highly specific noncovalent inhibitors of the essential bacterial thymidylate synthase ThyX in a cellular context. We used biochemical experiments and molecular dynamics simulations to elucidate the selective inhibition mechanism of NQ inhibitors of ThyX from Mycobacterium tuberculosis (Mtb). Free energy simulations rationalized how ThyX recognizes the natural substrate dUMP in the N3-ionized form using an arginine, Arg199, in Mtb. The results further demonstrated that 2-OH-NQ, similar to dUMP, binds to ThyX in the ionized form, and the strong and selective binding of 2-OH-NQ to ThyX is also explained by electrostatic interactions with Arg199. The stronger binding of the close analog 5F-dUMP to ThyX and its inhibitory properties compared with dUMP were explained by the stronger acidity of the uracil N3 atom. Our results, therefore, revealed that the ionization of 2-OH-NQs drives their biological activities by mimicking the interactions with the natural substrate. Our observations encourage the rational design of optimized ThyX inhibitors that ultimately may serve as antibiotics.

A Molecular Mechanics Model for Flavins.

Flavin containing molecules form a group of important cofactors that assist a wide range of enzymatic reactions. Flavins use the redox-active isoalloxazine system, which is capable of one- and two-electron transfer reactions and can exist in several protonation states. In this work, molecular mechanics force field parameters compatible with the CHARMM36 all-atom additive force field were derived for biologically important flavins, including riboflavin, flavin mononucleotide, and flavin adenine dinucleotide. The model was developed for important protonation and redox states of the isoalloxazine group. The partial charges were derived using the CHARMM force field parametrization strategy, where quantum mechanics water-solute interactions are used to target optimization. In addition to monohydrate energies and geometries, electrostatic potential around the compound was used to provide additional restraints during the charge optimization. Taking into account the importance of flavin-containing molecules special attention was given to the quality of bonded terms. All bonded terms, including stiff terms and torsion angle parameters, were parametrized using exhaustive potential energy surface scans. In particular, the model reproduces well the butterfly motion of isoalloxazine in the oxidized and reduced forms as predicted by quantum mechanics in gas phase. The model quality is illustrated by simulations of four flavoproteins. Overall, the presented molecular mechanics model will be of utility to model flavin cofactors in different redox states. © 2019 Wiley Periodicals, Inc.

The structure of an E. coli tRNAfMet A1-U72 variant shows an unusual conformation of the A1-U72 base pair.

Translation initiation in eukaryotes and archaea involves a methionylated initiator tRNA delivered to the ribosome in a ternary complex with e/aIF2 and GTP. Eukaryotic and archaeal initiator tRNAs contain a highly conserved A1-U72 base pair at the top of the acceptor stem. The importance of this base pair to discriminate initiator tRNAs from elongator tRNAs has been established previously using genetics and biochemistry. However, no structural data illustrating how the A1-U72 base pair participates in the accurate selection of the initiator tRNAs by the translation initiation systems are available. Here, we describe the crystal structure of a mutant E. coli initiator tRNAfMetA1-U72, aminoacylated with methionine, in which the C1:A72 mismatch at the end of the tRNA acceptor stem has been changed to an A1-U72 base pair. Sequence alignments show that the mutant E. coli tRNA is a good mimic of archaeal initiator tRNAs. The crystal structure, determined at 2.8 Å resolution, shows that the A1-U72 pair adopts an unusual arrangement. A1 is in a syn conformation and forms a single H-bond interaction with U72 This interaction requires protonation of the N1 atom of A1 Moreover, the 5' phosphoryl group folds back into the major groove of the acceptor stem and interacts with the N7 atom of G2 A possible role of this unusual geometry of the A1-U72 pair in the recognition of the initiator tRNA by its partners during eukaryotic and archaeal translation initiation is discussed.

Combining the polarizable Drude force field with a continuum electrostatic Poisson-Boltzmann implicit solvation model.

In this work, we have combined the polarizable force field based on the classical Drude oscillator with a continuum Poisson-Boltzmann/solvent-accessible surface area (PB/SASA) model. In practice, the positions of the Drude particles experiencing the solvent reaction field arising from the fixed charges and induced polarization of the solute must be optimized in a self-consistent manner. Here, we parameterized the model to reproduce experimental solvation free energies of a set of small molecules. The model reproduces well-experimental solvation free energies of 70 molecules, yielding a root mean square difference of 0.8 kcal/mol versus 2.5 kcal/mol for the CHARMM36 additive force field. The polarization work associated with the solute transfer from the gas-phase to the polar solvent, a term neglected in the framework of additive force fields, was found to make a large contribution to the total solvation free energy, comparable to the polar solute-solvent solvation contribution. The Drude PB/SASA also reproduces well the electronic polarization from the explicit solvent simulations of a small protein, BPTI. Model validation was based on comparisons with the experimental relative binding free energies of 371 single alanine mutations. With the Drude PB/SASA model the root mean square deviation between the predicted and experimental relative binding free energies is 3.35 kcal/mol, lower than 5.11 kcal/mol computed with the CHARMM36 additive force field. Overall, the results indicate that the main limitation of the Drude PB/SASA model is the inability of the SASA term to accurately capture non-polar solvation effects. © 2018 Wiley Periodicals, Inc.

Identification of a second GTP-bound magnesium ion in archaeal initiation factor 2.

Eukaryotic and archaeal translation initiation processes involve a heterotrimeric GTPase e/aIF2 crucial for accuracy of start codon selection. In eukaryotes, the GTPase activity of eIF2 is assisted by a GTPase-activating protein (GAP), eIF5. In archaea, orthologs of eIF5 are not found and aIF2 GTPase activity is thought to be non-assisted. However, no in vitro GTPase activity of the archaeal factor has been reported to date. Here, we show that aIF2 significantly hydrolyses GTP in vitro. Within aIF2γ, H97, corresponding to the catalytic histidine found in other translational GTPases, and D19, from the GKT loop, both participate in this activity. Several high-resolution crystal structures were determined to get insight into GTP hydrolysis by aIF2γ. In particular, a crystal structure of the H97A mutant was obtained in the presence of non-hydrolyzed GTP. This structure reveals the presence of a second magnesium ion bound to GTP and D19. Quantum chemical/molecular mechanical simulations support the idea that the second magnesium ion may assist GTP hydrolysis by helping to neutralize the developing negative charge in the transition state. These results are discussed in light of the absence of an identified GAP in archaea to assist GTP hydrolysis on aIF2.

Insights into molecular plasticity in protein complexes from Trm9-Trm112 tRNA modifying enzyme crystal structure.

Most of the factors involved in translation (tRNA, rRNA and proteins) are subject to post-transcriptional and post-translational modifications, which participate in the fine-tuning and tight control of ribosome and protein synthesis processes. In eukaryotes, Trm112 acts as an obligate activating platform for at least four methyltransferases (MTase) involved in the modification of 18S rRNA (Bud23), tRNA (Trm9 and Trm11) and translation termination factor eRF1 (Mtq2). Trm112 is then at a nexus between ribosome synthesis and function. Here, we present a structure-function analysis of the Trm9-Trm112 complex, which is involved in the 5-methoxycarbonylmethyluridine (mcm(5)U) modification of the tRNA anticodon wobble position and hence promotes translational fidelity. We also compare the known crystal structures of various Trm112-MTase complexes, highlighting the structural plasticity allowing Trm112 to interact through a very similar mode with its MTase partners, although those share less than 20% sequence identity.

Electrostatic free energies in translational GTPases: Classic allostery and the rest.

GTPases typically switch between an inactive, OFF conformation and an active, ON conformation when a GDP ligand is replaced by GTP. Their ON/OFF populations and activity thus depend on the stabilities of four protein complexes, two apo-protein forms, and GTP/GDP in solution. A complete characterization is usually not possible experimentally and poses major challenges for simulations. We review the most important methodological challenges and we review thermodynamic data for two GTPases involved in translation of the genetic code: archaeal Initiation Factors 2 and 5B (aIF2, aIF5B). One main challenge is the multiplicity of states and conformations, including those of GTP/GDP in solution. Another is force field accuracy, especially for interactions of GTP/GDP with co-bound divalent Mg(2+) ions. The calculation of electrostatic free energies also poses specific challenges, and requires careful protocols. For aIF2, experiments and earlier simulations showed that it is a "classic" GTPase, with distinct ON/OFF conformations that prefer to bind GTP and GDP, respectively. For aIF5B, we recently proposed a non-classic mechanism, where the ON/OFF states differ only in the protonation state of Glu81 in the nucleotide binding pocket. This model is characterized here using free energy simulations. The methodological analysis should help future studies, while the aIF2, aIF5B examples illustrate the diversity of ATPase/GTPase mechanisms. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Additive CHARMM force field for naturally occurring modified ribonucleotides.

More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all-atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs.