Global Fold Research Articles

On the role of physics and evolution in dictating protein structure and function.

How many of the structural and functional properties of proteins are inherent? Computer simulations provide a powerful tool to address this question. A series of studies on QS, quasi-spherical, compact polypeptides which lack any secondary structure; ART, artificial, proteins comprised of compact homopolypeptides with protein-like secondary structure; and PDB, native, single domain proteins shows that essentially all native global folds, pockets and protein-protein interfaces are in the ART library. This suggests that many protein properties are inherent and that evolution is involved in fine-tuning. The completeness of the space of ligand binding pockets and protein-protein interfaces suggests that promiscuous interactions are intrinsic to proteins and that the capacity to perform the biochemistry of life at low level does not require evolution. If so, this has profound consequences for the origin of life.

An approach to sequential NMR assignments of proteins: application to chemical shift restraint-based structure prediction.

A procedure for the simultaneous acquisition of {HNCOCANH & HCCCONH} chemical shift correlation spectra employing sequential [Formula: see text] data acquisition for moderately sized proteins is presented. The suitability of the approach for obtaining sequential resonance assignments, including complete [Formula: see text] and [Formula: see text] chemical shift information, is demonstrated experimentally for a [Formula: see text] and [Formula: see text] labelled sample of the C-terminal winged helix (WH) domain of the minichromosome maintenance (MCM) complex of Sulfolobus solfataricus. The chemical shift information obtained was used to calculate the global fold of this winged helix domain via CS-Rosetta. This demonstrates that our procedure provides a reliable and straight-forward protocol for a quick global fold determination of moderately-sized proteins.

From local structure to a global framework: recognition of protein folds.

Protein folding has been a major area of research for many years. Nonetheless, the mechanisms leading to the formation of an active biological fold are still not fully apprehended. The huge amount of available sequence and structural information provides hints to identify the putative fold for a given sequence. Indeed, protein structures prefer a limited number of local backbone conformations, some being characterized by preferences for certain amino acids. These preferences largely depend on the local structural environment. The prediction of local backbone conformations has become an important factor to correctly identifying the global protein fold. Here, we review the developments in the field of local structure prediction and especially their implication in protein fold recognition.

The S229L Colon Tumor-associated Variant of DNA Polymerase β Induces Cellular Transformation as a Result of Decreased Polymerization Efficiency

DNA polymerase β (Pol β) plays a key role in base excision repair (BER) by filling in small gaps that are generated after base adducts are excised from the DNA. Pol β is mutated in a large number of colorectal tumors, and these mutations may drive carcinogenesis. In the present study, we wished to determine whether the S229L somatic Pol β variant identified in a stage 3 colorectal tumor is a driver of carcinogenesis. We show that S229L does not possess any defects in binding to either DNA or nucleotides compared with the WT enzyme, but exhibits a significant loss of polymerization efficiency, largely due to an 8-fold decrease in the polymerization rate. S229L participates in BER, but due to its lower catalytic rate, does so more slowly than WT. Expression of S229L in mammalian cells induces the accumulation of BER intermediate substrates, chromosomal aberrations, and cellular transformation. Our results are consistent with the interpretation that S229L is a driver of carcinogenesis, likely as a consequence of its slow polymerization activity during BER in vivo.

Photosensitive Polyamines for High-Performance Photocontrol of DNA Higher-Order Structure

Polyamines are small, ubiquitous, positively charged molecules that play an essential role in numerous biological processes such as DNA packaging, gene regulation, neuron activity, and cell proliferation. Here, we synthesize the first series of photosensitive polyamines (PPAs) and demonstrate their ability to photoreversibly control nanoscale DNA higher-order structure with high efficiency. We show with fluorescence microscopy imaging that the efficiency of the PPAs as DNA-compacting agents is directly correlated to their molecular charge. Micromolar concentration of the most efficient molecule described here, a PPA containing three charges at neutral pH, compacts DNA molecules from a few kilobase pairs to a few hundred kilobase pairs, while subsequent 3 min UV illuminations at 365 nm triggers complete unfolding of DNA molecules. Additional application of blue light (440 nm for 3 min) induces the refolding of DNA into the compact state. Atomic force microscopy reveals that the compaction involves a global folding of the whole DNA molecule, whereas UV-induced unfolding is a modification initiated from the periphery of the compacted DNA, resulting in the occurrence of intermediate flower-like structures prior to the fully unfolded state.

A Molecular Interpretation of 2D IR Protein Folding Experiments with Markov State Models

The folding mechanism of the N-terminal domain of ribosomal protein L9 (NTL91–39) is studied using temperature-jump (T-jump) amide I′ two-dimensional infrared (2D IR) spectroscopy in combination with spectral simulations based on a Markov state model (MSM) built from millisecond-long molecular dynamics trajectories. The results provide evidence for a compact well-structured folded state and a heterogeneous fast-exchanging denatured state ensemble exhibiting residual secondary structure. The folding rate of 26.4 μs−1 (at 80°C), extracted from the T-jump response of NTL91–39, compares favorably with the 18 μs−1 obtained from the MSM. Structural decomposition of the MSM and analysis along the folding coordinate indicates that helix-formation nucleates the global folding. Simulated difference spectra, corresponding to the global folding transition of the MSM, are in qualitative agreement with measured T-jump 2D IR spectra. The experiments demonstrate the use of T-jump 2D IR spectroscopy as a valuable tool for studying protein folding, with direct connections to simulations. The results suggest that in addition to predicting the correct native structure and folding time constant, molecular dynamics simulations carried out with modern force fields provide an accurate description of folding mechanisms in small proteins.

The k-junction motif in RNA structure

The k-junction is a structural motif in RNA comprising a three-way helical junction based upon kink turn (k-turn) architecture. A computer program written to examine relative helical orientation identified the three-way junction of the Arabidopsis TPP riboswitch as an elaborated k-turn. The Escherichia coli TPP riboswitch contains a related k-junction, and analysis of >11 000 sequences shows that the structure is common to these riboswitches. The k-junction exhibits all the key features of an N1-class k-turn, including the standard cross-strand hydrogen bonds. The third helix of the junction is coaxially aligned with the C (canonical) helix, while the k-turn loop forms the turn into the NC (non-canonical) helix. Analysis of ligand binding by ITC and global folding by gel electrophoresis demonstrates the importance of the k-turn nucleotides. Clearly the basic elements of k-turn structure are structurally well suited to generate a three-way helical junction, retaining all the key features and interactions of the k-turn.

The binding mechanisms of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins are very common and instrumental for cellular signaling. Recently, a number of studies have investigated the kinetic binding mechanisms of IDPs and IDRs. These results allow us to draw conclusions about the energy landscape for the coupled binding and folding of disordered proteins. The association rate constants of IDPs cover a wide range (10(5)-10(9) M(-1) s(-1)) and are largely governed by long-range charge-charge interactions, similarly to interactions between well-folded proteins. Off-rate constants also differ significantly among IDPs (with half-lives of up to several minutes) but are usually around 0.1-1000 s(-1), allowing for rapid dissociation of complexes. Likewise, affinities span from pM to μM suggesting that the low-affinity high-specificity concept for IDPs is not straightforward. Overall, it appears that binding precedes global folding although secondary structure elements such as helices may form before the protein-protein interaction. Short IDPs bind in apparent two-state reactions whereas larger IDPs often display complex multi-step binding reactions. While the two extreme cases of two-step binding (conformational selection and induced fit) or their combination into a square mechanism is an attractive model in theory, it is too simplistic in practice. Experiment and simulation suggest a more complex energy landscape in which IDPs bind targets through a combination of conformational selection before binding (e.g., secondary structure formation) and induced fit after binding (global folding and formation of short-range intermolecular interactions).

Protein Stability in Living Cells

Theories concerning the effects of macromolecular crowding assert that the biophysical properties of proteins and nucleic acids can be significantly altered in native cellular environments relative to buffer alone. Despite a growing number of studies probing equilibrium thermodynamic protein stability in cells, there remains a lack of quantitative information, especially regarding residue-level stability under non-perturbing conditions. We have measured the in-cell stability of the 56-amino acid B1 domain of protein G (GB1) at the residue level without using destabilizing solutes or thermal modification by using NMR-detected hydrogen-deuterium exchange of quenched cell lysates. Comparison to dilute solution (pH 7.6 and 37 °C) shows that residues are stabilized in Escherichia coli cells by as much as 1.1 ±0.1 kcal/mol (Figure 1). We have also identified the residues most important for global folding of GB1 in cells. We discuss the implications of our findings with respect to structural models gleaned from studies conducted in buffer alone.View Large Image | View Hi-Res Image | Download PowerPoint Slide

Transient growth arrest in Escherichia coli induced by chromosome condensation.

MukB is a bacterial SMC (structural maintenance of chromosome) protein that regulates the global folding of the Escherichia coli chromosome by bringing distant DNA segments together. We report that moderate overproduction of MukB may lead, depending on strain and growth conditions, to transient growth arrest. In DH5α cells, overproduction of MukB or MukBEF using pBAD expression system triggered growth arrest 2.5 h after induction. The exit from growth arrest was accompanied by the loss of the overproducing plasmid and a decline in the abundance of MukBEF. The arrested cells showed a compound gene expression profile which can be characterized by the following features: (i) a broad and deep downregulation of ribosomal proteins (up to 80-fold); (ii) downregulation of groups of genes encoding enzymes involved in nucleotide metabolism, respiration, and central metabolism; (iii) upregulation of some of the genes responsive to general stress; and (iv) degradation of the patterns of spatial correlations in the transcriptional activity of the chromosome. The transcriptional state of the MukB induced arrest is most similar to stationary cells and cells recovered from stationary phase into a nutrient deprived medium, to amino acid starved cells and to the cells shifting from glucose to acetate. The mukB++ state is dissimilar from all examined transcriptional states generated by protein overexpression with the possible exception of RpoE and RpoH overexpression. Thus, the transcription profile of MukB-arrested cells can be described as a combination of responses typical for other growth-arrested cells and those for overproducers of DNA binding proteins with a particularly deep down-regulation of ribosomal genes.

Porphyrin Assemblies and Their Scaffolds

The chlorophyll and heme molecules of chloroplasts and mitochondria are brought to life by "the global fold of the protein scaffolds". Proteins in hydrophobic cell regions touch the dye platelets from both sides, pushing and orienting them according to their life-spending activities in light and electron transfers. The conjugated π-electron systems or planarity of the porphyrin macrocycles are never disturbed. Most artificial porphyrin assemblies contain meso-tetraphenylporphyrins (TPPs), because the four phenyl groups rotate freely and carry their substituents above or below the macrocycle. A single porphyrin molecule can, for example, be attached to an anionic surface with ammonium groups on its 2,3-carbons, be located within a hydrophobic membrane with its alkyl chains on the 4-position, and then fixate a cationic polymer with 4,5-sulphonates. Charged TPPs also show unique spectroscopic changes at different pH values and a reversible loss of the macrocycle's planarity. On smooth silicate, graphite, or gold scaffolds TPPs have been attached irreversibly as single molecules, as extended non-covalent H or J aggregates as well as acetylene or thiophene-linked polymers. Soft, mobile porphyrin ladders conduct excited electrons ("excitons") better than rigid porphyrin wires ("polarons").

Architecture of the hepatitis C virus E1 glycoprotein transmembrane domain studied by NMR

Oligomerization of hepatitis C viral envelope proteins E1 and E2 is essential to virus fusion and assembly. Although interactions within the transmembrane (TM) domains of these glycoproteins have proven contributions to the E1/E2 heterodimerization process and consequent infectivity, there is little structural information on this entry mechanism. Here, as a first step towards our long-term goal of understanding the interaction between E1 and E2 TM-domains, we have expressed, purified and characterized E1-TM using structural biomolecular NMR methods. An MBP-fusion expression system yielded sufficient quantities of pure E1-TM, which was solubilized in two membrane-mimicking environments, SDS- and LPPG-micelles, affording samples amenable to NMR studies. Triple resonance assignment experiments and relaxation measurements provided information on the secondary structure and global fold of E1-TM in these environments. In SDS micelles E1-TM adopts a helical conformation, with helical stretches at residues 354–363 and 371–379 separated by a more flexible segment of residues 364–370. In LPPG micelles a helical conformation was observed for residues 354–377 with greater flexibility in the 366–367 dyad, suggesting LPPG provides a more native environment for the peptide. Replacement of key positively charged residue K370 with an alanine did not affect the secondary structure of E1-TM but did change the relative positioning within the micelle of the two helices. These results lay the foundation for structure determination of E1-TM and a molecular understanding of how E1-TM flexibility enhances its interaction with E2-TM during heterodimerization and membrane fusion.

Global fold of human cannabinoid type 2 receptor probed by solid‐state 13C‐, 15N‐MAS NMR and molecular dynamics simulations

ABSTRACTThe global fold of human cannabinoid type 2 (CB2) receptor in the agonist‐bound active state in lipid bilayers was investigated by solid‐state 13C‐ and 15N magic‐angle spinning (MAS) NMR, in combination with chemical‐shift prediction from a structural model of the receptor obtained by microsecond‐long molecular dynamics (MD) simulations. Uniformly 13C‐ and 15N‐labeled CB2 receptor was expressed in milligram quantities by bacterial fermentation, purified, and functionally reconstituted into liposomes. 13C MAS NMR spectra were recorded without sensitivity enhancement for direct comparison of Cα, Cβ, and CO bands of superimposed resonances with predictions from protein structures generated by MD. The experimental NMR spectra matched the calculated spectra reasonably well indicating agreement of the global fold of the protein between experiment and simulations. In particular, the 13C chemical shift distribution of Cα resonances was shown to be very sensitive to both the primary amino acid sequence and the secondary structure of CB2. Thus the shape of the Cα band can be used as an indicator of CB2 global fold. The prediction from MD simulations indicated that upon receptor activation a rather limited number of amino acid residues, mainly located in the extracellular Loop 2 and the second half of intracellular Loop 3, change their chemical shifts significantly (≥1.5 ppm for carbons and ≥5.0 ppm for nitrogens). Simulated two‐dimensional 13Cα(i)13CO(i) and 13CO(i)15NH(i + 1) dipolar‐interaction correlation spectra provide guidance for selective amino acid labeling and signal assignment schemes to study the molecular mechanism of activation of CB2 by solid‐state MAS NMR. Proteins 2014; 82:452–465. © 2013 Wiley Periodicals, Inc.

RNA does the folding dance of twist, turn, stack

RNA molecules have long been known to adopt unusual conformations. The first crystal structure of an RNA was that of yeast tRNAphe (1, 2), which provided a portent of what was to come. With more structures available now from X-ray crystallography and NMR, that promise has been met beyond expectations. The RNAs that are folded in molecular dynamics simulations by Chen and Garcia (3) are tetraloops (RNA hairpin loops of only four nucleotides) that are among the most unusual small structures known today. Although or because they are unusual with atypical thermal stability, these loops are also ubiquitous in longer RNA strands [rUUCG (4); rGNRA and rCUUG (5)]; they were first identified as the most abundant hairpins in prokaryotic rRNA. The function of these loops varies: the GNRA tetraloops (where N is any base and R is a purine) is often found in RNA:RNA tertiary interactions to anchor a global fold. The UUCG tetraloop is solitary, with no known interactions with RNA or proteins: it is thought to nucleate folding and so direct in vivo RNA folding patterns. These little hairpins have been the objects of considerable in vitro experimental structural characterization by NMR and X-ray crystallography, and are found to consistently form a single unique structure. The rUUCG tetraloop in particular has been subjected to many thermodynamic investigations to measure the sources of its unusual stability (6). The rCUUG hairpin is structurally the most flexible, because it tucks its first U into the minor groove but its second U is relatively unconstrained. Among the three tetraloops, CUUG has a propensity for significant conformational exchange (7). Although the three tetraloops used in the Chen and Garcia (3) computational folding experiments are all unusually stable, they achieve their stability by very different means. Their structural diversity makes them excellent candidates for folding studies, because no one mechanism dominates their ability to adopt their unique structures.

Folding Dynamics of the Trp-Cage Miniprotein: Evidence for a Native-Like Intermediate from Combined Time-Resolved Vibrational Spectroscopy and Molecular Dynamics Simulations

Trp-cage is a synthetic 20-residue miniprotein which folds rapidly and spontaneously to a well-defined globular structure more typical of larger proteins. Due to its small size and fast folding, it is an ideal model system for experimental and theoretical investigations of protein folding mechanisms. However, Trp-cage's exact folding mechanism is still a matter of debate. Here we investigate Trp-cage's relaxation dynamics in the amide I' spectral region (1530-1700 cm(-1)) using time-resolved infrared spectroscopy. Residue-specific information was obtained by incorporating an isotopic label ((13)C═(18)O) into the amide carbonyl group of residue Gly11, thereby spectrally isolating an individual 310-helical residue. The folding-unfolding equilibrium is perturbed using a nanosecond temperature-jump (T-jump), and the subsequent re-equilibration is probed by observing the time-dependent vibrational response in the amide I' region. We observe bimodal relaxation kinetics with time constants of 100 ± 10 and 770 ± 40 ns at 322 K, suggesting that the folding involves an intermediate state, the character of which can be determined from the time- and frequency-resolved data. We find that the relaxation dynamics close to the melting temperature involve fast fluctuations in the polyproline II region, whereas the slower process can be attributed to conformational rearrangements due to the global (un)folding transition of the protein. Combined analysis of our T-jump data and molecular dynamics simulations indicates that the formation of a well-defined α-helix precedes the rapid formation of the hydrophobic cage structure, implying a native-like folding intermediate, that mainly differs from the folded conformation in the orientation of the C-terminal polyproline II helix relative to the N-terminal part of the backbone. We find that the main free-energy barrier is positioned between the folding intermediate and the unfolded state ensemble, and that it involves the formation of the α-helix, the 310-helix, and the Asp9-Arg16 salt bridge. Our results suggest that at low temperature (T ≪ Tm) a folding path via formation of α-helical contacts followed by hydrophobic clustering becomes more important.

Ligand Binding Site Detection by Local Structure Alignment and Its Performance Complementarity

Accurate determination of potential ligand binding sites (BS) is a key step for protein function characterization and structure-based drug design. Despite promising results of template-based BS prediction methods using global structure alignment (GSA), there is room to improve the performance by properly incorporating local structure alignment (LSA) because BS are local structures and often similar for proteins with dissimilar global folds. We present a template-based ligand BS prediction method using G-LoSA, our LSA tool. A large benchmark set validation shows that G-LoSA predicts drug-like ligands' positions in single-chain protein targets more precisely than TM-align, a GSA-based method, while the overall success rate of TM-align is better. G-LoSA is particularly efficient for accurate detection of local structures conserved across proteins with diverse global topologies. Recognizing the performance complementarity of G-LoSA to TM-align and a nontemplate geometry-based method, fpocket, a robust consensus scoring method, CMCS-BSP (Complementary Methods and Consensus Scoring for ligand Binding Site Prediction), is developed and shows improvement on prediction accuracy.

The study of pH-dependent stability shows that the TPLH-mediated hydrogen-bonding network is important for the conformation and stability of human gankyrin.

Ankyrin repeat (AR) proteins possess a distinctive modular and repetitive architecture, and their global folds are maintained by short-range interactions in terms of the primary sequence. In this work, we extended our previous study on the highly conserved TPLH tetrapeptide and investigated the impact of a solvent-exposed histidine residue on the pH-dependent stability of gankyrin, providing further insight into the contribution of the TPLH motif to the tertiary fold of AR proteins. Consisting of seven ARs, gankyrin has five histidine residues in TPLH motifs or its variants, all of which adopt a H(ε2)-tautermeric form and are shielded from solvent. By truncating the C-terminal ankyrin repeat (AR7), we exposed H177 in the (174)TPLH(177) of AR6 (the second C-terminal AR) to an aqueous environment. We showed that this truncated gankyrin mutant, namely, Gank(1-201), was well-folded at a neutral pH with a slightly lower stability with respect to gankyrin wild type (WT). However, unlike gankyrin WT, the stability of Gank(1-201) was markedly decreased together with a loss of conformation at a pH slightly below 6.0. It was rationalized that the protonation of the H177 imidazole ring triggered the disruption of the TPLH-mediated hydrogen-bonding network, which in turn led to the loss of conformation and stability. These results together with the work on Q210H mutant nicely explain that the C-terminal AR7 has a (207)TPLQ(210) variant and are in support of the notion that a string of TPLH/variant, which may arguably act like a zip lock to the elongated AR proteins via intra-/inter-repeat hydrogen-bonding, is important in maintaining the tertiary fold. Additionally, we made rational mutagenesis to introduce extra surface charge on AR7 of gankyrin and demonstrated that G214E and I219D mutations increased the stability of gankyrin while the function remained intact. Taken together, our results indicate that the TPLH-mediated hydrogen-bonding network is important for the conformation and stability of human gankyrin, and the C-terminal AR contributes to the conformational stability of gankyrin (AR proteins in general) through shielding this TPLH network from solvent as well as making the surface area more accessible to solvent.

Structural Determinants of RGS-RhoGEF Signaling Critical to Entamoeba histolytica Pathogenesis

G protein signaling pathways, as key components of physiologic responsiveness and timing, are frequent targets for pharmacologic intervention. Here, we identify an effector for heterotrimeric G proteinα subunit (EhGα1) signaling from Entamoeba histolytica, the causative agent of amoebic colitis. EhGα1 interacts with this effector and guanosine triphosphatase-accelerating protein, EhRGS-RhoGEF, in a nucleotide state-selective fashion. Coexpression of EhRGS-RhoGEF with constitutively active EhGα1 and EhRacC leads to Rac-dependent spreading in Drosophila S2 cells. EhRGS-RhoGEF overexpression in E.histolytica trophozoites leads to reduced migration toward serum and lower cysteine protease activity, as well as reduced attachment to, and killingof, host cells. A 2.3Å crystal structure of thefull-length EhRGS-RhoGEF reveals a putative inhibitory helix engaging the Dbl homology domain Rho-binding surface and the pleckstrin homology domain. Mutational analysis of the EhGα1/EhRGS-RhoGEF interface confirms a canonical "regulator of G protein signaling" domain rather than a RhoGEF-RGS ("rgRGS") domain, suggesting a convergent evolution toward heterotrimeric and small G protein cross-talk.

Mutational analysis of MukE reveals its role in focal subcellular localization of MukBEF

Bacterial condensin MukBEF is essential for global folding of the Escherichia coli chromosome. MukB, a SMC (structural maintenance of chromosome) protein, comprises the core of this complex and is responsible for its ATP-modulated DNA binding and reshaping activities. MukF serves as a kleisin that modulates MukB-DNA interactions and links MukBs into macromolecular assemblies. Little is known about the function of MukE. Using random mutagenesis, we generated six loss-of-function point mutations in MukE. The surface mutations clustered in two places. One of them was at or close to the interface with MukF while the other was away from the known interactions of the protein. All loss-of-function mutations affected focal localization of MukBEF in live cells. In vitro, however, only some of them interfered with the assembly of MukBEF into a complex or the ability of MukEF to disrupt MukB-DNA interactions. Moreover, some MukE mutants were able to join intracellular foci formed by endogenous MukBEF and most of the mutants were efficiently incorporated into MukBEF even in the presence of endogenous MukE. These data reveal that focal localization of MukBEF involves other activities besides DNA binding and that MukE plays a central role in them.

Novel Surface-Smoothing Based Local Gyrification Index

The quantication of cortical surface folding is important in identifying and classifying many neurodegenerative diseases. Much work has been done to identify regional and global brain folding, and in this paper we review some of these methods, as well as propose a new method that has advantages over the existing state of art. Using our novel proposed method, we mapped the local gyrification index on the cortical surface for subjects with mild Alzheimer's dementia (n=20) , very mild dementia (n=23) and age-matched healthy subjects (n=52). In our experiments we find a consistent pattern of gyrification changes in the dementia subjects, with regions generally affected early on in the progression of Alzheimer pathology, including medial temporal lobe, and cingulate gyrus, having decreased gyrification. At the same time we observe increased gyrification in dementia subjects, in frontal, anterior temporal and posteriorly located regions. We speculate that in neurodegenerative diseases including Alzheimer Disease, the folding of the entire cortical mantle undergoes dynamic changes as regional atrophy begins and expands, with both decreases and increases in gyrification.