Phosphoryl transfer enzymes and hypervalent phosphorus chemistry

Synthesis and characterization of a pair of O-fac/O-mer 12-P-6 alkyloxaphosphates with a P-O-C-C four-membered ring

The synthesis and characterization of a stereoisomer pair of 12-P-6 alkyloxaphosphates synthesised from pentacoordinate phosphoranes.

Structurally characterized hexacoordinate organophosphorus compounds remain rare due to their highly reactive nature and thermal instability. Herein we report the first synthesis of a pair of O-facial and O-meridional hexacoordinate oxaphosphates (5B and 5D) obtained from the O-apical and O-equatorial β-hydroxyalkylphosphoranes 3 and 4. This was achieved by using the bulky C₂F₅-groups on the ortho-substituted aryl backbone. Calculations of the relative energies of possible isomers indicate 5B and 5D are thermodynamic products. Although the mechanisms of their formation and the determining factor of stereo-selectivity are still unclear, their isolation and structure conformation contributes to a formulation of a viable strategy for diastereoselective synthesis of heteroleptic hexacoordinate organophosphates.

Nucleophilic Substitution in Solution: Activation Strain Analysis of Weak and Strong Solvent Effects

We have quantum chemically studied the effect of various polar and apolar solvents on the shape of the potential energy surface (PES) of a diverse collection of archetypal nucleophilic substitution reactions at carbon, silicon, phosphorus, and arsenic by using density functional theory at the OLYP/TZ2P level. In the gas phase, all our model SN 2 reactions have single-well PESs, except for the nucleophilic substitution reaction at carbon (SN 2@C), which has a double-well energy profile. The presence of the solvent can have a significant effect on the shape of the PES and, thus, on the nature of the SN 2 process. Solvation energies, charges on the nucleophile or leaving group, and structural features are compared for the various SN 2 reactions in a spectrum of solvents. We demonstrate how solvation can change the shape of the PES, depending not only on the polarity of the solvent, but also on how the charge is distributed over the interacting molecular moieties during different stages of the reaction. In the case of a nucleophilic substitution at three-coordinate phosphorus, the reaction can be made to proceed through a single-well [no transition state (TS)], bimodal barrier (two TSs), and then through a unimodal transition state (one TS) simply by increasing the polarity of the solvent.

The isomeric structure of pentacoordinate chiral spirophosphoranes in solution by the combined use of NMR experiments and GIAO DFT calculations of NMR parameters

The interplay of NMR experiments and DFT calculations of NMR parameters is a reliable method for determining the relative configurations of pentacoordinate chiral spirophosphoranes bearing two six- or five-membered rings at the phosphorus atom in solution. The major product of the Betti based derivatives corresponds to the isomers with both substituents at chiral carbons being opposite to the P-H proton. The next populated product corresponds to the isomer with different chiralities at carbons. The least populated isomer is one with both substituents being at the same side of the heterocycle as the P-H bond.

Stereoselective PCO/POC-Rearrangement of P-C-Cage Phosphorane in the Reaction of 4,5-Dimethyl-2-(2-oxo-1,2-diphenyl)ethoxy-1,3,2-dioxaphospholane with Hexafluoroacetone

Interaction of 4,5-dimethyl-2-(2-oxo-1,2-diphenyl)ethoxy-1,3,2-dioxaphospholane, bearing a carboxyl group in the γ-position with respect to the phosphorus atom and obtained from d,l-butanediol, with hexafluoroacetone (CCl4, -40 °C) leads to the simultaneous formation of regio- and stereoisomeric cage-like phosphoranes with phosphorus-carbon and phosphorus-oxygen bonds with a high stereoselectivity (>95%), whose structure was determined by 1D and 2D NMR spectroscopy and XRD. When stored as a solution in dichloromethane for one month, the PCO-isomer rearranges into the thermodynamically more stable POC-isomer of the cage-like phosphorane. Mild hydrolysis of the PCO/POC-isomers proceeds with a high chemoselectivity and leads to the formation of P(IV)-dioxaphospholane derivatives. Acidic hydrolysis of the POC-isomer leads to the formation of an oxirane derivative with an unexpectedly high stereoselectivity (>95%). DFT calculations (using the PBE functional) allowed us to obtain structures and energies of the initial phospholane, reaction products (PCO/POC-isomers), and an intermediate P(V)-oxaphosphirane.

Inductive and External Electric Field Effects in Pentacoordinated Phosphorus Compounds

Pentacoordination at phosphorus is associated with a nucleophilic displacement reaction at tetracoordinated phosphorus compounds and shows a great variability in what respects their geometrical and energetic features. By means of a systematic theoretical study on a series of elementary model compounds, we have analyzed the bonding features. The pentacoordinated phosphorus compounds are held together by dative bonds, and the geometry and stability depends on the inductive effects originated by different substitutes at phosphorus. We show also that an external electric field can modify the geometrical features and the reactivity of the nucleophilic substitution reactions. This issue may have great interest in biological reactions involving pentacoordinated phosphorus where the electric field originated by the folded protein could influence the catalytic process. We report also additional calculations on the geometry and NMR spectra on three triphenyl phosphonium ylide derivatives, and our results compare well with the experimental data.

A Theoretical Study of Phosphoryl Transfers of Tyrosyl-DNA Phosphodiesterase I (Tdp1) and the Possibility of a "Dead-End" Phosphohistidine Intermediate

Tyrosyl-DNA phosphodiesterase I (Tdp1) is a DNA repair enzyme conserved across eukaryotes that catalyzes the hydrolysis of the phosphodiester bond between the tyrosine residue of topoisomerase I and the 3'-phosphate of DNA. Atomic level details of the mechanism of Tdp1 are proposed and analyzed using a fully quantum mechanical, geometrically constrained model. The structural basis for the computational model is the vanadate-inhibited crystal structure of human Tdp1 (hTdp1, Protein Data Bank entry 1RFF ). Density functional theory computations are used to acquire thermodynamic and kinetic data along the catalytic pathway, including the phosphoryl transfer and subsequent hydrolysis. Located transition states and intermediates along the reaction coordinate suggest an associative phosphoryl transfer mechanism with five-coordinate phosphorane intermediates. Similar to both theoretical and experimental results for phospholipase D, the proposed mechanism for hTdp1 also includes the thermodynamically favorable possibility of a four-coordinate phosphohistidine "dead-end" product.

Phosphoryl transfers of the phospholipase D superfamily: a quantum mechanical theoretical study

The HKD-containing Phospholipase D superfamily catalyzes the cleavage of the headgroup of phosphatidylcholine to produce phosphatidic acid and choline. The mechanism of this cleavage process is studied theoretically. The geometric basis of our models is the X-ray crystal structure of the five-coordinate phosphohistidine intermediate from Streptomyces sp . Strain PMF (PDB Code = 1V0Y ). Hybrid ONIOM QM:QM methodology with Density Functional Theory (DFT) and semiempirical PM6 (DFT:PM6) is used to acquire thermodynamic and kinetic data for the initial phosphoryl transfer, subsequent hydrolysis, and finally, the formation of the experimentally observed ″dead-end″ phosphohistidine product (PDB Code = 1V0W ). The model contains nineteen amino acid residues (including the two highly conserved HKD-motifs), four explicit water molecules, and the substrate. Via computations, the persistence of the short-lived five-coordinate phosphorane intermediate on the minutes times scale is rationalized. This five-coordinate phosphohistidine intermediate energetically exists between the hydrolysis event and ″substrate reorganization″ (the reorganization of the in vitro model substrate within the active site). Computations directly support the thermodynamic favorability of the in vitro four-coordinate phosphohistidine product. In vivo, the activation energy of substrate reorganization is too high, perhaps due to a combination of substrate immobility when embedded in the lipid bilayer, as well as its larger steric bulk compared to the compound used in the in vitro substrate soaks. On this longer time scale, the enzyme will migrate along the lipid membrane toward its next substrate target, rather than promote the formation of the dead-end product.

Synthesis and characterization of antiapicophilic arsoranes and related compounds

Utilizing bulky bidentate ligand systems with C(2)F(5) and n-C(3)F(7) groups, antiapicophilic arsoranes (5b and 5c, respectively) were synthesized. A kinetic study on the isomerization of these arsoranes to their more stable isomers showed that the barriers increased in the order of CF(3) < C(2)F(5) < n-C(3)F(7) in accord with their steric bulk. It was also revealed that the degree of freezing isomerization was larger for the change from CF(3) to C(2)F(5) than from C(2)F(5) to n-C(3)F(7), obvious from the differences in activation free energy at 363 K of 1.6 and 0.3 kcal mol(-1), respectively. X-ray structural analysis of several precursors of these two systems disclosed the unique structures of these compounds.

Theoretical investigation of the enzymatic phosphoryl transfer of β-phosphoglucomutase: revisiting both steps of the catalytic cycle

Enzyme catalyzed phosphate transfer is a part of almost all metabolic processes. Such reactions are of central importance for the energy balance in all organisms and play important roles in cellular control at all levels. Mutases transfer a phosphoryl group while nucleases cleave the phosphodiester linkages between two nucleotides. The subject of our present study is the Lactococcus lactis β-phosphoglucomutase (β-PGM), which effectively catalyzes the interconversion of β-D-glucose-1-phosphate (β-G1P) to β-D-glucose-6-phosphate (β-G6P) and vice versa via stabile intermediate β-D-glucose-1,6-(bis)phosphate (β-G1,6diP) in the presence of Mg(2+). In this paper we revisited the reaction mechanism of the phosphoryl transfer starting from the bisphosphate β-G1,6diP in both directions (toward β-G1P and β-G6P) combining docking techniques and QM/MM theoretical method at the DFT/PBE0 level of theory. In addition we performed NEB (nudged elastic band) and free energy calculations to optimize the path and to identify the transition states and the energies involved in the catalytic cycle. Our calculations reveal that both steps proceed via dissociative pentacoordinated phosphorane, which is not a stabile intermediate but rather a transition state. In addition to the Mg(2+) ion, Ser114 and Lys145 also play important roles in stabilizing the large negative charge on the phosphate through strong coordination with the phosphate oxygens and guiding the phosphate group throughout the catalytic process. The calculated energy barrier of the reaction for the β-G1P to β-G1,6diP step is only slightly higher than for the β-G1,6diP to β-G6P step (16.10 kcal mol(-1) versus 15.10 kcal mol(-1)) and is in excellent agreement with experimental findings (14.65 kcal mol(-1)).

Cyclin-dependent kinases: bridging their structure and function through computations

Cyclin-dependent kinases (CDKs) are one of the most promising target families for drug discovery for several diseases, such as cancer and neurodegenerative disorders. Over the years, structural insights on CDKs have demonstrated high protein plasticity, with several cases where two or more structures of the same protein adopt different conformations. This has generated a great deal of interest in understanding the relationship between CDK structure and function. Here, we highlight how computer simulations have recently contributed in characterizing some key rare and transient events in CDKs, such as the reaction transition state and activation loop movement. Although not yet fully defined, we can now portray the enzymatic mechanism and plasticity of CDKs at high spatial and temporal resolution. These theoretical studies bridge with experiments and highlight structural determinants that could help in designing specific CDK inhibitors.

Fragmentation of pentacoordinate spirobicyclic aminoacyl-phosphoranes (P-AAs) by electrospray ionization tandem mass spectrometry concerning P-O and P-N bond cleavage

The fragmentation pathways of both protonated and sodiated pentacoordinate spirobicyclic aminoacylphosphoranes (P-AAs) have been studied by electrospray ionization multi-stage mass spectrometry (ESI-MS(n)) in positive mode. The possible pathways and their mechanisms are elucidated through the combination of ESI-MS/MS, isotope ((15)N and (2)H) labeling and high-resolution Fourier transform ion cyclotron resonance (FTICR)-MS/MS. The relative Gibbs free energies (ΔG) of the product ions and possible fragmentation pathways are estimated at the B3LYP/6-31 G(d) level of theory. The theoretical calculations show that both protonated and sodiated P-AAs would quickly fragment before Berry pseudorotation. For protonated P-AAs, they have different tendencies to P-O or P-N bond cleavage. For sodiated P-AAs, the P-N bond is easier to cleave and produces the tetracoordinated phosphorus ion H. These results to some extent may give a clue to the chemistry of the active sites of phosphoryl transfer enzymes and will enrich the gas-phase ESI-MS ion chemistry of pentacoordinate phosphoranes.

First prebiotic generation of a ribonucleotide from adenine, D-ribose and trimetaphosphate

Adenosine monophosphate isomers are obtained by self-assembling of adenine, D-ribose and trimetaphosphate in aqueous solution in good yields. This generation of a ribonucleotide from its three molecular components occurs in a one-pot reaction at room temperature for about 30-40 days and with high chemio-, regio-, and stereo-selectivity. Similar results are obtained with guanine. A mechanism is also proposed.

Vibrational circular dichroism spectroscopy of chiral molecules

In this chapter, new developments and main applications of vibrational circular dichroism (VCD) spectroscopy reported in the last 5 years are described. This includes the determinations of absolute configurations of chiral molecules, understanding solvent effects and modeling solvent-solute explicit hydrogen bonding networks using induced solvent chirality, studies of transition metal complexes and their peculiar and enormous intensity enhancements in VCD spectra, investigations of conformational preference of chiral ligands bound to gold nano particles, and two new advances in applying matrix isolation VCD spectroscopy to flexible, multi-conformational chiral molecules and complexes, and in development of femtosecond laser based VCD instruments for transient VCD monitoring. A brief review of the experimental techniques and theoretical methods is also given. The purpose of this chapter is to provide an up-to-date perspective on the capability of VCD to solve significant problems about chiral molecules in solution, in thin film states, or on surfaces.