Bridged Epipolythiodiketopiperazines from Penicillium raciborskii, an Endophytic Fungus of Rhododendron tomentosum Harmaja.

Three new epithiodiketopiperazine natural products [outovirin A (1), outovirin B (2), and outovirin C (3)] resembling the antifungal natural product gliovirin have been identified in extracts of Penicillium raciborskii, an endophytic fungus isolated from Rhododendron tomentosum. The compounds are unusual for their class in that they possess sulfide bridges between α- and ß-carbons rather than the typical α-α bridging. To our knowledge, outovirin A represents the first reported naturally produced epimonothiodiketopiperazine, and antifungal outovirin C is the first reported trisulfide gliovirin-like compound. This report describes the identification and structural elucidation of the compounds by LC-MS/MS and NMR.

Relative Binding Enthalpies from Molecular Dynamics Simulations Using a Direct Method.

The potential for reliably predicting relative binding enthalpies, ΔΔE, from a direct method utilizing molecular dynamics is examined for a system of three phosphotyrosyl peptides binding to a protein receptor, the Src SH2 domain. The binding enthalpies were calculated from the potential energy differences between the bound and the unbound end-states of each peptide from equilibrium simulations in explicit water. The statistical uncertainties in the ensemble-mean energy values from multiple, independent simulations were obtained using a bootstrap method. Simulations were initiated with different starting coordinates as well as different velocities. Statistical uncertainties in ΔΔE are 2 to 3 kcal/mol based on calculations from 40, 10 ns trajectories for each system (three SH2-peptide complexes or unbound peptides). Uncertainties in relative component energies, comprising solute-solute, solute-solvent and solvent-solvent interactions, are considerably larger. Energy values were estimated from an unweighted ensemble averaging of multiple trajectories with the a priori assumption that all trajectories are equally likely. Distributions in energy-rmsd space indicate that the trajectories sample the same basin and the difference in mean energy values between trajectories is due to sampling of alternative local regions of this superbasin. The direct estimate of relative binding enthalpies is concluded to be a reasonable approach for well-ordered systems with ΔΔE values greater than ∼3 kcal/mol, although the approach would benefit from future work to determine properly distributed starting points that would enable efficient sampling of conformational space using multiple trajectories.

Twinning and survivorship of captive common marmosets (Callithrix jacchus) and cotton-top tamarins (Saguinus oedipus).

Here we present the results of a demographic analysis of 25 y (1985 to 2010) of common marmoset (Callithrix jacchus) and cotton-top tamarin (Saguinus oedipus) records from the New England Primate Research Center. Summaries of longevity and survivorship are analyzed by birth-type category (including singletons, twins, triplets, and quadruplets) and sex. In addition, a brief evolutionary review is presented. Surrogates of hematopoietic chimerism, twinning, and reproductive output are explored in a large number of animals to help decipher the potential effects of chimerism on life history in marmosets and tamarins. In addition to exploring chimerism through demographic data, multiple-birth frequency and survivorship are compared between species. New World primates can make ideal translational models for disease and behavioral research across multiple disciplines. A better understanding of their reproductive success and survivorship in captivity helps develop these nonhuman primate models, their role in aging research, and understanding of their behavioral ecology. This mission is likely to only increase in its importance to biomedical research due to both the sequencing of the marmoset genome and the growing demand for alternatives to Old World primate models.

Large-scale polymorphism discovery in macaque G-protein coupled receptors.

BACKGROUND: G-protein coupled receptors (GPCRs) play an inordinately large role in human health. Variation in the genes that encode these receptors is associated with numerous disorders across the entire spectrum of disease. GPCRs also represent the single largest class of drug targets and associated pharmacogenetic effects are modulated, in part, by polymorphisms. Recently, non-human primate models have been developed focusing on naturally-occurring, functionally-parallel polymorphisms in candidate genes. This work aims to extend those studies broadly across the roughly 377 non-olfactory GPCRs. Initial efforts include resequencing 44 Indian-origin rhesus macaques (Macaca mulatta), 20 Chinese-origin rhesus macaques, and 32 cynomolgus macaques (M. fascicularis). RESULTS: Using the Agilent target enrichment system, capture baits were designed for GPCRs off the human and rhesus exonic sequence. Using next generation sequencing technologies, nearly 25,000 SNPs were identified in coding sequences including over 14,000 non-synonymous and more than 9,500 synonymous protein-coding SNPs. As expected, regions showing the least evolutionary constraint show greater rates of polymorphism and greater numbers of higher frequency polymorphisms. While the vast majority of these SNPs are singletons, roughly 1,750 non-synonymous and 2,900 synonymous SNPs were found in multiple individuals. CONCLUSIONS: In all three populations, polymorphism and divergence is highly concentrated in N-terminal and C-terminal domains and the third intracellular loop region of GPCRs, regions critical to ligand-binding and signaling. SNP frequencies in macaques follow a similar pattern of divergence from humans and new polymorphisms in primates have been identified that may parallel those seen in humans, helping to establish better non-human primate models of disease.

The resurgence and genetic implications of New World primates in biomedical research.

There has been a recent resurgence of interest in New World monkeys within the biomedical research community, driven by both the sequencing of the common marmoset (Callithrix jacchus) genome and a growing demand for alternatives to Old World primates. New World monkeys offer attractive advantages over Old World species, including cheaper and simpler husbandry, while still maintaining a greater evolutionary proximity to humans compared with other animal models. Although numerous commonalities across primate species exist, there are also important genetic and reproductive differences that can and should play a critical role in selecting appropriate animal models. Common marmosets in particular have significantly reduced diversity at the major histocompatibility complex (MHC) loci and are born as hematopoietic chimeras. New World primates can make ideal translational models for research, but scientists must necessarily incorporate complete understandings of their genetic and phenotypic differences from humans and other model organisms.

Very large residual dipolar couplings from deuterated ubiquitin.

Main-chain (1)H(N)-(15)N residual dipolar couplings (RDCs) ranging from approximately -200 to 200 Hz have been measured for ubiquitin under strong alignment conditions in Pf1 phage. This represents a ten-fold increase in the degree of alignment over the typical weakly aligned samples. The measurements are made possible by extensive proton-dilution of the sample, achieved by deuteration of the protein with partial back-substitution of labile protons from 25 % H(2)O / 75 % D(2)O buffer. The spectral quality is further improved by application of deuterium decoupling. Since standard experiments using fixed-delay INEPT elements cannot accommodate a broad range of couplings, the measurements were conducted using J-resolved and J-modulated versions of the HSQC and TROSY sequences. Due to unusually large variations in dipolar couplings, the trosy (sharp) and anti-trosy (broad) signals are often found to be interchanged in the TROSY spectra. To distinguish between the two, we have relied on their respective (15)N linewidths. This strategy ultimately allowed us to determine the signs of RDCs. The fitting of the measured RDC values to the crystallographic coordinates of ubiquitin yields the quality factor Q = 0.16, which confirms the perturbation-free character of the Pf1 alignment. Our results demonstrate that RDC data can be successfully acquired not only in dilute liquid crystals, but also in more concentrated ones. As a general rule, the increase in liquid crystal concentration improves the stability of alignment media and makes them more tolerant to variations in sample conditions. The technical ability to measure RDCs under moderately strong alignment conditions may open the door for development of alternative alignment media, including new types of media that mimic biologically relevant systems.

Microsecond time-scale conformational exchange in proteins: using long molecular dynamics trajectory to simulate NMR relaxation dispersion data.

With the advent of ultra-long MD simulations it becomes possible to model microsecond time-scale protein dynamics and, in particular, the exchange broadening effects (R(ex)) as probed by NMR relaxation dispersion measurements. This new approach allows one to identify the exchanging species, including the elusive "excited states". It further helps to map out the exchange network, which is potentially far more complex than the commonly assumed 2- or 3-site schemes. Under fast exchange conditions, this method can be useful for separating the populations of exchanging species from their respective chemical shift differences, thus paving the way for structural analyses. In this study, recent millisecond-long MD trajectory of protein BPTI (Shaw et al. Science 2010, 330, 341) is employed to simulate the time variation of amide (15)N chemical shifts. The results are used to predict the exchange broadening of (15)N lines and, more generally, the outcome of the relaxation dispersion measurements using Carr-Purcell-Meiboom-Gill sequence. The simulated R(ex) effect stems from the fast (~10-100 µs) isomerization of the C14-C38 disulfide bond, in agreement with the prior experimental findings (Grey et al. J. Am. Chem. Soc. 2003, 125, 14324).

Constraining binding hot spots: NMR and molecular dynamics simulations provide a structural explanation for enthalpy-entropy compensation in SH2-ligand binding.

NMR spectroscopy and molecular dynamics (MD) simulations were used to probe the structure and dynamics of complexes of three phosphotyrosine-derived peptides with the Src SH2 domain in an effort to uncover a structural explanation for enthalpy-entropy compensation observed in the binding thermodynamics. The series of phosphotyrosine peptide derivatives comprises the natural pYEEI Src SH2 ligand, a constrained mimic, in which the phosphotyrosine (pY) residue is preorganized in the bound conformation for the purpose of gaining an entropic advantage to binding, and a flexible analogue of the constrained mimic. The expected gain in binding entropy of the constrained mimic was realized; however, a balancing loss in binding enthalpy was also observed that could not be rationalized from the crystallographic structures. We examined protein dynamics to evaluate whether the observed enthalpic penalty might be the result of effects arising from altered motions in the complex. (15)N-relaxation studies and positional fluctuations from molecular dynamics indicate that the main-chain dynamics of the protein show little variation among the three complexes. Root mean squared (rms) coordinate deviations vary by less than 1.5 A for all non-hydrogen atoms for the crystal structures and in the ensemble average structures calculated from the simulations. In contrast to this striking similarity in the structures and dynamics, there are a number of large chemical shift differences from residues across the binding interface, but particularly from key Src SH2 residues that interact with pY, the "hot spot" residue, which contributes about one-half of the binding free energy. Rank-order correlations between chemical shifts and ligand binding enthalpy for several pY-binding residues, coupled with available mutagenesis and calorimetric data, suggest that subtle structural perturbations (<1 A) from the conformational constraint of the pY residue sufficiently alter the geometry of enthalpically critical interactions in the binding pocket to cause the loss of binding enthalpy, leading to the observed enthalpy-entropy compensation. We find no evidence to support the premise that enthalpy-entropy compensation is an inherent property and conclude that preorganization of Src SH2 ligand residues involved in binding hot spots may eventuate in suboptimal interactions with the domain. We propose that introducing constraints elsewhere in the ligand could minimize enthalpy-entropy compensation effects. The results illustrate the utility of the NMR chemical shift to highlight small, but energetically significant, perturbations in structure that might otherwise go unnoticed in an apparently rigid protein.

Dynamical transition in proteins and non-Gaussian behavior of low-frequency modes in self-consistent normal mode analysis.

Self-consistent normal mode analysis (SCNMA) is applied to heme c type cytochrome f to study temperature-dependent protein motion. Classical normal mode analysis assumes harmonic behavior and the protein mean-square displacement has a linear dependence on temperature. This is only consistent with low-temperature experimental results. To connect the protein vibrational motions between low and physiological temperatures, we have incorporated a fitted set of anharmonic potentials into SCNMA. In addition, quantum harmonic-oscillator theory has been used to calculate the displacement distribution for individual vibrational modes. We find that the modes involving soft bonds exhibit significant non-Gaussian dynamics at physiological temperature, which suggests that it may be the cause of the non-Gaussian behavior of the protein motions probed by elastic incoherent neutron scattering. The combined theory displays a dynamical transition caused by the softening of few "torsional" modes in the low-frequency regime ( <50 cm(-1) or <6 meV or >0.6 ps). These modes change from Gaussian to a classical distribution upon heating. Our theory provides an alternative way to understand the microscopic origin of the protein dynamical transition.