Elucidation of the disulfide-folding pathway of hirudin by a topology-based approach

Diselenide crosslinks for enhanced and simplified oxidative protein folding

The in vitro oxidative folding of proteins has been studied for over sixty years, providing critical insight into protein folding mechanisms. Hirudin, the most potent natural inhibitor of thrombin, is a 65-residue protein with three disulfide bonds, and is viewed as a folding model for a wide range of disulfide-rich proteins. Hirudin's folding pathway is notorious for its highly heterogeneous intermediates and scrambled isomers, limiting its folding rate and yield in vitro. Aiming to overcome these limitations, we undertake systematic investigation of diselenide bridges at native and non-native positions and investigate their effect on hirudin's folding, structure and activity. Our studies demonstrate that, regardless of the specific positions of these substitutions, the diselenide crosslinks enhanced the folding rate and yield of the corresponding hirudin analogues, while reducing the complexity and heterogeneity of the process. Moreover, crystal structure analysis confirms that the diselenide substitutions maintained the overall three-dimensional structure of the protein and left its function virtually unchanged. The choice of hirudin as a study model has implications beyond its specific folding mechanism, demonstrating the high potential of diselenide substitutions in the design, preparation and characterization of disulfide-rich proteins.

Performance assessment of different constraining potentials in computational structure prediction for disulfide-bridged proteins

The presence of disulfide bonds in proteins has very important implications on the three-dimensional structure and folding of proteins. An adequate treatment of disulfide bonds in de-novo protein simulations is therefore very important. Here we present a computational study of a set of small disulfide-bridged proteins using an all-atom stochastic search approach and including various constraining potentials to describe the disulfide bonds. The proposed potentials can easily be implemented in any code based on all-atom force fields and employed in simulations to achieve an improved prediction of protein structure. Exploring different potential parameters and comparing the structures to those from unconstrained simulations and to experimental structures by means of a scoring function we demonstrate that the inclusion of constraining potentials improves the quality of final structures significantly. For some proteins (1KVG and 1PG1) the native conformation is visited only in simulations in presence of constraints. Overall, we found that the Morse potential has optimal performance, in particular for the β-sheet proteins.

Knotting of random ring polymers in confined spaces

Stochastic simulations are used to characterize the knotting distributions of random ring polymers confined in spheres of various radii. The approach is based on the use of multiple Markov chains and reweighting techniques, combined with effective strategies for simplifying the geometrical complexity of ring conformations without altering their knot type. By these means we extend previous studies and characterize in detail how the probability to form a given prime or composite knot behaves in terms of the number of ring segments N and confining radius R. For 50 < or =N < or =450 we show that the probability of forming a composite knot rises significantly with the confinement, while the occurrence probability of prime knots are, in general, nonmonotonic functions of 1R. The dependence of other geometrical indicators, such as writhe and chirality, in terms of R and N is also characterized. It is found that the writhe distribution broadens as the confining sphere narrows.