Role of a nonnative interaction in the folding of the protein G B1 domain as inferred from the conformational analysis of the alpha-helix fragment

A three-state mechanism for trifluoroethanol denaturation of an intrinsically disordered protein (IDP)

Relating the amino acid composition and sequence to chain folding and binding preferences of intrinsically disordered proteins (IDPs) has emerged as a huge challenge. While globular proteins have respective 3D structures that are unique to their individual functions, IDPs violate this structure-function paradigm because rather than having a well-defined structure an ensemble of rapidly interconverting disordered structures characterize an IDP. This work measures 2,2,2-trifluoroethanol (TFE)-induced equilibrium transitions of an IDP called AtPP16-1 (Arabidopsis thaliana phloem protein type 16-1) by using fluorescence, circular dichroism, infrared and nuclear magnetic resonance (NMR) methods at pH 4, 298 K. Low TFE reversibly removes the tertiary structure to produce an ensemble of obligate intermediate ($\mathrm{I}$) retaining the native-state ($\mathrm{N}$) secondary structure. The intermediate $\mathrm{I}$ is preceded by a non-obligate tryptophan-specific intermediate ${\mathrm{I}}_{\mathrm{w}}$ whose population is detectable for AtPP16-1 specifically. Accumulation of such non-obligate intermediates is discriminated according to the sequence composition of the protein. In all cases, however, a tertiary structure-unfolded general obligate intermediate $\mathrm{I}$ is indispensable. The $\mathrm{I}$ ensemble has higher helical propensity conducive to the acquisition of an exceedingly large level of α-helices by a reversible denaturation transition of $\mathrm{I}$ to the denatured state $\mathrm{D}$ as the TFE level is increased. Strikingly, it is the same $\mathrm{N}\rightleftharpoons \mathrm{I}\rightleftharpoons \mathrm{D}$ scheme typifying the TFE transitions of globular proteins. The high-energy state $\mathrm{I}$ characterized by increased helical propensity is called a universal intermediate encountered in both genera of globular and disordered proteins. Neither $\mathrm{I}$ nor $\mathrm{D}$ strictly show molten globule (MG)-like properties, dismissing the belief that TFE promotes MGs.

A surprisingly simple three-state generic process for reversible protein denaturation by trifluoroethanol

Despite the rich knowledge of the influence of 2,2,2-trifluoroethanol (TFE) on the structure and conformation of peptides and proteins, the mode(s) of TFE-protein interactions and the mechanism by which TFE reversibly denatures a globular protein remain elusive. This study systematically examines TFE-induced equilibrium transition curves for six paradigmatic globular proteins by using basic fluorescence and circular dichroism measurements under neutral pH conditions. The results are remarkably simple. Low TFE invariably unfolds the tertiary structure of all proteins to produce the obligate intermediate (I) which retains nearly all of native-state secondary structure, but enables the formation of extra α-helices as the level of TFE is raised higher. Inspection of the transitions at once reveals that the tertiary structure unfolding is always a distinct process, necessitating the inclusion of at least one obligate intermediate in the TFE-induced protein denaturation. It appears that the intermediate in the minimal unfolding mechanism N⇌I⇌D somehow acquires higher α-helical propensity to generate α-helices in excess of that in the native state to produce the denatured state (D), also called the TFE state. The low TFE-populated intermediate I may be called a universal intermediate by virtue of its α-helical propensity. Contrary to many earlier suggestions, this study dismisses molten globule (MG)-like attribute of I or D.

Dried Protein Structure Revealed at the Residue Level by Liquid-Observed Vapor Exchange NMR

Water is key to protein structure and stability, yet the relationship between protein-water interactions and structure is poorly understood, in part because there are few techniques that permit the study of dehydrated protein structure at high resolution. Here, we describe liquid-observed vapor exchange (LOVE) NMR, a solution NMR-based method that provides residue-level information about the structure of dehydrated proteins. Using the model protein GB1, we show that LOVE NMR measurements reflect the fraction of the dried protein population trapped in a conformation where a given residue is protected from exchange with D₂O vapor. Comparisons to solution hydrogen-deuterium exchange data affirm that the dried protein structure is strongly influenced by local solution stability and that the mechanism of dehydration protection exerted by the widely used protectant trehalose differs from its mechanism of stabilization in solution. Our results highlight the need for refined models of cosolute-mediated dehydration protection and demonstrate the ability of LOVE NMR to inform such models.

Characterizing protein G B1 orientation and its effect on immunoglobulin G antibody binding using XPS, ToF-SIMS, and quartz crystal microbalance with dissipation monitoring

Controlling how proteins are immobilized (e.g., controlling their orientation and conformation) is essential for developing and optimizing the performance of in vitro protein-binding devices, such as enzyme-linked immunosorbent assays. Characterizing the identity, orientation, etc., of proteins in complex mixtures of immobilized proteins requires a multitechnique approach. The focus of this work was to control and characterize the orientation of protein G B1, an immunoglobulin G (IgG) antibody-binding domain of protein G, on well-defined surfaces and to measure the effect of protein G B1 orientation on IgG antibody binding. The surface sensitivity of time-of-flight secondary ion mass spectrometry (ToF-SIMS) was used to distinguish between different proteins and their orientation on both flat and nanoparticle gold surfaces by monitoring intensity changes of characteristic amino acid mass fragments. Amino acids distributed asymmetrically were used to calculate peak intensity ratios from ToF-SIMS data to determine the orientation of protein G B1 cysteine mutants covalently attached to a maleimide surface. To study the effect of protein orientation on antibody binding, multilayer protein films on flat gold surfaces were formed by binding IgG to the immobilized protein G B1 films. Quartz crystal microbalance with dissipation monitoring and x-ray photoelectron spectroscopy analysis revealed that coverage and orientation affected the antibody-binding process. At high protein G B1 coverage, the cysteine mutant immobilized in an end-on orientation with the C-terminus exposed bound 443 ng/cm² of whole IgG (H + L) antibodies. In comparison, the high coverage cysteine mutant immobilized in an end-on orientation with the N-terminus exposed did not bind detectable amounts of whole IgG (H + L) antibodies.

Quinary interactions with an unfolded state ensemble

Anfinsen's thermodynamic hypothesis states that the native three-dimensional fold of a protein represents the structure with the lowest Gibbs free energy. Changes in the free energy of denaturation can arise from changes to the folded state, the unfolded state, or both. It has been recently recognized that quinary interactions, transient contacts that take place only in cells, can modulate protein stability through interactions involving the folded state. Here we show that the cellular environment can also remodel the unfolded state ensemble.

Thermodynamic and structural effect of urea and guanidine chloride on the helical and on a hairpin fragment of GB1 from molecular simulations

With the help of molecular-dynamics simulations, we studied the effect of urea and guanidine chloride on the thermodynamic and structural properties of the helical fragment of protein GB1, comparing them with those of its second beta hairpin. We showed that the helical fragment in different solvents populates an ensemble of states that is more complex than that of the hairpin, and thus the associated experimental observables (circular-dichroism spectra, secondary chemical shifts, m values), that we back-calculated from the simulations and compared with the actual data, are more difficult to interpret. We observed that in the case of both peptides, urea binds tightly to their backbone, while guanidine exerts its denaturing effect in a more subtle way, strongly affecting the electrostatic properties of the solution. This difference can have consequences in the way denaturation experiments are interpreted. Proteins 2017; 85:753-763. © 2016 Wiley Periodicals, Inc.

The low populated folding intermediate of a mutant of the Fyn SH3 domain identified by a simple model

The low populated on-pathway folding intermediate of the A39V/N53P/V55L Fyn SH3 domain is captured by a native-centric model augmented by sequence-dependent nonnative hydrophobic interactions.

Radical Formation Initiates Solvent-Dependent Unfolding and β-sheet Formation in a Model Helical Peptide

We examined the effects of Cα-centered radical formation on the stability of a model helical peptide, N-Ac-KK(AL)10KK-NH2. Three, 100 ns molecular dynamics simulations using the OPLS-AA force field were carried out on each α-helical peptide in six distinct binary TIP4P water/2,2,2-trifluoroethanol (TFE) mixtures. The α-helicity was at a maximum in 20% TFE, which was inversely proportional to the number of H-bonds between water molecules and the peptide backbone. The radial distribution of TFE around the peptide backbone was highest in 20% TFE, which enhanced helix stability. The Cα-centered radical initiated the formation of a turn within 5 ns, which was a smaller kink at high TFE concentrations, and a loop at lower TFE concentrations. The highest helicity of the peptide radical was measured in 100% TFE. The formation of hydrogen bonds between the peptide backbone and water destabilized the helix, whereas the clustering of TFE molecules around the radical center stabilized the helix. Following radical termination, the once helical structure converted to a β-sheet rich state in 100% water only, and this transition did not occur in the nonradical control peptide. This study gives evidence on how the formation of peptide radicals can initiate α-helical to β-sheet transitions under oxidative stress conditions.

Simulation of chaperonin effect on protein folding: a shift from nucleation-condensation to framework mechanism

The iterative annealing mechanism (IAM) of chaperonin-assisted protein folding is explored in a framework of a well-established coarse-grained protein modeling tool, which enables the study of protein dynamics in a time-scale well beyond classical all-atom molecular mechanics. The chaperonin mechanism of action is simulated for two paradigm systems of protein folding, B domain of protein A (BdpA) and B1 domain of protein G (GB1), and compared to chaperonin-free simulations presented here for BdpA and recently published for GB1. The prediction of the BdpA transition state ensemble (TSE) is in perfect agreement with experimental findings. It is shown that periodic distortion of the polypeptide chains by hydrophobic chaperonin interactions can promote rapid folding and leads to a decrease in folding temperature. It is also demonstrated how chaperonin action prevents kinetically trapped conformations and modulates the observed folding mechanisms from nucleation-condensation to a more framework-like.

Nonnative interactions in coupled folding and binding processes of intrinsically disordered proteins

Proteins function by interacting with other molecules, where both native and nonnative interactions play important roles. Native interactions contribute to the stability and specificity of a complex, whereas nonnative interactions mainly perturb the binding kinetics. For intrinsically disordered proteins (IDPs), which do not adopt rigid structures when being free in solution, the role of nonnative interactions may be more prominent in binding processes due to their high flexibilities. In this work, we investigated the effect of nonnative hydrophobic interactions on the coupled folding and binding processes of IDPs and its interplay with chain flexibility by conducting molecular dynamics simulations. Our results showed that the free-energy profiles became rugged, and intermediate states occurred when nonnative hydrophobic interactions were introduced. The binding rate was initially accelerated and subsequently dramatically decreased as the strength of the nonnative hydrophobic interactions increased. Both thermodynamic and kinetic analysis showed that disordered systems were more readily affected by nonnative interactions than ordered systems. Furthermore, it was demonstrated that the kinetic advantage of IDPs (“fly-casting” mechanism) was enhanced by nonnative hydrophobic interactions. The relationship between chain flexibility and protein aggregation is also discussed.

beta-Strand interactions at the domain interface critical for the stability of human lens gammaD-crystallin

Human age-onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens gammaD-crystallin protein and its isolated domains. A partially unfolded folding intermediate of gammaD-crystallin is detected in simulations with its C-terminal domain (C-td) folded and N-terminal domain (N-td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N-td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C-td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non-native surface salt-bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non-native Glu135-Arg142 salt-bridge causes significant destabilization to the folding core of the isolated C-td, which, in turn, induces unfolding of the N-td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of gammaD-crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the gammaD-crystallin surface might lead to unfolding and subsequent aggregation.

Folding pathway of the b1 domain of protein G explored by multiscale modeling

The understanding of the folding mechanisms of single-domain proteins is an essential step in the understanding of protein folding in general. Recently, we developed a mesoscopic CA-CB side-chain protein model, which was successfully applied in protein structure prediction, studies of protein thermodynamics, and modeling of protein complexes. In this research, this model is employed in a detailed characterization of the folding process of a simple globular protein, the B1 domain of IgG-binding protein G (GB1). There is a vast body of experimental facts and theoretical findings for this protein. Performing unbiased, ab initio simulations, we demonstrated that the GB1 folding proceeds via the formation of an extended folding nucleus, followed by slow structure fine-tuning. Remarkably, a subset of native interactions drives the folding from the very beginning. The emerging comprehensive picture of GB1 folding perfectly matches and extends the previous experimental and theoretical studies.

Characterization of protein-folding pathways by reduced-space modeling

Ab initio simulations of the folding pathways are currently limited to very small proteins. For larger proteins, some approximations or simplifications in protein models need to be introduced. Protein folding and unfolding are among the basic processes in the cell and are very difficult to characterize in detail by experiment or simulation. Chymotrypsin inhibitor 2 (CI2) and barnase are probably the best characterized experimentally in this respect. For these model systems, initial folding stages were simulated by using CA-CB-side chain (CABS), a reduced-space protein-modeling tool. CABS employs knowledge-based potentials that proved to be very successful in protein structure prediction. With the use of isothermal Monte Carlo (MC) dynamics, initiation sites with a residual structure and weak tertiary interactions were identified. Such structures are essential for the initiation of the folding process through a sequential reduction of the protein conformational space, overcoming the Levinthal paradox in this manner. Furthermore, nucleation sites that initiate a tertiary interactions network were located. The MC simulations correspond perfectly to the results of experimental and theoretical research and bring insights into CI2 folding mechanism: unambiguous sequence of folding events was reported as well as cooperative substructures compatible with those obtained in recent molecular dynamics unfolding studies. The correspondence between the simulation and experiment shows that knowledge-based potentials are not only useful in protein structure predictions but are also capable of reproducing the folding pathways. Thus, the results of this work significantly extend the applicability range of reduced models in the theoretical study of proteins.

β-Hairpin formation in aqueous solution and in the presence of trifluoroethanol: A1H and13C nuclear magnetic resonance conformational study of designed peptides

In order to check our current knowledge on the principles involved in beta-hairpin formation, we have modified the sequence of a 3:5 beta-hairpin forming peptide with two different purposes, first to increase the stability of the formed 3:5 beta-hairpin, and second to convert the 3:5 beta-hairpin into a 2:2 beta-hairpin. The conformational behavior of the designed peptides was investigated in aqueous solution and in 30% trifluoroethanol (TFE) by analysis of the following nuclear magnetic resonance (NMR) parameters: nuclear Overhauser effect (NOE) data, and C(alpha)H, (13)C(alpha), and (13)C(beta) conformational shifts. From the differences in the ability to adopt beta-hairpin structures in these peptides, we have arrived to the following conclusions: (i) beta-Hairpin population increases with the statistical propensity of residues to occupy each turn position. (ii) The loop length, and in turn, the beta-hairpin type, can be modified as a function of the type of turn favored by the loop sequence. These two conclusions reinforce previous results about the importance of beta-turn sequence in beta-hairpin folding. (iii) Side-chain packing on each face of the beta-sheet may play a major role in beta-hairpin stability; hence simplified analysis in terms of isolated pair interactions and intrinsic beta-sheet propensities is insufficient. (iv) Contributions to beta-hairpin stability of turn and strand sequences are not completely independent. (v) The burial of hydrophobic surface upon beta-hairpin formation that, in turn, depends on side-chain packing also contributes to beta-hairpin stability. (vi) As previously observed, TFE stabilizes beta-hairpin structures, but the extent of the contribution of different factors to beta-hairpin formation is sometimes different in aqueous solution and in 30% TFE.

Matching simulation and experiment: a new simplified model for simulating protein folding

Simulations of simplified protein folding models have provided much insight into solving the protein folding problem. We propose here a new off-lattice bead model, capable of simulating several different fold classes of small proteins. We present the sequence for an alpha/beta protein resembling the IgG-binding proteins L and G. The thermodynamics of the folding process for this model are characterized using the multiple multihistogram method combined with constant-temperature Langevin simulations. The folding is shown to be highly cooperative, with chain collapse nearly accompanying folding. Two parallel folding pathways are shown to exist on the folding free energy landscape. One pathway contains an intermediate--similar to experiments on protein G, and one pathway contains no intermediates-similar to experiments on protein L. The folding kinetics are characterized by tabulating mean-first passage times, and we show that the onset of glasslike kinetics occurs at much lower temperatures than the folding temperature. This model is expected to be useful in many future contexts: investigating questions of the role of local versus nonlocal interactions in various fold classes, addressing the effect of sequence mutations affecting secondary structure propensities, and providing a computationally feasible model for studying the role of solvation forces in protein folding.

Synthesis, folding, and structure of the beta-turn mimic modified B1 domain of streptococcal protein G

The mechanism of beta-sheet formation remains a fundamental issue in our understanding of the protein folding process, but is hampered by the often encountered kinetic competition between folding and aggregation. The role of local versus nonlocal interactions has been probed traditionally by mutagenesis of both turn and strand residues. Recently, rigid organic molecules that impose a correct chain reversal have been introduced in several small peptides to isolate the importance of the long-range interactions. Here, we present the incorporation of a well-studied beta-turn mimic, designated as the dibenzofuran-based (DBF) amino acid, in the B1 domain of streptococcal protein G (B1G), and compare our results with those obtained upon insertion of the same mimic into the N-terminal beta-hairpin of B1G (O Melnyk et al., 1998, Lett Pept Sci 5:147-150). The DBF-B1G domain conserves the structure and the functional and thermodynamical properties of the native protein, whereas the modified peptide does not adopt a native-like conformation. The nature of the DBF flanking residues in the modified B1G domain prevents the beta-turn mimic from acting as a strong beta-sheet nucleator, which reinforces the idea that the native beta-hairpin formation is not driven by the beta-turn formation, but by tertiary interactions.

THEFASTPROTEINFOLDINGPROBLEM

During protein folding, many of the events leading to secondary and tertiary structure occur in milliseconds or faster. Modern nuclear magnetic resonance and laser detection techniques, coupled with fast initiation of the folding reaction, are probing these events in great detail. Theory, ranging from analytical models to molecular dynamics calculations, is beginning to match up with experiment. As a result, timescales, from such elementary steps as the addition of a residue to a helix to strange kinetics of collapsing protein backbones, can now be measured and interpreted.

A tale of two secondary structure elements: when a beta-hairpin becomes an alpha-helix

In this work, we have analyzed the relative importance of secondary versus tertiary interactions in stabilizing and guiding protein folding. For this purpose, we have designed four different mutants to replace the alpha-helix of the GB1 domain by a sequence with strong beta-hairpin propensity in isolation. In particular, we have chosen the sequence of the second beta-hairpin of the GB1 domain, which populates the native conformation in aqueous solution to a significant extent. The resulting protein has roughly 30 % of its sequence duplicated and maintains the 3D-structure of the wild-type protein, but with lower stability (up to -5 kcal/mol). The loss of intrinsic helix stability accounts for about 80 % of the decrease in free energy, illustrating the importance of local interactions in protein stability. Interestingly enough, all the mutant proteins, included the one with the duplicated beta-hairpin sequence, fold with similar rates as the GB1 domain. Essentially, it is the nature of the rate-limiting step in the folding reaction that determines whether a particular interaction will speed up, or not, the folding rates. While local contacts are important in determining protein stability, residues involved in tertiary contacts in combination with the topology of the native fold, seem to be responsible for the specificity of protein structures. Proteins with non-native secondary structure tendencies can adopt stable folds and be as efficient in folding as those proteins with native-like propensities.

Conformations of peptide fragments from the FK506 binding protein: comparison with the native and urea-unfolded states

The helix-forming tendency of seven peptide fragments corresponding with the entire sequence of the FK506 binding protein (FKBP) has been investigated in aqueous buffer and in 2,2,2-trifluoroethanol (TFE) using CD and NMR spectroscopy. All fragments exhibited random coil conformations in aqueous buffer, whereas the amount of helix induced in the peptide fragments by TFE varied. The fragment with the highest degree of helicity in TFE corresponded with the single (alpha-helix in native FKBP. Fragments corresponding with beta-strands 2 and 3 also exhibited strong propensity towards helix formation. In contrast, the fragment corresponding with beta-strand 1 did not form helix in TFE. The inherent helix-forming tendencies are interpreted in light of the native structure to suggest possible folding nucleation sites. Conformational sampling in each peptide fragment was also compared with that observed in urea-denatured FKBP. With the exception of the fragment corresponding with beta-strand 2, the formation of helical structures in the peptide fragments in TFE was correlated with the observation of turn and/or helix conformers in urea-unfolded FKBP. Surprisingly, peptide fragments in aqueous solution were less structured than the corresponding regions in urea-denatured FKBP. The conformational differences between the peptide fragments and unfolded FKBP were not due to the urea buffer or to differences in their rotational correlation times. We conclude that local amino acid interactions are not generally sufficient to account for the formation of non-random conformations in unfolded FKBP. Formation of non-random structures in unfolded FKBP may require stabilization of incipient turn or helical conformations through transient contact with non-local non-polar residues.

Molecular dynamics as a tool to detect protein foldability. A mutant of domain B1 of protein G with non-native secondary structure propensities

The usefulness of molecular dynamics to assess the structural integrity of mutants containing several mutations has been investigated. Our goal was to determine whether molecular dynamics would be able to discriminate mutants of a protein having a close-to-wild-type fold, from those that are not folded under the same conditions. We used as a model the B1 domain of protein G in which we replaced the unique central alpha-helix by the sequence of the second beta-hairpin, which has a strong intrinsic propensity to form this secondary structure in solution. In the resulting protein, one-third of the secondary structure has been replaced by a non-native one. Models of the mutants were built based on the three-dimensional structure of the wild-type GB1 domain. During 2 ns of molecular dynamics simulations on these models, mutants containing up to 10 mutations in the helix retained the native fold, while another mutant with an additional mutation unfolded. This result is in agreement with our circular dichroism and NMR experiments, which indicated that the former mutants fold into a structure similar to the wild-type, as opposed to the latter mutant which is partly unfolded. Additionally, a mutant containing six mutations scattered through the surface of the domain, and which is unfolded, was also detected by the simulation. This study suggests that molecular dynamics calculations could be performed on molecular models of mutants of a protein to evaluate their foldability, prior to a mutagenesis experiment.

Assembly of a ternary complex by the predicted minimal coiled-coil-forming domains of syntaxin, SNAP-25, and synaptobrevin. A circular dichroism study

The assembly of target (t-SNARE) and vesicle-associated SNAP receptor (v-SNARE) proteins is a critical step for the docking of synaptic vesicles to the plasma membrane. Syntaxin-1A, SNAP-25, and synaptobrevin-2 (also known as vesicle-associated membrane protein, or VAMP-2) bind to each other with high affinity, and their binding regions are predicted to form a trimeric coiled-coil. Here, we have designed three peptides, which correspond to sequences located in the syntaxin-1A H3 domain, the C-terminal domain of SNAP-25, and a conserved central domain of synaptobrevin-2, that exhibit a high propensity to form a minimal trimeric coiled-coil. The peptides were synthesized by solid phase methods, and their interactions were studied by CD spectroscopy. In aqueous solution, the peptides were unstructured and showed no interactions with each other. In contrast, upon the addition of moderate amounts of trifluoroethanol (30%), the peptides adopted an alpha-helical structure and displayed both homomeric and heteromeric interactions. The interactions observed in ternary mixtures induce a stabilization of peptide structure that is greater than that predicted from individual binary interactions, suggesting the formation of a higher order structure compatible with the assembly of a trimeric coiled-coil.

Elucidating the folding problem of alpha-helices: local motifs, long-range electrostatics, ionic-strength dependence and prediction of NMR parameters

The information about the conformational behavior of monomeric helical peptides in solution, as well as the alpha-helix stability in proteins, has been previously utilized to derive a database with the energy contributions for various interactions taking place in an alpha-helix: intrinsic helical propensities, side-chain-side-chain interactions, main-chain-main-chain hydrogen bonds, and capping effects. This database was implemented in an algorithm based on the helix/coil transition theory (AGADIR). Here, we have modified this algorithm to include previously described local motifs: hydrophobic staple, Schellman motif and Pro-capping motif, new variants of these, and newly described side-chain-side-chain interactions. Based on recent experimental data we have introduced a position dependence of the helical propensities for some of the 20 amino acid residues. A new electrostatic model that takes into consideration all electrostatic interactions up to 12 residues in distance in the helix and random-coil conformations, as well as the effect of ionic strength, has been implemented. We have synthesized and analyzed several peptides, and used data from peptides already analysed by other groups, to test the validity of our electrostatic model. The modified algorithm predicts, with an overall standard deviation value of 6.6 (maximum helix is 100%), the helical, content of 778 peptides of which 223 correspond to wild-type and modified protein fragments. To improve the prediction potential of the algorithm and to have a direct comparison with nuclear magnetic resonance data, the algorithm now predicts the conformational shift of the CalphaH protons, 13Calpha and 3JalphaN values. We have found that for those peptides correctly predicted from the point of view of circular dichroism, the prediction of the NMR parameters is very good.

Structure of human cyclin-dependent kinase inhibitor p19INK4d: comparison to known ankyrin-repeat-containing structures and implications for the dysfunction of tumor suppressor p16INK4a

BACKGROUND

The four members of the INK4 gene family (p16(INK4a), p15(INK4b), p18(INK4c) and p19(INK4d)) inhibit the closely related cyclin-dependent kinases CDK4 and CDK6 as part of the regulation of the G1-->S transition in the cell-division cycle. Loss of INK4 gene product function, particularly that of p16(INK4a), is found in 10-60% of human tumors, suggesting that broadly applicable anticancer therapies might be based on restoration of p16(INK4a) CDK inhibitory function. Although much less frequent, defects of p19(INK4d) have also been associated with human cancer (osteosarcomas). The protein structures of some INK4 family members, determined by nuclear magnetic resonance (NMR) spectroscopy and X-ray techniques, have begun to clarify the functional role of p16(INK4a) and the dysfunction introduced by the mutations associated with human tumors.

RESULTS

The crystal structure of human p19(INK4d) has been determined at 1.8 A resolution using multiple isomorphous replacement methods. The fold of p19(INK4d) produces an oblong molecule comprising five approximately 32-residue ankyrin-like repeats. The architecture of the protein demonstrates the high structural similarity within the INK4 family. Comparisons to other ankyrin-repeat-containing proteins (GABPbeta, 53BP2 and myotrophin) show similar structures with comparable hydrogen-bonding patterns and hydrophobic interactions. Such comparisons highlight the splayed beta-loop geometry that is specific to INK4 inhibitors. This geometry is the result of a modified ankyrin structure in the second repeat.

CONCLUSIONS

Among the INK4 inhibitors, the highest amino acid sequence conservation is found in the helical stacks; this conservation creates a conserved beta-loop geometry specific to INK4 inhibitors. Therefore, in addition to models which predict that the conserved helix alpha6 is responsible for CDK inhibition, a binding mode whereby the loops of INK4 proteins bind to the CDKs should also be considered. A similar loop-based interaction is seen in the complex formed between the ankyrin-repeat-containing protein GABPbeta and_GABPalpha. This mode of binding would be consistent with the observation that p16(INK4a) is sensitive to deleterious mutations found throughout this tumor suppressor protein; these mutations probably destabilize the three-dimensional structure.

Calculations on folding of segment B1 of streptococcal protein G

We present an investigation of the folding thermodynamics and mechanism of segment B1 of streptococcal protein G. Molecular dynamics simulations of the fully solvated protein are used to probe thermodynamically significant states at different stages of folding. We performed several unfolding simulations to generate a database of initial conditions. The database is analyzed and clustered. The cluster centers extracted from this database were then used as starting points for umbrella sampling of the folding free energy landscape under folding conditions. The resulting sampling was combined with the weighted histogram analysis method. One and two-dimensional free energy surfaces were constructed along several order parameters and used to analyze the folding process. Our findings indicate that an initial collapse precedes the formation of significant native structure. Elements of local structure originate in the regions of the protein shown to have higher H/2H exchange protection factors in early stages of folding. A non-native contact, observed experimentally at the N terminus of the alpha-helix in a peptide excised from the protein, is seen to pre-organize the chain in early stages of folding. Collapse and early structure formation yields a compact globule with a significant number of water molecules present. Desolvation of the protein core is coincident with the final stages of folding from the compact state.

Characterisation of urea-denatured states of an immunoglobulin superfamily domain by heteronuclear NMR

The structural and dynamic properties of an immunoglobulin superfamily domain (IgSF), Ig 18', have been characterised by NMR at 285 K, in the presence of 4.2 M and 6.0 M urea, respectively. Analysis of chemical shift deviations, 3JHNHalpha coupling constants, sequential NOE pattern, and 15N relaxation data reveals that although the two urea-denatured states are highly disordered, some local turn-like residual structures do exist. Moreover, some distinct differences between the properties of the two denatured states are observed. In 4.2 M urea-denatured Ig 18', regions 80-83 and 86-92 adopt turn-like conformations, furthermore, region 84-93 is involved in slow exchange processes that occur on a micro- to millisecond time-scale. In the 6.0 M urea-denatured state, these turn-like conformations are less occupied, and chemical exchange processes in region 84-93 are largely reduced. In contrast, region 32-36 has persistent turn-like structures in both urea-denatured states. Some correlation between the spectral density function at 0 frequency, Jeff(0), for the urea-denatured states and the secondary structure elements of the folded state have been observed. Except for the terminal regions, residues corresponding to beta-strands have higher Jeff(0) values compared to residues corresponding to loops. The characterisation and comparison of the two urea-denatured states highlight residues that possess properties that may be crucial for the initiation of folding of this domain.

Molecular picture of folding of a small alpha/beta protein

We characterize, at the atomic level, the mechanism and thermodynamics of folding of a small alpha/beta protein. The thermodynamically significant states of segment B1 of streptococcal protein G (GB1) are probed by using the statistical mechanical methods of importance sampling and molecular dynamics. From a thermodynamic standpoint, folding commences with overall collapse, accompanied by formation of approximately 35% of the native structure. Specific contacts form at the loci experimentally inferred to be structured early in folding kinetics studies. Our study reveals that these initially structured regions are not spatially adjacent. As folding progresses, fluid-like nonlocal native contacts form, with many contacts forming and breaking as the structure searches for the native conformation. Although the alpha-helix forms early, the beta-sheet forms concomitantly with the overall topology. Water is present in the protein core up to a late stage of folding, lubricating conformational transitions during the search process. Once 80% of the native contacts have formed, water is squeezed from the protein interior and the structure descends into the native manifold. Examination of the onset of side-chain mobility within our model indicates side-chain motion is most closely linked to the overall volume of the protein and no sharp order-disorder transition appears to occur. Exploration of models for hydrogen deuterium exchange show qualitative agreement with equilibrium measurement of hydrogen/deuterium protection factors.