Search

Stimulus-responsive hydrogel materials that stabilize and control protein dynamics have the potential to enable a range of applications that take advantage of the inherent specificity and catalytic efficiencies of proteins. Here we describe the modular construction of a hydrogel using an engineered calmodulin (CaM) within a poly(ethylene glycol) (PEG) matrix that involves the reversible tethering of proteins through an engineered CaM-binding sequence. For these measurements, maltose binding protein (MBP) was isotopically labeled with (13)C and (15)N, permitting dynamic structural measurements using TROSY-HSQC NMR spectroscopy. The protein dynamics is suppressed upon initial formation of hydrogels, with a concomitant increase in protein stability. Relaxation of the hydrogel matrix following transient heating results in enhanced protein dynamics and resolution of substrate-induced large-amplitude domain rearrangements.

The underlying mechanisms that lead to dramatic differences between closely related pathogens are not always readily apparent. For example, the genomes of Yersinia pestis (YP) the causative agent of plague with a high mortality rate and Yersinia pseudotuberculosis (YPT) an enteric pathogen with a modest mortality rate are highly similar with some species specific differences; however the molecular causes of their distinct clinical outcomes remain poorly understood. In this study, a temporal multi-omic analysis of YP and YPT at physiologically relevant temperatures was performed to gain insights into how an acute and highly lethal bacterial pathogen, YP, differs from its less virulent progenitor, YPT. This analysis revealed higher gene and protein expression levels of conserved major virulence factors in YP relative to YPT, including the Yop virulon and the pH6 antigen. This suggests that adaptation in the regulatory architecture, in addition to the presence of unique genetic material, may contribute to the increased pathogenenicity of YP relative to YPT. Additionally, global transcriptome and proteome responses of YP and YPT revealed conserved post-transcriptional control of metabolism and the translational machinery including the modulation of glutamate levels in Yersiniae. Finally, the omics data was coupled with a computational network analysis, allowing an efficient prediction of novel Yersinia virulence factors based on gene and protein expression patterns.

Formamidopyrimidine-DNA glycosylase (Fpg) is a base excision repair (BER) protein that removes oxidative DNA lesions. Recent crystal structures of Fpg bound to DNA revealed residues involved in damage recognition and enzyme catalysis, but failed to shed light on the dynamic nature of the processes. To examine the structural and dynamic changes that occur in solution when Fpg binds DNA, NMR spectroscopy was used to study Escherichia coli Fpg free in solution and bound to a double-stranded DNA oligomer containing 1,3-propanediol (13-PD), a non-hydrolyzable abasic-site analogue. Only 209 out of a possible 251 (83%) free-precession 15N/1H HSQC cross peaks were observed and 180 of these were assignable, indicating that approximately 30% of the residues undergo intermediate motion on the NMR timescale, broadening the resonances beyond detection or making them intractable in backbone assignment experiments. The majority of these affected residues were in the polypeptide linker region and the interface between the N- and C-terminal domains. DNA titration experiments revealed line broadening and chemical shift perturbations for backbone amides nearby and distant from the DNA binding surface, but failed to quench the intermediate timescale motion observed for free Fpg, including those residues directly involved in DNA binding, notwithstanding a nanomolar dissociation constant for 13-PD binding. Indeed, after binding to 13-PD, at least approximately 40% of the Fpg residues undergo intermediate timescale motion even though all other residues exhibit tight DNA binding characteristic of slow exchange. CPMG-HSQC experiments revealed millisecond to microsecond motion for the backbone amides of D91 and H92 that were quenched upon binding 13-PD. In free Fpg, heteronuclear 1H-15N NOE experiments detected picosecond timescale backbone motion in the alphaF-beta9 loop, the region primarily responsible for chemically discriminating 8-oxoguanine (8-oxoG) over normal guanine, that was quenched after binding 13-PD. Collectively, these observations reveal that, in solution, Fpg is a very dynamic molecule even after binding damaged DNA. Such motion, especially at the DNA binding surface, may be key to its processive search for DNA damage and its catalytic functions once it recognizes damaged DNA.

CylR2 is one of two regulatory proteins associated with the quorum-sensing-dependent synthesis of cytolysin in the common pathogen Enterococcus faecalis. The protein was expressed with a C-terminal six-histidine tag and purified to homogeneity with a cobalt-affinity column followed by size-exclusion chromatography. Both native and SeMet proteins were produced and crystallized. Complete X-ray diffraction data sets were collected from a native crystal, which diffracted to 2.3 angstroms resolution, and a SeMet crystal, which diffracted to 2.1 angstroms. The crystals were tetragonal, belonging to space group P4(1) or P4(3), with unit-cell parameters a = b = 66.2, c = 40.9 angstroms, alpha = beta = gamma = 90 degrees. Based on the calculated Matthews coefficient of 2.6 angstroms3 Da(-1) as well as analysis of anomalous difference Patterson maps, the asymmetric unit most likely contains two molecules of CylR2.

Restrained molecular dynamics is widely used to calculate DNA structures from NMR data. Here, results of an in silico experiment show that the force field can be significant compared to the NMR restraints in driving the final structures to converge. Specifically, we observed that i) the influence of the force field leads to artificially tight convergence within final families of structures and ii) the precision and character of resulting structures depend on the choice of force field used in the calculations. A canonical B-DNA model was used as a target structure. Distances, dihedral angles, and simulated residual dipolar couplings were measured in the target structure and used as restraints. X-PLOR and Discover, which use force fields developed for CHARMM and AMBER programs, respectively, were tested and found to produce different final structures despite the use of identical distance and dihedral restraints. Incorporation of residual dipolar coupling restraints in X-PLOR improves convergence with the target structure and between families of structures indicating that the force field dependence can potentially be overcome if residual dipolar coupling restraints are employed.

Recent NMR-based, chemical shift mapping experiments with the minimal DNA-binding domain of XPA (XPA-MBD: M98–F219) suggest that a basic cleft located in the loop-rich subdomain plays a role in DNA-binding. Here, XPA–DNA interactions are further characterized by NMR spectroscopy from the vantage point of the DNA using a single-stranded DNA nonamer, dCCAATAACC (d9). Up to 2.5 molar equivalents of XPA-MBD was titrated into a solution of d9. A subset of 31P resonances of d9 were observed to broaden and/or shift providing direct evidence that XPA-MBD binds d9 by a mechanism that perturbs the phosphodiester backbone of d9. The interior five residues of d9 broadened and/or shifted before 31P resonances of phosphate groups at the termini, suggesting that when d9 is bound to XPA-MBD the internal residues assume a correlation time that is characteristic of the molecular weight of the complex while the residues at the termini undergo a fraying motion away from the surface of the protein on a timescale such that the line widths are more characteristic of the molecular weight of ssDNA. A molecular model of the XPA-MBD complex with d9 was calculated based on the 15N (XPA-MBD) and 31P (d9) chemical shift mapping studies and on the assumption that electrostatic interactions drive the complex formation. The model shows that a nine residue DNA oligomer fully covers the DNA-binding surface of XPA and that there may be an energetic advantage to binding DNA in the 3′→5′ direction rather than in the 5′→3′ direction (relative to XPA-MBD α-helix-3).

Stereo-selectivedeuteration has been explored as an approach for improving the accuracy of NMR-derived, three-bond vicinal proton-proton coupling constants in the 12-base-pair DNA Dickerson sequence [d(CGCGAATTCGCG)(2)]. The coupling constants are useful for DNA structure determination in restrained molecular dynamics calculations. Specifically, the A5 and A6 residues were prepared with the H2" proton stereo-selectively replaced with a deuteron. Deuteration of the H2" leads to a 42-fold reduction in the transverse cross-relaxation rate of the H2' spin, effectively negating the contribution of transverse cross relaxation to the cross peak frequencies and phases. Calculated linewidth and polarization transfer functions indicated that the reduced dipolar interaction is also expected to result in a significant increase in intensity for all cross peaks involving the H1', H2', or H3' spin. The spectral complexity is also reduced by selective deuteration. Time-shared homonuclear decoupling of passive spins during acquisition was implemented, reducing the spin system, in some cases, to an effectively isolated two-spin system. This enables the use of a 90 degrees mixing pulse instead of the 35 degrees pulse commonly used in standard P.E.COSY experiments, leading to an additional 75% increase in signal intensity. Selective excitation pulses were used to reduce the number of increments required in the indirect dimension by as much as a factor of 4. The cumulative improvement in sensitivity is striking, approaching three orders of magnitude per unit time. Separate experiments, referred to as Stripe-COSY and Superstripe-COSY, were optimized for each coupling constant measured. Finally, J-doubling was used to obtain the most accurate peak separations. This comprehensive approach shows promise as an effective method for extracting highly accurate homonuclear vicinal coupling constants in DNA.

The solution structures of a duplex DNA dodecamer containing a cis-syn cyclobutane thymine dimer d(GCACGAAT[ cs ]TAAG)·d(CTTAATTCG TGC) and its native parent sequence were determined using NMR data collected at 750 MHz. The dodecamer sequence corresponds to the section of a site-specific cis-syn dimer containing 49-mer that was found to be the binding site for the dimer-specific T4 denV endonuclease V repair enzyme by chemical and enzymatic footprinting experiments. Structures of both sequences were derived from NOE restrained molecular dynamics/simulated annealing calculations using a fixed outer layer of water and an inner dynamic layer of water with sodium counterions. The resulting structures reveal a subtle distortion to the phosphodiester backbone in the dimer-containing sequence which includes a B II phosphate at the T9pA10 junction immediately 3′ to the dimer. The B II phosphate, established experimentally by analysis of the 31 P chemical shifts and interpretation of the 3 J P-H3′ values using an optimized Karplus relationship, enables the DNA helix to accommodate the dimer by destacking the base 3′ to the dimer. Furthermore, the structures provide explanations for the unusually shifted T8-N3H imino, A16-H2 and T8-Me proton resonances and T9pA10 31 P NMR resonance and are consistent with bending, unwinding, and thermodynamic data. The implications of the structural data for the mechanism by which cis-syn dimers are recognized by repair enzymes and bypassed by DNA polymerases are also discussed.

113Cd NMR spectroscopy in both the solution and solid state has been used to investigate the role of the metal ion and the proximal histidine on metalloporphyrin reorientation in myoglobin. Heme disorder has been known to exist for many years but understanding its mechanism has proved difficult due to the short-lived nature of the minor porphyrin isomer in native myoglobin. Cadmium-substituted myoglobin can be generated in one form which contains different insertion isomers or in another form which contains predominantly only one of these species. This allows for direct investigation of heme disorder in metal-substituted myoglobin.

Base dynamics, heretofore observed only at TpA steps in DNA, were investigated as a function of sequence context by NMR spectroscopy. The large amplitude conformational dynamics have been previously observed in TnAn segments where n > or = 2. In order to determine whether the dynamic characteristics occur in more general sequence contexts, we examined four self-complementary DNA sequences, [d(CTTTA-NATNTAAAG)2] (where N = A, C, T, G and N = complement of N). The anomalous broadening of the TpA adenine H2 resonance which is indicative of large amplitude base motion was observed in all nine unique four nucleotide contexts. Furthermore, all the adenine H2 resonances experienced a linewidth maximum as a function of temperature, which is a characteristic of the dynamic process. Interestingly, the temperature of the linewidth maximum varied with sequence indicating that the thermodynamics of TpA base dynamics are also sequence dependent. In one example, neither a T preceding nor an A trailing the TpA step was required for base dynamics. These results show that base dynamics, heretofore observed in only a few isolated sequences, occurs at all TpA steps which are either preceded or followed by a thymine or adenine, respectively, and may be characteristic of all TpA steps in DNA notwithstanding sequence context.

Salmonella and Yersinia are two distantly related genera containing species with wide host-range specificity and pathogenic capacity. The metabolic complexity of these organisms facilitates robust lifestyles both outside of and within animal hosts. Using a pathogen-centric systems biology approach, we are combining a multi-omics (transcriptomics, proteomics, metabolomics) strategy to define properties of these pathogens under a variety of conditions including those that mimic the environments encountered during pathogenesis. These high-dimensional omics datasets are being integrated in selected ways to improve genome annotations, discover novel virulence-related factors, and model growth under infectious states. We will review the evolving technological approaches toward understanding complex microbial life through multi-omic measurements and integration, while highlighting some of our most recent successes in this area.

Background Yersinia pestis is a gram-negative bacterium that causes plague, a disease linked historically to the Black Death in Europe during the Middle Ages and to several outbreaks during the modern era. Metabolism in Y. pestis displays remarkable flexibility and robustness, allowing the bacterium to proliferate in both warm-blooded mammalian hosts and cold-blooded insect vectors such as fleas. Results Here we report a genome-scale reconstruction and mathematical model of metabolism for Y. pestis CO92 and supporting experimental growth and metabolite measurements. The model contains 815 genes, 678 proteins, 963 unique metabolites and 1678 reactions, accurately simulates growth on a range of carbon sources both qualitatively and quantitatively, and identifies gaps in several key biosynthetic pathways and suggests how those gaps might be filled. Furthermore, our model presents hypotheses to explain certain known nutritional requirements characteristic of this strain. Conclusions Y. pestis continues to be a dangerous threat to human health during modern times. The Y. pestis genome-scale metabolic reconstruction presented here, which has been benchmarked against experimental data and correctly reproduces known phenotypes, provides an in silico platform with which to investigate the metabolism of this important human pathogen.

High-throughput (HTP) technologies offer the capability to evaluate the genome, proteome, and metabolome of an organism at a global scale. This opens up new opportunities to define complex signatures of disease that involve signals from multiple types of biomolecules. However, integrating these data types is difficult due to the heterogeneity of the data. We present a Bayesian approach to integration that uses posterior probabilities to assign class memberships to samples using individual and multiple data sources; these probabilities are based on lower-level likelihood functions derived from standard statistical learning algorithms. We demonstrate this approach on microbial infections of mice, where the bronchial alveolar lavage fluid was analyzed by three HTP technologies, two proteomic and one metabolomic. We demonstrate that integration of the three datasets improves classification accuracy to approximately 89% from the best individual dataset at approximately 83%. In addition, we present a new visualization tool called Visual Integration for Bayesian Evaluation (VIBE) that allows the user to observe classification accuracies at the class level and evaluate classification accuracies on any subset of available data types based on the posterior probability models defined for the individual and integrated data.

DNA curvature but no single model has been able to explain all the experimental data. An intriguing observation is that AnTn segments ligated in phase exhibit retarded migration in polyacrylamide gel electrophoresis (PAGE) but TnAn segments do not. We have determined a high-resolution structure of a 16 base-pair DNA oligomer [d(CGAGGTTTAAACCTCG)2] containing a T3A3 tract. The refinement included residual dipolar coupling (RDC) restraints. A radius of curvature, Rc, analysis was used to measure the overall bending in the molecule. A plot of the helical axis reference points showed a sigmoidal shape indicating a discontinuity at the central TpA step in the overall curvature. Because of the length of the 16mer it was possible to accurately measure Rc for the two halves of the molecule centered about the TpA step. The Rc for the overall molecule (Rc=128 nm) is large, i.e. a low overall magnitude of global bending, whereas the Rc for the two halves of the molecule (Rc=72 nm) is small indicating a much larger magnitude of localized bending. However the direction of bending in the two halves is in partial opposition leading to cancellation of the overall bending. This indicates that TnAn-containing DNA sequences are strongly and multiply bentmore » locally, but are nearly straight globally which is consistent with PAGE results. The RDC refined structure lacked anomalous features present in NOE-only structures indicating the RDC and NOE measurements have a different sensitivity to conformational dynamics at the central TpA step. Because of its increased length and refinement using RDC restraints, the structure of the 16mer reported here provides new insight into the structural origins of the enigmatic PAGE behavior of AnTn and TnAn tracts and the large amplitude, slow base dynamics observed at TpA steps.« less

NMR evidence is presented indicating that the exceptional conformational dynamics found at TpA steps in DNA is general to all immediate sequence contexts. One easily tractable NMR parameter that is sensitive to TpA base dynamics is the resonance linewidth of the TpA adenine H2 proton. This resonance experiences a temperature-dependent broadening due to conformational dynamics. Unusual dynamics at TpA steps were originally observed in the sequence context (T)pTpTpApAp(A). We have since shown that the evidence for TpA dynamics persists when either the thymine preceding the TpA step or the adenine following the TpA step is preserved [McAteer et al., Nucleic Acids Res. 23, 3962-3966 (1995)]. Here, in order establish whether or not exceptional TpA dynamics occurs in all DNA sequence contexts, we investigated a series of DNA sequences of the form GCNaTANbNbTANaGC, where N=A,T,C,G. In this family of sequences, all 16 possible immediate sequence context environments of the form NaTANb were examined using 10 DNA sequences. Our NMR results show that the TpA adenine H2 resonance contains a temperature dependent excess linewidth indicative of dynamics in all 16 sequence context environments. By studying a complete set of sequence contexts, it was possible to recognize trends relating resonance parameters and sequence environment. For example, the magnitude of the maximum linewidth is largely determined by the identity of the nucleotide following the TpA step and the magnitude of the linewidth maximum is moderately correlated (r=0.56) with the temperature of the linewidth maximum. The physical basis for these correlations is discussed.

Current approaches in the study of protein function following immobilization on solid supports do not address how the properties of the support material affect protein structure and dynamics. In order to more efficiently create appropriate immobilization systems, general principles regarding the interaction between the material properties and broad classes of proteins need to be understood. Our goal is to gain some insight as to what material properties are optimal in the functional immobilization of proteins through the study of maltose binding protein (MBP) immobilized in a calmodulin (CaM) hydrogel. Maltose binding protein undergoes a ∼35° rigid body domain reorientation upon binding of sugar that can be detected by solution NMR. We hypothesize that through the study of the interactions of MBP with the CaM hydrogel (through a range of pHs) we can form some general principles regarding material:protein interfaces that will be applicable in the functional immobilization of broad classes of proteins. In this study we tagged MBP with a peptide target of calmodulin which immobilized MBP to a micro-sized spherical hydrogel composed of Cys-mutated CaM cross-linked to poly(ethylene glycol) diacrylate (PEGDA). Immobilization of MBP was confirmed by fluorescence and NuPAGE protein gels, and gel optimization for maximum immobilization is underway. Various NMR experiments have been and will be done in order to confirm function and to analyze the interactions between the immobilized protein and the hydrogel.