SEARCH · PubMed Health

Results for “conformational ensemble”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46Linked to original sources

PubMed2000

Change in backbone torsion angle distribution on protein folding.

Understanding protein folding requires the determination of the configurational space accessible to the protein at different stages in folding. Here, computer simulation analysis of small angle neutron scattering results is used to probe the change in the distribution of configurations on strong denaturation of a globular protein, phosphoglycerate kinase, in 4 M guanidine hydrochloride solution. To do this atomic-detail ensembles of the unfolded protein chain are modeled and their scattering profiles compared with the experiment. The local conformational statistics are found to strongly influence the experimental intensity at scattering vectors between 0.05 and 0.3 A(-1). Denaturation leads to a reduction in the protein atom-pair distance distribution function over the approximately 3-15 A region that is associated with a quantifiable shift in the backbone torsional angle (phi, psi) distribution toward the beta region of the Ramachandran plot.

Computer Simulation↗

PubMed2006

On the characterization of protein native state ensembles.

Describing and understanding the biological function of a protein requires a detailed structural and thermodynamic description of the protein's native state ensemble. Obtaining such a description often involves characterizing equilibrium fluctuations that occur beyond the nanosecond timescale. Capturing such fluctuations remains nontrivial even for very long molecular dynamics and Monte Carlo simulations. We propose a novel multiscale computational method to exhaustively characterize, in atomistic detail, the protein conformations constituting the native state with no inherent timescale limitations. Applications of this method to proteins of various folds and sizes show that thermodynamic observables measured as averages over the native state ensembles obtained by the method agree remarkably well with nuclear magnetic resonance data that span multiple timescales. By characterizing equilibrium fluctuations at atomistic detail over a broad range of timescales, from picoseconds to milliseconds, our method offers to complement current simulation techniques and wet-lab experiments and can impact our understanding and description of the relationship between protein flexibility and function.

Amino Acids↗

PubMed2005

Determination of the folding transition states of barnase by using PhiI-value-restrained simulations validated by double mutant PhiIJ-values.

The protein barnase folds from the denatured state into its native conformation via a high-energy intermediate. Using PhiI-values determined experimentally from single-point mutations as restraints in all-atom molecular dynamics simulations, we have determined ensembles of structures corresponding to the transition states for the formation of the folding intermediate and its conversion into the native state. We have also introduced a stringent validation of the approach used to calculate such structures by predicting interaction PhiIJ-values determined experimentally from a series of double-mutant cycles. The ensembles that we have obtained illustrate the heterogeneity in the nucleation-condensation process by which barnase folds. Obligatory steps of this process include the sequential formation of two folding nuclei, which correspond to the two main hydrophobic cores of the protein. Nonobligatory steps include the elongation of the strand beta1 and the formation of the helix alpha2. The results confirm that the use of experimental observables such as PhiI-values as restraints in molecular dynamics simulations is a powerful general strategy to characterize the relatively heterogeneous structural ensembles that populate nonnative regions of the energy landscapes of proteins.

Bacterial Proteins↗

PubMed1993

Hevein: NMR assignment and assessment of solution-state folding for the agglutinin-toxin motif.

The first high-resolution solution-state structure of a member of the toxin-agglutinin folding motif with the WGA disulfide linkage is presented. The 1H NMR spectrum of hevein has been 100% assigned from residue 2 through residue 43, the C-terminus, using two-dimensional correlation and NOE spectroscopy. During the course of the NOESY analysis, the three-dimensional structural features of hevein were derived, using nonstereospecific distance constraints (with tight bounds) for XPLOR simulated annealing followed by unconstrained relaxation in the CHARMm force field, at two levels of long-range constraint density. In addition, a large number of low-bound-only constraints, corresponding to unobserved NOE's, were used in both refinements. The first structure elucidation employed a total of 180 distance constraints (60 of which were medium or long range, i/i+n with n < or = 2). The second refinement employed 244 (101 medium or long range) constraints: some conformation-insensitive intraresidue constraints were deleted, two misassigned long-range constraints were corrected, and 41 new i/i+n (n > or = 2) constraints were added. The average bounds precisions of the two refinements were comparable (+/- 0.44 A) and significantly tighter than those that result when a universal low bound corresponding to the sum of the van der Waals radii was used. (The more conservative treatment of NOE's gave the same final structure but required a higher constraint density before assignment errors would stand out during the refinement.) Constraint density also has a significant influence on convergence and accuracy using tight constraints. The study demonstrates that convergence within an ensemble of solution structures is not a dependable criterion for either the accuracy or precision of the derived structure. The best fitting conformers from the refinement at the higher constraint density bear a greater similarity to the solid-state structure of the domains of wheat germ agglutinin (0.95 A rmsd over residues 2-32) than to the recently reported 2.8-A X-ray structure of hevein (1.25 A rmsd over residues 2-32, 2.83 A rmsd over residues 2-42). The consensus conformer from the solution data is defined to a backbone rmsd of < 0.6 A over the full sequence for which NMR data could be collected.(ABSTRACT TRUNCATED AT 400 WORDS)

Amino Acid Sequence↗

PubMedDate not supplied

Ensemble docking into flexible active sites. Critical evaluation of FlexE against JNK-3 and beta-secretase.

One of the main complicating factors in structure-based drug design is the conformational rearrangement of the receptor upon ligand binding implicating protein flexibility as a crucial component in virtual screening. The FlexE approach allows flexibility through discrete alternative conformations of varying parts of the protein taken from structures having similar backbone traces. Here the performance of FlexE was tested against that of FlexX and FlexX-Pharm, by carrying out virtual screening experiments on two sets of structurally distinct complexes, for the enzymes beta-secretase (BACE), and c-jun N-terminal kinase 3 (JNK-3). A large number of incompatible instances occurred between structural elements of the proteins thus loop movements could not be studied in JNK-3 as well as in BACE. The investigation of the side-chain flexibility revealed that at the most FlexE could achieve the enrichment yielded by FlexX in JNK-3 but not in BACE. Although limited side-chain variations (e.g. different protonation states) can be treated by FlexE, docking into protein ensembles remains a practical tool that decreases the average run time for a ligand.

Amyloid Precursor Protein Secretases↗

PubMed1996

H-bond stability in the tRNA(Asp) anticodon hairpin: 3 ns of multiple molecular dynamics simulations.

Multiple molecular dynamics trajectories of the solvated and neutralized 17-residue tRNA(Asp) anticodon hairpin were generated for a total of 3 ns. Explicit treatment of all long-ranged electrostatic interactions by the particle mesh Ewald algorithm, as implemented in the AMBER MD software package, effected a degree of structural stabilization not previously achieved by use of a long 16-A solvent interaction truncation scheme. The increased stability of this multiple molecular dynamics set was appropriate for an in-depth analysis of the six 500-ps-long trajectories and allowed the characterization of a number of key structural interactions. The dynamical behavior of the standard Watson-Crick base pairs, the noncanonical G30-U40 "wobble" base pair, and the psi 32-C38 pseudo-base pair is presented as well as that of two C--H... O hydrogen bonds found to contribute to the array of tertiary interactions that stabilize the seven-nucleotide native loop conformation. The least mobile residue in the loop is U33, which forms the U-turn motif and which participates in several hydrogen-bonding interactions, whereas the most mobile residue is the apical residue G34 at the wobble position, a factor undoubtedly important in its biological function. The set of multiple molecular dynamics trajectories obtained does not converge on a 500-ps time scale to a unique dynamical model but instead describes an ensemble of structural microstates accessible to the system under the present simulation protocol, which is the result of local structural heterogeneity rather than of global conformational changes.

Anticodon↗

PubMed2006

Complementarity of ensemble and single-molecule measures of protein motion: a relaxation dispersion NMR study of an enzyme complex.

Single-molecule fluorescence experiments have shown that the conformation of the complex between Escherichia coli general NAD(P)H:flavin oxidoreductase (FRE) and flavin adenine dinucleotide (FAD) fluctuates over a range of timescales between 10(-4) and 1 s. Here we use (15)N and (13)C relaxation dispersion NMR methods to study millisecond-timescale dynamics in the complex. In this time regime, the protein is extremely flexible, with residues that undergo conformational exchange located throughout the molecule. Three distinct regions of dynamics are quantified, with two of them involving residues making contact to the donor (Tyr-35) and acceptor (FAD) sites that participate in the electron transfer reaction monitored in single-molecule experiments. Modulation of the donor-acceptor distance through these conformational exchange processes, occurring with rates of approximately 400 and 1,200 s(-1) (22 degrees C), affects the rate of electron transfer and partially accounts for the range of the observed dynamics monitored in the fluorescence experiments.

Biochemistry↗

PubMed2006

Unfolding transition state and intermediates of the tumor suppressor p16INK4a investigated by molecular dynamics simulations.

The ankyrin repeat is one of the most common protein motifs and is involved in protein-protein interactions. It consists of 33 residues that assume a beta-hairpin helix-loop-helix fold. Mutagenesis and kinetic experiments (Phi-value analysis of the folding transition state) have shown that the tumor suppressor p16(INK4a), a four-repeat protein, unfolds sequentially starting from the two N-terminal repeats. Here, the flexibility of p16(INK4a) at room temperature and its unfolding mechanism at high temperature have been investigated by multiple molecular dynamics runs in explicit water for a total simulation time of 0.65 micros. The transition state ensemble (TSE) of p16(INK4a) was identified by monitoring both the deviation from the experimental Phi values and sudden conformational changes along the unfolding trajectories. Conformations in the TSE have a mainly unstructured second repeat whereas the other repeats are almost completely folded. A rigid-body displacement of the first repeat involving both a rotation and translation is observed in all molecular dynamics simulations at high temperature. The Trp(15), Pro(75), and Ala(76) side-chains are more buried in the TSE than the native state. The sequential unfolding starting at the second repeat is in agreement with the mutagenesis studies whereas the displacement of the first repeat and the presence of nonnative interactions at the TSE are simulation results which supplement the experimental data. Furthermore, the unfolding trajectories reveal the presence of two on-pathway intermediates with partial alpha-helical structure. Finally, on the basis of the available experimental and simulation results we suggest that in modular proteins the shift of the folding TSE toward the native structure upon reduction of the number of tandem repeats is consistent with the Hammond effect.

Binding Sites↗

PubMed2005

A structural model for the adduct between cytochrome c and cytochrome c oxidase.

An ensemble of structural models of the adduct between cytochrome c and cytochrome c oxidase from Paracoccus denitrificans has been calculated based on the experimental data from site-directed mutagenesis and NMR experiments that have accumulated over the last years of research on this system. The residues from each protein that are at the protein-protein interface have been identified by the above experimental work, and this information has been converted in a series of restraints explicitly used in calculations. It is found that a single static structural model cannot satisfy all experimental data simultaneously. Therefore, it is proposed that the adduct exists as a dynamic ensemble of different orientations in equilibrium, and may be represented by a combination or average of the various limiting conformations calculated here. The equilibrium involves both conformations that are competent for electron transfer and conformations that are not. Long-range recognition of the partners is driven by non-specific electrostatic interactions, while at shorter distances hydrophobic contacts tune the reciprocal orientation. Electron transfer from cytochrome bc (1) to cytochrome c oxidase is mediated through cytochrome c experiencing multiple encounters with both of its partners, only part of which are productive. The number of encounters, and thus the electron transfer rate, may be increased by the formation of a cytochrome bc (1)-cytochrome c oxidase supercomplex and/or (in human) by increasing the concentration of the two enzymes in the membrane space.

Cytochromes c↗

PubMed2003

Role of the [65-72] disulfide bond in oxidative folding of bovine pancreatic ribonuclease A.

To assess the role of the [65-72] disulfide bond in the oxidative folding of RNase A, use has been made of [C65S, C72S], a three-disulfide-containing mutant of RNase A which regenerates from its two-disulfide precursor in an oxidation and conformational folding-coupled rate-determining step. The distribution of disulfide bonds in the one-disulfide-containing ensemble of this mutant has been characterized. In general, the disulfide-bond distribution in its 1S ensemble agrees relatively well with the corresponding distribution in wt-RNase A and with distributions based on calculations of loop entropy, except for the absence of the [65-72] disulfide bond. There is no bias (over the entropic influence) for the three native disulfide bonds, [26-84], [40-95], and [58-110]. Previous oxidative folding results for wt-RNase A indicated the predominance of the des [40-95] intermediate over des [65-72] after the rate-determining step in the regeneration process. Considering that there is no preferential distribution of disulfides in the 1S ensemble of [C65S, C72S], in contrast to the preferential population of the [65-72] disulfide bond in wt-RNase A, these results indicate a critical role for the [65-72] disulfide bond in the regeneration of wt-RNase A. Furthermore, analysis of the disulfide distribution of the 1S intermediates of [C65S, C72S] compared to that of wt-RNase A lends support for a physicochemical basis for the previously observed slow folding rate of this mutant, compared to its analogue (des [65-72]) of wt-RNase A.

Amino Acid Substitution↗

PubMed1977

Model of protein folding: incorporation of a one-dimensional short-range (Ising) model into a three-dimensional model.

In this paper, we have incorporated a one-dimensional short-range model into a three-dimensional model for protein folding. It has been applied, by extending the concept of the three-step mechanism for protein folding proposed in our previous paper, to simulate the folding of bovine pancreatic trypsin inhibitor, using a Monte Carlo procedure in all three steps, A, B, and C. The statistical mechanical ensemble treatment of the short-range model serves as a constraint on the Monte Carlo procedure, in which conformational transitions are introduced. The preliminary results of 10 independent Monte Carlo trials indicate that, while folding is achieved, improvements are required in order to account for the correct three-dimensional structure of a globular protein.

Methods↗

PubMed2001

The pair-functional method for direct solution of molecular structures. I. Statistical principles.

The pair-functional principle shows how to construct a unique statistical ensemble of strongly interacting atoms that corresponds to any feasible measured set of X-ray intensities. The ensemble and all its distribution functions are strictly periodic in the crystal lattice, so that each unit cell has exactly the same arrangement of atoms at all times. The mean particle density in the cell is uniform because the ensemble has undefined phases and the origin is not fixed. The atoms in this maximum-entropy ensemble interact through pairwise additive periodic statistical forces within the unit cell. The ensemble average pair-correlation function is matched to the observed originless Patterson function of the crystal. The derived pairing force then becomes approximately proportional to the Ornstein-Zernicke direct correlation function of the ensemble. The atoms have a many-body Boltzmann distribution and the logarithm of the likelihood of any particular conformation is related to its total pairing potential. The pairing potential of a group of atoms acts like a local field in the cell. This property is used in the pair-functional method. Molecular structures can be solved by a direct search in real space for clusters of atoms with high pair potentials. During a successful search, the atoms move from their original random positions to form larger and larger clusters of correctly formed fragments. Finally, every atom belongs to a single cluster, which is the correct solution.

Journal Article↗

PubMed2000

A rapid computational method for lead evolution: description and application to alpha(1)-adrenergic antagonists.

The high failure rate of drugs in the development phase requires a strategy to reduce risks by generating lead candidates from different chemical classes. We describe a new three-dimensional computational approach for lead evolution, based on multiple pharmacophore hypotheses. Using full conformational models for both active and inactive compounds, a large number of pharmacophore hypotheses are analyzed to select the set or "ensemble" of hypotheses that, when combined, is most able to discriminate between active and inactive molecules. The ensemble hypothesis is then used to search virtual chemical libraries to identify compounds for synthesis. This method is very rapid, allowing very large virtual libraries on the order of a million compounds to be filtered efficiently. In applying this method to alpha(1)-adrenergic receptor ligands, we have demonstrated lead evolution from heterocyclic alpha(1)-adrenergic receptor ligands to highly dissimilar active N-substituted glycine compounds. Our results also show that the active N-substituted glycines are part of our smaller filtered library and thus could have been identified by synthesizing only a portion of the N-substituted glycine library.

Adrenergic alpha-Antagonists↗

PubMed2000

Interhelical ion pairing in coiled coils: solution structure of a heterodimeric leucine zipper and determination of pKa values of Glu side chains.

Residues of opposite charge often populate heptad positions g (heptad i on chain 1) and e' (heptad i + 1 on chain 2) in dimeric coiled coils and may stabilize the dimer by formation of interchain ion pairs. To investigate the contribution to stability of such electrostatic interactions we have designed a disulfide-linked heterodimeric zipper (AB zipper) consisting of the acidic chain Ac-E-VAQLEKE-VAQAEAE-NYQLEQE-VAQLEHE-CG-NH(2) and the basic chain Ac-E-VQALKKR-VQALKAR-NYAAKQK-VQALRHK-CG-NH(2) in which all e and g positions are occupied by either E or K/R to form a maximum of seven interhelical salt bridges. Temperature-induced denaturation experiments monitored by circular dichroism reveal a stable coiled coil conformation below 50 degrees C and in the pH range 1.2-8.0. Stability is highest at pH approximately 4.0 [DeltaG(U) (37 degrees C) = 5.18 +/- 0.51 kcal mol(-)(1)]. The solution structure of the AB zipper at pH 5.65 has been elucidated on the basis of homonuclear (1)H NMR data collected at 800 MHz [heavy atom rmsd's for the ensemble of 50 calculated structures are 0.47 +/- 0.13 A (backbone) and 0.95 +/- 0.16 A (all)]. Both chains of the AB zipper are almost entirely in alpha-helical conformation and form a superhelix with a left-handed twist. Overhauser connectivities reveal close contacts between g position residues (heptad i on chain 1) and residues d/f (heptad i on chain 1), residues a/d (heptad i + 1 on chain 1), and residue a' (heptad i + 1 on chain 2). Residues in position e (heptad i on chain 1) are in contact with residues a/b/d/f (heptad i on chain 1) and residue d' (heptad i on chain 2). These connectivities hint at a relatively defined alignment of the side chains across the helix interface. Partial H-bond formation between the functional groups of residues g and e'(+1) is observed in the calculated structures. NMR pH titration experiments disclose pK(a) values for Glu delta-carboxylate groups: 4.14 +/- 0.02 (E(1)), 4.82 +/- 0.07 (E(6)), 4.52 +/- 0.01 (E(8)), 4.37 +/- 0.03 (E(13)), 4.11 +/- 0.02 (E(15)), 4.41 +/- 0.07 (E(20)), 4.82 +/- 0.03 (E(22)), 4.65 +/- 0.04 (E(27)), 4.63 +/- 0.03 (E(29)), 4.22 +/- 0.02 (E(1)(')). By comparison with pK(a) of Glu in unfolded peptides ( approximately 4. 3 +/- 0.1), our pK(a) data suggest marginal or even unfavorable contribution of charged Glu to the stability of the AB zipper. The electrostatic energy gained from interhelical ion pairs is likely to be surpassed by hydrophobic energy terms upon protonation of Glu, due to increased hydrophobicity of uncharged Glu and, thus, better packing against apolar residues at the chain interface.

Amino Acid Sequence↗

PubMed2006

Visualization of transient encounter complexes in protein-protein association.

Kinetic data on a number of protein-protein associations have provided evidence for the initial formation of a pre-equilibrium encounter complex that subsequently relaxes to the final stereospecific complex. Site-directed mutagenesis and brownian dynamics simulations have suggested that the rate of association can be modulated by perturbations in charge distribution outside the direct interaction surfaces. Furthermore, rate enhancement through non-specific binding may occur by either a reduction in dimensionality or the presence of a short-range, non-specific attractive potential. Here, using paramagnetic relaxation enhancement, we directly demonstrate the existence and visualize the distribution of an ensemble of transient, non-specific encounter complexes under equilibrium conditions for a relatively weak protein-protein complex between the amino-terminal domain of enzyme I and the phosphocarrier protein HPr. Neither the stereospecific complex alone nor any single alternative conformation can account fully for the intermolecular paramagnetic relaxation enhancement data. Restrained rigid-body simulated annealing refinement against the paramagnetic relaxation enhancement data enables us to obtain an atomic probability distribution map of the non-specific encounter complex ensemble that qualitatively correlates with the electrostatic surface potentials on the interacting proteins. Qualitatively similar results are presented for two other protein-protein complexes.

Bacterial Proteins↗

PubMed2000

Conformational dynamics in mixed alpha/beta-oligonucleotides containing polarity reversals: a molecular dynamics study using time-averaged restraints.

Nucleic acid duplexes featuring a single alpha-anomeric thymidine inserted into each DNA strand via 3'-3' and 5'-5' phosphodiester linkages exhibit local conformational dynamics that are not adequately depicted by conventional restrained molecular dynamics (rMD) methods. We have used molecular dynamics with time-averaged NMR restraints (MDtar) to explore its applicability to describing the conformational dynamics of two alpha-containing duplexes--d(GCGAAT-3'-3'-alphaT-5'-5'-CGC)2 and d(ATGG-3'-3'-alphaT-5'-5'-GCTC) x r(gagcaccau). In contrast to rMD, enforcing NOE-based distance restraints over a period of time in MDtar rather than instantaneously results in better agreement with the experimental NOE and J-data. This conclusion is based on the dramatic decreases in average distance and coupling constant violations (delta d(av), J(rms), and delta J(av)) and improvements in sixth-root R-factors (R(X)). In both duplexes, the deoxyribose ring puckering behavior predicted independently by pseudorotation analysis is portrayed remarkably well using this approach compared to rMD. This indicates that the local dynamic behavior is encoded within the NOE data, although this is not obvious from the local R(X) values. In both systems, the backbone torsion angles comprising the 3'-3' linkage as well as the (high S-) sugars of the alpha-nucleotide and preceding residue (alpha - 1) are relatively static, while the conformations of the 5'-5' linkage and the sugar in the neighboring beta-nucleotide (alpha + 1) show enhanced flexibility. To reduce the large ensembles generated by MDtar to more manageable clusters we utilized the PDQPRO program. The resulting PDQPRO clusters (in both cases, 13 structures and associated probabilities extracted from a pool of 300 structures) adequately represent the structural and dynamic characteristics predicted by the experimental data.

Base Pair Mismatch↗

PubMed2001

Carrier protein-modulated presentation and recognition of an N-glycan: observations on the interactions of Man(8) glycoform of ribonuclease B with conglutinin.

Conglutinin is a serum lectin of the innate immune system, which binds high mannose N-glycans when these are appropriately presented on proteins. Here we use the conglutinin-ribonuclease B (RNaseB)-recognition system as a model to investigate the structural basis of selective recognition of protein-bound oligosaccharides by this carbohydrate-binding receptor. Conglutinin shows little binding to the isolated RNaseB-Man(8 )glycoform, and no binding to Man(5-6) glycoforms. In contrast, when the protein moiety is reduced and denatured we observe that conglutinin binds strongly to the isolated RNaseB-Man(8) glycoform and weakly to the Man(5-6) glycoforms. These results are in accord with observations on the binding to the N-glycans in the absence of carrier protein. NMR analyses of native RNaseB-Man(8) and -Man(5-6) glycoforms reveal that the three-dimensional structure of the protein moiety is essentially identical to that of non-glycosylated RNase (RNaseA). Thus there are no perceptible differences between the RNase protein forms that could account for differential availability of the N-glycan for conglutinin-binding. After reduction and denaturation, the NMR spectrum became typical of a non-structured polypeptide, although the conformational preferences of the N-glycosidic linkage were unchanged, and most importantly, the Man(8 )oligosaccharide retained the average conformational behavior of the free oligosaccharide irrespective of the carrier protein fold. This conformational freedom is clearly not translated into full availability of the oligosaccharide for the carbohydrate-recognition protein. We propose, therefore, that the differing bioactivity of the N-glycan is a reflection of the existence of different geometries of presentation of the carbohydrate determinant in relation to the protein surface within the glycan:carrier protein ensemble.

Carbohydrate Conformation↗

PubMed2000

Charge-charge interactions influence the denatured state ensemble and contribute to protein stability.

Several recent studies have shown that it is possible to increase protein stability by improving electrostatic interactions among charged groups on the surface of the folded protein. However, the stability increases are considerably smaller than predicted by a simple Coulomb's law calculation, and in some cases, a charge reversal on the surface leads to a decrease in stability when an increase was predicted. These results suggest that favorable charge-charge interactions are important in determining the denatured state ensemble, and that the free energy of the denatured state may be decreased more than that of the native state by reversing the charge of a side chain. We suggest that when the hydrophobic and hydrogen bonding interactions that stabilize the folded state are disrupted, the unfolded polypeptide chain rearranges to compact conformations with favorable long-range electrostatic interactions. These charge-charge interactions in the denatured state will reduce the net contribution of electrostatic interactions to protein stability and will help determine the denatured state ensemble. To support this idea, we show that the denatured state ensemble of ribonuclease Sa is considerably more compact at pH 7 where favorable charge-charge interactions are possible than at pH 3, where unfavorable electrostatic repulsion among the positive charges causes an expansion of the denatured state ensemble. Further support is provided by studies of the ionic strength dependence of the stability of charge-reversal mutants of ribonuclease Sa. These results may have important implications for the mechanism of protein folding.

Enzyme Stability↗