Elsa Sanchez-Garcia

Enzyme-Directed Mutasynthesis: A Combined Experimental and Theoretical Approach to Substrate Recognition of a Polyketide Synthase

Acyltransferase domains control the extender unit recognition in Polyketide Synthases (PKS) and thereby the side-chain diversity of the resulting natural products. The enzyme engineering strategy presented here allows the alteration of the acyltransferase substrate profile to enable an engineered biosynthesis of natural product derivatives through the incorporation of a synthetic malonic acid thioester. Experimental sequence-function correlations combined with computational modeling revealed the origins of substrate recognition in these PKS domains and enabled a targeted mutagenesis. We show how a single point mutation was able to direct the incorporation of a malonic acid building block with a non-native functional group into erythromycin. This approach, introduced here as enzyme-directed mutasynthesis, opens a new field of possibilities beyond the state of the art for the combination of organic chemistry and biosynthesis toward natural product analogues.

Switching the Spin State of Diphenylcarbene via Halogen Bonding.

The interactions between diphenylcarbene DPC and the halogen bond donors CF3I and CF3Br were investigated using matrix isolation spectroscopy (IR, UV-vis, and EPR) in combination with QM and QM/MM calculations. Both halogen bond donors CF3X form very strong complexes with the singlet state of DPC, but only weakly interact with triplet DPC. This results in a switching of the spin state of DPC, the singlet complexes becoming more stable than the triplet complexes. CF3I forms a second complex (type II) with DPC that is thermodynamically slightly more stable. Calculations predict that in this second complex the DPC···I distance is shorter than the F3C···I distance, whereas in the first (type I) complex the DPC···I distance is, as expected, longer. CF3Br only forms the type I complex. Upon irradiation I or Br, respectively, are transferred to the DPC carbene center and radical pairs are formed. Finally, on annealing, the formal C-X insertion product of DPC is observed. Thus, halogen bonding is a powerful new principle to control the spin state of reactive carbenes.

QM/MM Study of the Absorption Spectra of DsRed.M1 Chromophores

We report geometries and vertical excitation energies for the red and green chromophores of the DsRed.M1 protein in the gas phase and in the solvated protein environment. Geometries are optimized using density functional theory (DFT, B3LYP functional) for the isolated chromophores and combined quantum mechanical/molecular mechanical (QM/MM) methods for the protein (B3LYP/MM). Vertical excitation energies are computed using DFT/MRCI, OM2/MRCI, and TDDFT as QM methods. In the case of the red chromophore, there is a general blue shift in the excitation energies when going from the isolated chromophore to the protein, which is caused both by structural changes and by electrostatic interactions with the environment. For the lowest ππ* transition, these two factors contribute to a similar extent to the overall DFT/MRCI shift of 0.4 eV. An enlargement of the QM region to include active‐site residues does not change the DFT/MRCI excitation energies much. The DFT/MRCI results are closest to experiment for both chromophores. OM2/MRCI and TDDFT overestimate the first vertical excitation energy by 0.3–0.5 and 0.2–0.4 eV, respectively, relative to the experimental or DFT/MRCI values. The experimental gap of 0.35 eV between the lowest ππ* excitation energies of the red (cis‐acylimine) and green (trans‐peptide) forms is well reproduced by DFT/MRCI and TDDFT (0.32 and 0.37 eV, respectively). A histogram spectrum for an equal mixture of the two forms, generated by OM2/MRCI calculations on 450 snapshots along molecular dynamics trajectories, matches the experimental spectrum quite well, with a gap of 0.23 eV and an overall blue shift of about 0.3 eV. DFT/MRCI appears as an attractive choice for calculating excitation energies in fluorescent proteins, without the shortcomings of TDDFT and computationally more affordable than CASSCF‐based approaches.

A molecular tweezer antagonizes seminal amyloids and HIV infection

Semen is the main vector for HIV transmission and contains amyloid fibrils that enhance viral infection. Available microbicides that target viral components have proven largely ineffective in preventing sexual virus transmission. In this study, we establish that CLR01, a ‘molecular tweezer’ specific for lysine and arginine residues, inhibits the formation of infectivity-enhancing seminal amyloids and remodels preformed fibrils. Moreover, CLR01 abrogates semen-mediated enhancement of viral infection by preventing the formation of virion–amyloid complexes and by directly disrupting the membrane integrity of HIV and other enveloped viruses. We establish that CLR01 acts by binding to the target lysine and arginine residues rather than by a non-specific, colloidal mechanism. CLR01 counteracts both host factors that may be important for HIV transmission and the pathogen itself. These combined anti-amyloid and antiviral activities make CLR01 a promising topical microbicide for blocking infection by HIV and other sexually transmitted viruses. DOI: http://dx.doi.org/10.7554/eLife.05397.001

Borazine and Benzene Homo- and Heterodimers

The homodimers of benzene and borazine as well as a heterodimer consisting of one benzene (bz) and one borazine (bor) molecule are investigated using MP2, SCS-MP2, and CCSD(T) theories in conjunction with basis sets of up to quadruple-zeta quality. Dimer geometries were completely optimized using the resolution of the identity approximation of MP2 with a QZVPP basis set and characterized by computation of harmonic vibrational frequencies using triple-zeta basis sets. While significant higher order correlation effects beyond MP2 are important for the benzene dimer, these are very small for the borazine dimer and intermediate for the heterodimer. The spin-component scaling (SCS) correction of MP2 produces binding energies for the borazine dimer that are too low but yields very good agreement with CCSD(T) for the heterodimer. The decrease in the intermolecular distance in the sandwich (S) configurations from bz(2) via bz-bor to bor(2) is accompanied by an increased binding energy and a change from second-order stationary points to a minimum for bor(2). The T isomer is less stable than the S configuration for bor(2), but it is preferred over the S and a parallel-displaced (PD) arrangement in the heterodimer. The following order of stability is obtained for the minima at the extrapolated CCSD(T) level: T(bz-bor) > S(bor(2)) > PD(bz-bor) > PD(bor(2)) > T(bor(2)) > PD(bz(2)). The most stable isomer at all levels of theory, T(bz-bor), features a NH...pi interaction.

Interaction and Reaction of the Phenyl Radical with Water: A Source of OH Radicals

Thats radical! A photochemical reaction between the phenyl radical and water results in the abstraction of a hydrogen atom from water and the formation of a hydroxyl radical. The hydroxyl radical forms an OHpi hydrogen bond with benzene (see picture) and does not react with benzene thermally under the conditions of matrix isolation.

Predicted incorporation of non-native substrates by a polyketide synthase yields bioactive natural product derivatives

The polyether ionophore monensin is biosynthesized by a polyketide synthase that delivers a mixture of monensins A and B by the incorporation of ethyl‐ or methyl‐malonyl‐CoA at its fifth module. Here we present the first computational model of the fifth acyltransferase domain (AT5mon) of this polyketide synthase, thus affording an investigation of the basis of the relaxed specificity in AT5mon, insights into the activation for the nucleophilic attack on the substrate, and prediction of the incorporation of synthetic malonic acid building blocks by this enzyme. Our predictions are supported by experimental studies, including the isolation of a predicted derivative of the monensin precursor premonensin. The incorporation of non‐native building blocks was found to alter the ratio of premonensins A and B. The bioactivity of the natural product derivatives was investigated and revealed binding to prenyl‐binding protein. We thus show the potential of engineered biosynthetic polyketides as a source of ligands for biological macromolecules.

Quantum refinement of protein structures: implementation and application to the red fluorescent protein DsRed.M1.

Quantum refinement is an improvement upon the molecular mechanics (MM)-based crystallographic refinement. In the latter, X-ray data are supplemented with additional chemical information through MM force fields, whereas quantum refinement describes crucial regions of interest in the macromolecule by quantum mechanics (QM) instead of MM. In this paper, we report the implementation of quantum refinement in the ChemShell QM/MM framework and its application in an investigation of the chromophore structure of the red fluorescent protein DsRed.M1. Both mechanical and electrostatic QM/MM embedding schemes are implemented and tested. In the quantum refinement of DsRed.M1, the anionic red acylimine chromophore adopts a nearly orthogonal arrangement (rather than a cis or trans form), and the bond lengths in the acylimine moiety are more consistent with a phenolate (rather than a quinoid) structure. These findings are in contrast to the structure deduced from a standard crystallographic refinement (PDB: 2VAD), but in agreement with our earlier results from a purely theoretical QM/MM study. On the other hand, the quantum refinement of the anionic acylimine form of DsRed.M1 yields a hydrogen bonding network around the chromophore, especially with regard to the arrangement of the water molecules and the Glu148 residue, that is closer to the 2VAD structure than to the previously optimized QM/MM structure. In our earlier study the initial classical molecular dynamics (MD) simulations during QM/MM setup apparently exaggerated the mobility of the water molecules around the chromophore. On the basis of the present results, it seems likely that the Glu148 residue is protonated in the DsRed.M1 protein. The calculation of electronic excitation energies allows for further assessment of the proposed structures, especially in the chromophore region. Using a combination of density functional theory and multireference configuration interaction (DFT/MRCI), we find excellent agreement between experiment and theory only for the structures obtained from quantum refinement and from QM/MM optimization, but not for the 2VAD structure. The present case study on DsRed.M1 thus demonstrates the merits of combining reliable theoretical and experimental information in the determination of protein structures.

The Phenoxyl Radical–Water Complex—A Matrix Isolation and Computational Study

The phenoxyl radical 1 was generated in high yields by flash vacuum pyrolysis of allyl phenyl ether 2 with subsequent trapping of the products in argon at 3 K. In water-doped argon matrices, an OH···O complex between 1 and water is formed that could be characterized by IR spectroscopy. Several isotopomers of the complex were generated, and the IR spectra compared to results of density functional theory calculations. Other dimers between 1 and water were not found under these conditions. QM/MM calculations in simulated argon matrices reveal that an OH···π complex is unstable even at a time scale of picoseconds. This finding has implications on the related interaction between the tyrosyl radical and the water in biological systems.

Hybrid Quantum Mechanics/Molecular Mechanics/Coarse Grained Modeling: A Triple-Resolution Approach for Biomolecular Systems.

We present a hybrid quantum mechanics/molecular mechanics/coarse-grained (QM/MM/CG) multiresolution approach for solvated biomolecular systems. The chemically important active-site region is treated at the QM level. The biomolecular environment is described by an atomistic MM force field, and the solvent is modeled with the CG Martini force field using standard or polarizable (pol-CG) water. Interactions within the QM, MM, and CG regions, and between the QM and MM regions, are treated in the usual manner, whereas the CG-MM and CG-QM interactions are evaluated using the virtual sites approach. The accuracy and efficiency of our implementation is tested for two enzymes, chorismate mutase (CM) and p-hydroxybenzoate hydroxylase (PHBH). In CM, the QM/MM/CG potential energy scans along the reaction coordinate yield reaction energies that are too large, both for the standard and polarizable Martini CG water models, which can be attributed to adverse effects of using large CG water beads. The inclusion of an atomistic MM water layer (10 Å for uncharged CG water and 5 Å for polarizable CG water) around the QM region improves the energy profiles compared to the reference QM/MM calculations. In analogous QM/MM/CG calculations on PHBH, the use of the pol-CG description for the outer water does not affect the stabilization of the highly charged FADHOOH-pOHB transition state compared to the fully atomistic QM/MM calculations. Detailed performance analysis in a glycine-water model system indicates that computation times for QM energy and gradient evaluations at the density functional level are typically reduced by 40-70% for QM/MM/CG relative to fully atomistic QM/MM calculations.