Structure, Folding, and Energetics of Cooperative Interactions between the β-Strands of a de Novo Designed Three-Stranded Antiparallel β-Sheet Peptide

Tau Inclusions in Alzheimer's, Chronic Traumatic Encephalopathy and Pick's Disease. A Speculation on How Differences in Backbone Polarization Underlie Divergent Pathways of Tau Aggregation

Tau-related dementias appear to involve specific to each disease aggregation pathways and morphologies of filamentous tau assemblies. To understand etiology of these differences, here we elucidate molecular mechanism of formation of tau PHFs based on the PMO theory of misfolding and aggregation of pleiomorphic proteins associated with neurodegenerative diseases. In this model, fibrillization of tau is initiated by the coupled binding and folding of the MTB domains that yields antiparallel homodimers, in analogy to folding of split inteins. The free energy of binding is minimized when the antiparallel alignment brings about backbone-backbone H-bonding between the MTBD segments of similar "strand" propensities. To assess these propensities, a function of the NMR shielding tensors of the Cα atoms is introduced as the folding potential function FP i ; the Cα tensors are obtained by the quantum mechanical modeling of protein secondary structure (GIAO//B3LYP/D95**). The calculated FP i plots show that the "strand" propensities of the MBTD segments, and hence the homodimer's register, can be affected by the relatively small changes in the environment's pH, as a result of protonation of MBTD's conserved histidines. The assembly of the antiparallel tau dimers into granular aggregates and their subsequent conversion into the parallel cross-β structure of paired helical filaments is expected to follow the same path as the previously described fibrillization of Aβ. Consequently, the core structure of the nascent tau fibril is determined by the register of the tau homodimer. This model accounts for the reported differences in (i) fibril-core structure of in vivo and in vitro filaments, (ii) cross-seeding of isoforms, (iii) effects of reducing/non-reducing conditions, (iv) effects of PHF6 mutations, and (v) homologs' aggregation properties. The proposed model also suggests that in contrast to Alzheimer's and chronic traumatic encephalopathy disease, the assembly of tau prions in Pick's disease would be facilitated by a moderate drop in pH that accompanies e.g., transit in the endosomal system, inflammation response or an ischemic injury.

Isotopically Site-Selected Dynamics of a Three-Stranded β-Sheet Peptide Detected with Temperature-Jump Infrared-Spectroscopy

Infrared detected temperature-jump (T-jump) spectroscopy and site-specific isotopic labeling were applied to study a model three-stranded β-sheet peptide with the goal of individually probing the dynamics of strand and turn structural elements. This peptide had two ^DPro-Gly (pG) turn sequences to stabilize the two component hairpins, which were labeled with ¹³C═O on each of the Gly residues to resolve them spectroscopically. Labeling the second turn on the amide preceding the ^DPro (Xxx-^DPro amide) provided an alternate turn label as a control. Placing ¹³C═O labels on specific in-strand residues gave shifted modes that overlap the Xxx-^DPro amide I' modes. Their impact could be separated from the turn dynamics by a novel difference transient analysis approach. Fourier-transform infrared spectra were modeled with density functional theory-computations which showed the local, isotope-selected vibrations were effectively uncoupled from the other amide I modes. Our T-jump dynamics results, combined with nuclear magnetic resonance structures and equilibrium spectral measurements, showed the first turn to be most stable and best formed with the slowest dynamics, whereas the second turn and first strand (N-terminus) had similar dynamics, and the third strand (C-terminus) had the fastest dynamics and was the least structured. The relative dynamics of the strands, Xxx-^DPro amides, and ¹³C-labeled Gly residues on the turns also qualitatively corresponded to molecular dynamics (MD) simulations of turn and strand fluctuations. MD trajectories indicated the turns to be bistable, with the first turn being Type I' and the second turn flipping from I' to II'. The differences in relaxation times for each turn and the separate strands revealed that the folding process of this turn-stabilized β-sheet structure proceeds in a multistep process.

Conformational stabilization of a β-hairpin through a triazole–tryptophan interaction

Triazole and indole rings stabilize a β-hairpin conformation through an aromatic–aromatic interaction.

Protein folding, misfolding and aggregation: The importance of two-electron stabilizing interactions

Proteins associated with neurodegenerative diseases are highly pleiomorphic and may adopt an all-α-helical fold in one environment, assemble into all-β-sheet or collapse into a coil in another, and rapidly polymerize in yet another one via divergent aggregation pathways that yield broad diversity of aggregates’ morphology. A thorough understanding of this behaviour may be necessary to develop a treatment for Alzheimer’s and related disorders. Unfortunately, our present comprehension of folding and misfolding is limited for want of a physicochemical theory of protein secondary and tertiary structure. Here we demonstrate that electronic configuration and hyperconjugation of the peptide amide bonds ought to be taken into account to advance such a theory. To capture the effect of polarization of peptide linkages on conformational and H-bonding propensity of the polypeptide backbone, we introduce a function of shielding tensors of the C^α atoms. Carrying no information about side chain-side chain interactions, this function nonetheless identifies basic features of the secondary and tertiary structure, establishes sequence correlates of the metamorphic and pH-driven equilibria, relates binding affinities and folding rate constants to secondary structure preferences, and manifests common patterns of backbone density distribution in amyloidogenic regions of Alzheimer’s amyloid β and tau, Parkinson’s α-synuclein and prions. Based on those findings, a split-intein like mechanism of molecular recognition is proposed to underlie dimerization of Aβ, tau, αS and PrP^C, and divergent pathways for subsequent association of dimers are outlined; a related mechanism is proposed to underlie formation of PrP^Sc fibrils. The model does account for: (i) structural features of paranuclei, off-pathway oligomers, non-fibrillar aggregates and fibrils; (ii) effects of incubation conditions, point mutations, isoform lengths, small-molecule assembly modulators and chirality of solid-liquid interface on the rate and morphology of aggregation; (iii) fibril-surface catalysis of secondary nucleation; and (iv) self-propagation of infectious strains of mammalian prions.

Direct Determination of Site-Specific Noncovalent Interaction Strengths of Proteins from NMR-Derived Fast Side Chain Motional Parameters

A novel approach to accurately determine residue-specific noncovalent interaction strengths (ξ) of proteins from NMR-measured fast side chain motional parameters (O_axis²) is presented. By probing the environmental sensitivity of side chain conformational energy surfaces of individual residues of a diverse set of proteins, the microscopic connections between ξ, O_axis², conformational entropy (S_conf), conformational barriers, and rotamer stabilities established here are found to be universal among proteins. The results reveal that side chain flexibility and conformational entropy of each residue decrease with increasing ξ and that for each residue type there exists a critical range of ξ, determined primarily by the mean side chain conformational barriers, within which flexibility of any residue can be reversibly tuned from highly flexible (with O_axis² ∼ 0) to highly restricted (with O_axis² ∼ 1) by increasing ξ by ∼3 kcal/mol. Beyond this critical range of ξ, both side chain flexibility and conformational entropy are insensitive to ξ. The interrelationships between conformational dynamics, conformational entropy, and noncovalent interactions of protein side chains established here open up new avenues to probe perturbation-induced (for example, ligand-binding, temperature, pressure) changes in fast side chain dynamics and thermodynamics of proteins by comparing their conformational energy surfaces in the native and perturbed states.

Automated protein design: Landmarks and operational principles

Protein design has an eventful history spanning over three decades, with handful of success stories reported, and numerous failures not reported. Design practices have benefited tremendously from improvements in computer hardware and advances in scientific algorithms. Though protein folding problem still remains unsolved, the possibility of having multiple sequence solutions for a single fold makes protein design a more tractable problem than protein folding. One of the most significant advancement in this area is the implementation of automated design algorithms on pre-defined templates or completely new folds, optimized through deterministic and heuristic search algorithms. This progress report provides a succinct presentation of important landmarks in automated design attempts, followed by brief account of operational principles in automated design methods.

Density Functionals for Noncovalent Interaction Energies of Biological Importance

Forty density functionals and one wavefunction method are assessed against a recently published database of accurate noncovalent interaction energies of biological importance. The comparison shows that two newly developed density functional theory (DFT) methods, PWB6K and M05-2X, give the best performance for this benchmark database of 22 noncovalent complexes, including both hydrogen-bonding and dispersion-dominated complexes. In contrast, the more popular B3LYP and PBEh functionals fail to describe the interactions in the dispersion-dominated complexes. The local spin density approximation and BHandH functionals give good performance for dispersion-dominated interactions at the expense of a large error for hydrogen bonding. PWB6K and M05-2X constitute a new generation of DFT methods based on simultaneously optimized exchange and correlation functionals that include kinetic energy density in both the exchange and correlation functional, and the present study confirms that they have greatly improved performance for noncovalent interactions as compared to previous DFT methods. We interpret this as being due to an improved treatment of medium-range correlation effects by the exchange-correlation functional. We recommend the PWB6K and M05-2X methods for investigating large biological systems and soft materials.

Design of monomeric water-soluble β-hairpin and β-sheet peptides

Since the first report in 1993 (JACS 115, 5887-5888) of a peptide able to form a monomeric β-hairpin structure in aqueous solution, the design of peptides forming either β-hairpins (two-stranded antiparallel β-sheets) or three-stranded antiparallel β-sheets has become a field of growing interest and activity. These studies have yielded great insights into the principles governing the stability and folding of β-hairpins and antiparallel β-sheets. This chapter provides an overview of the reported β-hairpin/β-sheet peptides focussed on the applied design criteria, reviews briefly the factors contributing to β-hairpin/β-sheet stability, and describes a protocol for the de novo design of β-sheet-forming peptides based on them. Guidelines to select appropriate turn and strand residues and to avoid self-association are provided. The methods employed to check the success of new designed peptides are also summarized. Since NMR is the best technique to that end, NOEs and chemical shifts characteristic of β-hairpins and three-stranded antiparallel β-sheets are given.

Structure of a simplified β-hairpin and its ATP complex

The capacity of three designed duodecamer peptides with the low diversity sequence: H1ϕ2I3K4I5D6G7K8ϕ9I10K11H12 where ϕ is His, Phe or Trp, to adopt a β-hairpin conformation was studied using NMR spectroscopy. Whereas KIAβH, the variant with His at positions two and nine, is disordered, KIAβF, the peptide with Phe at these positions, adopts a small population of β-hairpin. A high population of β-hairpin structure was detected for KIAβW, the variant with Trp. Utilizing NMR data, the structure of KIAβW was solved and it reveals a β-hairpin stabilized by hydrophobic interactions between Ile residues on one face and Trp-Trp and cation-π interactions on the opposite face. Upon adding ATP, these peptides show chemical shift changes indicative of ATP binding. The binding of ATP to KIAβW shows a KD ≈ 20 μM at pH 5, 5 °C and has a 1:1 stoichiometry. The KIAβW-ATP complex was determined using NMR spectroscopy and reveals the adenine ring sandwiched between the two Trp indole rings and that ATP binding induces important conformational changes in His1, Trp2, Lys4, Trp9 and Lys11 in the β-hairpin. The implications of these results for the hypothetic presence of β-hairpins and amyloids alongside RNAs on the prebiotic Earth are discussed.

A model for the controlled assembly of semiconductor peptides

The self-assembly of small molecules provides a potentially powerful method to create functional nanomaterials for many applications ranging from optoelectronics to oncology. However, the design of well-defined nanostructures via molecular assembly is a highly empirical process, which severely hampers efforts to create functional nanostructures using this method. In this review, we describe a simple strategy to control the assembly of functionalized peptides by balancing attractive hydrophobic effects that drive assembly with opposing electrostatic repulsions. Extended π-π contacts are created in the nanostructures when assembly is driven by π-stacking interactions among chromophores that are appended to the peptide. The formation of insoluble β-sheet aggregates are mitigated by incorporating charged side-chains capable of attenuating the assembly process. Although the application of this approach to the assembly of organic semiconductors is described, we expect this strategy to be effective for many other functional organic materials.

Constrained peptides as miniature protein structures

This paper discusses the recent developments of protein engineering using both covalent and noncovalent bonds to constrain peptides, forcing them into designed protein secondary structures. These constrained peptides subsequently can be used as peptidomimetics for biological functions such as regulations of protein-protein interactions.

Modeling Structural Flexibility of Proteins with Go-Models

Structure-based models are an efficient tool for folding studies of proteins since by construction their energy landscape is only minimal frustrated. However, their intrinsic drawback is a lack of structural flexibility as usually only one target structure is employed to construct the potentials. Hence, a Go-model may not capture differences in mutation-induced protein dynamics, if - as in the case of the disease-related A629P mutant of the Menkes protein ATP7A - the structural differences between mutant and wild type are small. In this work, we introduced three implementations of Go-models that take into account the flexibility of proteins in the NMR ensemble. Comparing the wild type and the mutant A629P of the 75-residue large 6th domain Menkes protein, we find that these new Go-potentials lead to broader distributions than Go-models relying on a single member of the NMR-ensemble. This allows us to detect the transient unfolding of a loosely formed β1β4-sheet in the mutant protein. Our results are consistent with previous simulations using physical force field and an explicit solvent, and suggests a mechanism by which these mutations cause Menkes disease. In addition, the improved Go-models suggest differences in the folding pathway between wild type and mutant, an observation that was not accessible to simulations of this 75-residue protein with a physical all-atom force field and explicit solvent.

β-Hairpin stabilization through an interstrand triazole bridge

β-Hairpin peptides were conformationally stabilized through a 1,4 disubstituted 1,2,3-triazole interstrand linkage. A NMR conformational analysis revealed that the β-hairpin content depends on the number and position of substituent methylene units of the 1,2,3-triazole ring. These results will allow the design of metabolically stable peptidomimetic analogs of bioactive β-hairpin peptides.

Contributions of aromatic pairs to the folding and stability of long-lived human γD-crystallin

Human γD-crystallin (HγD-Crys) is a highly stable protein that remains folded in the eye lens for the majority of an individual's lifetime. HγD-Crys exhibits two homologous crystallin domains, each containing two Greek key motifs with eight β-strands. Six aromatic pairs (four Tyr/Tyr, one Tyr/Phe and one Phe/Phe) are present in the β-hairpin sequences of the Greek keys. Ultraviolet damage to the aromatic residues in lens crystallins may contribute to the genesis of cataract. Mutant proteins with these aromatic residues substituted with alanines were constructed and expressed in E. coli. All mutant proteins except F115A and F117A had lower thermal stability than the WT protein. In equilibrium experiments in guanidine hydrochloride (GuHCl), all mutant proteins had lower thermodynamic stability than the WT protein. N-terminal domain (N-td) substitutions shifted the N-td transition to lower GuHCl concentration, but the C-terminal domain (C-td) transition remained unaffected. C-td substitutions led to a more cooperative unfolding/refolding process, with both the N-td and C-td transitions shifted to lower GuHCl concentration. The aromatic pairs conserved for each Greek key motif (Greek key pairs) had larger contributions to both thermal stability and thermodynamic stability than the other pairs. Aromatic-aromatic interaction was estimated as 1.5-2.0 kcal/mol. In kinetic experiments, N-td substitutions accelerated the early phase of unfolding, while C-td substitutions accelerated the late phase, suggesting independent domain unfolding. Only substitutions of the second Greek key pair of each crystallin domain slowed refolding. The second Greek keys may provide nucleation sites during the folding of the double-Greek-key crystallin domains.

Redesign of a WW domain peptide for selective recognition of single-stranded DNA

A β-sheet miniprotein based on the FBP11 WW1 domain sequence has been redesigned for the molecular recognition of ssDNA. A previous report showed that a β-hairpin peptide dimer, (WKWK)(2), binds ssDNA with low micromolar affinity but with little selectivity over duplex DNA. This report extends those studies to a three-stranded β-sheet miniprotein designed to mimic the OB-fold. The new peptide binds ssDNA with low micromolar affinity and shows about 10-fold selectivity for ssDNA over duplex DNA. The redesigned peptide no longer binds its native ligand, the polyproline helix, confirming that the peptide has been redesigned for the function of binding ssDNA. Structural studies provide evidence that this peptide consists of a well-structured β-hairpin made of strands 2 and 3 with a less structured first strand that provides affinity for ssDNA but does not improve the stability of the full peptide. These studies provide insight into protein-DNA interactions as well as a novel example of protein redesign.

Influence of the solvent on the self-assembly of a modified amyloid beta peptide fragment. II. NMR and computer simulation investigation

The conformation of a model peptide AAKLVFF based on a fragment of the amyloid beta peptide Abeta16-20, KLVFF, is investigated in methanol and water via solution NMR experiments and molecular dynamics computer simulations. In previous work, we have shown that AAKLVFF forms peptide nanotubes in methanol and twisted fibrils in water. Chemical shift measurements were used to investigate the solubility of the peptide as a function of concentration in methanol and water. This enabled the determination of critical aggregation concentrations. The solubility was lower in water. In dilute solution, diffusion coefficients revealed the presence of intermediate aggregates in concentrated solution, coexisting with NMR-silent larger aggregates, presumed to be beta-sheets. In water, diffusion coefficients did not change appreciably with concentration, indicating the presence mainly of monomers, coexisting with larger aggregates in more concentrated solution. Concentration-dependent chemical shift measurements indicated a folded conformation for the monomers/intermediate aggregates in dilute methanol, with unfolding at higher concentration. In water, an antiparallel arrangement of strands was indicated by certain ROESY peak correlations. The temperature-dependent solubility of AAKLVFF in methanol was well described by a van't Hoff analysis, providing a solubilization enthalpy and entropy. This pointed to the importance of solvophobic interactions in the self-assembly process. Molecular dynamics simulations constrained by NOE values from NMR suggested disordered reverse turn structures for the monomer, with an antiparallel twisted conformation for dimers. To model the beta-sheet structures formed at higher concentration, possible model arrangements of strands into beta-sheets with parallel and antiparallel configurations and different stacking sequences were used as the basis for MD simulations; two particular arrangements of antiparallel beta-sheets were found to be stable, one being linear and twisted and the other twisted in two directions. These structures were used to simulate circular dichroism spectra. The roles of aromatic stacking interactions and charge transfer effects were also examined. Simulated spectra were found to be similar to those observed experimentally (in water or methanol) which show a maximum at 215 or 218 nm due to pi-pi* interactions, when allowance is made for a 15-18 nm red-shift that may be due to light scattering effects.

Understanding the mechanism of beta-sheet folding from a chemical and biological perspective

Perturbing the structure of the Pin1 WW domain, a 34-residue protein comprised of three beta-strands and two intervening loops has provided significant insight into the structural and energetic basis of beta-sheet folding. We will review our current perspective on how structure acquisition is influenced by the sequence, which determines local conformational propensities and mediates the hydrophobic effect, hydrogen bonding, and analogous intramolecular interactions. We have utilized both traditional site-directed mutagenesis and backbone mutagenesis approaches to alter the primary structure of this beta-sheet protein. Traditional site-directed mutagenesis experiments are excellent for altering side-chain structure, whereas amide-to-ester backbone mutagenesis experiments modify backbone-backbone hydrogen bonding capacity. The transition state structure associated with the folding of the Pin1 WW domain features a partially H-bonded, near-native reverse turn secondary structure in loop 1 that has little influence on thermodynamic stability. The thermodynamic stability of the Pin1 WW domain is largely determined by the formation of a small hydrophobic core and by the formation of desolvated backbone-backbone H-bonds enveloped by this hydrophobic core. Loop 1 engineering to the consensus five-residue beta-bulge-turn found in most WW domains or a four-residue beta-turn found in most beta-hairpins accelerates folding substantially relative to the six-residue turn found in the wild type Pin1 WW domain. Furthermore, the more efficient five- and four-residue reverse turns now contribute to the stability of the three-stranded beta-sheet. These insights have allowed the design of Pin1 WW domains that fold at rates that approach the theoretical speed limit of folding.

Transformation of an alpha-helix peptide into a beta-hairpin induced by addition of a fragment results in creation of a coexisting state

Intrinsic rules of determining the tertiary structure of a protein have been unknown partly because physicochemical factors that contribute to stabilization of a protein structure cannot be represented as a linear combination of local interactions. To clarify the rules on the nonlinear term caused by nonlocal interaction in a protein, we tried to transform a peptide that has a fully helical structure (Target Peptide or TP) into a peptide that has a beta-hairpin structure (Designed Peptide or DP) by adding seven residues to the C terminus of TP. According to analyses of nuclear magnetic resonance measurements, while the beta-hairpin structure is stabilized in some DPs, it is evident that the helical structure observed in TP is also persistent and even extended throughout the length of the molecule. As a result, we have produced a peptide molecule that contains both the alpha-helix and beta-hairpin conformation at an almost equally populated level. The helical structures contained in these DPs were more stable than the helix in TP, suggesting that stabilizing one conformation does not result in destabilizing the other conformation. These DPs can thus be regarded as an isolated peptide version of the chameleon sequence, which has the capability of changing the secondary structure depending on the context of the surrounding environment in a protein structure. The fact that the transformation of one secondary structure caused stabilization of both the original and the induced structure would shed light on the mechanism of protein folding.

Polyisocyanides Derived from Tripeptides of Alanine

Helical polymers of isocyanotripeptides derived from alanine have been synthesized and their architectures studied in detail. The helical conformation of the polyisocyanotripeptides is stabilized by internal hydrogen-bonding arrays between the tripeptide side chains. The possibility of extending the well-defined hydrogen-bonded array, from dipeptide to tripeptide side chains, depends strongly on the stereochemistry of the constituent alanine amino acids, as has been shown by circular dichroism and IR studies. In polymers containing a weaker hydrogen-bonding array adjacent to the polymer backbone, due to steric interactions between the alanine methyl groups, a stronger second hydrogen-bonding array was present between the peptide bonds furthest away from the main chain, which is probably a result of stretching/compression of the helical-polymer conformation.

Measuring cooperativity in the formation of a three-stranded beta sheet (double hairpin)

Recently validated chemical shift measures of hairpin structuring have been applied to a series of turn mutants of the Schenck-Gellman three-strand beta-sheet model with the aim of measuring the entropic advantage associated with aligning an additional strand onto an existing hairpin versus aligning the same two strands in an initial hairpin formation. In a four-state analysis (unfolded, 2 single hairpins, and the double hairpin fold in equilibrium) a cooperativity index can be defined as the factor by which the equilibrium constant for hairpin formation is improved when one strand is prestructured. This cooperativity index is 2.7 +/- 0.7 for hairpin formation about a stable D-Pro-Gly turn locus and increases to 7.6 +/- 1.2 for an Asn-Gly turn locus. The latter corresponds to a cooperativity induced DeltaDeltaG increment of 4.9 kJ/mol for the folding of a hairpin. Although larger than previous experimental measures of folding cooperativity in three-stranded sheets, the magnitude of this effect (which is considerably less than the TDeltaDeltaS expectation for prestructuring three or more beta-strand residue sites) likely reflects the intrinsic preference of these designed sequences for extended conformations. If similar or larger effects apply to protein beta-sheet folding, it is not surprising that particularly favorable hairpin alignments serve as nucleation sites in protein folding pathways.

Design and NMR conformational study of a beta-sheet peptide based on Betanova and WW domains

A good approach to test our current knowledge on formation of protein beta-sheets is de novo protein design. To obtain a three-stranded beta-sheet mini-protein, we have built a series of chimeric peptides by taking as a template a previously designed beta-sheet peptide, Betanova-LLM, and incorporating N- and/or C-terminal extensions taken from WW domains, the smallest natural beta-sheet domain that is stable in absence of disulfide bridges. Some Betanova-LLM strand residues were also substituted by those of a prototype WW domain. The designed peptides were cloned and expressed in Escherichia coli. The ability of the purified peptides to adopt beta-sheet structures was examined by circular dichroism (CD). Then, the peptide showing the highest beta-sheet population according to the CD spectra, named 3SBWW-2, was further investigated by 1H and 13C NMR. Based on NOE and chemical shift data, peptide 3SBWW-2 adopts a well defined three-stranded antiparallel beta-sheet structure with a disordered C-terminal tail. To discern between the contributions to beta-sheet stability of strand residues and the C-terminal extension, the structural behavior of a control peptide with the same strand residues as 3SBWW-2 but lacking the C-terminal extension, named Betanova-LYYL, was also investigated. beta-Sheet stability in these two peptides, in the parent Betanova-LLM and in WW-P, a prototype WW domain, decreased in the order WW-P > 3SBWW-2 > Betanova-LYYL > Betanova-LLM. Conclusions about the contributions to beta-sheet stability were drawn by comparing structural properties of these four peptides.

A multibody, whole-residue potential for protein structures, with testing by Monte Carlo simulated annealing

A new multibody, whole-residue potential for protein tertiary structure is described. The potential is based on the local environment surrounding each main-chain alpha carbon (CA), defined as the set of all residues whose CA coordinates lie within a spherical volume of set radius in 3-dimensional (3D) space surrounding that position. It is shown that the relative positions of the CAs in these local environments belong to a set of preferred templates. The templates are derived by cluster analysis of the presently available database of over 3000 protein chains (750,000 residues) having not more than 30% sequence similarity. For each template is derived also a set of residue propensities for each topological position in the template. Using lookup tables of these derived templates, it is then possible to calculate an energy for any conformation of a given protein sequence. The application of the potential to ab initio protein tertiary structure prediction is evaluated by performing Monte Carlo simulated annealing on test protein sequences.

Engineering Diverse Changes in β-Turn Propensities in the N-Terminal β-Hairpin of Ubiquitin Reveals Significant Effects on Stability and Kinetics but a Robust Folding Transition State

Using the N-terminal 17-residue beta-hairpin of ubiquitin as a "host" for mutational studies, we have investigated the influence of the beta-turn sequence on protein stability and folding kinetics by replacing the native G-bulged turn (TLTGK) with more flexible analogues (TG3K and TG5K) and a series of four-residue type I' beta-turn sequences, commonly found in beta-hairpins. Although a statistical analysis of type I' turns demonstrates residue preferences at specific sites, the frequency of occurrence appears to only broadly correlate with experimentally determined protein stabilities. The subsequent engineering of context-dependent non-native tertiary contacts involving turn residues is shown to produce large changes in stability. Relatively few point mutations have been described that probe secondary structure formation in ubiquitin in a manner that is independent of tertiary contacts. To this end, we have used the more rigorous rate-equilibrium free energy relationship (Leffler analysis), rather than the two-point phi value analysis, to show for a family of engineered beta-turn mutants that stability (range of approximately 20 kJ/mol) and folding kinetics (190-fold variation in refolding rate) are linearly correlated (alpha(f) = 0.74 +/- 0.08). The data are consistent with a transition state that is robust with regard to a wide range of statistically favored and disfavored beta-turn mutations and implicate a loosely assembled beta-hairpin as a key template in transition state stabilization with the beta-turn playing a central role.

Cooperativity in the self-assembly of porphyrin ladders

Cooperativity is a general feature of intermolecular interactions in biomolecular systems, but there are many different facets of the phenomenon that are not well understood. Positive cooperativity stabilizes a system as progressively more interactions are added, and the origin of the beneficial free energy may be entropic or enthalpic in origin. An "enthalpic chelate effect" has been proposed to operate through structural tightening that improves all of the functional group interactions in a complex, when it is more strongly bound. Here, we present direct calorimetric evidence that no such enthalpic effects exist in the cooperative assembly of supramolecular ladder complexes composed of metalloporphyrin oligomers coordinated to bipyridine ligands. The enthalpic contributions of the individual coordination interactions are 35 kJ.mol(-1) and constant over a range of free energies of self-assembly of -35 to -111 kJ.mol(-1). In rigid well defined systems of this type, the enthalpies of individual interactions are additive, and no enthalpic cooperative effects are apparent. The implication is that in more flexible, less well defined systems such as biomolecular assemblies, the enthalpy contributions available from specific functional group interactions are well defined and constant parameters.

Chemical shifts provide fold populations and register of beta hairpins and beta sheets

A detailed analysis of peptide backbone amide (H(N)) and H alpha chemical shifts reveals a consistent pattern for beta hairpins and three-stranded beta sheets. The H alpha's at non-hydrogen-bonded strand positions are inwardly directed and shifted downfield approximately 1 ppm due largely to an anisotropy contribution from the cross-strand amide function. The secondary structure associated H alpha shift deviations for the H-bonded strand positions are also positive but much smaller (0.1-0.3 ppm) and the turn residues display negative H alpha chemical shift deviations (CSDs). The pattern of (H(N)) shift deviations is an even better indicator of both hairpin formation and register, with the cross-strand H-bonded sites shifted downfield (also by approximately 1 ppm) and with diagnostic values for the first turn residue and the first strand position following the turn. These empirical observations, initially made for [2:2]/[2:4]-type-I' and -II' hairpins, are rationalized and can be extended to the analysis of other turns, hairpin classes ([3:5], [4:4]/[4:6]), and three-stranded peptide beta-sheet models. The H alpha's at non-hydrogen-bonded sites and (H(N))'s in the intervening H-bonded sites provide the largest and most dependable measures of hairpin structuring and can be used for melting studies; however the intrinsic temperature dependence of (H(N)) shifts deviations needs to reflect the extent of solvent sequestration in the folded state. Several observations made in the course of this study provide insights into beta-sheet folding mechanisms: (1) The magnitude of the (H(N)) shifts suggests that the cross-strand H-bonds in peptide hairpins are as short as those in protein beta sheets. (2) Even L-Pro-Gly turns, which are frequently used in unfolded controls for hairpin peptides, can support hairpin populations in aqueous fluoroalcohol media. (3) The good correlation between hairpin population estimates from cross-strand H-bonded (H(N)) shift deviations, H alpha shift deviations, and structuring shifts at the turn locus implies that hairpin folding transitions approximate two-state behavior.

Folding cooperativity in a three-stranded beta-sheet model

The thermodynamic behavior of a previously designed three-stranded beta-sheet was studied via several microseconds of standard and replica exchange molecular dynamics simulations. The system is shown to populate at least four thermodynamic minima, including two partially folded states in which only a single hairpin is formed. Simulated melting curves show different profiles for the C and N-terminal hairpins, consistent with differences in secondary structure content in published NMR and CD/FTIR measurements, which probed different regions of the chain. Individual beta-hairpins that comprise the three-stranded beta-sheet are observed to form cooperatively. Partial folding cooperativity between the component hairpins is observed, and good agreement between calculated and experimental values quantifying this cooperativity is obtained when similar analysis techniques are used. However, the structural detail in the ensemble of conformations sampled in the simulations permits a more direct analysis of this cooperativity than has been performed on the basis of experimental data. The results indicate the actual folding cooperativity perpendicular to strand direction is significantly larger than the lower bound obtained previously.

Thermodynamic basis for antibody binding to Z-DNA: comparison of a monoclonal antibody and its recombinant derivatives

Antibody engineering represents a promising area in biotechnology. Recombinant antibodies can be easily manipulated generating new ligand and effector activities that can be used as prototype magic bullets. On the other hand, an extensive knowledge of recombinant antibody binding and stability features are essential for an efficient substitution. In this study, we compared the stability and protein binding properties of two recombinant antibody fragments with their parental monoclonal antibody. The recombinant fragments were a monomeric scFv and a dimeric one, harboring human IgG1 CH2-CH3 domains. We have used fluorescence titration quenching to determine the thermodynamics of the interaction between an anti-Z-DNA monoclonal antibody and its recombinant antibody fragments with Z-DNA. All the antibody fragments seemed to bind DNA similarly, in peculiar two-affinity states. Enthalpy-entropy compensation was observed for both affinity states, but a marked entropy difference was observed for the monomeric scFv antibody fragment, mainly for the high affinity binding. In addition, we compared the stability of the dimeric antibody fragment and found differences favoring the monoclonal antibody. These differences seem to derive from the heterologous expression system used.

Hairpin Folding Dynamics: The Cold-Denatured State Is Predisposed for Rapid Refolding

Cold denaturation is a general phenomenon in globular proteins, and the associated cold-denatured states of proteins have important fundamental and practical significance. Here, we have characterized the cold-denatured state of a beta-hairpin forming peptide, MrH3a, in 8% hexafluoro-2-propanol (HFIP) and the dynamics of its refolding following a laser-induced T-jump. Beta-hairpins constitute an important class of protein structural elements, yet their folding mechanisms are not fully understood. Characterization of MrH3a using NMR, CD, and IR spectroscopies reveals residual structure in the cold-denatured state, in contrast with the highly disordered heat-denatured state. The residual structure in the cold-denatured state comprises relatively compact and solvent protected conformations. Furthermore, we find a substantial acceleration in the rate of folding from the cold-denatured state compared to that of the heat-denatured state. In addition, the cold-denatured state is not populated in 20% HFIP; folding occurs only from the fully unfolded state and is significantly slower. We interpret the acceleration of the folding rate of MrH3a in 8% HFIP as a direct consequence of the collapsed conformations of the cold-denatured state. Finally, there may be some reduction of the loop search cost when starting from the cold-denatured state, since this state may have some of the stabilizing cross-strand interactions already formed.

Multicoefficient Extrapolated Density Functional Theory Studies of π···π Interactions: The Benzene Dimer

We report tests of new (2005) and established (1999-2003) multilevel methods against essentially converged benchmark results for nonbonded interactions in benzene dimers. We found that the newly developed multicoefficient extrapolated density functional theory (DFT) methods (which combine DFT with correlated wave function methods) give better performance than multilevel methods such as G3SX, G3SX(MP3), and CBS-QB3 that are based purely on wave function theory (WFT); furthermore, they have a lower computational cost. We conclude that our empirical approach for combining WFT methods with DFT methods is a very efficient and effective way for describing not only covalent interactions (as shown previously) but also nonbonded interactions.

Exploring amino-acid radical chemistry: protein engineering and de novo design

Amino-acid radical enzymes are often highly complex structures containing multiple protein subunits and cofactors. These properties have in many cases hampered the detailed characterization of their amino-acid redox cofactors. To address this problem, a range of approaches has recently been developed in which a common strategy is to reduce the complexity of the radical-containing system. This work will be reviewed and it includes the light-induced generation of aromatic radicals in small-molecule and peptide systems. Natural redox proteins, including the blue copper protein azurin and a bacterial photosynthetic reaction center, have been engineered to introduce amino-acid radical chemistry. The redesign strategies to achieve this remarkable change in the properties of these proteins will be described. An additional approach to gain insights into the properties of amino-acid radicals is to synthesize de novo designed model proteins in which the redox chemistry of these species can be studied. Here we describe the design, synthesis and characteristics of monomeric three-helix bundle and four-helix bundle proteins designed to study the redox chemistry of tryptophan and tyrosine. This work demonstrates that de novo protein design combined with structural, electrochemical and quantum chemical analyses can provide detailed information on how the protein matrix tunes the thermodynamic properties of tryptophan.

Structure and stability of β-pleated sheets

Beside alpha-helices, beta-sheets are the most common secondary structure elements of proteins. In this article, the question of structure and stability of parallel and antiparallel sheets of various lengths is addressed. All data obtained are compared to a selected set of protein structures. In antiparallel beta-sheets, one of the two possible H-bonded structures (containing 14 atoms in the H-bonded pseudoring) is energetically more favored and also more abundant in proteins than the other one (with 10 atoms involved in the pseudoring). Parallel beta-sheets and their subunits are energetically less stable and indeed found to occur more rarely in proteins. Antiparallel hairpins are disfavored compared to beta-sheets formed by sequentially separated strands. Agreement between theory and experimental data indicates that characterization of structural building blocks at an appropriately accurate level of theory is a useful tool to get insight into fundamentals of protein structure.

How well can new-generation density functional methods describe stacking interactions in biological systems?

We compare the performance of four recently developed DFT methods (MPW1B95, MPWB1K, PW6B95, and PWB6K) and two previous, generally successful DFT methods (B3LYP and B97-1) for the calculation of stacking interactions in six nucleic acid bases complexes and five amino acid pairs and for the calculation of hydrogen bonding interactions in two Watson-Crick type base pairs. We found that the four newly developed DFT methods give reasonable results for the stacking interactions in the DNA base pairs and amino acid pairs, whereas the previous DFT methods fail to describe interactions in these stacked complexes. We conclude that the new generation of DFT methods have greatly improved performance for stacking interaction as compared to previously available methods. We recommend the PWB6K method for investigating large DNA or protein systems where stacking plays an important role.