Towards a molecular understanding of phase separation in the lens: a comparison of the X-ray structures of two high Tc gamma-crystallins, gammaE and gammaF, with two low Tc gamma-crystallins, gammaB and gammaD

Fundamental Aspects of Phase-Separated Biomolecular Condensates

Biomolecular condensates, formed through phase separation, are upending our understanding in much of molecular, cell, and developmental biology. There is an urgent need to elucidate the physicochemical foundations of the behaviors and properties of biomolecular condensates. Here we aim to fill this need by writing a comprehensive, critical, and accessible review on the fundamental aspects of phase-separated biomolecular condensates. We introduce the relevant theoretical background, present the theoretical basis for the computation and experimental measurement of condensate properties, and give mechanistic interpretations of condensate behaviors and properties in terms of interactions at the molecular and residue levels.

Interaction of βL- and γ-Crystallin with Phospholipid Membrane Using Atomic Force Microscopy

Highly concentrated lens proteins, mostly β- and γ-crystallin, are responsible for maintaining the structure and refractivity of the eye lens. However, with aging and cataract formation, β- and γ-crystallin are associated with the lens membrane or other lens proteins forming high-molecular-weight proteins, which further associate with the lens membrane, leading to light scattering and cataract development. The mechanism by which β- and γ-crystallin are associated with the lens membrane is unknown. This work aims to study the interaction of β- and γ-crystallin with the phospholipid membrane with and without cholesterol (Chol) with the overall goal of understanding the role of phospholipid and Chol in β- and γ-crystallin association with the membrane. Small unilamellar vesicles made of Chol/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (Chol/POPC) membranes with varying Chol content were prepared using the rapid solvent exchange method followed by probe tip sonication and then dispensed on freshly cleaved mica disk to prepare a supported lipid membrane. The βL- and γ-crystallin from the cortex of the bovine lens was used to investigate the time-dependent association of βL- and γ-crystallin with the membrane by obtaining the topographical images using atomic force microscopy. Our study showed that βL-crystallin formed semi-transmembrane defects, whereas γ-crystallin formed transmembrane defects on the phospholipid membrane. The size of semi-transmembrane defects increases significantly with incubation time when βL-crystallin interacts with the membrane. In contrast, no significant increase in transmembrane defect size was observed in the case of γ-crystallin. Our result shows that Chol inhibits the formation of membrane defects when βL- and γ-crystallin interact with the Chol/POPC membrane, where the degree of inhibition depends upon the amount of Chol content in the membrane. At a Chol/POPC mixing ratio of 0.3, membrane defects were observed when both βL- and γ-crystallin interacted with the membrane. However, at a Chol/POPC mixing ratio of 1, no association of γ-crystallin with the membrane was observed, which resulted in a defect-free membrane, and the severity of the membrane defect was decreased when βL-crystallin interacted with the membrane. The semi-transmembrane or transmembrane defects formed by the interaction of βL- and γ-crystallin on phospholipid membrane might be responsible for light scattering and cataract formation. However, Chol suppressed the formation of such defects in the membrane, likely maintaining lens membrane homeostasis and protecting against cataract formation.

Atomistic modeling of liquid-liquid phase equilibrium explains dependence of critical temperature on γ-crystallin sequence

Liquid-liquid phase separation of protein solutions has regained heightened attention for its biological importance and pathogenic relevance. Coarse-grained models are limited when explaining residue-level effects on phase equilibrium. Here we report phase diagrams for γ-crystallins using atomistic modeling. The calculations were made possible by combining our FMAP method for computing chemical potentials and Brownian dynamics simulations for configurational sampling of dense protein solutions, yielding the binodal and critic temperature (Tc). We obtain a higher Tc for a known high-Tc γ-crystallin, γF, than for a low-Tc paralog, γB. The difference in Tc is corroborated by a gap in second virial coefficient. Decomposition of inter-protein interactions reveals one amino-acid substitution between γB and γF, from Ser to Trp at position 130, as the major contributor to the difference in Tc. This type of analysis enables us to link phase equilibrium to amino-acid sequence and to design mutations for altering phase equilibrium.

Atomistic Modeling of Liquid-Liquid Phase Equilibrium Explains Dependence of Critical Temperature on γ-Crystallin Sequence

Liquid-liquid phase separation of protein solutions has regained heightened attention for its biological importance and pathogenic relevance. Coarse-grained models are limited when explaining residue-level effects on phase equilibrium. Here we report phase diagrams for γ-crystallins using atomistic modeling. The calculations were made possible by combining our FMAP method for computing chemical potentials and Brownian dynamics simulations for configurational sampling of dense protein solutions, yielding the binodal and critic temperature ( T c ). We obtain a higher T c for a known high- T c γ-crystallin, γF, than for a low- T c paralog, γB. The difference in T c is corroborated by a gap in second virial coefficient. Decomposition of inter-protein interactions reveals one amino-acid substitution between γB and γF, from Ser to Trp at position 130, as the major contributor to the difference in T c . This type of analysis enables us to link phase equilibrium to amino-acid sequence and to design mutations for altering phase equilibrium.

Controlling Liquid-Liquid Phase Separation of Cold-Adapted Crystallin Proteins from the Antarctic Toothfish

Liquid-liquid phase separation (LLPS) of proteins is important to a variety of biological processes both functional and deleterious, including the formation of membraneless organelles, molecular condensations that sequester or release molecules in response to stimuli, and the early stages of disease-related protein aggregation. In the protein-rich, crowded environment of the eye lens, LLPS manifests as cold cataract. We characterize the LLPS behavior of six structural γ-crystallins from the eye lens of the Antarctic toothfish Dissostichus mawsoni, whose intact lenses resist cold cataract in subzero waters. Phase separation of these proteins is not strongly correlated with thermal stability, aggregation propensity, or cross-species chaperone protection from heat denaturation. Instead, LLPS is driven by protein-protein interactions involving charged residues. The critical temperature of the phase transition can be tuned over a wide temperature range by selective substitution of surface residues, suggesting general principles for controlling this phenomenon, even in compactly folded proteins.

Crystal Structure of Chicken γS-Crystallin Reveals Lattice Contacts with Implications for Function in the Lens and the Evolution of the βγ-Crystallins

Previous attempts to crystallize mammalian γS-crystallin were unsuccessful. Native L16 chicken γS crystallized avidly while the Q16 mutant did not. The x-ray structure for chicken γS at 2.3Å resolution shows the canonical structure of the superfamily plus a well-ordered N-arm aligned with a β-sheet of a neighboring N-domain. L16 is also in a lattice contact, partially shielded from solvent. Unexpectedly, the major lattice contact matches a conserved interface (QR) in the multimeric β-crystallins. QR shows little conservation of residue contacts, except for one between symmetry-related tyrosines, but molecular dipoles for the proteins with QR show striking similarities while other γ-crystallins differ. In γS, QR has few hydrophobic contacts and features a thin layer of tightly bound water. The free energy of QR is slightly repulsive and AUC confirms no dimerization in solution. The lattice contacts suggest how γcrystallins allow close packing without aggregation in the crowded environment of the lens.

Direct Determination of Site-Specific Noncovalent Interaction Strengths of Proteins from NMR-Derived Fast Side Chain Motional Parameters

A novel approach to accurately determine residue-specific noncovalent interaction strengths (ξ) of proteins from NMR-measured fast side chain motional parameters (O_axis²) is presented. By probing the environmental sensitivity of side chain conformational energy surfaces of individual residues of a diverse set of proteins, the microscopic connections between ξ, O_axis², conformational entropy (S_conf), conformational barriers, and rotamer stabilities established here are found to be universal among proteins. The results reveal that side chain flexibility and conformational entropy of each residue decrease with increasing ξ and that for each residue type there exists a critical range of ξ, determined primarily by the mean side chain conformational barriers, within which flexibility of any residue can be reversibly tuned from highly flexible (with O_axis² ∼ 0) to highly restricted (with O_axis² ∼ 1) by increasing ξ by ∼3 kcal/mol. Beyond this critical range of ξ, both side chain flexibility and conformational entropy are insensitive to ξ. The interrelationships between conformational dynamics, conformational entropy, and noncovalent interactions of protein side chains established here open up new avenues to probe perturbation-induced (for example, ligand-binding, temperature, pressure) changes in fast side chain dynamics and thermodynamics of proteins by comparing their conformational energy surfaces in the native and perturbed states.

Fast Method for Computing Chemical Potentials and Liquid-Liquid Phase Equilibria of Macromolecular Solutions

Chemical potential is a fundamental property for determining thermodynamic equilibria involving exchange of molecules, such as between two phases of molecular systems. Previously, we developed the fast Fourier transform (FFT)-based method for Modeling Atomistic Protein-crowder interactions (FMAP) to calculate excess chemical potentials according to the Widom insertion. Intermolecular interaction energies were expressed as correlation functions and evaluated via FFT. Here, we extend this method to calculate liquid-liquid phase equilibria of macromolecular solutions. Chemical potentials are calculated by FMAP over a wide range of molecular densities, and the condition for coexistence of low- and high-density phases is determined by the Maxwell equal-area rule. When benchmarked on Lennard-Jones fluids, our method produces an accurate phase diagram at 18% of the computational cost of the current best method. Importantly, the gain in computational speed increases dramatically as the molecules become more complex, leading to many orders of magnitude in speed up for atomistically represented proteins. We demonstrate the power of FMAP by reporting the first results for the liquid-liquid coexistence curve of γII-crystallin represented at the all-atom level. Our method may thus open the door to accurate determination of phase equilibria for macromolecular mixtures such as protein-protein mixtures and protein-RNA mixtures, that are known to undergo liquid-liquid phase separation, both in vitro and in vivo.

The joys and perils of flexible fitting

While performing their functions, biological macromolecules often form large, dynamically changing macromolecular assemblies. Only a relatively small number of such assemblies have been accessible to the atomic-resolution techniques X-ray crystallography and NMR. Electron microscopy in conjunction with image reconstruction has become the preferred alternative for revealing the structures of such macromolecular complexes. However, for most assemblies the achievable resolution is too low to allow accurate atomic modeling directly from the data. Yet, useful models often can be obtained by fitting atomic models of individual components into a low-resolution reconstruction of the entire assembly. Several algorithms for achieving optimal fits in this context were developed recently, many allowing considerable degrees of flexibility to account for binding-induced conformational changes of the assembly components. This chapter describes the advantages and potential pitfalls of these methods and puts them into perspective with alternative approaches such as iterative modular fitting of rigid-body domains.

Solution properties of γ-crystallins: hydration of fish and mammal γ-crystallins

Lens γ crystallins are found at the highest protein concentration of any tissue, ranging from 300 mg/mL in some mammals to over 1000 mg/mL in fish. Such high concentrations are necessary for the refraction of light, but impose extreme requirements for protein stability and solubility. γ-crystallins, small stable monomeric proteins, are particularly associated with the lowest hydration regions of the lens. Here, we examine the solvation of selected γ-crystallins from mammals (human γD and mouse γS) and fish (zebrafish γM2b and γM7). The thermodynamic water binding coefficient B₁ could be probed by sucrose expulsion, and the hydrodynamic hydration shell of tightly bound water was probed by translational diffusion and structure-based hydrodynamic boundary element modeling. While the amount of tightly bound water of human γD was consistent with that of average proteins, the water binding of mouse γS was found to be relatively low. γM2b and γM7 crystallins were found to exhibit extremely low degrees hydration, consistent with their role in the fish lens. γM crystallins have a very high methionine content, in some species up to 15%. Structure-based modeling of hydration in γM7 crystallin suggests low hydration is associated with the large number of surface methionine residues, likely in adaptation to the extremely high concentration and low hydration environment in fish lenses. Overall, the degree of hydration appears to balance stability and tissue density requirements required to produce and maintain the optical properties of the lens in different vertebrate species.

Evolution of crystallins for a role in the vertebrate eye lens

The camera eye lens of vertebrates is a classic example of the re-engineering of existing protein components to fashion a new device. The bulk of the lens is formed from proteins belonging to two superfamilies, the α-crystallins and the βγ-crystallins. Tracing their ancestry may throw light on the origin of the optics of the lens. The α-crystallins belong to the ubiquitous small heat shock proteins family that plays a protective role in cellular homeostasis. They form enormous polydisperse oligomers that challenge modern biophysical methods to uncover the molecular basis of their assembly structure and chaperone-like protein binding function. It is argued that a molecular phenotype of a dynamic assembly suits a chaperone function as well as a structural role in the eye lens where the constraint of preventing protein condensation is paramount. The main cellular partners of α-crystallins, the β- and γ-crystallins, have largely been lost from the animal kingdom but the superfamily is hugely expanded in the vertebrate eye lens. Their structures show how a simple Greek key motif can evolve rapidly to form a complex array of monomers and oligomers. Apart from remaining transparent, a major role of the partnership of α-crystallins with β- and γ-crystallins in the lens is to form a refractive index gradient. Here, we show some of the structural and genetic features of these two protein superfamilies that enable the rapid creation of different assembly states, to match the rapidly changing optical needs among the various vertebrates.

Characterization of a transient unfolding intermediate in a core mutant of γS-crystallin

In many age-related and neurological diseases, formerly native proteins aggregate via formation of a partially unfolded intermediate. γS-Crystallin is a highly stable structural protein of the eye lens. In the mouse Opj cataract, a non-conservative F9S mutation in the N-terminal domain core of γS allows the adoption of a native fold but renders the protein susceptible to temperature- and concentration-dependent aggregation, including fibril formation. Hydrogen/deuterium exchange and denaturant unfolding studies of this mutant protein (Opj) have suggested the existence of a partially unfolded intermediate in its aggregation pathway. Here, we used NMR and fluorescence spectroscopy to obtain evidence for this intermediate. In 3.5 M urea, Opj forms a stable and partially unfolded entity that is characterized by an unstructured N-terminal domain and a largely intact C-terminal domain. Under physiologically relevant conditions, Carr-Purcell-Meiboom-Gill T(2)-relaxation dispersion experiments showed that the N-terminal domain residues were in conformational exchange with a loosely structured intermediate with a population of 1-2%, which increased with temperature. This provides direct evidence for a model in which proteins of native fold can explore an intermediate state with an increased propensity for formation of aggregates, such as fibrils. For the crystallins, this shows how inherited sequence variants or environmentally induced modifications can destabilize a well-folded protein, allowing the formation of intermediates able to act as nucleation sites for aggregation and the accumulation of light-scattering centers in the cataractous lens.

A single destabilizing mutation (F9S) promotes concerted unfolding of an entire globular domain in gammaS-crystallin

Conformational change and aggregation of native proteins are associated with many serious age-related and neurological diseases. gammaS-Crystallin is a highly stable, abundant structural component of vertebrate eye lens. A single F9S mutation in the N-terminal domain of mouse gammaS-crystallin causes the severe Opj cataract, with disruption of cellular organization and appearance of fibrillar structures in the lens. Although the mutant protein has a near-native fold at room temperature, significant increases in hydrogen/deuterium exchange rates were observed by NMR for all the well-protected beta-sheet core residues throughout the entire N-terminal domain of the mutant protein, resulting in up to a 3.5-kcal/mol reduction in the free energy of the folding/unfolding equilibrium. No difference was detected for the C-terminal domain. At a higher temperature, this effect further increases to allow for a much more uniform exchange rate among the N-terminal core residues and those of the least well-structured surface loops. This suggests a concerted unfolding intermediate of the N-terminal domain, while the C-terminal domain stays intact. Increasing concentrations of guanidinium chloride produced two transitions for the Opj mutant, with an unfolding intermediate at approximately 1 M guanidinium chloride. The consequence of this partial unfolding, whether by elevated temperature or by denaturant, is the formation of thioflavin T staining aggregates, which demonstrated fibril-like morphology by atomic force microscopy. Seeding with the already unfolded protein enhanced the formation of fibrils. The Opj mutant protein provides a model for stress-related unfolding of an essentially normally folded protein and production of aggregates with some of the characteristics of amyloid fibrils.

The subunit interfaces of weakly associated homodimeric proteins

We analyzed subunit interfaces in 315 homodimers with an X-ray structure in the Protein Data Bank, validated by checking the literature for data that indicate that the proteins are dimeric in solution and that, in the case of the "weak" dimers, the homodimer is in equilibrium with the monomer. The interfaces of the 42 weak dimers, which are smaller by a factor of 2.4 on average than in the remainder of the set, are comparable in size with antibody-antigen or protease-inhibitor interfaces. Nevertheless, they are more hydrophobic than in the average transient protein-protein complex and similar in amino acid composition to the other homodimer interfaces. The mean numbers of interface hydrogen bonds and hydration water molecules per unit area are also similar in homodimers and transient complexes. Parameters related to the atomic packing suggest that many of the weak dimer interfaces are loosely packed, and we suggest that this contributes to their low stability. To evaluate the evolutionary selection pressure on interface residues, we calculated the Shannon entropy of homologous amino acid sequences at 60% sequence identity. In 93% of the homodimers, the interface residues are better conserved than the residues on the protein surface. The weak dimers display the same high degree of interface conservation as other homodimers, but their homologs may be heterodimers as well as homodimers. Their interfaces may be good models in terms of their size, composition, and evolutionary conservation for the labile subunit contacts that allow protein assemblies to share and exchange components, allosteric proteins to undergo quaternary structure transitions, and molecular machines to operate in the cell.

Confidence intervals for fitting of atomic models into low-resolution densities

The fitting of high-resolution structures into low-resolution densities obtained from techniques such as electron microscopy or small-angle X-ray scattering can yield powerful new insights. While several algorithms for achieving optimal fits have recently been developed, relatively little effort has been devoted to developing objective measures for judging the quality of the resulting fits, in particular with regard to the danger of overfitting. Here, a general method is presented for obtaining confidence intervals for atomic coordinates resulting from fitting of atomic resolution domain structures into low-resolution densities using well established statistical tools. It is demonstrated that the resulting confidence intervals are sufficiently accurate to allow meaningful statistical tests and to provide tools for detecting potential overfitting.

Mechanism of the very efficient quenching of tryptophan fluorescence in human gamma D- and gamma S-crystallins: the gamma-crystallin fold may have evolved to protect tryptophan residues from ultraviolet photodamage

Proteins exposed to UV radiation are subject to irreversible photodamage through covalent modification of tryptophans (Trps) and other UV-absorbing amino acids. Crystallins, the major protein components of the vertebrate eye lens that maintain lens transparency, are exposed to ambient UV radiation throughout life. The duplicated beta-sheet Greek key domains of beta- and gamma-crystallins in humans and all other vertebrates each have two conserved buried Trps. Experiments and computation showed that the fluorescence of these Trps in human gammaD-crystallin is very efficiently quenched in the native state by electrostatically enabled electron transfer to a backbone amide [Chen et al. (2006) Biochemistry 45, 11552-11563]. This dispersal of the excited state energy would be expected to minimize protein damage from covalent scission of the excited Trp ring. We report here both experiments and computation showing that the same fast electron transfer mechanism is operating in a different crystallin, human gammaS-crystallin. Examination of solved structures of other crystallins reveals that the Trp conformation, as well as favorably oriented bound waters, and the proximity of the backbone carbonyl oxygen of the n - 3 residues before the quenched Trps (residue n), are conserved in most crystallins. These results indicate that fast charge transfer quenching is an evolved property of this protein fold, probably protecting it from UV-induced photodamage. This UV resistance may have contributed to the selection of the Greek key fold as the major lens protein in all vertebrates.

Molecular evolution of the betagamma lens crystallin superfamily: evidence for a retained ancestral function in gamma N crystallins?

Within the vertebrate eye, betagamma crystallins are extremely stable lens proteins that are uniquely adapted to increase refractory power while maintaining transparency. Unlike alpha crystallins, which are well-characterized, multifunctional proteins that have important functions both in and out of the lens, betagamma lens crystallins are a diverse group of proteins with no clear ancestral or contemporary nonlens role. We carried out phylogenetic and molecular evolutionary analyses of the betagamma-crystallin superfamily in order to study the evolutionary history of the gamma N crystallins, a recently discovered, biochemically atypical family suggested to possess a divergent or ancestral function. By including nonlens, betagamma-motif-containing sequences in our analysis as outgroups, we confirmed the phylogenetic position of the gamma N family as sister to other gamma crystallins. Using maximum likelihood codon models to estimate lineage-specific nonsynonymous-to-synonymous rate ratios revealed strong positive selection in all of the early lineages within the betagamma family, with the striking exception of the lineage leading to the gamma N crystallins which was characterized by strong purifying selection. Branch-site analysis, used to identify candidate sites involved in functional divergence between gamma N crystallins and its sister clade containing all other gamma crystallins, identified several positively selected changes at sites of known functional importance in the betagamma crystallin protein structure. Further analyses of a fish-specific gamma N crystallin gene duplication revealed a more recent episode of positive selection in only one of the two descendant lineages (gamma N2). Finally, from the guppy, Poecilia reticulata, we isolated complete gamma N1 and gamma N2 coding sequence data from cDNA and partial coding sequence data from genomic DNA in order to confirm the presence of a novel gamma N2 intron, discovered through data mining of two pufferfish genomes. We conclude that the function of the gamma N family likely resembles the ancestral vertebrate betagamma crystallin more than other betagamma families. Furthermore, owing to the presence of an additional intron in some fish gamma N2 crystallins, and the inferred action of positive selection following the fish-specific gamma N duplication, we suggest that further study of fish gamma N crystallins will be critical in further elucidating possible ancestral functions of gamma N crystallins and any nonstructural role they may have.

Biophysical properties of gammaC-crystallin in human and mouse eye lens: the role of molecular dipoles

The eye lens is packed with soluble crystallin proteins, providing a lifetime of transparency and light refraction. gamma-Crystallins are major components of the dense, high refractive index central regions of the lens and generally have high solubility, high stability and high levels of cysteine residues. Human gammaC belongs to a group of gamma-crystallins with a pair of cysteine residues at positions 78 and 79. Unlike other gamma-crystallins it has relatively low solubility, whereas mouse gammaC, which has the exposed C79 replaced with arginine, and a novel mouse splice variant, gammaCins, are both highly soluble. Furthermore, human gammaC is extremely stable, while the mouse orthologs are less stable. Evolutionary pressure may have favoured stability over solubility for human gammaC and the reverse for the orthologs in the mouse. Mutation of C79 to R79, in human gammaC, greatly increased solubility, however, neither form produced crystals. Remarkably, when the human gammaD R36S crystallization cataract mutation was mimicked in human gammaC-crystallin, the solubility of gammaC was dramatically increased, although it still did not crystallize. The highly soluble mouse gammaC-crystallin did crystallize. Its X-ray structure was solved and used in homology modelling of human gammaC, and its mutants C79R and R36S. The human gammaD R36S mutant was also modelled from human gammaD coordinates. Molecular dynamics simulation of the six molecules in the solution state showed that the human gammaCs differed from gammaDs in domain pairing, behaviour that correlates with interface sequence changes. When the fluctuations of the calculated molecular dipoles, for the six structures, over time were analysed, characteristic patterns for soluble gammaC and gammaD proteins were observed. Individual sequence changes that increase or decrease solubility correlated well with changes in the magnitude and direction of these dipoles. It is suggested that changes in surface residues have allowed adaptation for the differing needs of human and mouse lenses.

A reference dataset for circular dichroism spectroscopy tailored for the βγ-crystallin lens proteins

Circular dichroism (CD) spectroscopy is a powerful solution technique for the study of protein secondary structure. As hierarchical euclidean clustering analyses of high quality crystallin synchrotron radiation circular dichroism (SRCD) spectral data can be separated into structural groups based solely on spectral information, the technique can potentially be improved to more accurately determine secondary structures and monitor conformational changes in crystallins. Secondary structure estimates can be determined through use of reference datasets of circular dichroism spectra from proteins with determined crystal structures. As with any empirical method, the accuracies of the analyses are dependent upon how closely the reference dataset characteristics match those of the protein to be studied. To date, crystallin proteins have not been well analysed by CD because existing reference datasets do not contain good representations of their structural characteristics. This work describes a betagamma-crystallin specific reference dataset, CRYST175, which was created solely for the study of betagamma-crystallin secondary structures. Prediction accuracy was assessed for the new dataset using several deconvolution algorithms and it was found to substantially outperform existing more general reference datasets.

Solution structure of (gamma)S-crystallin by molecular fragment replacement NMR

The solution structure of murine gammaS-crystallin (gammaS) has been determined by multidimensional triple resonance NMR spectroscopy, using restraints derived from two sets of dipolar couplings, recorded in different alignment media, and supplemented by a small number of NOE distance restraints. gammaS consists of two topologically similar domains, arranged with an approximate twofold symmetry, and each domain shows close structural homology to closely related (approximately 50% sequence identity) domains found in other members of the gamma-crystallin family. Each domain consists of two four-strand "Greek key" beta-sheets. Although the domains are tightly anchored to one another by the hydrophobic surfaces of the two inner Greek key motifs, the N-arm, the interdomain linker and several turn regions show unexpected flexibility and disorder in solution. This may contribute entropic stabilization to the protein in solution, but may also indicate nucleation sites for unfolding or other structural transitions. The method used for solving the gammaS structure relies on the recently introduced molecular fragment replacement method, which capitalizes on the large database of protein structures previously solved by X-ray crystallography and NMR.

Cold-stable eye lens crystallins of the Antarctic nototheniid toothfish Dissostichus mawsoni Norman

The eye lenses of the Antarctic nototheniid fishes that inhabit the perennially freezing Antarctic seawater are transparent at -2 degrees C, whereas the cold-sensitive mammalian and tropical fish lenses display cold-induced cataract at 20 degrees C and 7 degrees C, respectively. No cold-cataract occurs in the giant Antarctic toothfish Dissostichus mawsoni lens when cooled to temperatures as low as -12 degrees C, indicating highly cold-stable lens proteins. To investigate this cold stability, we characterised the lens crystallin proteins of the Antarctic toothfish, in parallel with those of the sub-tropical bigeye tuna Thunnus obesus and the endothermic cow Bos taurus, representing three disparate thermal climes (-2 degrees C, 18 degrees C and 37 degrees C, respectively). Sizing chromatography resolved their lens crystallins into three groups, alpha/betaH, beta and gamma, with gamma crystallins being the most abundant (>40%) lens proteins in fish, in contrast to the cow lens where they comprise only 19%. The upper thermal stability of these crystallin components correlated with the body temperature of the species. In vitro chaperone assays showed that fish alpha crystallin can protect same-species gamma crystallins from heat denaturation, as well as lysozyme from DTT-induced unfolding, and therefore are small Heat Shock Proteins (sHSP) like their mammalian counterparts. Dynamic light scattering measured an increase in size of alphagamma crystallin mixtures upon heating, which supports formation of the alphagamma complex as an integral part of the chaperone process. Surprisingly, in cross-species chaperone assays, tuna alpha crystallins only partly protected toothfish gamma crystallins, while cow alpha crystallins completely failed to protect, indicating partial and no alphagamma interaction, respectively. Toothfish gamma was likely to be the component that failed to interact, as the supernatant from a cow alpha plus toothfish gamma incubation could chaperone cow gamma crystallins in a subsequent heat incubation, indicating the presence of uncomplexed cow alpha. This suggests that the inability of toothfish gamma crystallins to fully complex with tuna alpha, and not at all with the cow alpha crystallins, may have its basis in adaptive changes in the protein that relate to the extreme cold-stability of the toothfish lens.

Sulfur in human crystallins

Molecular models of human gamma-crystallins and the 'alpha-crystallin domain' of human alphaA-crystallin have been built based on available related X-ray crystal structures. The accessibilities of the component cysteine, methionine and tryptophan side chains in the crystallin models have been calculated. The reactivities of these cysteines, which are oxidised in cataract, are assessed based on their known modifications and within the context of their location within the 3D models.

Ageing and vision: structure, stability and function of lens crystallins

The alpha-, beta- and gamma-crystallins are the major protein components of the vertebrate eye lens, alpha-crystallin as a molecular chaperone as well as a structural protein, beta- and gamma-crystallins as structural proteins. For the lens to be able to retain life-long transparency in the absence of protein turnover, the crystallins must meet not only the requirement of solubility associated with high cellular concentration but that of longevity as well. For proteins, longevity is commonly assumed to be correlated with long-term retention of native structure, which in turn can be due to inherent thermodynamic stability, efficient capture and refolding of non-native protein by chaperones, or a combination of both. Understanding how the specific interactions that confer intrinsic stability of the protein fold are combined with the stabilizing effect of protein assembly, and how the non-specific interactions and associations of the assemblies enable the generation of highly concentrated solutions, is thus of importance to understand the loss of transparency of the lens with age. Post-translational modification can have a major effect on protein stability but an emerging theme of the few studies of the effect of post-translational modification of the crystallins is one of solubility and assembly. Here we review the structure, assembly, interactions, stability and post-translational modifications of the crystallins, not only in isolation but also as part of a multi-component system. The available data are discussed in the context of the establishment, the maintenance and finally, with age, the loss of transparency of the lens. Understanding the structural basis of protein stability and interactions in the healthy eye lens is the route to solve the enormous medical and economical problem of cataract.

High-resolution X-ray crystal structures of human gammaD crystallin (1.25 A) and the R58H mutant (1.15 A) associated with aculeiform cataract

Several human cataracts have been linked to mutations in the gamma crystallin gene. One of these is the aculeiform cataract, which is caused by an R58H mutation in gammaD crystallin. We have shown previously that this cataract is caused by crystallization of the mutant protein, which is an order of magnitude less soluble than the wild-type. Here, we report the very high-resolution crystal structures of the mutant and wild-type proteins. Both proteins crystallize in the same space group and lattice. Thus, a strict comparison of the protein-protein and protein-water intermolecular interactions in the two crystal lattices is possible. Overall, the differences between the mutant and wild-type structures are small. At position 58, the mutant protein loses the direct ion-pair intermolecular interaction present in the wild-type, due to the differences between histidine and arginine at the atomic level; the interaction in the mutant is mediated by water molecules. Away from the mutation site, the mutant and wild-type lattice structures differ in the identity of side-chains that occupy alternate conformations. Since the interactions in the crystal phase are very similar for the two proteins, we conclude that the reduction in the solubility of the mutant is mainly due to the effect of the R58H mutation in the solution phase. The results presented here are also important as they are the first high-resolution X-ray structures of human gamma crystallins.

In vitro unfolding, refolding, and polymerization of human gammaD crystallin, a protein involved in cataract formation

Human gammaD crystallin (HgammaD-Crys), a major protein of the human eye lens, is a primary component of cataracts. This 174-residue primarily beta-sheet protein is made up of four Greek keys separated into two domains. Mutations in the human gene sequence encoding HgammaD-Crys are implicated in early-onset cataracts in children, and the mutant protein expressed in Escherichia coli exhibits properties that reflect the in vivo pathology. We have characterized the unfolding, refolding, and competing aggregation of human wild-type HgammaD-Crys as a function of guanidinium hydrochloride (GuHCl) concentration at neutral pH and 37 degrees C, using intrinsic tryptophan fluorescence to monitor in vitro folding. Wild-type HgammaD-Crys exhibited reversible refolding above 1.0 M GuHCl. The GuHCl unfolded protein was more fluorescent than its native counterpart despite the absence of metal or ion-tryptophan interactions. Aggregation of refolding intermediates of HgammaD-Crys was observed in both equilibrium and kinetic refolding processes. The aggregation pathway competed with productive refolding at denaturant concentrations below 1.0 M GuHCl, beyond the major conformational transition region. Atomic force microscopy of samples under aggregating conditions revealed the sequential appearance of small nuclei, thin protofibrils, and fiber bundles. The HgammaD-Crys fibrous aggregate species bound bisANS appreciably, indicating the presence of exposed hydrophobic pockets. The mechanism of HgammaD-Crys aggregation may provide clues to understanding age-onset cataract formation in vivo.

Homology models of human gamma-crystallins: structural study of the extensive charge network in gamma-crystallins

The lens is composed of highly stable and long-lived proteins, the crystallins which are divided into alpha-, beta-, and gamma-crystallins. Human gamma-crystallins belong to the betagamma superfamily. A large number of gamma-crystallins have been sequenced and have been found to share remarkable sequence homology with each other. Some of the gamma-crystallins from various sources have also been elucidated structurally by X-ray crystallographic or NMR spectroscopic experiments. Their three-dimensional structures are also similar having consisted of two domains each possessing two Greek key motifs. In this study we have constructed the comparative or homology models of the four major human gamma-crystallins, gammaA-,gammaB-, gammaC-, and gammaD-crystallins and studied the charge network in these crystallins. Despite an overall structural similarity between these crystallins, differences in the ion pair formation do exist which is partly due to the differences in their primary sequence and partly due to the structural orientation of the neighboring amino acids. In this study, we present an elaborate analysis of these charged interactions and their formation or loss with respect to the structural changes.

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.

Enhanced crystallization of the Cys18 to Ser mutant of bovine gammaB crystallin

The cysteine residues of the gamma crystallins, a family of ocular lens proteins, are involved in the aggregation and phase separation of these proteins. Both these phenomena are implicated in cataract formation. We have used bovine gammaB crystallin as a model system to study the role of the individual cysteine residues in the aggregation and phase separation of the gamma crystallins. Here, we compare the thermodynamic and kinetic behavior of the recombinant wild-type protein (WT) and the Cys18 to Ser (C18S) mutant. We find that the solubilities of the two proteins are similar. The kinetics of crystallization, however, are different. The WT crystallizes slowly enough for the metastable liquid-liquid coexistence to be easily observed. C18S, on the other hand, crystallizes rapidly; the metastable coexisting liquid phases of the pure mutant do not form. Nevertheless, the coexistence curve of C18S can be determined provided that crystallization is kinetically suppressed. In this way we found that the coexistence curve coincides with that of the WT. Despite the difference in the kinetics of crystallization, the two proteins were found to have the same crystal forms and almost identical X-ray structures. Our results demonstrate that even conservative point mutations can bring about dramatic changes in the kinetics of crystallization. The implications of our findings for cataract formation and protein crystallization are discussed.

Association behaviour of human betaB1-crystallin and its truncated forms

betaB1-crystallin plays an important role in the assembly of betaH-crystallin yet is known to be subject to N-terminal sequence truncations during human lens development and ageing. Here we have over-expressed human betaB1-crystallin, and various truncated forms in Escherichia coli and used mass spectrometry to monitor the monomer molecular weight. Gel permeation chromatography and laser light scattering have been used to estimate the assembly size of the various polypeptides as a function of protein concentration. The full-length betaB1-crystallin behaves as a dimer, like recombinant human betaB2-crystallin, but undergoes further self-association at high protein concentrations, unlike the betaB2-crystallin. Major truncations from the N-terminal extension lead to anomalous behaviour on gel permeation chromatography indicative of altered interactions with the column matrix, whereas light scattering indicated dimers at low protein concentration that self-associate as a function of protein concentration. Loss of 41 residues from the N-terminus, equivalent to an in vivo truncation site, resulted in temperature-dependent phase separation behaviour of the shortened betaB1-crystallin. Good crystals have been grown of a truncated version of human betaB1-crystallin using an in vitro cleavage protocol.

Lens crystallins and oxidation: the special case of gammaS

Among lens crystallins, gamma-crystallins are particularly sensitive to oxidation, because of their high amount of Cys and Met residues. They have the reputation to induce, upon ageing, lens structural modifications leading to opacities. A combination of small angle X-ray scattering and chromatography was used to study the oxidation of gamma-crystallins. At pH 7.0, all the gamma-crystallins under study were checked to have the same structure in solution. Under gentle oxidation conditions at pH 8.0, human gammaS (hgammaS) and bovine gammaS (bgammaS) formed disulfide-linked dimers, whereas the other bgamma-crystallins did not. Cys20 was shown to be responsible for dimer formation since the C20S mutant only formed monomers. The hgammaS dimers were stable for weeks and did not form higher oligomers. In contrast, monomeric gammaS-crystallins freshly prepared at pH 8.0, and submitted to more drastic oxidation by X-ray induced free radicals, were rapidly transformed into higher oligomers. So, only extensive oxidation causing partial unfolding could be detrimental to the lens and linked to cataract formation. The gammaS-crystallins lack the temperature-induced opacification observed with the other gamma-crystallins and known as cold cataract. The oxidation-induced associative behaviour and cold cataract are therefore demonstrated to be uncoupled.

Structure of the crystallins

The lens is formed from two protein superfamilies, the alpha- and beta gamma-crystallins. Representative three-dimensional structures show they both have a basic 2-beta-sheet domain fold, with the beta gamma-domain being made from two intercalating Greek keys. X-ray structures of monomeric gamma-crystallins and simple oligomeric beta-crystallins show how multiple gene duplications can give rise to highly symmetrical assemblies based on paired domains. These protein folds have been engineered by directed mutagenesis to investigate the roles of the critical region in domain pairing and assembly. Inherited human cataracts have been described that are associated with representatives of each of the crystallin protein families. Mutations to certain beta- and gamma-crystallin genes cause expression of truncated polypeptides that would not be expected to fold properly; instead they would randomly aggregate causing light scattering. As crystallin proteins are not renewed, age-related cataract is a gradual accumulation of small changes to pre-existing normal proteins. The precise sites of post-translational modifications are now being mapped to the various crystallins.

Thermodynamic and kinetic characterization of calf lens gammaF-crystallin

Gamma-crystallin is reported to be conformationally stable because of its internal structural symmetry, and gammaF (gammaIVa) is the most stable among the various gamma-crystallin gene products. However, there is no detailed report on its thermodynamic and kinetic stability. In the present study, detailed unfolding of gammaF-crystallin was investigated by equilibrium and kinetics methods with fluorescence and far-UV CD spectroscopic measurements. The GdnHCl-induced unfolding curves probed by Trp emission maximum and intensity showed a sharp single-step transition. Upon widening the unfolding transition with the use of urea in 1.5 M GdnHCl, a more proper fit for thermodynamic analysis was obtained. GammaF-Crystallin underwent a straightforward two-state process (N <==> U) without showing any measurable amount of intermediate. The conformational stability, as measured by deltaG(D)H2O (approximately 9 kcal/mol), indicates that gammaF-crystallin is a very stable protein. The high activation energy deltaG++H2O (approximately 24 kcal/mol), calculated from unfolding kinetics monitored by far-UV CD at 218 nm, also indicates that the native and unfolded states are separated by a high activation energy barrier.

Three murine cataract mutants (Cat2) are defective in different gamma-crystallin genes

A number of murine cataract mutations have been localized to chromosome 1 close to the gamma-crystallin gene cluster (Cryg) (Everett et al., 1994, Genomics 20: 429-434; Löster et al., 1994, Genomics 23: 240-242). Based on the size of the mapping or allelism tests they have not been shown to be genetically distinct and have been assigned to locus symbol Cat2. Here we assign three mutations to the respective gamma-crystallin gene. Using a systematic candidate gene approach to analyze the entire Cryg cluster, an A-->G transition was found in exon 2 of Cryga for the ENU-436 mutation and is designated Cryga1Neu. The mutant allele Crygbnop (formerly Cat2(nop)) is caused by a replacement of 11 bp by 4 bp in the third exon of Crygb, while a C-->G transversion in exon 3 of Cryge has been found for the Cryget (formerly Cat2(t)) mutation. For the mutation Cryga1Neu, an Asp-->Gly exchange is deduced, whereas the mutations Crygbnop and Cryget lead to the formation of in-frame stop codons and give rise to truncated proteins of 144 and 143 amino acids, respectively. The effects of the mutations upon gamma-crystallin structure are likely to be quite different. The Cryga1Neu mutation is expected to affect the link between Greek-key motifs 2 and 3, whereas both Crygbnop and Cryget mutations are supposed to truncate the fourth Greek-key motif. All three mutations are predicted to alter protein folding of the gamma-crystallins and result in lens cataract, but the phenotype for each is quite distinctive.

X-ray structures of three interface mutants of gammaB-crystallin from bovine eye lens

GammaB-crystallin consists of two domains each comprising two "Greek key" motifs. Both domains fold independently, and domain interactions contribute significantly to the stability of the C-terminal domain. In a previous study (Palme S et al., 1996, Protein Sci 6:1529-1636) it was shown that Phe56 from the N-terminal domain, a residue involved in forming a hydrophobic core at the domain interface, effects the interaction of the two domains, and therefore, the stability of the C-terminal domain. Ala or Asp at position 56 drastically decreased the stability of the C-terminal domain, whereas Trp had a more moderate effect. In this article we present the X-ray structures of these interface mutants and correlate them with the stability data. The mutations do not effect the overall structure of the molecule. No structural changes are observed in the vicinity of the replaced residue, suggesting that the local structure is too rigid to allow compensations for the amino acid replacements. In the mutants gammaB-F56A and -F56D, a solvent-filled groove accessible to the bulk solvent is created by the replacement of the bulky Phe side chain. In gammaB-F56W, the pyrrole moiety of the indole ring replaces the phenyl side chain of the wild type. With the exception of gammaB-F56W, there is a good correlation between the hydrophobicity of the amino acid at position 56 according to the octanol scale and the stability of the C-terminal domain. In gammaB-F56W, the C-terminal domain is less stable than estimated from the hydrophobicity, presumably because the ring nitrogen (Nepsilon1) has no partner to form hydrogen bonds. The data suggest that the packing of hydrophobic residues in the interface core is important for domain interactions and the stability of gammaB-crystallin. Apparently, for protein stability, the same principles apply for hydrophobic cores within domains and at domain interfaces.