Electrostatic interactions play an essential role in the binding of oleic acid with α‐lactalbumin in the HAMLET‐like complex: A study using charge‐specific chemical modifications

H uman α ‐lactalbumin m ade le thal to t umor cells (HAMLET) and its analogs are partially unfolded protein‐oleic acid (OA) complexes that exhibit selective tumoricidal activity normally absent in the native protein itself. To understand the nature of the interaction between protein and OA moieties, charge‐specific chemical modifications of lysine side chains involving citraconylation, acetylation, and guanidination were employed and the biophysical and biological properties were probed. Upon converting the original positively‐charged lysine residues to negatively‐charged citraconyl or neutral acetyl groups, the binding of OA to protein was eliminated, as were any cytotoxic activities towards osteosarcoma cells. Retention of the positive charges by converting lysine residues to homoarginine groups (guanidination); however, yielded unchanged binding of OA to protein and identical tumoricidal activity to that displayed by the wild‐type α‐lactalbumin‐oleic acid complex. With the addition of OA, the wild‐type and guanidinated α‐lactalbumin proteins underwent substantial conformational changes, such as partial unfolding, loss of tertiary structure, but retention of secondary structure. In contrast, no significant conformational changes were observed in the citraconylated and acetylated α‐lactalbumins, most likely because of the absence of OA binding. These results suggest that electrostatic interactions between the positively‐charged basic groups on α‐lactalbumin and the negatively‐charged carboxylate groups on OA molecules play an essential role in the binding of OA to α‐lactalbumin and that these interactions appear to be as important as hydrophobic interactions. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Cysteine‐free rop: A four‐helix bundle core mutant has wild‐type stability and structure but dramatically different unfolding kinetics

Cysteine residues can complicate the folding and storage of proteins due to improper formation of disulfide bonds or oxidation of residues that are natively reduced. Wild‐type Rop is a homodimeric four‐helix bundle protein and an important model for protein design in the understanding of protein stability, structure and folding kinetics. In the native state, Rop has two buried, reduced cysteine residues in its core, but these are prone to oxidation in destabilized variants, particularly upon extended storage. To circumvent this problem, we designed and characterized a Cys‐free variant of Rop, including solving the 2.3 Å X‐ray crystal structure. We show that the C38A C52V variant has similar structure, stability and in vivo activity to wild‐type Rop, but that it has dramatically faster unfolding kinetics like virtually every other mutant of Rop that has been characterized. This cysteine‐free Rop has already proven useful for studies on solution topology and on the relationship of core mutations to stability. It also suggests a general strategy for removal of pairs of Cys residues in proteins, both to make them more experimentally tractable and to improve their storage properties for therapeutic or industrial purposes.     

Analysis of control elements for position‐independent expression of human α‐lactalbumin YAC

A major problem in the production of transgenic animal bioreactors using microinjections is the low production rate of high‐expressing transgenic animals due to the position effect. We previously reported that transgenic rats carrying the 210 kb yeast artificial chromosome (YAC) including the human α‐lactalbumin gene express the transgene in a position‐independent manner. The 210 kb YAC was thought to have all the elements necessary for position‐independent expression. In this paper, we constructed fragmented YAC clones and a cosmid clone, and produced transgenic rats to analyze these elements. Transgenic rats with both the 50 kb upstream and downstream regions of the α‐lactalbumin gene had position‐independent expression. Transgenic rats with the 20 kb upstream and downstream regions, however, had position‐dependent expression. Therefore, all the elements necessary for position‐independent expression are thought to be located in the 50 kb upstream to 50 kb downstream region of the α‐lactalbumin gene. Furthermore, we replaced the human α‐lactalbumin promoter with the bovine αS1‐casein promoter in the 210 kb YAC and produced transgenic rats. Position‐dependent expression was observed. The elements required for position‐independent expression of the bovine αS1‐casein gene are different from those required for the human α‐lactalbumin gene, despite the fact that the two genes have the same tissue and developmental specificity. Mol. Reprod. Dev. 54:17–23, 1999. © 1999 Wiley‐Liss, Inc.     

Small‐angle X‐ray scattering of BAMLET at pH 12: A complex of α‐lactalbumin and oleic acid

BAMLET (Bovine Alpha‐lactalbumin Made LEthal to Tumors) is a member of the family of the HAMLET‐like complexes, a novel class of protein‐based anti‐cancer complexes that incorporate oleic acid and deliver it to cancer cells. Small angle X‐ray scattering (SAXS) was performed on the complex at pH 12, examining the high pH structure as a function of oleic acid added. The SAXS data for BAMLET species prepared with a range of oleic acid concentrations indicate extended, irregular, partially unfolded protein conformations that vary with the oleic acid concentration. Increases in oleic acid concentration correlate with increasing radius of gyration without an increase in maximum particle dimension, indicating decreasing protein density. The models for the highest oleic acid content BAMLET indicate an unusual coiled elongated structure that contrasts with apo‐α‐lactalbumin at pH 12, which is an elongated globular molecule, suggesting that oleic acid inhibits the folding or collapse of the protein component of BAMLET to the globular form. Circular dichroism of BAMLET and apo‐α‐lactalbumin was performed and the results suggest that α‐lactalbumin and BAMLET unfold in a continuum of increasing degree of unfolded states. Taken together, these results support a model in which BAMLET retains oleic acid by non‐specific association in the core of partially unfolded protein, and represent a new type of lipoprotein structure. Proteins 2014; 82:1400–1408. © 2014 Wiley Periodicals, Inc.     

Comprehensive secondary‐structure analysis of disulfide variants of lysozyme by synchrotron‐radiation vacuum‐ultraviolet circular dichroism

To elucidate the effects of specific disulfide bridges (Cys6‐Cys127, Cys30‐Cys115, Cys64‐Cys80, and Cys76‐Cys94) on the secondary structure of hen lysozyme, the vacuum‐ultraviolet circular dichroism (VUVCD) spectra of 13 species of disulfide‐deficient variants in which Cys residues were replaced with Ala or Ser residues were measured down to 170 nm at pH 2.9 and 25°C using a synchrotron‐radiation VUVCD spectrophotometer. Each variant exhibited a VUVCD spectrum characteristic of a considerable amount of residual secondary structures depending on the positions and numbers of deleted disulfide bridges. The contents of α‐helices, β‐strands, turns, and unordered structures were estimated with the SELCON3 program using the VUVCD spectra and PDB data of 31 reference proteins. The numbers of α‐helix and β‐strand segments were also estimated from the VUVCD data. In general, the secondary structures were more effectively stabilized through entropic forces as the number of disulfide bridges increased and as they were formed over larger distances in the primary structure. The structures of three‐disulfide variants were similar to that of the wild type, but other variants exhibited diminished α‐helices with a border between the ordered and disordered structures around the two‐disulfide variants. The sequences of the secondary structures were predicted for all the variants by combining VUVCD data with a neural‐network method. These results revealed the characteristic role of each disulfide bridge in the formation of secondary structures. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Denaturation induced aggregation in α‐crystallin: differential action of chaotropes

α‐Crystallin is a member of small heat shock proteins and is believed to play an exceptional role in the stability of eye lens proteins. The disruption or denaturation of the protein arrangement or solubility of the crystallin proteins can lead to vision problems including cataract. In the present study, we have examined the effect of chemical denaturants urea and guanidine hydrochloride (GdnHCl) on α‐crystallin aggregation, with special emphasis on protein conformational changes, unfolding, and amyloid fibril formation. GdnHCl (4 M) induced a 16 nm red shift in the intrinsic fluorescence of α‐crystallin, compared with 4 nm shift by 8 M urea suggesting a major change in α‐crystallin structure. Circular dichroism analysis showed marked increase in the ellipticity of α‐crystallin at 216 nm, suggesting gain in β‐sheet structure in the presence of GdnHCl (0.5–1 M) followed by unfolding at higher concentration (2–6 M). However, only minor changes in the secondary structure of α‐crystallin were observed in the presence of urea. Moreover, 8‐anilinonaphthalene‐1‐sulfonic acid fluorescence measurement in the presence of GdnHCl and urea showed changes in the hydrophobicity of α‐crystallin. Amyloid studies using thioflavin T fluorescence and congo red absorbance showed that GdnHCl induced amyloid formation in α‐crystallin, whereas urea induced aggregation in this protein. Electron microscopy studies further confirmed amyloid formation of α‐crystallin in the presence of GdnHCl, whereas only aggregate‐like structures were observed in α‐crystallin treated with urea. Our results suggest that α‐crystallin is susceptible to unfolding in the presence of chaotropic agents like urea and GdnHCl. The destabilized protein has increased likelihood to fibrillate. Copyright © 2016 John Wiley & Sons, Ltd.     

Nanotubular structures developed from whey‐based α‐lactalbumin fractions for food applications

Whey proteins have high nutritional value providing use in dietary purposes and improvement of technological properties in processed foods. Functionality of the whey‐based α‐lactalbumin (α‐La) may be increased when assembled in the form of nanotubes, promising novel potential applications subject to investigation. The purpose of this study was to extract highly pure α‐La from whey protein isolate (WPI) and whey powder (WP) and to construct protein nanotubes from them for industrial applications. For protein fractionation, WPI was directly fed to chromatography, however, WP was first subjected to membrane filtration and the retentate fraction, whey protein concentrate (WPC), was obtained and then used for chromatographic separation. α‐La and, additionally β‐Lg, were purified at the same batches with the purities in the range of 95%–99%. After enzymatic hydrolysis, WPI‐based α‐La produced chain‐like and long nanotubules with ～20 nm width while WPC‐based α‐La produced thinner, miscellaneous, and fibril‐like nanostructures by self‐assembly. Raman and FT‐IR spectroscopies revealed that α‐La fractions, obtained from both sources and the nanostructures, developed using both fractions have some structural differences due to conformation of secondary structure elements. Nanotube formation induced gelation and nanotubular gel network entrapped a colorant uniformly with a transparent appearance. Dairy‐based α‐La protein nanotubules could be served as alternative gelling agents and the carriers of natural colorants in various food processes. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 30:1301–1310, 2014     

Unfolding study of a trimeric membrane protein AcrB

The folding of a multi‐domain trimeric α‐helical membrane protein, Escherichia coli inner membrane protein AcrB, was investigated. AcrB contains both a transmembrane domain and a large periplasmic domain. Protein unfolding in sodium dodecyl sulfate (SDS) and urea was monitored using the intrinsic fluorescence and circular dichroism spectroscopy. The SDS denaturation curve displayed a sigmoidal profile, which could be fitted with a two‐state unfolding model. To investigate the unfolding of separate domains, a triple mutant was created, in which all three Trp residues in the transmembrane domain were replaced with Phe. The SDS unfolding profile of the mutant was comparable to that of the wild type AcrB, suggesting that the observed signal change was largely originated from the unfolding of the soluble domain. Strengthening of trimer association through the introduction of an inter‐subunit disulfide bond had little effect on the unfolding profile, suggesting that trimer dissociation was not the rate‐limiting step in unfolding monitored by fluorescence emission. Under our experimental condition, AcrB unfolding was not reversible. Furthermore, we experimented with the refolding of a monomeric mutant, AcrB_Δloop, from the SDS unfolded state. The CD spectrum of the refolded AcrB_Δloop superimposed well onto the spectra of the original folded protein, while the fluorescence spectrum was not fully recovered. In summary, our results suggested that the unfolding of the trimeric AcrB started with a local structural rearrangement. While the refolding of secondary structure in individual monomers could be achieved, the re‐association of the trimer might be the limiting factor to obtain folded wild‐type AcrB.     

Oxidative folding and preparation of α‐conotoxins for use in high‐throughput structure–activity relationship studies

α‐Conotoxins are peptide neurotoxins that selectively inhibit various subtypes of nicotinic acetylcholine receptors. They are important research tools for studying numerous pharmacological disorders, with profound potential for developing drug leads for treating pain, tobacco addiction, and other conditions. They are characterized by the presence of two disulfide bonds connected in a globular arrangement, which stabilizes a bioactive helical conformation. Despite extensive structure–activity relationship studies that have produced α‐conotoxin analogs with increased potency and selectivity towards specific nicotinic acetylcholine receptor subtypes, the efficient production of diversity‐oriented α‐conotoxin combinatorial libraries has been limited by inefficient folding and purification procedures. We have investigated the optimized conditions for the reliable folding of α‐conotoxins using simplified oxidation procedures for use in the accelerated production of synthetic combinatorial libraries of α‐conotoxins. To this end, the effect of co‐solvent, redox reagents, pH, and temperature on the proportion of disulfide bond isomers was determined for α‐conotoxins exhibiting commonly known Cys loop spacing frameworks. In addition, we have developed high‐throughput ‘semi‐purification’ methods for the quick and efficient parallel preparation of α‐conotoxin libraries for use in accelerated structure–activity relationship studies. Our simplified procedures represent an effective strategy for the preparation of large arrays of correctly folded α‐conotoxin analogs and permit the rapid identification of active hits directly from high‐throughput pharmacological screening assays. Copyright © 2012 European Peptide Society and John Wiley & Sons, Ltd.     

α‐chymotrypsin in water–acetone and water–dimethyl sulfoxide mixtures: Effect of preferential solvation and hydration

We investigated water/organic solvent sorption and residual enzyme activity to simultaneously monitor preferential solvation/hydration of protein macromolecules in the entire range of water content at 25°C. We applied this approach to estimate protein destabilization/stabilization due to the preferential interactions of bovine pancreatic α‐chymotrypsin with water‐acetone (moderate‐strength H‐bond acceptor) and water‐DMSO (strong H‐bond acceptor) mixtures. There are three concentration regimes for the dried α‐chymotrypsin. α‐Chymotrypsin is preferentially hydrated at high water content. The residual enzyme activity values are close to 100%. At intermediate water content, the dehydrated α‐chymotrypsin has a higher affinity for acetone/DMSO than for water. Residual enzyme activity is minimal in this concentration range. The acetone/DMSO molecules are preferentially excluded from the protein surface at the lowest water content, resulting in preferential hydration. The residual catalytic activity in the water‐poor acetone is ～80%, compared with that observed after incubation in pure water. This effect is very small for the water‐poor DMSO. Two different schemes are operative for the hydrated enzyme. At high and intermediate water content, α‐chymotrypsin exhibits preferential hydration. However, at intermediate water content, in contrast to the dried enzyme, the initially hydrated α‐chymotrypsin possesses increased preferential hydration parameters. At low water content, no residual enzyme activity was observed. Preferential binding of DMSO/acetone to α‐chymotrypsin was detected. Our data clearly demonstrate that the hydrogen bond accepting ability of organic solvents and the protein hydration level constitute key factors in determining the stability of protein–water–organic solvent systems.     

Analysis of peptides and proteins in their binding to GroEL

The GroEL–GroES is an essential molecular chaperon system that assists protein folding in cell. Binding of various substrate proteins to GroEL is one of the key aspects in GroEL‐assisted protein folding. Small peptides may mimic segments of the substrate proteins in contact with GroEL and allow detailed structural analysis of the interactions. A model peptide SBP has been shown to bind to a region in GroEL that is important for binding of substrate proteins. Here, we investigated whether the observed GroEL–SBP interaction represented those of GroEL–substrate proteins, and whether SBP was able to mimic various aspects of substrate proteins in GroE‐assisted protein folding cycle. We found that SBP competed with substrate proteins, including α‐lactalbumin, rhodanese, and malate dehydrogenase, in binding to GroEL. SBP stimulated GroEL ATP hydrolysis rate in a manner similar to that of α‐lactalbumin. SBP did not prevent GroES from binding to GroEL, and GroES association reduced the ATPase rates of GroEL/SBP and GroEL/α‐lactalbumin to a comparable extent. Binding of both SBP and α‐lactalbumin to apo GroEL was dominated by hydrophobic interaction. Interestingly, association of α‐lactalbumin to GroEL/GroES was thermodynamically distinct from that to GroEL with reduced affinity and decreased contribution from hydrophobic interaction. However, SBP did not display such differential binding behaviors to apo GroEL and GroEL/GroES, likely due to the lack of a contiguous polypeptide chain that links all of the bound peptide fragments. Nevertheless, studies using peptides provide valuable information on the nature of GroEL–substrate protein interaction, which is central to understand the mechanism of GroEL‐assisted protein folding. Copyright © 2010 European Peptide Society and John Wiley & Sons, Ltd.     

In vitro folding of methionine‐arginine human lyspro‐proinsulin S‐sulfonate—Disulfide formation pathways and factors controlling yield

We investigated the in vitro folding of an oxidized proinsulin (methionine‐arginine human lyspro‐proinsulin S‐sulfonate), using cysteine as a reducing agent at 5°C and high pH (10.5–11). Folding intermediates were detected and characterized by means of matrix‐assisted laser desorption ionization mass spectrometry (MALDI‐MS), reversed‐phase chromatography (RPC), size‐exclusion chromatography, and gel electrophoresis. The folding kinetics and yield depended on the protein and cysteine concentrations. RPC coupled with MALDI‐MS analyses indicated a sequential formation of intermediates with one, two, and three disulfide bonds. The MALDI‐MS analysis of Glu‐C digested, purified intermediates indicated that an intra‐A‐chain disulfide bond formed first among A6, A7, and A11. Various non‐native intra‐A (A20 with A6, A7, or A11), intra‐B (between B7 and B19), and inter‐A‐B disulfide bonds were observed in the intermediates with two disulfide bonds. The intermediates with three disulfide bonds had mainly the non‐native intra‐A and intra‐B bonds. At a cysteine‐to‐proinsulin‐SH ratio of 3.5, all intermediates with the non‐native disulfide bonds were converted to properly folded proinsulin via disulfide bond reshuffling, which was the slowest step. Aggregation via the formation of intermolecular disulfide bonds of early intermediates was the major cause of yield loss. At a higher cysteine‐to‐proinsulin‐SH ratio, some intermediates and folded MR‐KPB‐hPI were reduced to proteins with thiolate anions, which caused unfolding and even more yield loss than what resulted from aggregation of the early intermediates. Reducing protein concentration, while keeping an optimal cysteine‐to‐protein ratio, can improve folding yield significantly. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Thermodynamic analysis of ANS binding to partially unfolded α‐lactalbumin: correlation of endothermic to exothermic changeover with formation of authentic molten globules

A fluorescent reporter, 8‐anilino‐1‐naphthalene sulfonic acid (ANS), can serve as a reference molecule for conformational transition of a protein because its aromatic carbons have strong affinity with hydrophobic cores of partially unfolded molten globules. Using a typical calcium‐binding protein, bovine α‐lactalbumin (BLA), as a model protein, we compared the ANS binding thermodynamics to the decalcified (10 mM EDTA treated) apo‐BLA at two representative temperatures: 20 and 40 °C. This is because the authentic molten globule is known to form more heavily at an elevated temperature such as 40 °C. Isothermal titration calorimetry experiments revealed that the BLA–ANS interactions at both temperatures were entropy‐driven, and the dissociation constants were similar on the order of 10^?4 M, but there was a dramatic changeover in the binding thermodynamics from endothermic at 20 °C to exothermic at 40 °C. We believe that the higher subpopulation of authentic molten globules at 40 °C than 20 °C would be responsible for the results, which also indicate that weak binding is sufficient to alter the ANS binding mechanisms. We expect that the thermodynamic properties obtained from this study would serve as a useful reference for investigating the binding of other hydrophobic ligands such as oleic acid to apo‐BLA, because oleic acid is known to have tumor‐selective cytotoxicity when complexed with partially unfolded α‐lactalbumin. Copyright © 2016 John Wiley & Sons, Ltd.     

Dramatic stabilization of the native state of human carbonic anhydrase II by an engineered disulfide bond

To find a disulfide pair that could stabilize the enzyme human carbonic anhydrase II (HCA II), we grafted the disulfide bridge from the related and unusually stable carbonic anhydrase form from Neisseria gonorrhoeae (NGCA) into the human enzyme. Thus, the two Cys residues at positions 23 and 203 were engineered into a pseudo-wild-type form of HCA II (C206S), giving the mutant C206S/A23C/L203C. The disulfide bond was not formed spontaneously. The native state of the reduced form of the mutant was markedly destabilized (2.9 kcal/mol) compared to that of HCA II. Formation of a disulfide bridge was achieved by treatment by oxidized glutathione. This led to a significant stabilization of the native conformation. Compared to HCA II the unfolding midpoint for the variant was increased from 0.9 to 1.7 M guanidine HCl, corresponding to a stabilization of 3.7 kcal/mol. This makes the human enzyme almost as stable as the model protein NGCA, for which the unfolding of the native state has a midpoint at 2.1 M guanidine HCl. The stabilized protein underwent, contrary to all other investigated variants of HCA II, an apparent two-state unfolding transition, as judged from intrinsic Trp fluorescence measurements. A molten-globule intermediate is nevertheless formed but is suppressed because of the high denaturant pressure it faces upon rupture of the native state.     

Identification of a α‐helical molten globule intermediate and structural characterization of β‐cardiotoxin,an all β‐sheet protein isolated from the venom of Ophiophagus hannah (king cobra)

β‐Cardiotoxin is a novel member of the snake venom three‐finger toxin (3FTX) family. This is the first exogenous protein to antagonize β‐adrenergic receptors and thereby causing reduction in heart rates (bradycardia) when administered into animals, unlike the conventional cardiotoxins as reported earlier. 3FTXs are stable all β‐sheet peptides with 60–80 amino acid residues. Here, we describe the three‐dimensional crystal structure of β‐cardiotoxin together with the identification of a molten globule intermediate in the unfolding pathway of this protein. In spite of the overall structural similarity of this protein with conventional cardiotoxins, there are notable differences observed at the loop region and in the charge distribution on the surface, which are known to be critical for cytolytic activity of cardiotoxins. The molten globule intermediate state present in the thermal unfolding pathway of β‐cardiotoxin was however not observed during the chemical denaturation of the protein. Interestingly, circular dichroism (CD) and NMR studies revealed the presence of α‐helical secondary structure in the molten globule intermediate. These results point to substantial conformational plasticity of β‐cardiotoxin, which might aid the protein in responding to the sometimes conflicting demands of structure, stability, and function during its biological lifetime.     

Spectroscopic studies of the unfolding of a multimeric protein α‐crystallin

α‐Crystallin is a multimeric eye lens protein having molecular chaperone‐like function which is crucial for lens transparency. The stability and unfolding‐refolding properties of α‐crystallin plays important roles for its function. We undertook a multi probe based fluorescence spectroscopic approach to explore the changes in the various levels of organization of this protein at different urea concentration. Steady state fluorescence studies reveal that at 0.6M urea a compact structural intermediate is formed which has a native‐like secondary structure with enhanced surface exposure of hydrophobic groups. At 2.8M urea the tertiary interactions are largely collapsed with partial collapse of secondary and quaternary structure. The surface solvation probed by picosecond time resolved fluorescence of acrylodan labeled α‐crystallin revealed dry native‐like core of α‐crystallin at 0.6M urea compared to enhanced water penetration at 2.8M urea and extensive solvation at 6M urea. Activation energy for the subunit exchange decreased by 22 kJ mol^?1 on changing urea concentration from 0 to 0.6M compared with over 75 kJ mol^?1 on changing urea concentration from 0 to 2.8M. Light scattering and analytical ultracentrifugation techniques were used to determine size and oligomerization of the unfolding intermediates. The data indicated swelling but no oligomer breakdown at 0.6M urea. At 2.8M urea the oligomeric size is considerably reduced and a monomer is produced at 6M urea. The data clearly reveals that structural breakdown of α‐crystallin does not follow hierarchical sequence as tertiary structure dissolution takes place before complete oligomeric dissociation. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 549–560, 2014.     

Association of bovine β‐casein protein variant I with milk production and milk protein composition

The aim of this study was to detect new polymorphisms in the bovine β‐casein (β‐CN) gene and to evaluate association of (new) β‐CN protein variants with milk production traits and milk protein composition. Screening of the β‐CN gene in genomic DNA from 72 Holstein Friesian (HF) bulls resulted in detection of 19 polymorphisms and revealed the presence of β‐CN protein variant I in the Dutch HF population. Studies of association of β‐CN protein variants with milk composition usually do not discriminate protein variant I from variant A2. Association of β‐CN protein variants with milk composition was studied in 1857 first‐lactation HF cows and showed that associations of protein variants A2 and I were quite different for several traits. β‐CN protein variant I was significantly associated with protein percentage and protein yield, and with α_s1‐casein (α_s1‐CN), α_s2‐casein (α_s2‐CN), κ‐casein (κ‐CN), α‐lactalbumin (α‐LA), β‐lactoglobulin (β‐LG), casein index and casein yield. Inferring β‐κ‐CN haplotypes showed that β‐CN protein variant I occurred only with κ‐CN variant B. Consequently, associations of β‐κ‐CN haplotype IB with protein percentage, κ‐CN, α‐LA, β‐LG and casein index are likely resulting from associations of κ‐CN protein variant B, while associations of β‐κ‐CN haplotype IB with α_s1‐CN and α_s2‐CN seem to be resulting from associations of β‐CN variant I.     

The effect of additional disulfide bonds on the stability and folding of ribonuclease A

The significant contribution of disulfide bonds to the conformational stability of proteins is generally considered to result from an entropic destabilization of the unfolded state causing a faster escape of the molecules to the native state. However, the introduction of extra disulfide bonds into proteins as a general approach to protein stabilization yields rather inconsistent results. By modeling studies, we selected positions to introduce additional disulfide bonds into ribonuclease A at regions that had proven to be crucial for the initiation of the folding or unfolding process, respectively. However, only two out of the six variants proved to be more stable than unmodified ribonuclease A. The comparison of the thermodynamic and kinetic data disclosed a more pronounced effect on the unfolding reaction for all variants regardless of the position of the extra disulfide bond. Native-state proteolysis indicated a perturbation of the native state of the destabilized variants that obviously counterbalances the stability gain by the extra disulfide bond.     

2,2′‐Dipyridyl diselenide: A chemoselective tool for cysteine deprotection and disulfide bond formation

There are many examples of bioactive, disulfide‐rich peptides and proteins whose biological activity relies on proper disulfide connectivity. Regioselective disulfide bond formation is a strategy for the synthesis of these bioactive peptides, but many of these methods suffer from a lack of orthogonality between pairs of protected cysteine (Cys) residues, efficiency, and high yields. Here, we show the utilization of 2,2′‐dipyridyl diselenide (PySeSePy) as a chemical tool for the removal of Cys‐protecting groups and regioselective formation of disulfide bonds in peptides. We found that peptides containing either Cys(Mob) or Cys(Acm) groups treated with PySeSePy in trifluoroacetic acid (TFA) (with or without triisopropylsilane (TIS) were converted to Cys‐S–SePy adducts at 37 °C and various incubation times. This novel Cys‐S–SePy adduct is able to be chemoselectively reduced by five‐fold excess ascorbate at pH 4.5, a condition that should spare already installed peptide disulfide bonds from reduction. This chemoselective reduction by ascorbate will undoubtedly find utility in numerous biotechnological applications. We applied our new chemistry to the iodine‐free synthesis of the human intestinal hormone guanylin, which contains two disulfide bonds. While we originally envisioned using ascorbate to chemoselectively reduce one of the formed Cys‐S–SePy adducts to catalyze disulfide bond formation, we found that when pairs of Cys(Acm) residues were treated with PySeSePy in TFA, the second disulfide bond formed spontaneously. Spontaneous formation of the second disulfide is most likely driven by the formation of the thermodynamically favored diselenide (PySeSePy) from the two Cys‐S–SePy adducts. Thus, we have developed a one‐pot method for concomitant deprotection and disulfide bond formation of Cys(Acm) pairs in the presence of an existing disulfide bond.     

Temperature‐induced unfolding behavior of proteins studied by tensorial elastic network model

Motivated by single molecule experiments and recent molecular dynamics (MD) studies, we propose a simple and computationally efficient method based on a tensorial elastic network model to investigate the unfolding pathways of proteins under temperature variation. The tensorial elastic network model, which relies on the native state topology of a protein, combines the anisotropic network model, the bond bending elasticity, and the backbone twist elasticity to successfully predicts both the isotropic and anisotropic fluctuations in a manner similar to the Gaussian network model and anisotropic network model. The unfolding process is modeled by breaking the native contacts between residues one by one, and by assuming a threshold value for strain fluctuation. Using this method, we simulated the unfolding processes of four well‐characterized proteins: chymotrypsin inhibitor, barnase, ubiquitein, and adenalyate kinase. We found that this step‐wise process leads to two or more cooperative, first‐order‐like transitions between partial denaturation states. The sequence of unfolding events obtained using this method is consistent with experimental and MD studies. The results also imply that the native topology of proteins “encrypts” information regarding their unfolding process. Proteins 2016; 84:1767–1775. © 2016 Wiley Periodicals, Inc.