Dynamic and conformational switching in proteins.

Using bioinformatic methods for treating protein dynamics, developed in earlier work, we study the relationship between sequence mobility and dynamics in proteins. It is shown that sequence mobility drives a transition between two dynamic regimes in proteins, and that the specific details of this transition differ qualitatively between α-helical proteins and those in other structural classes. We examine the possibility that conformational switching is related to dynamic switching, by considering a specific system of sequences which exhibit the switching phenomenon. It is shown that a relationship between dynamic and conformational switching is entirely plausible.

Sequence-specific dynamic information in proteins.

We examine the local and global properties of the average B-factor, 〈B〉, as a residue-specific indicator of protein dynamic characteristics. It has been shown that values of 〈B〉 for the 20 amino acids differ in a statistically significant manner, and that, while strongly determined by the static physical properties of amino acids, they also encode averaged information about the influence of global fold on single-residue dynamics. Therefore, complete sequences of amino acids also encode fold-related global dynamic information, in addition to the local information that arises from static physical properties. We show that the relative magnitudes of these two contributions can be determined using Fourier methods, which represent the global properties of the sequences. It has also been shown that the behavior of Fourier components of 〈B〉 differs, with very high statistical significance, between structural groups, and that this information is not available from a comparable analysis of static amino acid properties.

Tracking the mechanism of fibril assembly by simulated two-dimensional ultraviolet spectroscopy.

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by the accumulation of plaque deposits in the human brain. The main component of these plaques consists of highly ordered structures called amyloid fibrils, formed by the amyloid ß-peptide (Aß). The mechanism connecting Aß and AD is yet undetermined. In a previous study, a coarse-grained united-residue model and molecular dynamics simulations were used to model the growth mechanism of Aß amyloid fibrils. On the basis of these simulations, a dock/lock mechanism was proposed, in which Aß fibrils grow by adding monomers at either end of an amyloid fibril template. To examine the structures in the early time-scale formation and growth of amyloid fibrils, simulated two-dimensional ultraviolet spectroscopy is used. These early structures are monitored in the far ultraviolet regime (λ = 190-250 nm) in which the computed signals originate from the backbone nπ* and ππ* transitions. These signals show distinct cross-peak patterns that can be used, in combination with molecular dynamics, to monitor local dynamics and conformational changes in the secondary structure of Aß-peptides. The protein geometry-correlated chiral xxxy signal and the non-chiral combined signal xyxy-xyyx were found to be sensitive to, and in agreement with, a dock/lock pathway.

Dynamic Formation and Breaking of Disulfide Bonds in Molecular Dynamics Simulations with the UNRES Force Field.

Many proteins contain disulfide bonds that are usually essential for maintaining function and a stable structure. Several algorithms attempt to predict the arrangement of disulfide bonds in the context of protein structure prediction, but none can simulate the entire process of oxidative folding, including dynamic formation and breaking of disulfide bonds. In this work, a potential function developed to model disulfide bonds is coupled with the united-residue (UNRES) force field, and used in both canonical and replica exchange molecular dynamics simulations to produce complete oxidative folding pathways. The potential function is obtained by introducing a transition barrier that separates the bonded and nonbonded states of the half-cystine residues. Tests on several helical proteins show that improved predictions are obtained when dynamic disulfide-bond formation and breaking are considered. The effect of the disulfide bonds on the folding kinetics is also investigated, particularly their role in stabilizing folding intermediates, resulting in slower folding.

A third blind test of crystal structure prediction.

Following the interest generated by two previous blind tests of crystal structure prediction (CSP1999 and CSP2001), a third such collaborative project (CSP2004) was hosted by the Cambridge Crystallographic Data Centre. A range of methodologies used in searching for and ranking the likelihood of predicted crystal structures is represented amongst the 18 participating research groups, although most are based on the global minimization of the lattice energy. Initially the participants were given molecular diagrams of three molecules and asked to submit three predictions for the most likely crystal structure of each. Unlike earlier blind tests, no restriction was placed on the possible space group of the target crystal structures. Furthermore, Z' = 2 structures were allowed. Part-way through the test, a partial structure report was discovered for one of the molecules, which could no longer be considered a blind test. Hence, a second molecule from the same category (small, rigid with common atom types) was offered to the participants as a replacement. Success rates within the three submitted predictions were lower than in the previous tests - there was only one successful prediction for any of the three ;blind' molecules. For the ;simplest' rigid molecule, this lack of success is partly due to the observed structure crystallizing with two molecules in the asymmetric unit. As in the 2001 blind test, there was no success in predicting the structure of the flexible molecule. The results highlight the necessity for better energy models, capable of simultaneously describing conformational and packing energies with high accuracy. There is also a need for improvements in search procedures for crystals with more than one independent molecule, as well as for molecules with conformational flexibility. These are necessary requirements for the prediction of possible thermodynamically favoured polymorphs. Which of these are actually realised is also influenced by as yet insufficiently understood processes of nucleation and crystal growth.

Combined solid-phase/solution synthesis of large ribonuclease A C-terminal peptides containing a non-natural proline analog.

Three large peptides corresponding to the 65-124 (60-mer), 72-124 (53-mer), and 77-124 (48-mer) sequence of bovine pancreatic ribonuclease A (RNase A) were assembled from either two or three shorter protected peptide fragments by chemical coupling in solution. The fragments were synthesized manually by 9-fluorenylmethyloxycarbonyl (Fmoc)-based solid-phase peptide chemistry in plastic syringes, and subsequently purified by normal-phase high-performance liquid chromatography on a silica gel column. The main aim of this work was to incorporate sterically hindered l-5,5-dimethylproline (dmP) as a substitute for Pro(93) into the sequence of RNase A in order to constrain the -Tyr(92)-Pro(93)- peptide group to a single cis-conformation.

Physics-based protein-structure prediction using a hierarchical protocol based on the UNRES force field: assessment in two blind tests.

Recent improvements in the protein-structure prediction method developed in our laboratory, based on the thermodynamic hypothesis, are described. The conformational space is searched extensively at the united-residue level by using our physics-based UNRES energy function and the conformational space annealing method of global optimization. The lowest-energy coarse-grained structures are then converted to an all-atom representation and energy-minimized with the ECEPP/3 force field. The procedure was assessed in two recent blind tests of protein-structure prediction. During the first blind test, we predicted large fragments of alpha and alpha+beta proteins [60-70 residues with C(alpha) rms deviation (rmsd) <6 A]. However, for alpha+beta proteins, significant topological errors occurred despite low rmsd values. In the second exercise, we predicted whole structures of five proteins (two alpha and three alpha+beta, with sizes of 53-235 residues) with remarkably good accuracy. In particular, for the genomic target TM0487 (a 102-residue alpha+beta protein from Thermotoga maritima), we predicted the complete, topologically correct structure with 7.3-A C(alpha) rmsd. So far this protein is the largest alpha+beta protein predicted based solely on the amino acid sequence and a physics-based potential-energy function and search procedure. For target T0198, a phosphate transport system regulator PhoU from T. maritima (a 235-residue mainly alpha-helical protein), we predicted the topology of the whole six-helix bundle correctly within 8 A rmsd, except the 32 C-terminal residues, most of which form a beta-hairpin. These and other examples described in this work demonstrate significant progress in physics-based protein-structure prediction.

Impact of an easily reducible disulfide bond on the oxidative folding rate of multi-disulfide-containing proteins.

The burial of native disulfide bonds, formed within stable structure in the regeneration of multi-disulfide-containing proteins from their fully reduced states, is a key step in the folding process, as the burial greatly accelerates the oxidative folding rate of the protein by sequestering the native disulfide bonds from thiol-disulfide exchange reactions. Nevertheless, several proteins retain solvent-exposed disulfide bonds in their native structures. Here, we have examined the impact of an easily reducible native disulfide bond on the oxidative folding rate of a protein. Our studies reveal that the susceptibility of the (40-95) disulfide bond of Y92G bovine pancreatic ribonuclease A (RNase A) to reduction results in a reduced rate of oxidative regeneration, compared with wild-type RNase A. In the native state of RNase A, Tyr 92 lies atop its (40-95) disulfide bond, effectively shielding this bond from the reducing agent, thereby promoting protein oxidative regeneration. Our work sheds light on the unique contribution of a local structural element in promoting the oxidative folding of a multi-disulfide-containing protein.

The protein folding problem: global optimization of the force fields.

The evolutionary development of a theoretical approach to the protein folding problem, in our laboratory, is traced. The theoretical foundations and the development of a suitable empirical all-atom potential energy function and a global optimization search are examined. Whereas the all-atom approach has thus far succeeded for relatively small molecules and for alpha-helical proteins containing up to 46 residues, it has been necessary to develop a hierarchical approach to treat larger proteins. In the hierarchical approach to single- and multiple-chain proteins, global optimization is carried out for a simplified united residue (UNRES) description of a polypeptide chain to locate the region in which the global minimum lies. Conversion of the UNRES structures in this region to all-atom structures is followed by a local search in this region. The performance of this approach in successive CASP blind tests for predicting protein structure by an ab initio physics-based method is described. Finally, a recent attempt to compute a folding pathway is discussed.

A convenient incorporation of conformationally constrained 5,5-dimethylproline into the ribonuclease A 89-124 sequence by condensation of synthetic peptide fragments.

The presence of l-5,5-dimethylproline (dmP) within an amino acid sequence results in the formation of an X-dmP peptide bond predominantly locked in a cis conformation. However, the common use of this unnatural amino acid has been hampered by the difficulty of the economical incorporation of the dmP residue into longer peptide segments due to the steric hindrance imposed by the dimethyl moieties. Here, we describe synthesis of the C-terminal 36-residue peptide, corresponding to the 89-124 sequence of bovine pancreatic ribonuclease A (RNase A), in which dmP is incorporated as a substitute for Pro93. The peptide was assembled by condensation of protected 5- and 31-residue peptide fragments, which were synthesized by solid-phase peptide methodology using fluorenylmethyloxycarbonyl chemistry. We focused on optimizing the synthesis of the Fmoc-Ser(tBu)-Ser(tBu)-Lys(Boc)-Tyr(tBu)-dmP-OH pentapeptide (residues 89-93) with efficient acylation of the sterically hindered dmP residue. In a comparative study, the reagent O-(7-azabenzotriazol-1-yl)-1,1,3,3-tetramethyluronium hexafluorophosphate was found to be superior to bromo-tris-pyrrolidino-phosphonium hexafluorophosphate and tetramethylfluoroformamidinium hexafluorophosphate for the synthesis of the -Tyr(tBu)-dmP- peptide bond in solution as well as on a resin.

Ab initio folding of multiple-chain proteins.

Our previous methodology for ab initio prediction of protein structure is extended here to treat multiple-chain proteins. This involved modification of our united-residue (UNRES) force field and our Conformational Space Annealing (CSA) Global Optimization procedure. Good results have been obtained for both a four- and a three-helix protein from the CASP3 exercise.

Folding of a disulfide-bonded protein species with free thiol(s): competition between conformational folding and disulfide reshuffling in an intermediate of bovine pancreatic ribonuclease A.

The conformational folding of the nativelike intermediate des-[40-95] on the major oxidative folding pathway of bovine pancreatic ribonuclease A (RNase A) has been examined at various pHs and temperatures in the absence of a redox reagent. Des-[40-95] has three of the four disulfide bonds of native RNase A and lacks the bond between Cys40 and Cys95. This three-disulfide species was unfolded at low pH to inhibit any disulfide reshuffling and was refolded at higher pH, allowing both conformational folding and disulfide-reshuffling reactions to take place. As a result of this competition, 15-85% of des-[40-95], depending on the experimental conditions, undergoes intramolecular disulfide-reshuffling reactions. That portion of the des-[40-95] population which has native isomers of essential proline residues appears to fold faster than the disulfide reaction can occur. However, when the folding is retarded, conceivably by the presence of non-native isomers of essential proline residues, des-[40-95] may reshuffle before completing the conformational folding process. These results enable us to distinguish among current models for the critical structure-forming step in oxidative folding and reveal a new model for coupling proline isomerization to disulfide-bond formation. These experiments also demonstrate that the reshuffling-folding competition assay is a useful tool for detecting structured populations in conformational folding intermediates.

Structure of the interferon-receptor complex determined by distance constraints from double-mutant cycles and flexible docking.

The pleiotropic activity of type I interferons has been attributed to the specific interaction of IFN with the cell-surface receptor components ifnar1 and ifnar2. To date, the structure of IFN has been solved, but not that of the receptor or the complex. In this study, the structure of the IFN-alpha 2-ifnar2 complex was generated with a docking procedure, using nuclear Overhauser effect-like distance constraints obtained from double-mutant cycle experiments. The interaction free energy between 13 residues of the ligand and 11 of the receptor was measured by double-mutant cycles. Of the 100 pairwise interactions probed, five pairs of residues were found to interact. These five interactions were incorporated as distance constraints into the flexible docking program prodock by using fixed and movable energy-gradient grids attached to the receptor and ligand, respectively. Multistart minimization and Monte Carlo minimization docking of IFN-alpha 2 onto ifnar2 converged to a well-defined average structure, with the five distance constraints being satisfied. Furthermore, no structural artifacts or intraloop energy strain were observed. The mutual binding sites on IFN-alpha 2 and ifnar2 predicted from the model showed an almost complete superposition with the ones determined from mutagenesis studies. Based on this structure, differences in IFN-alpha 2 versus IFN-beta binding are discussed.

Conformation-family Monte Carlo: a new method for crystal structure prediction.

A new global optimization method, Conformation-family Monte Carlo, has been developed recently for searching the conformational space of macromolecules. In the present paper, we adapted this method for prediction of crystal structures of organic molecules without assuming any symmetry constraints except the number of molecules in the unit cell. This method maintains a database of low energy structures that are clustered into families. The structures in this database are improved iteratively by a Metropolis-type Monte Carlo procedure together with energy minimization, in which the search is biased toward the regions of the lowest energy families. The Conformation-family Monte Carlo method is applied to a set of nine rigid and flexible organic molecules by using two popular force fields, AMBER and W99. The method performed well for the rigid molecules and reasonably well for the molecules with torsional degrees of freedom.

Coupling of conformational folding and disulfide-bond reactions in oxidative folding of proteins.

The oxidative folding of proteins consists of conformational folding and disulfide-bond reactions. These two processes are coupled significantly in folding-coupled regeneration steps, in which a single chemical reaction (the "forward" reaction) converts a conformationally unstable precursor species into a conformationally stable, disulfide-protected successor species. Two limiting-case mechanisms for folding-coupled regeneration steps are described. In the folded-precursor mechanism, the precursor species is preferentially folded at the moment of the forward reaction. The (transient) native structure increases the effective concentrations of the reactive thiol and disulfide groups, thus favoring the forward reaction. By contrast, in the quasi-stochastic mechanism, the forward reaction occurs quasi-stochastically in an unfolded precursor; i.e., reactive groups encounter each other with a probability determined primarily by loop entropy, albeit modified by conformational biases in the unfolded state. The resulting successor species is initially unfolded, and its folding competes with backward chemical reactions to the unfolded precursors. The folded-precursor and quasi-stochastic mechanisms may be distinguished experimentally by the dependence of their kinetics on factors affecting the rates of thiol--disulfide exchange and conformational (un)folding. Experimental data and structural and biochemical arguments suggest that the quasi-stochastic mechanism is more plausible than the folded-precursor mechanism for most proteins.

Effect of mutation of proline 93 on redox unfolding/folding of bovine pancreatic ribonuclease A.

Both the reductive unfolding and oxidative regeneration of a P93A mutant and wild-type RNase A have been studied at 15 degrees C and pH 8.0. The rate of reduction of the 40--95 disulfide bond is accelerated about 120-fold by the P93A mutation, while the reduction of the 65--72 disulfide bond is not accelerated by this mutation (within the experimental error). Moreover, the reduction of native P93A to des[40--95] is about 10 times faster than the further reduction of the same des[40--95] species. These results demonstrate that the reduction of the mutant proceeds through a local unfolding event and provides strong support for our model in which the reduction of wild-type RNase A to the des species proceeds through two independent local conformational unfolding events. The oxidative regeneration rate of the P93A mutant is comparable to that of wild-type RNase A, suggesting that a cis 92--93 peptide group that is present in native wild-type RNase A and in native des[40--95], is not obligatory for the formation of the third (final) native disulfide bond of des[40--95] by reshuffling from an unstructured 3S precursor. Thus, the trans to cis isomerization of the Tyr92-Pro93 peptide group during the regeneration of wild-type RNase A may occur after the formation of the third native disulfide bond.

A role for intermolecular disulfide bonds in prion diseases?

The key event in prion diseases seems to be the conversion of the prion protein PrP from its normal cellular isoform (PrP(C)) to an aberrant "scrapie" isoform (PrP(Sc)). Earlier studies have detected no covalent modification in the scrapie isoform and have concluded that the PrP(C) --> PrP(Sc) conversion is a purely conformational transition involving no chemical reactions. However, a reexamination of the available biochemical data suggests that the PrP(C) --> PrP(Sc) conversion also involves a covalent reaction of the (sole) intramolecular disulfide bond of PrP(C). Specifically, the data are consistent with the hypothesis that infectious prions are composed of PrP(Sc) polymers linked by intermolecular disulfide bonds. Thus, the PrP(C) --> PrP(Sc) conversion may involve not only a conformational transition but also a thiol/disulfide exchange reaction between the terminal thiolate of such a PrP(Sc) polymer and the disulfide bond of a PrP(C) monomer. This hypothesis seems to account for several unusual features of prion diseases.