Search PubMed⌕ Search

Biomedical subjects

H A Scheraga

Publications and source records attributed to H A Scheraga.

At least 19 recordsLinked to original sources
PubMed2005

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.

Algorithms↗
PubMed2005

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.

Amino Acid Sequence↗
PubMed2005

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.

Animals↗
PubMed2005

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.

Peptide Fragments↗
PubMed2004

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.

Algorithms↗
PubMed2003

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.

Amino Acids↗
PubMed2002

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.

Dimerization↗
PubMed2001

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.

Animals↗
PubMed2001

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.

Interferon-alpha↗
PubMed2001

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.

Crystallization↗
PubMed2001

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.

Disulfides↗
PubMed2001

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.

Alanine↗
PubMed2001

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.

Disulfides↗
PubMed2001

Structural determinants of oxidative folding in proteins.

A method for determining the kinetic fate of structured disulfide species (i.e., whether they are preferentially oxidized or reshuffle back to an unstructured disulfide species) is introduced. The method relies on the sensitivity of unstructured disulfide species to low concentrations of reducing agents. Because a structured des species that preferentially reshuffles generally first rearranges to an unstructured species, a small concentration of reduced DTT (e.g., 260 microM) suffices to distinguish on-pathway intermediates from dead-end species. We apply this method to the oxidative folding of bovine pancreatic ribonuclease A (RNase A) and show that des[40-95] and des[65-72] are productive intermediates, whereas des[26-84] and des[58-110] are metastable dead-end species that preferentially reshuffle. The key factor in determining the kinetic fate of these des species is the relative accessibility of both their thiol groups and disulfide bonds. Productive intermediates tend to be disulfide-secure, meaning that their structural fluctuations preferentially expose their thiol groups, while keeping their disulfide bonds buried. By contrast, dead-end species tend to be disulfide-insecure, in that their structural fluctuations expose their disulfide bonds in concert with their thiol groups. This distinction leads to four generic types of oxidative folding pathways. We combine these results with those of earlier studies to suggest a general three-stage model of oxidative folding of RNase A and other single-domain proteins with multiple disulfide bonds.

Animals↗
PubMed2001

Recent improvements in prediction of protein structure by global optimization of a potential energy function.

Recent improvements of a hierarchical ab initio or de novo approach for predicting both alpha and beta structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic alpha-helical protein of 70 residues was predicted to within 4.3 A alpha-carbon (C(alpha)) rms deviation (rmsd) whereas, for other alpha-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 A C(alpha) rmsd. Whereas beta structures can now be predicted with the new procedure, the success rate for alpha/beta- and beta-proteins is lower than that for alpha-proteins at present. For the beta portions of alpha/beta structures, the C(alpha) rmsd's are less than 6.0 A for contiguous fragments of 30-40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a C(alpha) rmsd less than 6.0 A. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Models, Molecular↗
PubMed2001

Local and long-range interactions in the thermal unfolding transition of bovine pancreatic ribonuclease A.

This research was undertaken to distinguish between local and global unfolding in the reversible thermal denaturation of bovine pancreatic ribonclease A (RNase A). Local unfolding was monitored by steady-state and time-resolved fluorescence of nine mutants in each of which a single tryptophan was substituted for a wild-type residue. Global unfolding was monitored by far-UV circular dichroism and UV absorbance. All the mutants (except F8W and D38W) exhibited high specific enzymatic activity, and their far-UV CD spectra were very close to that of wild-type RNase A, indicating that the tryptophan substitutions did not affect the structure of any of the mutants (excluding K1W and Y92W) under folding conditions at 20 degrees C. Like wild-type RNase A, the various mutants exhibited reversible cooperative thermal unfolding transitions at pH 5, with transition temperatures 2.5-11 degrees C lower than that of the wild-type transition, as detected by far-UV CD or UV absorbance. Even at 80 degrees C, well above the cooperative transition of all the RNase A mutants, a considerable amount of secondary and tertiary structure was maintained. These studies suggest the following two-stage mechanism for the thermal unfolding transition of RNase A as the temperature is increased. First, at temperatures lower than those of the main cooperative transition, long-range interactions within the major hydrophobic core are weakened, e.g., those involving residues Phe-8 (in the N-terminal helix) and Lys-104 and Tyr-115 (in the C-terminal beta-hairpin motif). The structure of the chain-reversal loop (residues 91-95) relaxes in the same temperature range. Second, the subsequent higher-temperature cooperative unfolding transition is associated with a loss of secondary structure and additional changes in the tertiary contacts of the major hydrophobic core, e.g., those involving residues Tyr-73, Tyr-76, and Asp-38 on the other side of the molecule. The hydrophobic interactions of the C-terminal loop of the protein are enhanced by high temperature, and perhaps are responsible for the preservation of the local structural environment of Trp-124 at temperatures slightly above the major cooperative transition. The results shed new light on the thermal unfolding transitions, generally supporting the thermal unfolding hypothesis of Burgess and Scheraga, as modified by Matheson and Scheraga.

Amino Acid Substitution↗
PubMed2001

Distributions of intramolecular distances in the reduced and denatured states of bovine pancreatic ribonuclease A. Folding initiation structures in the C-terminal portions of the reduced protein.

The purpose of this investigation is to characterize the reduced state of RNase A (r-RNase A) in terms of (i) intramolecular distances, (ii) the sequence of formation of stable loops in the initial stages of folding, and (iii) the unfolding transitions induced by GdnHCl. This is accomplished by identifying specific subdomain structures and local and long-range interactions that direct the folding process of this protein and lead to the native fold and formation of the disulfide bonds. Eleven pairs of dispersed sites in the RNase A molecule were labeled with fluorescent donor and acceptor probes, and the distributions of intramolecular distances (IDDs) were determined by means of time-resolved dynamic nonradiative excitation energy transfer (TR-FRET) measurements. The mutants were designed to search for (a) a possible nonrandom fold of the backbone in the collapsed state and (b) possible loops stabilized by long-range interactions. It was found that, under folding conditions, (i) the labeled mutants of r-RNase A in refolding buffer (the R(N) state) exhibit features of specific (nonrandom) compact but very dispersed subdomain structures (indicated by short mean distances, broad IDDs, and a weak dependence of the mean distances on segment length), (ii) the backbone fold in the C-terminal beta-like portion of the molecule appears to adopt a native-like overall fold, (iii) the N-terminal alpha-like portion of the chain is separated from the C-terminal core by very large intramolecular distances, larger than those in the crystal structure, and (iv) perturbations by addition of GdnHCl reveal several conformational transitions in different sections of the chain. Addition of GdnHCl to the native disulfide-intact protein provided a reference state for the extent of expansion of intramolecular distances under denaturing conditions. In conclusion, r-RNase A under folding conditions (the R(N) state) is poised for the final folding step(s) with a native-like trace of the chain fold but a large separation between the two subdomains which is then decreased upon introduction of three of the four native disulfide cross-links.

Animals↗