Arginine Kinase Activates Arginine for Phosphorylation by Pyramidalization and Polarization.

Arginine phosphorylation plays numerous roles throughout biology. Arginine kinase (AK) catalyzes the delivery of an anionic phosphoryl group (PO3-) from ATP to a planar, trigonal nitrogen in a guanidinium cation. Density functional theory (DFT) calculations have yielded a model of the transition state (TS) for the AK-catalyzed reaction. They reveal a network of over 50 hydrogen bonds that delivers unprecedented pyramidalization and out-of-plane polarization of the arginine guanidinium nitrogen (Nη2) and aligns the electron density on Nη2 with the scissile P-O bond, leading to in-line phosphoryl transfer via an associative mechanism. In the reverse reaction, the hydrogen-bonding network enforces the conformational distortion of a bound phosphoarginine substrate to increase the basicity of Nη2. This enables Nη2 protonation, which triggers PO3- migration to generate ATP. This polarization-pyramidalization of nitrogen in the arginine side chain is likely a general phenomenon that is exploited by many classes of enzymes mediating the post-translational modification of arginine.

Octahedral Trifluoromagnesate, an Anomalous Metal Fluoride Species, Stabilizes the Transition State in a Biological Motor.

Isoelectronic metal fluoride transition state analogue (TSA) complexes, MgF3 - and AlF4 -, have proven to be immensely useful in understanding mechanisms of biological motors utilizing phosphoryl transfer. Here we report a previously unobserved octahedral TSA complex, MgF3(H2O)-, in a 1.5 Å resolution Zika virus NS3 helicase crystal structure. 19F NMR provided independent validation and also the direct observation of conformational tightening resulting from ssRNA binding in solution. The TSA stabilizes the two conformations of motif V of the helicase that link ATP hydrolysis with mechanical work. DFT analysis further validated the MgF3(H2O)- species, indicating the significance of this TSA for studies of biological motors.

A GAP-GTPase-GDP-Pi Intermediate Crystal Structure Analyzed by DFT Shows GTP Hydrolysis Involves Serial Proton Transfers.

Cell signaling by small G proteins uses an ON to OFF signal based on conformational changes following the hydrolysis of GTP to GDP and release of dihydrogen phosphate (Pi ). The catalytic mechanism of GTP hydrolysis by RhoA is strongly accelerated by a GAP protein and is now well defined, but timing of inorganic phosphate release and signal change remains unresolved. We have generated a quaternary complex for RhoA-GAP-GDP-Pi . Its 1.75 Šcrystal structure shows geometry for ionic and hydrogen bond coordination of GDP and Pi in an intermediate state. It enables the selection of a QM core for DFT exploration of a 20 H-bonded network. This identifies serial locations of the two mobile protons from the original nucleophilic water molecule, showing how they move in three rational steps to form a stable quaternary complex. It also suggests how two additional proton transfer steps can facilitate Pi release.

"Skinny" and "Fat" DNA: Two New Double Helices.

According to the iconic model, the Watson-Crick double helix exploits nucleobase pairs that are both size complementary (big purines pair with small pyrimidines) and hydrogen bond complementary (hydrogen bond donors pair with hydrogen bond acceptors). Using a synthetic biology strategy, we report here the discovery of two new DNA-like systems that appear to support molecular recognition with the same proficiency as standard Watson-Crick DNA. However, these both violate size complementarity (big pairs with small), retaining hydrogen bond complementarity (donors pair with acceptors) as their only specificity principle. They exclude mismatches as well as standard Watson-Crick DNA excludes mismatches. In crystal structures, these "skinny" and "fat" systems form the expected hydrogen bonds, while conferring novel minor groove properties to the resultant duplex regions of the DNA oligonucleotides. Further, computational tools, previously tested primarily on natural DNA, appear to work well for these two new molecular recognition systems, offering a validation of the power of modern computational biology. These new molecular recognition systems may have application in materials science and synthetic biology, and in developing our understanding of alternative ways that genetic information might be stored and transmitted.

Structure and Biophysics for a Six Letter DNA Alphabet that Includes Imidazo[1,2-a]-1,3,5-triazine-2(8H)-4(3H)-dione (X) and 2,4-Diaminopyrimidine (K).

A goal of synthetic biology is to develop new nucleobases that retain the desirable properties of natural nucleobases at the same time as expanding the genetic alphabet. The nonstandard Watson-Crick pair between imidazo[1,2-a]-1,3,5-triazine-2(8H)-4(3H)-dione (X) and 2,4-diaminopyrimidine (K) does exactly this, pairing via complementary arrangements of hydrogen bonding in these two nucleobases, which do not complement any natural nucleobase. Here, we report the crystal structure of a duplex DNA oligonucleotide in B-form including two consecutive X:K pairs in GATCXK DNA determined as a host-guest complex at 1.75 Å resolution. X:K pairs have significant propeller twist angles, similar to those observed for A:T pairs, and a calculated hydrogen bonding pairing energy that is weaker than that of A:T. Thus, although inclusion of X:K pairs results in a duplex DNA structure that is globally similar to that of an analogous G:C structure, the X:K pairs locally and energetically more closely resemble A:T pairs.

Assessing the Influence of Mutation on GTPase Transition States by Using X-ray Crystallography, 19 F†NMR, and DFT Approaches.

We report X-ray crystallographic and 19 F NMR studies of the G-protein RhoA complexed with MgF3- , GDP, and RhoGAP, which has the mutation Arg85'Ala. When combined with DFT calculations, these data permit the identification of changes in transition state (TS) properties. The X-ray data show how Tyr34 maintains solvent exclusion and the core H-bond network in the active site by relocating to replace the missing Arg85' sidechain. The 19 F NMR data show deshielding effects that indicate the main function of Arg85' is electronic polarization of the transferring phosphoryl group, primarily mediated by H-bonding to O3G and thence to PG . DFT calculations identify electron-density redistribution and pinpoint why the TS for guanosine 5'-triphosphate (GTP) hydrolysis is higher in energy when RhoA is complexed with RhoGAPArg85'Ala relative to wild-type (WT) RhoGAP. This study demonstrates that 19 F NMR measurements, in combination with X-ray crystallography and DFT calculations, can reliably dissect the response of small GTPases to site-specific modifications.

Metal Fluorides: Tools for Structural and Computational Analysis of Phosphoryl Transfer Enzymes.

The phosphoryl group, PO3-, is the dynamic structural unit in the biological chemistry of phosphorus. Its transfer from a donor to an acceptor atom, with oxygen much more prevalent than nitrogen, carbon, or sulfur, is at the core of a great majority of enzyme-catalyzed reactions involving phosphate esters, anhydrides, amidates, and phosphorothioates. The serendipitous discovery that the phosphoryl group could be labeled by "nuclear mutation," by substitution of PO3- by MgF3- or AlF4-, has underpinned the application of metal fluoride (MF x ) complexes to mimic transition states for enzymatic phosphoryl transfer reactions, with sufficient stability for experimental analysis. Protein crystallography in the solid state and 19F NMR in solution have enabled direct observation of ternary and quaternary protein complexes embracing MF x transition state models with precision. These studies have underpinned a radically new mechanistic approach to enzyme catalysis for a huge range of phosphoryl transfer processes, as varied as kinases, phosphatases, phosphomutases, and phosphohydrolases. The results, without exception, have endorsed trigonal bipyramidal geometry (tbp) for concerted, "in-line" stereochemistry of phosphoryl transfer. QM computations have established the validity of tbp MF x complexes as reliable models for true transition states, delivering similar bond lengths, coordination to essential metal ions, and virtually identical hydrogen bond networks. The emergence of protein control of reactant orbital overlap between bond-forming species within enzyme transition states is a new challenging theme for wider exploration.

Consecutive non-natural PZ nucleobase pairs in DNA impact helical structure as seen in 50 µs molecular dynamics simulations.

Z: Little is known about the influence of multiple consecutive 'non-standard' ( , 6-amino-5-nitro-2(1H)-pyridone, and , 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) nucleobase pairs on the structural parameters of duplex DNA. nucleobase pairs follow standard rules for Watson-Crick base pairing but have rearranged hydrogen bonding donor and acceptor groups. Using the X-ray crystal structure as a starting point, we have modeled the motions of a DNA duplex built from a self-complementary oligonucleotide (5΄-CTTATPPPZZZATAAG-3΄) in water over a period of 50 µs and calculated DNA local parameters, step parameters, helix parameters, and major/minor groove widths to examine how the presence of multiple, consecutive nucleobase pairs might impact helical structure. In these simulations, the -containing DNA duplex exhibits a significantly wider major groove and greater average values of stagger, slide, rise, twist and h-rise than observed for a 'control' oligonucleotide in which nucleobase pairs are replaced by . The molecular origins of these structural changes are likely associated with at least two differences between and . First, the electrostatic properties of differ from in terms of density distribution and dipole moment. Second, differences are seen in the base stacking of pairs in dinucleotide steps, arising from energetically favorable stacking of the nitro group in with π-electrons of the adjacent base.

Gas phase RDX decomposition pathways using coupled cluster theory.

Electronic and free energy barriers for a series of gas-phase RDX decomposition mechanisms have been obtain using coupled cluster singles, doubles, and perturbative triples with complete basis set (CCSD(T)/CBS) electronic energies for MBPT(2)/cc-pVTZ structures. Importantly, we have located a well-defined transition state for NN homolysis, in the initial RDX decomposition step, thereby obtaining a true barrier for this reaction. These calculations support the view that HONO elimination is preferred at STP over other proposed mechanisms, including NN homolysis, "triple whammy" and NONO isomerization. Indeed, our calculated values of Arrhenius parameters are in agreement with experimental findings for gas phase RDX decomposition. We also investigate a number of new pathways leading to breakdown of the intermediate formed by the initial HONO elimination, and find that NN homolysis in this intermediate has an activation energy barrier comparable with that computed for HONO elimination.

(19)F†NMR and DFT Analysis Reveal Structural and Electronic Transition State Features for RhoA-Catalyzed GTP Hydrolysis.

Molecular details for RhoA/GAP catalysis of the hydrolysis of GTP to GDP are poorly understood. We use (19)F NMR chemical shifts in the MgF3(-) transition state analogue (TSA) complex as a spectroscopic reporter to indicate electron distribution for the γ-PO3(-) oxygens in the corresponding TS, implying that oxygen coordinated to Mg has the greatest electron density. This was validated by QM calculations giving a picture of the electronic properties of the transition state (TS) for nucleophilic attack of water on the γ-PO3(-) group based on the structure of a RhoA/GAP-GDP-MgF3(-) TSA complex. The TS model displays a network of 20 hydrogen bonds, including the GAP Arg85' side chain, but neither phosphate torsional strain nor general base catalysis is evident. The nucleophilic water occupies a reactive location different from that in multiple ground state complexes, arising from reorientation of the Gln-63 carboxamide by Arg85' to preclude direct hydrogen bonding from water to the target γ-PO3(-) group.

Facile C(sp(2))-C(sp(2)) bond cleavage in oxalic acid-derived radicals.

Oxalate decarboxylase (OxDC) catalyzes the Mn-dependent conversion of the oxalate monoanion into CO2 and formate. Many questions remain about the catalytic mechanism of OxDC although it has been proposed that the reaction proceeds via substrate-based radical intermediates. Using coupled cluster theory combined with implicit solvation models we have examined the effects of radical formation on the structure and reactivity of oxalic acid-derived radicals in aqueous solution. Our results show that the calculated solution-phase free-energy barrier for C-C bond cleavage to form CO2 is decreased from 34.2 kcal/mol for oxalic acid to only 9.3 kcal/mol and a maximum of 3.5 kcal/mol for the cationic and neutral oxalic acid-derived radicals, respectively. These studies also show that the C-C σ bonding orbital of the radical cation contains only a single electron, giving rise to an elongated C-C bond distance of 1.7 Å; a similar lengthening of the C-C bond is not observed for the neutral radical. This study provides new chemical insights into the structure and stability of plausible intermediates in the catalytic mechanism of OxDC, and suggests that removal of an electron to form a radical (with or without the concomitant loss of a proton) may be a general strategy for cleaving the unreactive C-C bonds between adjacent sp(2)-hybridized carbon atoms.

The great diversity of HMX conformers: probing the potential energy surface using CCSD(T).

The octahydro-1,3,5,7-tetranitro-1,3,5,7-tetraazocine (HMX) molecule is a very commonly studied system, in all 3 phases, because of its importance as an explosive; however, no one has ever attempted a systematic study of what all the major gas-phase conformers are. This is critical to a mechanistic study of the kinetics involved, as well as the viability of various crystalline polymorphs based on the gas-phase conformers. We have used existing knowledge of basic cyclooctane chemistry to survey all possible HMX conformers based on its fundamental ring structure. After studying what geometries are possible after second-order many-body perturbation theory (MBPT(2)) geometry optimization, we calculated the energetics using coupled cluster singles, doubles, and perturbative triples (CCSD(T))/cc-pVTZ. These highly accurate energies allow us to better calculate starting points for future mechanistic studies. Additionally, the plethora of structures are compared to existing experimental data of crystals. It is found that the crystal field effect is sometimes large and sometimes small for HMX.

Pragmatic ab initio prediction of enthalpies of formation for large molecules: accuracy of MP2 geometries and frequencies using CCSD(T) correlation energies.

We have addressed the accuracy of calculating the enthalpy of formation of an arbitrary single reference molecule using practical ab initio methodologies. It is known that MP2 geometries with a triple zeta basis set are almost as reliable as CCSD(T) geometries. It is also known that CCSD(T) correlation energies, with basis extrapolation, feature chemical accuracy for single-reference molecules. We investigate what accuracy one might expect in enthalpies of formation from a MP2 geometry, MP2 harmonic vibrational frequencies, a CCSD(T) correlation energy using triple zeta basis sets. It is far from obvious, a priori, as to which error source contributes most significantly. We observe that the accuracy in calculating enthalpies of formation of single-reference molecules with this protocol is 4 kcal mol(-1); our error analysis shows this comes almost exclusively from the correlation energy basis extrapolation, rather than errors intrinsic to MP2.

Conformers of CL-20 explosive and ab initio refinement using perturbation theory: implications to detonation mechanisms.

We have identified the major conformers of CL-20 explosive, otherwise known as 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane, more formally known as 2,4,6,8,10,12-hexanitrohexaazatetracyclo[5.5.0.0]-dodecane, via Monte Carlo search in conformational space through molecular mechanics and subsequent quantum mechanical refinement using perturbation theory. Our search produced enough conformers to account for all of the various forms of CL-20 found in crystals. This suggests that our methodology will be useful in studying the conformational landscape of other nitramines. The energy levels of the conformers found are all within 0.25 eV of one another based on MBPT(2)/6-311G(d,p); consequently, without further refinement from a method such as coupled cluster theory, all conformers may reasonably be populated at STP in the gas phase. We also report the harmonic vibrational frequencies of conformers, including the implications on the mechanism of detonation. In particular, we establish that the weakest N-N nitramine of CL-20 is the cyclohexane equatorial nitramine. This preliminary mapping of the conformers of CL-20 makes it possible to study the mechanism of detonation of this explosive rigorously in future work.

RDX geometries, excited states, and revised energy ordering of conformers via MP2 and CCSD(T) methodologies: insights into decomposition mechanism.

The geometries, harmonic frequencies, elec-tronic excitation levels, and energetic orderings of various conformers of RDX have been computed at the ab initio MP2 and CCSD(T) levels, providing more reliable results than have been previously obtained. We observe that the various local minimum-energy conformers are all competitive for being the absolute minimum and that, at reasonable temperatures, several conformers will appreciably contribute to the population of RDX. As a result, we have concluded that any mechanistic study to investigate thermal decomposition can reasonably begin from any one of the cyclohexane conformers of RDX. As such, it is necessary to consider the transition states for each RDX conformer to gauge what the activation energy is. Homolytic bond dissociation has long been speculated to be critical to detonation; we report here the most accurate estimates of homolytic BDEs yet calculated, likely to be accurate within 3 kcal mol(-1). The differences in energy for homolytic BDEs among all the possible RDR conformers are again small, such that most all of the conformers can reasonably be speculated as the next step in the mechanism starting from the RDR radical.