Unfolding crystallins: the destabilizing role of a beta-hairpin cysteine in betaB2-crystallin by simulation and experiment.

The thermodynamic and kinetic stabilities of the eye lens family of betagamma-crystallins are important factors in the etiology of senile cataract. They control the chance of proteins unfolding, which can lead to aggregation and loss of transparency. betaB2-Crystallin orthologs are of low stability and comprise two typical betagamma-crystallin domains, although, uniquely, the N-terminal domain has a cysteine in one of the conserved folded beta-hairpins. Using high-temperature (500 K) molecular dynamics simulations with explicit solvent on the N-terminal domain of rodent betaB2-crystallin, we have identified in silico local flexibility in this folded beta-hairpin. We have shown in vitro using two-domain human betaB2-crystallin that replacement of this cysteine with a more usual aromatic residue (phenylalanine) results in a gain in conformational stability and a reduction in the rate of unfolding. We have used principal components analysis to visualize and cluster the coordinates from eight separate simulated unfolding trajectories of both the wild-type and the C50F mutant N-terminal domains. These data, representing fluctuations around the native well, show that although the mutant and wild-type appear to behave similarly over the early time period, the wild type appears to explore a different region of conformational space. It is proposed that the advantage of having this low-stability cysteine may be correlated with a subunit-exchange mechanism that allows betaB2-crystallin to interact with a range of other beta-crystallin subunits.

Protein stability and plasticity of the hydrophobic cavity in wheat ns-LTP.

Plant ns-LTPs display an original structure with four helices and a flexible C-terminus, maintained together by four disulphide bridges and delineating an elongated central hydrophobic cavity. In order to relate these structural features to the protein stability and plasticity, combined molecular mechanics and simulated annealing calculations were undertaken on a wheat ns-LTP "mutant" with Cys-Ala replacement and with the application of core inter-residue restraints up to 2 A, reducing the cross-section size of the hydrophobic cavity. Analysis of the energy-minimized structures shows that removal of the disulphide bridges results in structures with a lower total energy and a smaller cavity volume. A 1-ns MD simulation at 300K in water, underlines that, despite the absence of a well-packed hydrophobic core, the native structure is extremely stable at room temperature and the cavity is not hydrated. This confirms that the disulphide bridges are essential for the existence of the cavity, whereas its plasticity depends both on the hydrophobic chain lining the cavity and on the C-terminal flexibility. A high temperature (500K) MD simulation confirms the stability of the secondary structure elements and the flexibility of the loops and of the C-terminal segment. Two important structural transitions during this simulation are discussed and possible routes for the insertion and release of hydrophobic ligands are suggested.

Simulations of human lysozyme: probing the conformations triggering amyloidosis.

A natural mutant of human lysozyme, D67H, causes hereditary systemic nonneuropathic amyloidosis, which can be fatal. In this disease, insoluble beta-stranded fibrils (amyloids) are found in tissues stemming from the aggregation of partially folded intermediates of the mutant. In this study, we specifically compare the conformation and properties of the structures adopted from the induced unfolding, at elevated temperature, using molecular dynamics. To increase the sampling of the unfolding conformational landscape, three 5 ns trajectories are performed for each of the wild-type and mutant D67H proteins resulting in a total of 30 ns simulation. Our results show that the mutant unfolds slightly faster than the wild-type with both wild-type and mutant proteins losing most of their native secondary structure within the first 2 ns. They both develop random transient beta-strands across the whole polypeptide chain. Clustering analysis of all the conformations shows that a high population of the mutant protein conformations have a distorted beta-domain. This is consistent with experimental results suggesting that this region is pivotal in the formation of conformations prone to act as "seeds" for amyloid fiber formation.

Histidine pK(a) values for the N-lobe of human transferrin: effect of substitution of binding site Asp by Ser (D63S).

The pK(a) values have been determined for eight of the nine histidine residues and the amino terminus of the N-lobe of human apo-transferrin (hTF/2N), and for seven of the nine histidine residues and the amino terminus of the protein Asp63Ser hTF/2N containing a mutation of the Fe(3+)-ligand Asp63 to Ser63. Calculations suggested that substitution of aspartate by serine would result in decreases of the pK(a) values of most of the histidine residues in the protein. This was found to be the case experimentally, and allowed assignment of the varepsilonCH resonance of His249. For the wild-type protein, the His residue with a pK(a) of 7.40 was assigned as His249, whereas for the mutant, no observable His residue had a pK(a) value higher than 6.9. The protonated form of His249 appears to be stabilised by interactions with Asp63, and the high pK(a) value may be critical for ensuring the release of iron at endosomal pH (5.5). The mutation lowered the apparent binding constant of hTF/2N for the synergistic anion oxalate from log K 4.0 to log K 3.3. (1)H NMR spectral changes induced by Ga(3+) binding to the mutant are compared to those observed for the wild-type protein.

The X-ray crystal structure of human gamma S-crystallin C-terminal domain.

gammaS-crystallin is a major human lens protein found in the outer region of the eye lens, where the refractive index is low. Because crystallins are not renewed they acquire post-translational modifications that may perturb stability and solubility. In common with other members of the betagamma-crystallin superfamily, gammaS-crystallin comprises two similar beta-sheet domains. The crystal structure of the C-terminal domain of human gammaS-crystallin has been solved at 2.4 A resolution. The structure shows that in the in vitro expressed protein, the buried cysteines remain reduced. The backbone conformation of the "tyrosine corner" differs from that of other betagamma-crystallins because of deviation from the consensus sequence. The two C-terminal domains in the asymmetric unit are organized about a slightly distorted 2-fold axis to form a dimer with similar geometry to full-length two-domain family members. Two glutamines found in lattice contacts may be important for short range interactions in the lens. An asparagine known to be deamidated in human cataract is located in a highly ordered structural region.