Inhibitors of orotidine monophosphate decarboxylase (ODCase) have applications in RNA viral, parasitic, and other infectious diseases. ODCase catalyzes the decarboxylation of orotidine monophosphate (OMP), producing uridine monophosphate (UMP). Novel inhibitors 6-amino-UMP and 6-cyano-UMP were designed on the basis of the substructure volumes in the substrate OMP and in an inhibitor of ODCase, barbituric acid monophosphate, BMP. A new enzyme assay method using isothermal titration calorimetry (ITC) was developed to investigate the inhibition kinetics of ODCase. The reaction rates were measured by monitoring the heat generated during the decarboxylation reaction of orotidine monophosphate. Kinetic parameters (k(cat) = 21 s(-1) and KM = 5 microM) and the molar enthalpy (DeltaH(app) = 5 kcal/mol) were determined for the decarboxylation of the substrate by ODCase. Competitive inhibition of the enzyme was observed and the inhibition constants (Ki) were determined to be 12.4 microM and 29 microM for 6-aza-UMP and 6-cyano-UMP, respectively. 6-Amino-UMP was found to be among the potent inhibitors of ODCase, having an inhibition constant of 840 nM. We reveal here the first inhibitors of ODCase designed by the principles of bioisosterism and a novel method of using isothermal calorimetry for enzyme inhibition studies.

Recently, we reported the development of sugar-assisted ligation (SAL), a novel peptide ligation method for the synthesis of glycopeptides. After screening a large number of glycoprotein sequences in a glycoprotein database, it became evident that a large proportion (approximately 53%) of O-glycosylation sites contain amino acid residues that will not undergo SAL reactions. To overcome these inherent limitations and broaden the scope of the method we report here the development of an extended SAL method. Glycopeptides containing up to six amino acid extensions N-terminal to the glycosylated residue were shown to facilitate ligation reactions with peptide thioesters, and these products were isolated in good yields. Kinetic analysis was used to show that as glycopeptides were extended by further amino acid residues, ligation reactions became slower. This finding was rationalized by molecular dynamics simulations using AMBER9. These studies suggested a general trend whereby the proximal distance between the reactive sites of the thioester intermediate (the N-terminal amine and the carbonyl carbon of the thioester) increased as glycopeptides were extended, thus slowing down the ligation rate. Each of the extended SAL methods showed broad tolerance to a number of different amino acid combinations at the ligation junction. Re-evaluation of the glycoprotein database suggested that 95% of the O-linked glycosylation sites can now be utilized to facilitate SAL or extended SAL reactions. As such, this method represents an extremely valuable tool for the synthesis of naturally occurring glycopeptides and glycoproteins. To demonstrate the applicability of the method, extended SAL was successfully implemented in the synthesis of the starting unit of the cancer-associated MUC1 glycoprotein.

Individual peptide groups in proteins must exhibit some variation in the chemical shift anisotropy (CSA) of their constituent atoms, but not much is known about the extent or origins of this dispersion. Direct spectroscopic measurement of CSA remains technically challenging, and theoretical methods can help to overcome these limitations by estimating shielding tensors for arbitrary structures. Here we use an automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach to compute 15N, 13C′ and 1H chemical shift tensors for human ubiquitin and the GB1 and GB3 fragments of staphylococcal protein G. The average and range of variation of the anisotropies is in good agreement with experimental estimates from solid-state NMR, and the variation among residues is somewhat smaller than that estimated from solution-state measurements. Hydrogen-bond effects account for much of the variation, both between helix and sheet regions, and within elements of secondary structure, but other effects (including variations in torsion angles) may play a role as well.

The effects of chemical shift anisotropy (CSA) are evident in line-shapes or side-band analysis in solid-state NMR, in the observed line positions in partially oriented samples, and in relaxation effects in liquid-state studies. In all of these cases, the effective shielding tensor is influenced by fast vibrational averaging in addition to larger-amplitude internal motions and to overall libration or rotation. Here we compute the contributions of vibrational averaging (including zero-point motions) to the CSA relaxation strengths for the nitrogen and carbonyl carbon in two simple peptide models, and for snapshots taken from a path-integral simulation of a small protein. Because the (15)N shielding tensor is determined by all the atoms of the peptide group, it is less influenced by vibrational motion than (for example) the N-H dipolar interaction, which is more sensitive to the motion of the light hydrogen atom. Computed order parameters for CSA averaging are hence much closer to unity than are N-H dipolar order parameters. This leads to a reduction by about 9% in the magnitude of the amide nitrogen CSA that is needed to fit liquid-state relaxation data. Similar considerations apply to the carbonyl carbon shielding tensor, but in this case the differences between dipolar and CSA averaging are smaller. These considerations will be important for making comparisons between CSA tensors extracted from various NMR experiments, and for comparisons to quantum chemical calculations carried out on static conformers.

The human severe acute respiratory syndrome coronavirus (SARS-CoV) and the NL63 coronaviruses are human respiratory pathogens for which no effective antiviral treatment exists. The papain-like cysteine proteases encoded by the coronavirus (SARS-CoV: PLpro; NL63: PLP1 and PLP2) represent potential targets for antiviral drug development. Three recent inhibitor-bound PLpro structures highlight the role of an extremely flexible six-residue loop in inhibitor binding. The high binding site plasticity is a major challenge in computational drug discovery/design efforts. From conventional molecular dynamics and accelerated molecular dynamics (aMD) simulations, we find that with conventional molecular dynamics simulation, PLpro translationally samples the open and closed conformation of BL2 loop on a picosecond–nanosecond timescale but does not reproduce the peptide bond inversion between loop residues Tyr269 and Gln270 that is observed on inhibitor GRL0617 binding. Only aMD simulation, starting from the closed loop conformation, reproduced the 180° ϕ–Ψ dihedral rotation back to the open loop state. The Tyr–Gln peptide bond inversion appears to involve a progressive conformational change of the full loop, starting at one side, and progressing to the other. We used the SARS-CoV apo X-ray structure to develop a model of the NL63-PLP2 catalytic site. Superimposition of the PLP2 model on the PLpro X-ray structure identifies binding site residues in PLP2 that contribute to the distinct substrate cleavage site specificities between the two proteases. The topological and electrostatic differences between the two protease binding sites also help explain the selectivity of non-covalent PLpro inhibitors.