Polymer principles and protein folding

ChSeq: A database of chameleon sequences

Chameleon sequences (ChSeqs) refer to sequence strings of identical amino acids that can adopt different conformations in protein structures. Researchers have detected and studied ChSeqs to understand the interplay between local and global interactions in protein structure formation. The different secondary structures adopted by one ChSeq challenge sequence-based secondary structure predictors. With increasing numbers of available Protein Data Bank structures, we here identify a large set of ChSeqs ranging from 6 to 10 residues in length. The homologous ChSeqs discovered highlight the structural plasticity involved in biological function. When compared with previous studies, the set of unrelated ChSeqs found represents an about 20-fold increase in the number of detected sequences, as well as an increase in the longest ChSeq length from 8 to 10 residues. We applied secondary structure predictors on our ChSeqs and found that methods based on a sequence profile outperformed methods based on a single sequence. For the unrelated ChSeqs, the evolutionary information provided by the sequence profile typically allows successful prediction of the prevailing secondary structure adopted in each protein family. Our dataset will facilitate future studies of ChSeqs, as well as interpretations of the interplay between local and nonlocal interactions. A user-friendly web interface for this ChSeq database is available at prodata.swmed.edu/chseq.

Steered quantum dynamics for energy minimization

We introduce a quantum optimal control algorithm for energy minimization that combines the diffeomorphic modulation under observable response preserving homotopy (D-MORPH) gradient and the Broyden Fletcher Goldfarb Shanno (BFGS) iterative scheme for nonlinear optimization. An extended set of controls defining the time-dependent mass, dipole moment, and external perturbational field are optimized to find an effective Hamiltonian that steers the dynamics of the system into the global minimum without getting trapped into local minima. The algorithm is illustrated as applied to energy minimization on rugged surfaces and golf potentials comparable to those previously explored for testing quantum annealing methodologies.

Sequence determines degree of knottedness in a coarse-grained protein model

Knots are abundant in globular homopolymers but rare in globular proteins. To shed new light on this long-standing conundrum, we study the influence of sequence on the formation of knots in proteins under native conditions within the framework of the hydrophobic-polar lattice protein model. By employing large-scale Wang-Landau simulations combined with suitable Monte Carlo trial moves we show that even though knots are still abundant on average, sequence introduces large variability in the degree of self-entanglements. Moreover, we are able to design sequences which are either almost always or almost never knotted. Our findings serve as proof of concept that the introduction of just one additional degree of freedom per monomer (in our case sequence) facilitates evolution towards a protein universe in which knots are rare.

NMR Spectroscopic Studies of the Conformational Ensembles of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are characterized by substantial conformational flexibility and thus not amenable to conventional structural biology techniques. Given their inherent structural flexibility NMR spectroscopy offers unique opportunities for structural and dynamic studies of IDPs. The past two decades have witnessed significant development of NMR spectroscopy that couples advances in spin physics and chemistry with a broad range of applications. This chapter will summarize key advances in NMR methodology. Despite the availability of efficient (multi-dimensional) NMR experiments for signal assignment of IDPs it is discussed that NMR of larger and more complex IDPs demands spectral simplification strategies capitalizing on specific isotope-labeling strategies. Prototypical applications of isotope labeling-strategies are described. Since IDP-ligand association and dissociation processes frequently occur on time scales that are amenable to NMR spectroscopy we describe in detail the application of CPMG relaxation dispersion techniques to studies of IDP protein binding. Finally, we demonstrate that the complementary usage of NMR and EPR data provide a more comprehensive picture about the conformational states of IDPs and can be employed to analyze the conformational ensembles of IDPs.

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.

Conformational flexibility in the binding surface of the potassium channel blocker ShK

ShK is a 35-residue peptide that binds with high affinity to human voltage-gated potassium channels through a conserved K-Y dyad. Here we have employed NMR measurements of backbone-amide (15)N spin-relaxation rates to investigate motions of the ShK backbone. Although ShK is rigid on the ps to ns timescale, increased linewidths observed for 11 backbone-amide (15)N resonances identify chemical or conformational exchange contributions to the spin relaxation. Relaxation dispersion profiles indicate that exchange between major and minor conformers occurs on the sub-millisecond timescale. Affected residues are mostly clustered around the central helix-kink-helix structure and the critical K22-Y23 motif. We suggest that the less structured minor conformer increases the exposure of Y23, known to contribute to binding affinity and selectivity, thereby facilitating its interaction with potassium channels. These findings have potential implications for the design of new channel blockers based on ShK.

Effect of single-site mutations on hydrophobic-polar lattice proteins

We developed a heuristic method for determining the ground-state degeneracy of hydrophobic-polar (HP) lattice proteins, based on Wang-Landau and multicanonical sampling. It is applied during comprehensive studies of single-site mutations in specific HP proteins with different sequences. The effects in which we are interested include structural changes in ground states, changes of ground-state energy, degeneracy, and thermodynamic properties of the system. With respect to mutations, both extremely sensitive and insensitive positions in the HP sequence have been found. That is, ground-state energies and degeneracies, as well as other thermodynamic and structural quantities, may be either largely unaffected or may change significantly due to mutation.

NMR contributions to structural dynamics studies of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are characterized by substantial conformational plasticity. Given their inherent structural flexibility X-ray crystallography is not applicable to study these proteins. In contrast, NMR spectroscopy offers unique opportunities for structural and dynamic studies of IDPs. The past two decades have witnessed significant development of NMR spectroscopy that couples advances in spin physics and chemistry with a broad range of applications. This article will summarize key advances in basic physical-chemistry and NMR methodology, outline their limitations and envision future R&D directions.

Physicochemical bases for protein folding, dynamics, and protein-ligand binding

Proteins are essential parts of living organisms and participate in virtually every process within cells. As the genomic sequences for increasing number of organisms are completed, research into how proteins can perform such a variety of functions has become much more intensive because the value of the genomic sequences relies on the accuracy of understanding the encoded gene products. Although the static three-dimensional structures of many proteins are known, the functions of proteins are ultimately governed by their dynamic characteristics, including the folding process, conformational fluctuations, molecular motions, and protein-ligand interactions. In this review, the physicochemical principles underlying these dynamic processes are discussed in depth based on the free energy landscape (FEL) theory. Questions of why and how proteins fold into their native conformational states, why proteins are inherently dynamic, and how their dynamic personalities govern protein functions are answered. This paper will contribute to the understanding of structure-function relationship of proteins in the post-genome era of life science research.

Computational and theoretical methods for protein folding

A computational approach is essential whenever the complexity of the process under study is such that direct theoretical or experimental approaches are not viable. This is the case for protein folding, for which a significant amount of data are being collected. This paper reports on the essential role of in silico methods and the unprecedented interplay of computational and theoretical approaches, which is a defining point of the interdisciplinary investigations of the protein folding process. Besides giving an overview of the available computational methods and tools, we argue that computation plays not merely an ancillary role but has a more constructive function in that computational work may precede theory and experiments. More precisely, computation can provide the primary conceptual clues to inspire subsequent theoretical and experimental work even in a case where no preexisting evidence or theoretical frameworks are available. This is cogently manifested in the application of machine learning methods to come to grips with the folding dynamics. These close relationships suggested complementing the review of computational methods within the appropriate theoretical context to provide a self-contained outlook of the basic concepts that have converged into a unified description of folding and have grown in a synergic relationship with their computational counterpart. Finally, the advantages and limitations of current computational methodologies are discussed to show how the smart analysis of large amounts of data and the development of more effective algorithms can improve our understanding of protein folding.

Conformational analysis of Gly-Ala-NHMe in D(2)O and DMSO solutions: a two-dimensional infrared spectroscopy study

A relevant number of experiments on short peptides has been performed in recent years. One of the major problems rises from the simultaneous presence of slightly different conformers at equilibrium in solution. In the present paper, the conformational characteristics of the Gly-l-Ala-Methyl amide dipeptide in D2O and DMSO solutions are investigated by nonlinear IR spectroscopy. The pump-probe scheme with ultrashort mid-infrared pulses, in the Amide I region, is used to determine the mutual orientation of the two C═O bonds and the dynamics due to solute-solvent interactions. The coupling between Amide I modes is evaluated from both linear and 2D spectra. The interconversion between the different conformations occurs on time scales longer than the vibrational lifetime, and the spectral diffusion observed in 2D spectra is attributed to the solvent dynamics. Quantum mechanical calculations and molecular dynamics simulations are performed to identify the most stable geometries. By comparing the experimental and the theoretical data, we establish the prevalence of β-like polar conformers in both water and DMSO solvents.

On the subject of rigor in the study of amyloid β-protein assembly

According to Thomas Kuhn, the success of 'normal science,' the science we all practice on a daily basis, depends on the adherence to, and practice of, a paradigm accepted by the scientific community. When great scientific upheavals occur, they involve the rejection of the current paradigm in favor of a new paradigm that better integrates the facts available and better predicts the behavior of a particular scientific system. In the field of Alzheimer's disease, a recent example of such a paradigm shift has been the apparent rejection of the 'amyloid cascade hypothesis,' promulgated by Hardy and Higgins in 1992 to explain the etiology of Alzheimer's disease, in favor of what has been referred to as the 'oligomer cascade hypothesis'. This paradigm shift has been breathtaking in its rapidity, its pervasiveness in the Alzheimer's disease field, and its adoption in an increasing number of other fields, including those of Parkinson's disease, Huntington's disease, amyotrophic lateral sclerosis, and the prionoses. However, these facts do not mean, a priori, that the experiments extant, and any re-interpretation of them, should be accepted by rote as support for the new paradigm. In the discussion that follows, I consider the foundational studies leading to the oligomer cascade hypothesis and evaluate the current state of the paradigm. I argue here that, more often than not, insufficient rigor has been applied in studies upon which this new paradigm has been based. Confusion, rather than clarity, has resulted. If the field is to make progress forward using as its paradigmatic basis amyloid β-protein oligomerization, then an epistemological re-evaluation of the amyloid β-protein oligomer system is required.

Protein-protein alternative binding modes do not overlap

Proteins often bind other proteins in more than one way. Thus alternative binding modes is an essential feature of protein interactions. Such binding modes may be detected by X-ray crystallography and thus reflected in Protein Data Bank. The alternative binding is often observed not for the protein itself but for its structural homolog. The results of this study based on the analysis of a comprehensive set of co-crystallized protein-protein complexes show that the alternative binding modes generally do not overlap, but are spatially separated. This effect is based on molecular recognition characteristics of the protein structures. The results are also in excellent agreement with the intermolecular energy funnel size estimates obtained previously by an independent methodology. The results provide an important insight into the principles of protein association, as well as potential guidelines for modeling of protein complexes and the design of protein interfaces.

Stabilization of a protein conferred by an increase in folded state entropy

Entropic stabilization of native protein structures typically relies on strategies that serve to decrease the entropy of the unfolded state. Here we report, using a combination of experimental and computational approaches, on enhanced thermodynamic stability conferred by an increase in the configurational entropy of the folded state. The enhanced stability is observed upon modifications of a loop region in the enzyme acylphosphatase and is achieved despite significant enthalpy losses. The modifications that lead to increased stability, as well as those that result in destabilization, however, strongly compromise enzymatic activity, rationalizing the preservation of the native loop structure even though it does not provide the protein with maximal stability or kinetic foldability.

A method to predict edge strands in beta-sheets from protein sequences

There is a need for rules allowing three-dimensional structure information to be derived from protein sequences. In this work, consideration of an elementary protein folding step allows protein sub-sequences which optimize folding to be derived for any given protein sequence. Classical mechanics applied to this system and the energy conservation law during the elementary folding step yields an equation whose solutions are taken over the field of rational numbers. This formalism is applied to beta-sheets containing two edge strands and at least two central strands. The number of protein sub-sequences optimized for folding per amino acid in beta-strands is shown in particular to predict edge strands from protein sequences. Topological information on beta-strands and loops connecting them is derived for protein sequences with a prediction accuracy of 75%. The statistical significance of the finding is given. Applications in protein structure prediction are envisioned such as for the quality assessment of protein structure models.

Exploiting conformational ensembles in modeling protein-protein interactions on the proteome scale

Cellular functions are performed through protein-protein interactions; therefore, identification of these interactions is crucial for understanding biological processes. Recent studies suggest that knowledge-based approaches are more useful than "blind" docking for modeling at large scales. However, a caveat of knowledge-based approaches is that they treat molecules as rigid structures. The Protein Data Bank (PDB) offers a wealth of conformations. Here, we exploited an ensemble of the conformations in predictions by a knowledge-based method, PRISM. We tested "difficult" cases in a docking-benchmark data set, where the unbound and bound protein forms are structurally different. Considering alternative conformations for each protein, the percentage of successfully predicted interactions increased from ~26 to 66%, and 57% of the interactions were successfully predicted in an "unbiased" scenario, in which data related to the bound forms were not utilized. If the appropriate conformation, or relevant template interface, is unavailable in the PDB, PRISM could not predict the interaction successfully. The pace of the growth of the PDB promises a rapid increase of ensemble conformations emphasizing the merit of such knowledge-based ensemble strategies for higher success rates in protein-protein interaction predictions on an interactome scale. We constructed the structural network of ERK interacting proteins as a case study.

Ensemble-based characterization of unbound and bound states on protein energy landscape

Physicochemical description of numerous cell processes is fundamentally based on the energy landscapes of protein molecules involved. Although the whole energy landscape is difficult to reconstruct, increased attention to particular targets has provided enough structures for mapping functionally important subspaces associated with the unbound and bound protein structures. The subspace mapping produces a discrete representation of the landscape, further called energy spectrum. We compiled and characterized ensembles of bound and unbound conformations of six small proteins and explored their spectra in implicit solvent. First, the analysis of the unbound-to-bound changes points to conformational selection as the binding mechanism for four proteins. Second, results show that bound and unbound spectra often significantly overlap. Moreover, the larger the overlap the smaller the root mean square deviation (RMSD) between the bound and unbound conformational ensembles. Third, the center of the unbound spectrum has a higher energy than the center of the corresponding bound spectrum of the dimeric and multimeric states for most of the proteins. This suggests that the unbound states often have larger entropy than the bound states. Fourth, the exhaustively long minimization, making small intrarotamer adjustments (all-atom RMSD ≤ 0.7 Å), dramatically reduces the distance between the centers of the bound and unbound spectra as well as the spectra extent. It condenses unbound and bound energy levels into a thin layer at the bottom of the energy landscape with the energy spacing that varies between 0.8-4.6 and 3.5-10.5 kcal/mol for the unbound and bound states correspondingly. Finally, the analysis of protein energy fluctuations showed that protein vibrations itself can excite the interstate transitions, including the unbound-to-bound ones.

Understanding the folding-function tradeoff in proteins

When an amino-acid sequence cannot be optimized for both folding and function, folding can get compromised in favor of function. To understand this tradeoff better, we devise a novel method for extracting the "function-less" folding-motif of a protein fold from a set of structurally similar but functionally diverse proteins. We then obtain the β-trefoil folding-motif, and study its folding using structure-based models and molecular dynamics simulations. CompariA protein sequence serves two purpson with the folding of wild-type β-trefoil proteins shows that function affects folding in two ways: In the slower folding interleukin-1β, binding sites make the fold more complex, increase contact order and slow folding. In the faster folding hisactophilin, residues which could have been part of the folding-motif are used for function. This reduces the density of native contacts in functional regions and increases folding rate. The folding-motif helps identify subtle structural deviations which perturb folding. These may then be used for functional annotation. Further, the folding-motif could potentially be used as a first step in the sequence design of function-less scaffold proteins. Desired function can then be engineered into these scaffolds.

Membrane protein native state discrimination by implicit membrane models

Four implicit membrane models [IMM1, generalized Born (GB)-surface area-implicit membrane (GBSAIM), GB with a simple switching (GBSW), and heterogeneous dielectric GB (HDGB)] were tested for their ability to discriminate the native conformation of five membrane proteins from 450 decoys generated by the Rosetta-Membrane program. The energy ranking of the native state and Z-scores were used to assess the performance of the models. The effect of membrane thickness was examined and was found to be substantial. Quite satisfactory discrimination was achieved with the all-atom IMM1 and GBSW models at 25.4 Å thickness and with the HDGB model at 28.5 Å thickness. The energy components by themselves were not discriminative. Both van der Waals and electrostatic interactions contributed to native state discrimination, to a different extent in each model. Computational efficiency of the models decreased in the order: extended-atom IMM1 > all-atom IMM1 > GBSAIM > GBSW > HDGB. These results encourage the further development and use of implicit membrane models for membrane protein structure prediction.

Protein dynamics and motions in relation to their functions: several case studies and the underlying mechanisms

Proteins are dynamic entities in cellular solution with functions governed essentially by their dynamic personalities. We review several dynamics studies on serine protease proteinase K and HIV-1 gp120 envelope glycoprotein to demonstrate the importance of investigating the dynamic behaviors and molecular motions for a complete understanding of their structure–function relationships. Using computer simulations and essential dynamic (ED) analysis approaches, the dynamics data obtained revealed that: (i) proteinase K has highly flexible substrate-binding site, thus supporting the induced-fit or conformational selection mechanism of substrate binding; (ii) Ca²⁺ removal from proteinase K increases the global conformational flexibility, decreases the local flexibility of substrate-binding region, and does not influence the thermal motion of catalytic triad, thus explaining the experimentally determined decreased thermal stability, reduced substrate affinity, and almost unchanged catalytic activity upon Ca²⁺ removal; (iii) substrate binding affects the large concerted motions of proteinase K, and the resulting dynamic pocket can be connected to substrate binding, orientation, and product release; (iv) amino acid mutations 375 S/W and 423 I/P of HIV-1 gp120 have distinct effects on molecular motions of gp120, facilitating 375 S/W mutant to assume the CD4-bound conformation, while 423 I/P mutant to prefer for CD4-unliganded state. The mechanisms underlying protein dynamics and protein–ligand binding, including the concept of the free energy landscape (FEL) of the protein–solvent system, how the ruggedness and variability of FEL determine protein's dynamics, and how the three ligand-binding models, the lock-and-key, induced-fit, and conformational selection are rationalized based on the FEL theory are discussed in depth.

Where soft matter meets living matter--protein structure, stability, and folding in the cell

A protein is a biopolymer that self-assembles through the process of protein folding. A cell is a crowded space where the surrounding macromolecules of a protein can limit the number of ways of folding. These crowding macromolecules can also affect the shape and the size of a physically malleable, or 'soft, squishy', protein with regulatory purposes. In this review, we focus on the in silico approaches of coarse-grained molecular simulations that enable the investigation of protein folding in a cell-like environment. When these simulation results were compared with experimentally measured properties of a protein, such joint effort has yielded new ideas on the specific function of a protein in cells. We also highlighted the recent developments of computer modeling and simulations that encompass the importance of the shape of a macromolecule, the interactions between macromolecules, and the hydrodynamic interactions on the kinetics and thermodynamics of a protein in a high concentration of protein solution and in cytoplasmic environments.

Folding of a heterogeneous β-hairpin peptide from temperature-jump 2D IR spectroscopy

We provide a time- and structure-resolved characterization of the folding of the heterogeneous β-hairpin peptide Tryptophan Zipper 2 (Trpzip2) using 2D IR spectroscopy. The amide I' vibrations of three Trpzip2 isotopologues are used as a local probe of the midstrand contacts, β-turn, and overall β-sheet content. Our experiments distinguish between a folded state with a type I' β-turn and a misfolded state with a bulged turn, providing evidence for distinct conformations of the peptide backbone. Transient 2D IR spectroscopy at 45 °C following a laser temperature jump tracks the nanosecond and microsecond kinetics of unfolding and the exchange between conformers. Hydrogen bonds to the peptide backbone are loosened rapidly compared with the 5-ns temperature jump. Subsequently, all relaxation kinetics are characterized by an observed 1.2 ± 0.2-μs exponential. Our time-dependent 2D IR spectra are explained in terms of folding of either native or nonnative contacts from a common compact disordered state. Conversion from the disordered state to the folded state is consistent with a zip-out folding mechanism.

Limiting assumptions in molecular modeling: electrostatics

Molecular mechanics attempts to represent intermolecular interactions in terms of classical physics. Initial efforts assumed a point charge located at the atom center and coulombic interactions. It is been recognized over multiple decades that simply representing electrostatics with a charge on each atom failed to reproduce the electrostatic potential surrounding a molecule as estimated by quantum mechanics. Molecular orbitals are not spherically symmetrical, an implicit assumption of monopole electrostatics. This perspective reviews recent evidence that requires use of multipole electrostatics and polarizability in molecular modeling.

Unusual biophysics of intrinsically disordered proteins

Research of a past decade and a half leaves no doubt that complete understanding of protein functionality requires close consideration of the fact that many functional proteins do not have well-folded structures. These intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered protein regions (IDPRs) are highly abundant in nature and play a number of crucial roles in a living cell. Their functions, which are typically associated with a wide range of intermolecular interactions where IDPs possess remarkable binding promiscuity, complement functional repertoire of ordered proteins. All this requires a close attention to the peculiarities of biophysics of these proteins. In this review, some key biophysical features of IDPs are covered. In addition to the peculiar sequence characteristics of IDPs these biophysical features include sequential, structural, and spatiotemporal heterogeneity of IDPs; their rough and relatively flat energy landscapes; their ability to undergo both induced folding and induced unfolding; the ability to interact specifically with structurally unrelated partners; the ability to gain different structures at binding to different partners; and the ability to keep essential amount of disorder even in the bound form. IDPs are also characterized by the "turned-out" response to the changes in their environment, where they gain some structure under conditions resulting in denaturation or even unfolding of ordered proteins. It is proposed that the heterogeneous spatiotemporal structure of IDPs/IDPRs can be described as a set of foldons, inducible foldons, semi-foldons, non-foldons, and unfoldons. They may lose their function when folded, and activation of some IDPs is associated with the awaking of the dormant disorder. It is possible that IDPs represent the "edge of chaos" systems which operate in a region between order and complete randomness or chaos, where the complexity is maximal. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

Hydrophobic enhancement of Dopa-mediated adhesion in a mussel foot protein

Dopa (3,4-dihydroxyphenylalanine) is recognized as a key chemical signature of mussel adhesion and has been adopted into diverse synthetic polymer systems. Dopa's notorious susceptibility to oxidation, however, poses significant challenges to the practical translation of mussel adhesion. Using a surface forces apparatus to investigate the adhesion of mussel foot protein 3 (Mfp3) "slow", a hydrophobic protein variant of the Mfp3 family in the plaque, we have discovered a subtle molecular strategy correlated with hydrophobicity that appears to compensate for Dopa instability. At pH 3, where Dopa is stable, Mfp3 slow, like Mfp3 "fast" adhesion to mica, is directly proportional to the mol % of Dopa present in the protein. At pH of 5.5 and 7.5, however, loss of adhesion in Mfp3 slow was less than half that occurring in Mfp3 fast, purportedly because Dopa in Mfp3 slow is less prone to oxidation. Indeed, cyclic voltammetry showed that the oxidation potential of Dopa in Mfp3 slow is significantly higher than in Mfp3 fast at pH of 7.5. A much greater difference between the two variants was revealed in the interaction energy of two symmetric Mfp3 slow films (E(ad) = -3 mJ/m(2)). This energy corresponds to the energy of protein cohesion which is notable for its reversibility and pH independence. Exploitation of aromatic hydrophobic sequences to protect Dopa against oxidation as well as to mediate hydrophobic and H-bonding interactions between proteins provides new insights for developing effective artificial underwater adhesives.

Levinthal's question revisited, and answered

Attempts to answer the Levinthal question "How proteins fold to give such a unique structure" are discussed. In the first part of this article, we focus on a few reasons as to why the solution to the protein-folding problem (PFP) has been elusive for a very long time. One is a result of the misinterpretation of Anfinsen's Thermodynamic hypothesis which led to the conclusion that the native structure of a protein must be at a global minimum of the Gibbs energy. The second is the result of the adherence to the hydrophobic paradigm, and at the same time ignoring a whole repertoire of hydrophilic effects. It is argued that switching from a target-based to a caused-based approach, and adopting the hydrophilic paradigm leads straightforwardly to a simple answer to Levinthal's question, as well as to a solution of the PFP.

Reducing the dimensionality of the protein-folding search problem

How does a folding protein negotiate a vast, featureless conformational landscape and adopt its native structure in biological real time? Motivated by this search problem, we developed a novel algorithm to compare protein structures. Procedures to identify structural analogs are typically conducted in three-dimensional space: the tertiary structure of a target protein is matched against each candidate in a database of structures, and goodness of fit is evaluated by a distance-based measure, such as the root-mean-square distance between target and candidate. This is an expensive approach because three-dimensional space is complex. Here, we transform the problem into a simpler one-dimensional procedure. Specifically, we identify and label the 11 most populated residue basins in a database of high-resolution protein structures. Using this 11-letter alphabet, any protein's three-dimensional structure can be transformed into a one-dimensional string by mapping each residue onto its corresponding basin. Similarity between the resultant basin strings can then be evaluated by conventional sequence-based comparison. The disorder → order folding transition is abridged on both sides. At the onset, folding conditions necessitate formation of hydrogen-bonded scaffold elements on which proteins are assembled, severely restricting the magnitude of accessible conformational space. Near the end, chain topology is established prior to emergence of the close-packed native state. At this latter stage of folding, the chain remains molten, and residues populate natural basins that are approximated by the 11 basins derived here. In essence, our algorithm reduces the protein-folding search problem to mapping the amino acid sequence onto a restricted basin string.

Molecular dynamics simulation of thionated hen egg white lysozyme

Understanding of the driving forces of protein folding is a complex challenge because different types of interactions play a varying role. To investigate the role of hydrogen bonding involving the backbone, the effect of thio substitutions in a protein, hen egg white lysozyme (HEWL), was investigated through molecular dynamics simulations of native as well as partly (only residues in loops) and fully thionated HEWL using the GROMOS 54A7 force field. The results of the three simulations show that the structural properties of fully thionated HEWL clearly differ from those of the native protein, while for partly thionated HEWL they only changed slightly compared with native HEWL. The analysis of the torsional-angle distributions and hydrogen bonds in the backbone suggests that the α-helical segments of native HEWL tend to show a propensity to convert to 3(10)-helical geometry in fully thionated HEWL. A comparison of the simulated quantities with experimental NMR data such as nuclear overhauser effect (NOE) atom-atom distance bounds and (3)J((H)(N)(H)(α))-couplings measured for native HEWL illustrates that the information content of these quantities with respect to the structural changes induced by thionation of the protein backbone is rather limited.

Use of SANS and biophysical techniques to reveal subtle conformational differences between native apo-calmodulin and its unfolded states

Apo-calmodulin, a small, mainly α, soluble protein is a calcium-dependent protein activator. It is made of two N- and C-terminal domains having a sequence homology of 70%, an identical folding but different stabilities, and is thus an interesting system for unfolding studies. The use of small angle neutron scattering (SANS) and other biophysical techniques has permitted to reveal conformational difference between native and thermal denatured states of apo-calmodulin. The results show that secondary and tertiary structures of apo-calmodulin evolve in a synchronous way, indicating the absence in the unfolding pathway of molten-globule state sufficiently stable to affect transition curves. From SANS experiments, at 85 °C, apo-calmodulin adopts a polymer chain conformation with some residual local structures. After cooling down, apo-calmodulin recovers a compact state, with a secondary structure close to the native one but with a higher radius of gyration and a different tyrosine environment. In fact on a timescale of few minutes, heat denaturation of apo-calmodulin is partially reversible, but on a time scale of hours (for SANS experiments), the long exposure to heat may lead to a non-reversibility due to some chemical perturbation of the protein. In fact, from Mass Spectrometry measurements, we got evidence of dehydration and deamidation of heated apo-calmodulin.

UV resonance Raman studies of the NaClO4 dependence of poly-L-lysine conformation and hydrogen exchange kinetics

We used 204 nm excitation UV Resonance Raman (UVRR) spectroscopy to examine the effects of NaClO(4) on the conformation of poly-L-lysine (PLL). The presence of NaClO(4) induces the formation of α-helix, π-helix/bulge, and turn conformations. The dependence of the AmIII(3) frequency on the peptide Ψ Ramachandran angle allows us to experimentally determine the conformational population distributions and the energy landscape of PLL along the Ramachandran Ψ angle. We also used UVRR to measure the NaClO(4) concentration dependence of PLL amide hydrogen exchange kinetics. Exchange rates were determined by fitting the D(2)O exchanging PLL UVRR AmII' band time evolution. Hydrogen exchange is slowed at high NaClO(4) concentrations. The PLL AmII' band exchange kinetics at 0.0, 0.2, and 0.35 M NaClO(4) can be fit by single exponentials, but the AmII' band kinetics of PLL at 0.8 M NaClO(4) requires a double exponential fit. The exchange rates for the extended conformations were monitored by measuring the C(α)-H band kinetics. These kinetics are identical to those of the AmII' band until 0.8 M NaClO(4) whereupon the extended conformation exchange becomes clearly faster than that of the α-helix-like conformations. Our results indicate that ClO(4)(-) binds to the PLL backbone to protect it from OH(-) exchange catalysis. In addition, ClO(4)(-) binding also slows the conformational exchange between the extended and α-helix-like conformations, probably by increasing the activation barriers for conformational interchanges.

Predicting three-dimensional structures of transmembrane domains of β-barrel membrane proteins

β-Barrel membrane proteins are found in the outer membrane of gram-negative bacteria, mitochondria, and chloroplasts. They are important for pore formation, membrane anchoring, and enzyme activity. These proteins are also often responsible for bacterial virulence. Due to difficulties in experimental structure determination, they are sparsely represented in the protein structure databank. We have developed a computational method for predicting structures of the transmembrane (TM) domains of β-barrel membrane proteins. Based on physical principles, our method can predict structures of the TM domain of β-barrel membrane proteins of novel topology, including those from eukaryotic mitochondria. Our method is based on a model of physical interactions, a discrete conformational state space, an empirical potential function, as well as a model to account for interstrand loop entropy. We are able to construct three-dimensional atomic structure of the TM domains from sequences for a set of 23 nonhomologous proteins (resolution 1.8-3.0 Å). The median rmsd of TM domains containing 75-222 residues between predicted and measured structures is 3.9 Å for main chain atoms. In addition, stability determinants and protein-protein interaction sites can be predicted. Such predictions on eukaryotic mitochondria outer membrane protein Tom40 and VDAC are confirmed by independent mutagenesis and chemical cross-linking studies. These results suggest that our model captures key components of the organization principles of β-barrel membrane protein assembly.

The different ways through which specificity works in orthosteric and allosteric drugs

Currently, there are two types of drugs on the market: orthosteric, which bind at the active site; and allosteric, which bind elsewhere on the protein surface, and allosterically change the conformation of the protein binding site. In this perspective we argue that the different mechanisms through which the two drug types affect protein activity and their potential pitfalls call for different considerations in drug design. The key problem facing orthosteric drugs is side effects which can occur by drug binding to homologous proteins sharing a similar binding site. Hence, orthosteric drugs should have very high affinity to the target; this would allow a low dosage to selectively achieve the goal of target-only binding. By contrast, allosteric drugs work by shifting the free energy landscape. Their binding to the protein surface perturbs the protein surface atoms, and the perturbation propagates like waves, finally reaching the binding site. Effective drugs should have atoms in good contact with the 'right' protein atoms; that is, the contacts should elicit propagation waves optimally reaching the protein binding site target. While affinity is important, the design should consider the protein conformational ensemble and the preferred propagation states. We provide examples from functional in vivo scenarios for both types of cases, and suggest how high potency can be achieved in allosteric drug development.

Exploration of the relationship between topology and designability of conformations

Protein structures are evolutionarily more conserved than sequences, and sequences with very low sequence identity frequently share the same fold. This leads to the concept of protein designability. Some folds are more designable and lots of sequences can assume that fold. Elucidating the relationship between protein sequence and the three-dimensional (3D) structure that the sequence folds into is an important problem in computational structural biology. Lattice models have been utilized in numerous studies to model protein folds and predict the designability of certain folds. In this study, all possible compact conformations within a set of two-dimensional and 3D lattice spaces are explored. Complementary interaction graphs are then generated for each conformation and are described using a set of graph features. The full HP sequence space for each lattice model is generated and contact energies are calculated by threading each sequence onto all the possible conformations. Unique conformation giving minimum energy is identified for each sequence and the number of sequences folding to each conformation (designability) is obtained. Machine learning algorithms are used to predict the designability of each conformation. We find that the highly designable structures can be distinguished from other non-designable conformations based on certain graphical geometric features of the interactions. This finding confirms the fact that the topology of a conformation is an important determinant of the extent of its designability and suggests that the interactions themselves are important for determining the designability.

Gene-specific transcription activation via long-range allosteric shape-shifting

How is specificity transmitted over long distances at the molecular level? REs (regulatory elements) are often far from transcription start sites. In the present review we discuss possible mechanisms to explain how information from specific REs is conveyed to the basal transcription machinery through TFs (transcription factors) and the Mediator complex. We hypothesize that this occurs through allosteric pathways: binding of a TF to a RE results in changes in the AD (activation domain) of the TF, which binds to Mediator and alters the distribution of the Mediator conformations, thereby affecting transcription initiation/activation. We argue that Mediator is formed by highly disordered proteins with large densely packed interfaces that make efficient long-range signal propagation possible. We suggest two possible general mechanisms for Mediator action: one in which Mediator influences PIC (pre-initiation complex) assembly and transcription initiation, and another in which Mediator exerts its effect on the already assembled but stalled transcription complex. We summarize (i) relevant information from the literature about Mediator composition, organization and structure; (ii) Mediator interaction partners and their effect on Mediator conformation, function and correlation to the RNA Pol II (polymerase II) CTD (C-terminal domain) phosphorylation; and (iii) propose that different allosteric signal propagation pathways in Mediator relate to PIC assembly and polymerase activation of the stalled transcription complex. The emerging picture provides for the first time a mechanistic view of allosteric signalling from the RE sequence to transcription activation, and an insight into how gene specificity and signal transmission can take place in transcription initiation.