Synthesis, Characterization, and Computational Modeling of N-(1-Ethoxyvinyl)pyridinium Triflates, an Unusual Class of Pyridinium Salts.

N-Substituted pyridinium salts constitute one of the most valuable reagent classes in organic synthesis, due to their versatility and ease of use. Herein we report a preliminary synthesis and detailed structural analysis of several N-(1-ethoxyvinyl)pyridinium triflates, an unusual class of pyridinium salts with potentially broad use as a reagent in organic synthesis. Treatment of pyridines with trifluoromethane sulfonic acid and ethoxyacetylene generates stable, isolable adducts which have been extensively characterized, due to their novelty. Three-dimensional structural stability is perpetuated by an array of C-H•••O hydrogen bonds involving oxygen atoms from the -SO3 groups of the triflate anion, and hydrogen atoms from the aromatic ring and vinyl group of the pyridinium cation. Predictions from density functional theory calculations of the energy landscape for rotation about the exocyclic C-N bond of 2-chloro-1-(1-ethoxyvinyl)pyridine-1-ium trifluoromethanesulfonate (7) and 1-(1-ethoxyvinyl)pyridine-1-ium trifluoromethanesulfonate (16) are also reported. Notably, the predicted global energy minimum of 7 was nearly identical to that found within the crystal structure.

Solution X-ray scattering combined with computational modeling reveals multiple conformations of covalently bound ubiquitin on PCNA.

PCNA ubiquitination in response to DNA damage leads to the recruitment of specialized translesion polymerases to the damage locus. This constitutes one of the initial steps in translesion synthesis (TLS)--a critical pathway for cell survival and for maintenance of genome stability. The recent crystal structure of ubiquitinated PCNA (Ub-PCNA) sheds light on the mode of association between the two proteins but also revealed that paradoxically, the ubiquitin surface engaged in PCNA interactions was the same as the surface implicated in translesion polymerase binding. This finding implied a degree of flexibility inherent in the Ub-PCNA complex that would allow it to transition into a conformation competent to bind the TLS polymerase. To address the issue of segmental flexibility, we combined multiscale computational modeling and small angle X-ray scattering. This combined strategy revealed alternative positions for ubiquitin to reside on the surface of the PCNA homotrimer, distinct from the position identified in the crystal structure. Two mutations originally identified in genetic screens and known to interfere with TLS are positioned directly beneath the bound ubiquitin in the alternative models. These computationally derived positions, in an ensemble with the crystallographic and flexible positions, provided the best fit to the solution scattering, indicating that ubiquitin dynamically associated with PCNA and is capable of transitioning between a few discrete sites on the PCNA surface. The finding of new docking sites and the positional equilibrium of PCNA-Ub occurring in solution provide unexpected insight into previously unexplained biological observations.

Conservation and variation of structural flexibility in protein families.

In this issue, Raimondi et al. (2010) obtained interesting insights concerning structural flexibilities in the Ras superfamily that are essential to both function retention and specialization by analyzing the deformation patterns from physical models of protein structure and from crystal structures of homologous proteins.

Role of secondary sialic acid binding sites in influenza N1 neuraminidase.

Within influenza viral particles, the intricate balance between host cell binding and sialic acid receptor destruction is carefully maintained by the hemagglutinin (HA) and neuraminidase (NA) glycoproteins, respectively. A major outstanding question in influenza biology is the function of a secondary sialic acid binding site on the NA enzyme. Through a series of Brownian dynamics (BD) simulations of the avian N1, human pandemic N2, and currently circulating pandemic (H1)N1 enzymes, we have probed the role of this secondary sialic acid binding site in the avian N1 subtype. Our results suggest that electrostatic interactions at the secondary and primary sites in avian NA may play a key role in the recognition process of the sialic acid receptors and catalytic efficiency of NA. This secondary site appears to facilitate the formation of complexes with the NA protein and the sialic acid receptors, as well as provide HA activity to a lesser extent. Moreover, this site is able to steer inhibitor binding as well, albeit with reduced capacity in N1, and may have potential implications for drug resistance or optimal inhibitor design. Although the secondary sialic acid binding site has previously been shown to be nonconserved in swine NA strains, our investigations of the currently circulating pandemic H1N1 strain of swine origin appears to have retained some of the key features of the secondary sialic acid binding site. Our results indicate possible lowered HA activity for this secondary sialic acid site, which may be an important event in the emergence of the current pandemic strain.

Protein structural variation in computational models and crystallographic data.

Normal mode analysis offers an efficient way of modeling the conformational flexibility of protein structures. We use anisotropic displacement parameters from crystallography to test the quality of prediction of both the magnitude and directionality of conformational flexibility. Normal modes from four simple elastic network model potentials and from the CHARMM force field are calculated for a data set of 83 diverse, ultrahigh-resolution crystal structures. While all five potentials provide good predictions of the magnitude of flexibility, all-atom potentials have a clear edge at prediction of directionality, and the CHARMM potential has the highest prediction quality. The low-frequency modes from different potentials are similar, but those computed from the CHARMM potential show the greatest difference from the elastic network models. The comprehensive evaluation demonstrates the costs and benefits of using normal mode potentials of varying complexity.

Interpreting correlated motions using normal mode analysis.

With the increased popularity of normal mode analyses in structural biology, it is important to carefully consider how to best utilize the results for gaining biological insights without over interpretation. The discussion in this article argues that for the purpose of identifying correlated motions in biomolecules, a case separate from concomitant conformational changes of structural motifs, it is generally important to use a large number of normal modes. This is illustrated through three increasingly complex examples. The simplest case includes two bilinearly coupled harmonic oscillators and serves as a straightforward problem where the important considerations are explicit and transparent. The argument is then generalized to include a system of N-coupled harmonic oscillators and finally to a realistic biomolecule. Although a small number of normal modes are useful for probing structural flexibility, it is clear that a much larger number of modes are required for properly investigating correlated motions in biomolecules.

Comparison of mode analyses at different resolutions applied to nucleic acid systems.

More than two decades of different types of mode analyses has shown that these techniques can be useful in describing large-scale motions in protein systems. A number of mode analyses are available and include quasiharmonics, classical normal mode, block normal mode, and the elastic network model. Each of these methods has been validated for protein systems and this variety allows researchers to choose the technique that gives the best compromise between computational cost and the level of detail in the calculation. These same techniques have not been systematically tested for nucleic acid systems, however. Given the differences in interactions and structural features between nucleic acid and protein systems, the validity of these techniques in the protein regime cannot be directly translated into validity in the nucleic acid realm. In this work, we investigate the usefulness of the above mode analyses as applied to two RNA systems, i.e., the hammerhead ribozyme and a guanine riboswitch. We show that classical normal-mode analysis can match the magnitude and direction of residue fluctuations from the more detailed, anharmonic technique, quasiharmonic analysis of a molecular dynamics trajectory. The block normal-mode approximation is shown to hold in the nucleic acid systems studied. Only the mode analysis at the lowest level of detail, the elastic network model, produced mixed results in our calculations. We present data that suggest that the elastic network model, with the popular parameterization, is not best suited for systems that do not have a close packed structure; this observation also hints at why the elastic network model has been found to be valid for many globular protein systems. The different behaviors of block normal-mode analysis and the elastic network model, which invoke similar degrees of coarse-graining to the dynamics but use different potentials, suggest the importance of applying a heterogeneous potential function in a robust analysis of the dynamics of biomolecules, especially those that are not closely packed. In addition to these comparisons, we briefly discuss insights into the conformational space available to the hammerhead ribozyme.

The exclusion of glycine betaine from anionic biopolymer surface: why glycine betaine is an effective osmoprotectant but also a compatible solute.

Paradoxically, glycine betaine (N,N,N-trimethyl glycine; GB) in vivo is both an effective osmoprotectant (efficient at increasing cytoplasmic osmolality and growth rate) and a compatible solute (without deleterious effects on biopolymer function, including stability and activity). For GB to be an effective osmoprotectant but not greatly affect biopolymer stability, we predict that it must interact very differently with folded protein surface than with that exposed in unfolding. To test this hypothesis, we quantify the preferential interaction of GB with the relatively uncharged surface exposed in unfolding the marginally stable lacI helix-turn-helix (HTH) DNA binding domain using circular dichroism and with the more highly charged surfaces of folded hen egg white lysozyme (HEWL) and bovine serum albumin (BSA) using all-gravimetric vapor pressure osmometry (VPO) and compare these results with results of VPO studies (Hong et al. (2004), Biochemistry, 43, 14744-14758) of the interaction of GB with polyanionic duplex DNA. For these four biopolymer surfaces, we observe that the extent of exclusion of GB per unit of biopolymer surface area increases strongly with increasing fraction of anionic oxygen (protein carboxylate or DNA phosphate) surface. In addition, GB is somewhat more excluded from the surface exposed in unfolding the lacI HTH and from the folded surface of HEWL than expected from their small fraction of anionic surface, consistent with moderate exclusion of GB from polar amide surface, as predicted by the osmophobic model of protein stability (Bolen and Baskakov (2001) J. Mol. Biol. 310, 955-963). Strong exclusion of GB from anionic surface explains how it can be both an effective osmoprotectant and a compatible solute; analysis of this exclusion yields a lower bound on the hydration of anionic protein carboxylate surface of two layers of water (>or=0.22 H(2)O A(-)(2)).