Validity of Gō models: comparison with a solvent-shielded empirical energy decomposition

Do Gō-type model potentials provide a valid approach for studying protein folding? They have been widely used for this purpose because of their simplicity and the speed of simulations based on their use. The essential assumption in such models is that only contact interactions existing in the native state determine the energy surface of a polypeptide chain, even for non-native configurations sampled along folding trajectories. Here we use an all-atom molecular mechanics energy function to investigate the adequacy of Gō-type potentials. We show that, although the contact approximation is accurate, non-native contributions to the energy can be significant. The assumed relation between residue-residue interaction energies and the number of contacts between them is found to be only approximate. By contrast, individual residue energies correlate very well with the number of contacts. The results demonstrate that models based on the latter should give meaningful results (e.g., as used to interpret phi values), whereas those that depend on the former are only qualitative, at best.     

Incidence and molecular epidemiology of Pseudomonas aeruginosa bacteremias in patients with acute leukemia: analysis by pulsed-field gel electrophoresis

The incidence and molecular epidemiology of P. aeruginosa bacteremias, were monitored in patients with acute leukemia to define mechanisms of possible nosocomial transmission. From September 1997 to March 2001 febrile episodes were examined and blood isolates of P. aeruginosa were studied employing Pulsed-Field gel Electrophoresis (PFGE). Evaluation of DNA correlation was performed according to Tenover criteria. A total of 309 febrile episodes occurred in 187 patients. Of 139 organisms isolated in 116 bacteremias, 48% were gram negative bacilli (GNB); P. aeruginosa bacteremias were recorded in 34 (51%) of GNB sepsis. Evaluation of DNA correlation showed 2 related in 1997, 7 related in 1998, 10 related in 1999, 6 related in 2000-2001 (mainly closely and possibly related); therefore isolates closely related among themselves were also possibly related with other strains. About 60% of patients with related strains were hospitalized in the same room or in different rooms but became infected in the same period. Our data suggest a horizontal spread among the patients even if other sources were possible. The study assessed the usefulness of PFGE in bacteriological epidemiology.     

Pulling geometry defines the mechanical resistance of a beta-sheet protein

Proteins show diverse responses when placed under mechanical stress. The molecular origins of their differing mechanical resistance are still unclear, although the orientation of secondary structural elements relative to the applied force vector is thought to have an important function. Here, by using a method of protein immobilization that allows force to be applied to the same all-beta protein, E2lip3, in two different directions, we show that the energy landscape for mechanical unfolding is markedly anisotropic. These results, in combination with molecular dynamics (MD) simulations, reveal that the unfolding pathway depends on the pulling geometry and is associated with unfolding forces that differ by an order of magnitude. Thus, the mechanical resistance of a protein is not dictated solely by amino acid sequence, topology or unfolding rate constant, but depends critically on the direction of the applied extension.     

The conformation of peptide thymosin α1 in solution and in a membrane-like environment by circular dichroism and NMR spectroscopy. a possible model for its interaction with the lymphocyte membrane

The 28-residue peptide thymosin α1 was studied by circular dichroism and two-dimensional NMR. Circular dichroism indicates that thymosin α1 in water solution does not assume a preferred conformation, while in the presence of small unilamellar vesicles of dimiristoylphosphatidylcholine and dimiristoylphosphatidic acid (10:1) and in sodium dodecyl sulphate, it assumes a partly structured conformation. Presence of zinc ions produces similar effects. In a more hydrophobic environment like a solution of a mixed solvent water-2,2,2 trifluoroethanol, it adopts a structured conformation. NMR spectra indicated that in this mixture as solvent, thymosin α1 has a structure characterized by two regions. A β-turn is present between residue 5 and residue 8, while the region between residues 17 and 24 shows an α helix conformation. These changes of conformation in different environments may be considered structural requirements in the steps of its interaction with the lymphocyte membrane. In fact, these conformational changes may correspond to the first event of the mechanism of lymphocyte activation in the immune response modulation by thymosin α1.     

The Role of High-Dimensional Diffusive Search,Stabilization, and Frustration in Protein Folding

Proteins are polymeric molecules with many degrees of conformational freedom whose internal energetic interactions are typically screened to small distances. Therefore, in the high-dimensional conformation space of a protein, the energy landscape is locally relatively flat, in contrast to low-dimensional representations, where, because of the induced entropic contribution to the full free energy, it appears funnel-like. Proteins explore the conformation space by searching these flat subspaces to find a narrow energetic alley that we call a hypergutter and then explore the next, lower-dimensional, subspace. Such a framework provides an effective representation of the energy landscape and folding kinetics that does justice to the essential characteristic of high-dimensionality of the search-space. It also illuminates the important role of nonnative interactions in defining folding pathways. This principle is here illustrated using a coarse-grained model of a family of three-helix bundle proteins whose conformations, once secondary structure has formed, can be defined by six rotational degrees of freedom. Two folding mechanisms are possible, one of which involves an intermediate. The stabilization of intermediate subspaces (or states in low-dimensional projection) in protein folding can either speed up or slow down the folding rate depending on the amount of native and nonnative contacts made in those subspaces. The folding rate increases due to reduced-dimension pathways arising from the mere presence of intermediate states, but decreases if the contacts in the intermediate are very stable and introduce sizeable topological or energetic frustration that needs to be overcome. Remarkably, the hypergutter framework, although depending on just a few physically meaningful parameters, can reproduce all the types of experimentally observed curvature in chevron plots for realizations of this fold.     

Determination of free and total S-phenylmercapturic acid by HPLC/MS/MS in the biological monitoring of benzene exposure

Urinary S-phenylmercapturic acid (SPMA) is a biomarker suggested by the American Conference of Governmental Industrial Hygienists (ACGIH) for assessing occupational exposure to benzene. A possible cause of the miscorrelation between environmental monitoring and biological monitoring for benzene exposure, which many authors complain about, is the existence of a urinary metabolite that turns into SPMA by acid hydrolysis. Forty urine samples were tested to determine which concentration value would correspond to the ACGIH Biological Exposure Index (BEI) of 25 µg g^-1 creatinine if exposure assessment was based on the determination of SPMA after quantitative hydrolysis of its precursor. An aliquot of each sample was hydrolysed with 9 M H₂SO₄, a second one was brought to pH 2 and a third one was used as it was (free SPMA). SPMA was determined by high-performance liquid chromatography/tandem mass spectrometric technique (HPLC/MS/MS) using an internal standard. The analytical method was validated in the range 0.5-50 µg l^-1. The average SPMA in pH 2 samples is 45-60% of the total, while free SPMA varies from 1% to 66%. The hydrolysis of pre-SPMA reduces the likelihood of variability in the results by reducing pH differences in urine samples and increasing the amount of measured SPMA. The BEI limit value would be about 50 µg g^-1 creatinine.     

Solution structure of ApaG from Xanthomonas axonopodis pv. citri reveals a fibronectin-3 fold

ApaG proteins are found in a wide variety of bacterial genomes but their function is as yet unknown. Some eukaryotic proteins involved in protein-protein interactions, such as the human polymerase delta-interacting protein (PDIP38) and the F Box A (FBA) proteins, contain ApaG homology domains. We have used NMR to determine the solution structure of ApaG protein from the plant pathogen Xanthomonas axonopodis pv. citri (ApaG(Xac)) with the aim to shed some light on its molecular function. ApaG(Xac) is characterized by seven antiparallel beta strands forming two beta sheets, one containing three strands (ABE) and the other four strands (GFCC'). Relaxation measurements indicate that the protein has a quite rigid structure. In spite of the presence of a putative GXGXXG pyrophosphate binding motif ApaG(Xac) does not bind ATP or GTP, in vitro. On the other hand, ApaG(Xac) adopts a fibronectin type III (Fn3) fold, which is consistent with the hypothesis that it is involved in mediating protein-protein interactions. The fact that the proteins of ApaG family do not display significant sequence similarity with the Fn3 domains found in other eukaryotic or bacterial proteins suggests that Fn3 domain may have arisen earlier in evolution than previously estimated.     

Prediction of stability changes upon mutation in an icosahedral capsid

Identifying the contributions to thermodynamic stability of capsids is of fundamental and practical importance. Here we use simulation to assess how mutations affect the stability of lumazine synthase from the hyperthermophile Aquifex aeolicus, a T = 1 icosahedral capsid; in the simulations the icosahedral symmetry of the capsid is preserved by simulating a single pentamer and imposing crystal symmetry, in effect simulating an infinite cubic lattice of icosahedral capsids. The stability is assessed by estimating the free energy of association using an empirical method previously proposed to identify biological units in crystal structures. We investigate the effect on capsid formation of seven mutations, for which it has been experimentally assessed whether they disrupt capsid formation or not. With one exception, our approach predicts the effect of the mutations on the capsid stability. The method allows the identification of interaction networks, which drive capsid assembly, and highlights the plasticity of the interfaces between subunits in the capsid. Proteins 2015; 83:1733–1741. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc