Molecular Mechanics: The Cross-conjugated Carbonyl Group in Heterocyclic Compounds. 2. Parameterisation (MM2) of the Adjacent C-O and C-N Bonds

The definition, in a previous paper, of a new type of carbon atom (type 44) in the MM2 force field, i.e. when it is cross-conjugated, implies the reparameterisation of the bonds adjacent to the carbonyl bond. By means of a statistical study based on data issued from the Cambridge Structural Data System, we propose here new σ dipole moments and new stretching parameters for the C(44)-O(41) and C(44)-N(40) bonds included in heterocyclic compounds. These stretching parameters, which are π-bond-order dependent, are quite satisfactory (mean of unsigned deviations: 0.023Å) within limited π-bond-order ranges: 0.25 to 0.32 for the C-O bond, 0.37 to 0.53, for the C-N bond. As for the C(44)=O(7) bond the parameters are generally inappropriate for compounds in which both atoms connected to the type 44 carbon are heteroatoms. The maximum deviation from the experimental bond length, for both bond types can reach ±0.03Å. Here also, part of the dispersion of the results could be attributed to the variation of the effective dielectric constant D from one crystal to another.     

Protein structure refinement: Streptomyces griseus serine protease A at 1.8 A resolution.

The crystal structure of the bacterial serine protease from Streptomyces griseus (SGPA) has been refined at 1.8 Å resolution by a restrained parameter least-squares procedure (Konnert, 1976) to a conventional R factor of 0.139 for 12662 statistically significant reflections [I > 3σ(I)]. The number of variable parameters in the final model was 5912 which included positional and individual thermal parameters of the enzyme, and positions, B factors and occupancies of 175 solvent molecules. The algorithm used for this refinement allows for the simultaneous restraint on bond distances and distances related to interbond angles, the coplanarity of atoms in planar groups, the conservation of chirality of asymmetric centres, non-bonded contact distances, conformational torsional angles and individual isotropic temperature factors.The refined structure of SGPA differs from ideal bond lengths by an overall root-mean-square deviation of 0.02 Å; the corresponding value for angle distances is 0.038 Å. Comparison of the phase angles for the shell of data, 8.0 to 2.8 Å, between the multiple isomorphous replacement phases (Brayer et al., 1978a) and the refined phases, indicates an overall difference (r.m.s.) of 56.6 °. The average conformational angle of the peptide bond (ω) is 179.7 ° (root-mean-square deviation ± 2.5 °) for the 180 peptide bonds of SGPA. Of the 175 solvent molecules included during the course of the refinement, 22 with occupancies ranging from 1.00 to 0.38 are located in the active site and the substrate binding region. It was not until these water molecules were included in the refinement process that the active Ser195 adopted its final conformation (χ¹ = ?77 °). The resulting distance from O^γ of Ser195 to N^ε2 of His57 is 3.1 Å, which, when taken with the observed distortion from linearity (50 °), indicates a rather weak interaction.     

The structure of sodium beta-fluoropyruvate: a gem-diol.

Sodium β-fluoropyruvate is shown by X-ray diffraction and infrared spectroscopy to exist as a gem-diol when crystallized from water in the monoclinic space group $\text{P2}{}_1\text{a}$. The unit cell dimensions are a = 8.693(1), b = 11.556(1), c = 5.252(1) Å, β = 94.85(8) °. Molecular dimensions of the gem-diol are presented. The carbon-carbon bond to the carboxyl group is long [1.551(2) Å], indicative of a bond that can break easily. The hydration of the carbonyl group is attributed to the high electronegativity of the fluorine atom. While reports exist in the literature of a crystalline form of β-fluoropyruvic acid that contains a carbonyl, rather than a gem-diol group, ¹³C nmr studies presented here demonstrate that, in aqueous solution, the gem-diol is the predominant (≥95%) form. The significance of these results in the study of the mechanisms of enzymes utilizing pyruvate as a substrate is discussed. If the carbonyl oxygen activity is important for the action of an enzyme, e.g., in the formation of a Schiff base, it is possible that fluoropyruvate may not be a substitute for pyruvate. In other cases the fluoropyruvate can behave as a substrate analog, although its greater bulk may slow down or prevent its entry into the active site.     

Molecular mechanics:the cross-conjugated carbonyl group in heterocyclic compounds. 3. Parameterisation (MM2) of the adjacent C=C bond: evaluation tests

In the MM2 force field, the definition of a new type of carbon (carbonyl atom, when it is cross-conjugated) has led to the reestimation of the mechanical parameters of the adjacent C(O)-O and C(O)-N bonds in fully-conjugated cyclic compouds: !-pyrones, '-lactones, and conjugated "lactams". New parameters, based on the study of 97 bonds, are presented here for the similar adjacent C(O)-C bond in the same compounds. Comparison of calculated bond lengths to experimental X-ray bond lengths shows that, statistically, the results are substantially improved but the dispersion remains large. Full optimisation of the molecules concerned shows that in some cases the errors accumulate on the C(O)-O bond which is more sensitive to errors in the evaluation of its own ? bond order. The origins of the discrepancies are discussed. Using caffeine as a test molecule, the MM2 method with the parameters proposed here, appears less acurate than the ab initio and DFT methods (both with 6-31G**basis) but still better than the semi-empirical methods (AM1-PM3).     

Refinement of human lysozyme at 1.5 Å resolution analysis of non-bonded and hydrogen-bond interactions

The structure of human lysozyme has been crystallographically refined at 1.5 Å resolution by difference map and restrained least-squares procedures to an R factor of 0.187. A comprehensive analysis of the non-bonded and hydrogen-bonded contacts in the lysozyme molecule, which were not restrained, revealed by the refinement has been carried out. The non-bonded CC contacts begin at ~3.45 Å, and the shorter contacts are dominated, as expected, by interactions between trigonal and tetrahedral carbon atoms. The CO contact distances have a “foot” at 3.05 Å. The CN distance plot shows a significant peak at 3.25 Å, which results from close contact between peptide NHs and carbonyl carbons involved in N_iC′_{i ? 2} interactions in α-helices and reverse turns. The distances involving sulphur atoms discriminate SC trigonal interactions at 3.4 to 3.6 Å from SC tetrahedral interactions at 3.7 Å. All these types of non-bonded interactions show minimum distances close to standard van der Waals' separations.Analysis of hydrogen-bond distances has been carried out by using standard geometry to place hydrogen atoms and measuring the XHO distances. On this basis, there are 130 intramolecular hydrogens: 111 NHO bonds, of which 69 are between main-chain atoms, 13 between side-chain atoms and 29 between mainchain and side-chain atoms. If a cluster of four well-defined internal water molecules is included in the protein structure, there is a total of 19 OHO hydrogen bonds. The mean NO, NHO distances and HN?O angles are 2.96 ± 0.17 Å, 2.05 ± 0.18 Å and 18.5 ± 9.6 °, and the mean OO, OHO distances and HÔO angles are 2.83 ± 0.19 Å, 1.98 ± 0.26 Å and 23.8 ± 13.4 °. The distances agree well with standard values, although the hydrogen bonds are consistently more non-linear than in equivalent small molecules. An analysis of the hydrogen-bond angles at the receptor atom indicates that the α-helix, β-sheet and reverse turn have characteristic angular values. A detailed analysis of the regularity of the α-helices and reverse turns shows small but consistent differences between the α-helices in lysozyme and the current standard model, which may now need revision. Of the 21 reverse turns that include a hydrogen bond, the conformations of 19 agree very closely with four of the five standard types. We conclude that the restrained least-squares method of refinement has been validated by these analyses.     

Hydrogen-bonding in enzyme catalysis. Fourier-transform infrared detection of ground-state electronic strain in acyl-chymotrypsins and analysis of the kinetic consequences.

I.r. difference spectra are presented for 3-(indol-3-yl)acryloyl-, cinnamoyl-, 3-(5-methylthien-2-yl)acryloyl-, dehydrocinnamoyl- and dihydrocinnamoyl-chymotrypsins at low pH, where the acyl-enzymes are catalytically inactive. At least two absorption bands are seen in each case in the ester carbonyl stretching region of the spectrum. Cinnamoyl-chymotrypsin substituted at the carbonyl carbon atom with 13C was prepared. A difference spectrum in which 13C-substituted acyl-enzyme was subtracted from [12C]acyl-enzyme shows two bands in the ester carbonyl region and thus confirms the assignment of the features to the single ester carbonyl group. The frequencies of the ester carbonyl bands are interpreted in terms of differential hydrogen-bonding. In each case a lower-frequency relatively narrow band is assigned to a productive potentially reactive binding mode in which the carbonyl oxygen atom is inserted in the oxyanion hole of the enzyme active centre. The higher-frequency band, which is broader, is assigned to a non-productive binding mode in each case, where a water molecule bridges from the carbonyl oxygen atom to His-57; this mode is equivalent to the crystallographically determined structure of 3-(indol-3-yl)acryloyl-chymotrypsin, i.e. the Henderson structure. A difference spectrum of dihydrocinnamoyl-chymotrypsin taken at higher pH shows resolution of a feature centred upon 1731 cm-1, which is assigned to a non-bonded conformer in which the carbonyl oxygen atom is not hydrogen-bonded. Perturbation of the protein spectrum in the presence of acyl groups is interpreted in terms of enhanced structural rigidity. It is reported that the ester carbonyl region of the difference spectrum of cinnamoyl-subtilisin is complicated by overlap of features that arise from protein perturbation. Measurements of carbonyl absorption frequencies in a number of solvents of the methyl esters of the acyl groups used to make acyl-enzymes have permitted determination of the apparent dielectric constants experienced by carbonyl groups in the enzyme active centre as well as a discussion of the effects of polarity. The ester carbonyl bond strengths of the various conformations were estimated by using simple harmonic oscillator theory and an empirical relation between the force constants and bond strengths. The fractional bond breaking induced by hydrogen-bonding was used to calculate rate enhancement factors by using absolute reaction rate theory.(ABSTRACT TRUNCATED AT 400 WORDS)     

Torsion angle dynamics: Reduced variable conformational sampling enhances crystallographic structure refinement

A reduced variable conformational sampling strategy for macromolecules based on molecular dynamics in torsion angle space is evaluated using crystallographic refinement as a prototypical search problem. Bae and Haug's algorithm for constrained dynamics [Bae, D. S., Haug, E. J. A recursive formulation for constrained mechanical system dynamics. Mech. Struct. Mach. 15:359–382, 1987], originally developed for robotics, was used. Their formulation solves the equations of motion exactly for arbitrary holonomic constraints, and hence differs from commonly used approximation algorithms. It uses gradients calculated in Cartesian coordinates, and thus also differs from internal coordinate formulations. Molecular dynamics can be carried out at significantly higher temperatures due to the elimination of the high frequency bond and angle vibrations. The sampling strategy presented here combines high temperature torsion angle dynamics with repeated trajectories using different initial velocities. The best solutions can be identified by the free R value, or the R value if experimental phase information is appropriately included in the refinement. Applications to crystallographic refinement show a significantly increased radius of convergence over conventional techniques. For a test system with diffraction data to 2 Å resolution, slow-cooling protocols fail to converge if the backbone atom root mean square (rms) coordinate deviation from the crystal structure is greater than 1.25 Å, but torsion angle refinement can correct backbone atom rms coordinate deviations up to approximately 1.7 Å. © 1994 Wiley-Liss, Inc.     

RNA loop structure prediction via bond scaling and relaxation

We have developed a method for predicting the structure of small RNA loops that can be used to augment already existing RNA modeling techniques. The method requires no input constraints on loop configuration other than end-to-end distance. Initial loop structures are generated by randomizing the torsion angles, beginning at one end of the polynucleotide chain and correlating each successive angle with the previous. The bond lengths of these structures are then scaled to fit within the known end constraints and the equilibrium bond lengths of the potential energy function are scaled accordingly. Through a series of rescaling and minimization steps the structures are allowed to relax to lower energy configurations with standard bond lengths and reduced van der Waals clashes. This algorithm has been tested on the variable loops of yeast tRNA-Asp and yeast tRNA-Phe, as well as the sarcin-ricin tetraloop and the anticodon loop of yeast tRNA-Phe. The results indicate good correlation between potential energy and the loop structure predictions that are closest to the variable loop crystal structures, but poorer correlation for the more isolated stem loops. The number of stacking interactions has proven to be a good objective measure of the best loop predictions. Selecting on the basis of energy and stacking, we obtain two structures with 0.65 and 0.75 Å all-atom rms deviations (RMSD) from the crystal structure for the tRNA-Asp variable loop. The best structure prediction for the tRNA-Phe variable loop has an all-atom RMSD of 2.2 Å and a backbone RMSD of 1.6 Å, with a single base responsible for most of the deviation. For the sarcin-ricin loop from 28S ribosomal RNA, the predicted structure's all-atom RMSD from the nmr structure is 1.0 Å. We obtain a 1.8 Å RMSD structure for the tRNA-Phe anticodon loop. © 1996 John Wiley & Sons, Inc.     

Packing analysis of carbohydrates and polysaccharides. V. Crystal structures of two polymorphs of pachyman triacetate

The crystal structures of two polymorphic forms of pachyman triacetate, the fully acetylated derivative of a naturally occuring β-(1 → 3)-D -glucan, were determined by a combination of stereochemical and x-ray diffraction analysis. The two polymorphs could be obtained depending on the temperature and the degree of stretching of film specimens of the substance: polymorph I resulted from stretching 25–50% at 125°C and polymorph II resulted from further stretching to 300% at 215°C. Both polymorphs had previously been shown to have sixfold helical chain conformations, but of unequal pitch. Subsequent detailed structure refinement performed with bond lengths, bond angles, conformational angles, and helix-packing parameters as refinement variables, and the simultaneous minimization of packing and conformational energy and the crystallographic R-factor as refinement criteria, resulted in a complete determination of the two crystal structures. Pachyman triacetate I was found to be a right-handed helix packing with antiparallel polarity and space group P2₁2₁2₁ symmetry (unit-cell parameters a = 11.0, b = 19.0, c (fiber repeat) = 22.38 Å). The acetate groups were nearly planar and the O(2) and O(4) acetates were oriented in such a fashion that the carbonyl double-bond nearly eclipsed the corresponding C—H bond of the ring. The O(6) was in the tg position and its acetate was oriented in such a fashion that the bond sequence C(6)—O(6)—C(6C)—C(6M) was nearly trans-planar, with the carbonyl double-bond bisecting the tetrahedral angle formed by C(6) and its two hydrogens. The final R = 0.221. Pachyman triacetate II was similarly found to be a right-handed helix, but packing as a 50:50 mixture of parallel and antiparallel polarities (unit-cell parameters a = 11.49, b = 20.13, c (fiber repeat) = 18.6 Å). The acetate positions in pachyman triacetate II were substantially the same as in pachyman triacetate I. The final R for the 50:50 mixture was 0.234. Probable reasons for the change in packing polarities are discussed, as are the difficulties encountered in the structure refinement of acetate derivatives.     

Computational studies of novel carbonyl-containing diazabicyclic ligands interacting with α4β2 nicotinic acetylcholine receptor (nAChR) reveal alternative binding modes

We have carried out computational studies on interactions of diazabicyclic amide analogs with α4β2 nAChR using homology modeling, docking and pharmacophore elucidation techniques. We have found alternative ligand binding modes in most cases. All these diverse poses exhibit the quintessential hydrogen-bonding interaction between the ligand basic nitrogen and the backbone carbonyl oxygen atom of the highly conserved Trp-149. This hydrogen bond was always found to be shorter than the one contracted by the ligand carbonyl group and a second hydrogen-bond made by the cationic center with Tyr-93 of the principal face of the protein. In most of the poses observed, cation–π interactions involved three aromatic residues located in the principal face of the protein: Trp-149, Tyr-190 and Tyr-197. The latter amino acid residue appears to often donate a hydrogen-bond to the ligand carbonyl oxygen atom. We also describe two rings of alternative receptor-based hydrogen-bond donor features equidistantly separated from the carbonyl oxygen of the highly conserved Trp-149 approximately by 5 and 8 Å, respectively. These findings could be exploited to design diverse and selective novel chemical libraries for the treatment of diseases and conditions where the α4β2 nAChR is disrupted, such as Alzheimer disease, Parkinson’s disease and l-dopa-induced dyskinesia (LID).     

A comparison of staphostatin B with standard mechanism serine protease inhibitors

Staphostatins are the endogenous, highly specific inhibitors of staphopains, the major secreted cysteine proteases from Staphylococcus aureus. We have previously shown that staphostatins A and B are competitive, active site-directed inhibitors that span the active site clefts of their target proteases in the same orientation as substrates. We now report the crystal structure of staphostatin B in complex with wild-type staphopain B at 1.9 A resolution. In the complex structure, the catalytic residues are found in exactly the positions that would be expected for uncomplexed papain-type proteases. There is robust, continuous density for the staphostatin B binding loop and no indication for cleavage of the peptide bond that comes closest to the active site cysteine of staphopain B. The carbonyl carbon atom C of this peptide bond is 4.1 A away from the active site cysteine sulfur Sgamma atom. The carbonyl oxygen atom O of this peptide bond points away from the putative oxyanion hole and lies almost on a line from the Sgamma atom to the C atom. The arrangement is strikingly similar to the "ionmolecule" arrangement for the complex of papain-type enzymes with their substrates but differs significantly from the arrangement conventionally assumed for the Michaelis complex of papain-type enzymes with their substrates and also from the arrangement that is crystallographically observed for complexes of standard mechanism inhibitors and their target serine proteases.     

Geometric characteristics of hydrogen bonds involving sulfur atoms in proteins

Sulfur atoms have been known to participate in hydrogen bonds (H‐bonds) and these sulfur‐containing H‐bonds (SCHBs) are suggested to play important roles in certain biological processes. This study aims to comprehensively characterize all the SCHBs in 500 high‐resolution protein structures (≤1.8 Å). We categorized SCHBs into six types according to donor/acceptor behaviors and used explicit hydrogen approach to distinguish SCHBs from those of nonhydrogen bonding interactions. It is revealed that sulfur atom is a very poor H‐bond acceptor, but a moderately good H‐bond donor. In α‐helix, considerable SCHBs were found between the sulphydryl group of cysteine residue i and the carbonyl oxygen of residue i‐4, and these SCHBs exert effects in stabilizing helices. Although for other SCHBs, they possess no specific secondary structural preference, their geometric characteristics in proteins and in free small compounds are significantly distinct, indicating the protein SCHBs are geometrically distorted. Interestingly, sulfur atom in the disulfide bond tends to form bifurcated H‐bond whereas in cysteine‐cysteine pairs prefer to form dual H‐bond. These special H‐bonds remarkably boost the interaction between H‐bond donor and acceptor. By oxidation/reduction manner, the mutual transformation between the dual H‐bonds and disulfide bonds for cysteine‐cysteine pairs can accurately adjust the structural stability and biological function of proteins in different environments. Furthermore, few loose H‐bonds were observed to form between the sulphydryl groups and aromatic rings, and in these cases the donor H is almost over against the rim rather than the center of the aromatic ring. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

A crystallographic study of the complex of phosphoramidon with thermolysin. A model for the presumed catalytic transition state and for the binding of extended substrates

The structure of the thermolysin inhibitor phosphoramidon (N-(α-l-rhamnopyranosyl-oxyhydroxyphosphinyl)-l-leucyl-l-tryptophan bound to the crystalline enzyme has been determined to a resolution of 2.3 Å by X-ray crystallography. The study shows that the complex of phosphoramidon with thermolysin resembles that of the presumed catalytic transition state inferred from the geometry of binding of dipeptide inhibitors. Also, the study reveals the mode of binding of thermolysin substrates extended on the imino side of the scissile peptide bond.The crystallographic results are consistent with a variety of other studies on the catalytic activity of thermolysin, and suggest a mechanism of action which is analogous to one of the two alternative mechanisms proposed by Lipscomb and co-workers (1968) for carboxypeptidase A. Key features of the proposed mechanism are that the substrate is initially bound to the enzyme with the carbonyl oxygen of the scissile peptide liganded to the zinc; that Glu143 promotes the nucleophilic attack of a buried water molecule on the carbonyl carbon, forming a tetrahedral intermediate; and that His231 acts as a proton donor. The observed binding of phosphoramidon to thermolysin provides further evidence supporting the mechanism in which Glu143 acts as a general base, promoting the attack of water on the carbonyl carbon, rather than the alternative mechanism in which Glu143 attacks the carbonyl carbon directly, forming an anhydride intermediate.     

Active site dynamics of acyl-chymotrypsin

The motions of water molecules, the acyl moiety, the catalytic triad, and the oxyanion binding site of acyl-chymotrypsin were studied by means of a stochastic boundary molecular dynamics simulation. A water molecule that could provide the nucleophilic OH^? for the deacylation stage of the catalysis was found to be trapped between the imidazole ring of His-57 and the carbonyl carbon of the acyl group. It makes a hydrogen bond with the N^ε2 of His-57 and is heldin place through a network of hydrogen-bonded water molecules in theactive site. The water molecule was found as close as 2.8 Å to the carbonyl carbon. This appears to be due to the constraints imposed by nonbonded interaction in the active site. Configurations were found in which one hydrogen of the trapped water shared a bifurcated hydrogen bond with His-57-N^ε2 and Ser-195-0^γ with the water oxygen very close to the carbonyl carbon. The existence of such a water molecule suggests that large movement of the His-57 imidazole ring between positions suitable for providing general-base catalyzed assistance and for providing general-acid catalyzed assistance may notbe required during the reaction. The simulation indicates that the side chains of residues involved in catalysis (i.e., His-57, Ser-195, and Asp-102) are significantly less flexible than other side chains in the protein. The 40% reduction in rms fluctuations is consistent with a comparable reduction calculated from the temperature factors obtained in the X-ray crystal-lographic data of γ-chymotrypsin. The greater rigidity of active site residues seems to result from interconnected hydrogen bonding networks among the residues and between the residues and the solvent water in the active site. © Wiley-Liss, Inc.     

Molecular Mechanics Studies of Molybdenum Disulphide Catalysts Parameterisation of Molybdenum and Sulphur

Abstract

A method for the parameterisation of molybdenum disulphide is presented which reproduces the crystal structure accurately. The method involves calculating parameters such that there is no net force contribution from any individual term of the potential on any atom. Ideal bond lengths and bond angles are taken from the X-ray crystal structure; stretching and bending force constants are calculated from a combination of spectroscopic data and quantum mechanics calculations, whereby the energy function with bond length or bond angle is calculated and fitted with an harmonic potential. For the non-bonded Lennard-Jones parameters, the dispersion coefficient C was calculated by an interpolation of existing published parameters using a multiple regression and then the crystal energy was minimised with respect to the van der Waals radius r₀ using a fixed crystal fragment.

These parameters were tested for various models of the hexagonal and rhombohedral forms of MoS₂. RMS fits between structures minimised with molecular mechanics and experimental models ranged from 0.006 Å to 0.012 Å.     

Benchmarking of copper(II) LFMM parameters for studying amyloid-β peptides

Ligand field molecular mechanics (LFMM) parameters have been benchmarked for copper (II) bound to the amyloid-β_1–16 peptide fragment. Several density functional theory (DFT) optimised small test models, representative of different possible copper coordination modes, have been used to test the accuracy of the LFMM copper bond lengths and angles, resulting in errors typically less than 0.1 Å and 5°. Ligand field molecular dynamics (LFMD) simulations have been carried out on the copper bound amyloid-β_1–16 peptide and snapshots extracted from the subsequent trajectory. Snapshots have been optimised using DFT and the semi-empirical PM7 method resulting in good agreement against the LFMM calculated geometry. Analysis of substructures within snapshots shows that the larger contribution of geometrical difference, as measured by RMSD, lies within the peptide backbone, arising from differences in DFT and AMBER, and the copper coordination sphere is reproduced well by LFMM. PM7 performs excellently against LFMM with an average RMSD of 0.2 Å over 21 tested snapshots. Further analysis of the LFMD trajectory shows that copper bond lengths and angles have only small deviations from average values, with the exception of a carbonyl moiety from the N-terminus, which can act as a weakly bound fifth ligand.     

Biosynthesis of chloramphenicol. Studies on the origin of the dichloroacetyl moiety

Chloramphenicol produced by cultures of Streptomyces species 3022a supplemented with sodium [1,2-13C]acetate was labelled with 13C exclusively in the dichloromethine (2.6 +/- 0.1%) and carbonyl (0.59 +/- 0.05% carbon atoms. Satellite signals from 13C-13C coupling between covalently bonded 13C-enriched carbon atoms were too intense to be attributed to random combination of labelled atoms at the average enrichments measured, but their intensity relative to those of the signals for uncoupled 13C atoms indicated that most of the precursor had been incorporated after 13C-13C bond fission. Since [2,3-13c]succinic acid enriched only the carbonyl carbon atom of chloramphenicol, these results suggest that neither acetate nor a Krebs cycle intermediate is a direct precursor of the dichloroacetyl group. Cultures supplemented with [2-3h]-or [2h2]-dichloroacetic acid incorporated negligible amounts of isotope into the antibiotic; on this evidence, the free acid is not an intermediate in chloramphenicol biosynthesis and the acylation step may precede chlorination.     

Alpha-ketoamides and alpha-ketocarbonyls: conformational analysis and development of all-atom OPLS force field

The molecular structures and barriers for the internal rotation around the OC-CO single bond in four alpha-ketoamides and eight alpha-ketocarbonyls have been determined from the MP3/aug-cc-pVDZ and MP2/aug-cc-pVDZ calculations. Alpha-ketocarbonyls with non-bulky substituents adopt planar conformations with two carbonyl oxygens in s-trans arrangement. The s-cis conformation is significantly less stable due to the electrostatic repulsion between the two carbonyl groups. Primary and secondary alpha-ketoamides are planar when the substituent at the carbonyl carbon is hydrogen or methyl group but tertiary alpha-ketoamides adopt a conformation where the OC-CO unit is significantly bent. Based on current ab initio structural data, a set of OPLS-AA force field parameters has been derived. These parameters can be used for the modeling of a variety of alpha-ketoamide or alpha-ketocarbonyl containing drugs such as novel protease inhibitors or neuroregenerative polyketides.     

Crystal structures of an intein from the split dnaE gene of Synechocystis sp. PCC6803 reveal the catalytic model without the penultimate histidine and the mechanism of zinc ion inhibition of protein splicing

The first naturally occurring split intein was found in the dnaE gene of Synechocystis sp. PCC6803 and belongs to a subclass of inteins without a penultimate histidine residue. We describe two high-resolution crystal structures, one derived from an excised Ssp DnaE intein and the second from a splicing-deficient precursor protein. The X-ray structures indicate that His147 in the conserved block F activates the side-chain N(delta) atom of the intein C-terminal Asn159, leading to a nucleophilic attack on the peptide bond carbonyl carbon atom at the C-terminal splice site. In this process, Arg73 appears to stabilize the transition state by interacting with the carbonyl oxygen atom of the scissile bond. Arg73 also seems to substitute for the conserved penultimate histidine residue in the formation of an oxyanion hole, as previously identified in other inteins. The finding that the precursor structure contains a zinc ion chelating the highly conserved Cys160 and Asp140 reveals the structural basis of Zn2+-mediated inhibition of protein splicing. Furthermore, it is of interest to observe that the carbonyl carbon atom of Asn159 and N(eta) of Arg73 are 2.6 angstroms apart in the free intein structure and 10.6 angstroms apart in the precursor structure. The orientation change of the aromatic ring of Tyr-1 following the initial acyl shift may be a key switching event contributing to the alignment of Arg73 and the C-terminal scissile bond, and may explain the sequential reaction property of the Ssp DnaE intein.     

Monomeric molybdenum(V) complexes. 4. The structure of tetraphenylarsonium oxodichloro(N-2-oxophenylsalicylideniminato)molybdate(V), [Ph4As] [MoOCl2(SalphO)]

The structure of [Ph₄As] [MoOCl₂(SalphO)], where SalphO is N-2-oxophenylsalicylideniminate dianion, has been determined by X-ray crystallography. The complex crystallizes in the monoclinic space group P2₁/n with a = 11.829(2), b = 16.149(3), c = 17.410(3) Å, β = 97.485(15)° and Z = 4. The calculated and observed densities and 1.566 and 1.573(10) g cm^?3, respectively. Block-diagonal least-squares refinement of the structure using 4722 independent reflections with I ? 3σ(I) converged at R = 0.0345 and R_w = 0.0484. The crystal contains [Ph₄As]⁺ cations and [MoOCl₂(SalphO)]^? anions. The Mo atom in the anion is in a distorted octahedral coordination environment. A planar terdentate Schiff base ligand occupies meridional positions with the N atom trans to the terminal oxo group (O_t). Two Cl atoms are cis to the O_t atom. The Mo atom is displaced by 0.33 Å from the equatorial plane toward the O_t atom. The MoO_t distance is 1.673(3) Å. The MoN bond trans to the O_t atom is 2.298(4) Å. The two MoCl bond lengths are 2.371(1) and 2.408(1) Å. The difference of 0.037 Å is significant (30 σ). Preparations of the title complex and the related complexes are also described.