Using Schematic Models to Understand the Microscopic Basis for Inverted Solubility in γD-Crystallin

Insights into the solubility of γ $$ \gamma $$ D-crystallin from multiscale atomistic simulations

The molecular basis underlying the rich phase behavior of globular proteins remains poorly understood. We use atomistic multiscale molecular simulations to model the solution-state conformational dynamics and interprotein interactions of γ $$ \gamma $$ D-crystallin and its P23T-R36S mutant, which drastically limits the protein solubility, at both infinite dilution and at a concentration where the mutant fluid phase and crystalline phase coexist. We find that while the mutant conserves the protein fold, changes to the surface exposure of residues in the neighborhood of residue-36 enhance protein-protein interactions and develop specific protein-protein contacts found in the protein crystal lattice.

Network Hamiltonian Models for Unstructured Protein Aggregates, with Application to γD-Crystallin

Network Hamiltonian models (NHMs) are a framework for topological coarse-graining of protein-protein interactions, in which each node corresponds to a protein, and edges are drawn between nodes representing proteins that are noncovalently bound. Here, this framework is applied to aggregates of γD-crystallin, a structural protein of the eye lens implicated in cataract disease. The NHMs in this study are generated from atomistic simulations of equilibrium distributions of wild-type and the cataract-causing variant W42R in solution, performed by Wong, E. K.; Prytkova, V.; Freites, J. A.; Butts, C. T.; Tobias, D. J. Molecular Mechanism of Aggregation of the Cataract-Related γD-Crystallin W42R Variant from Multiscale Atomistic Simulations. Biochemistry2019, 58 (35), 3691-3699. Network models are shown to successfully reproduce the aggregate size and structure observed in the atomistic simulation, and provide information about the transient protein-protein interactions therein. The system size is scaled from the original 375 monomers to a system of 10000 monomers, revealing a lowering of the upper tail of the aggregate size distribution of the W42R variant. Extrapolation to higher and lower concentrations is also performed. These results provide an example of the utility of NHMs for coarse-grained simulation of protein systems, as well as their ability to scale to large system sizes and high concentrations, reducing computational costs while retaining topological information about the system.

Surface Exposed Free Cysteine Suppresses Crystallization of Human γD-Crystallin

Human γD-crystallin (HGD) has remarkable stability against condensation in the human lens, sometimes over a whole lifetime. The native protein has a surface exposed free cysteine that forms dimers (Benedek, 1997; Ramkumar et al., 1864)^1,2 without specific biological function and leads to further protein association and/or aggregation, which creates a paradox for understanding its stability. Previous work has demonstrated that chemical modification of the protein at the free cysteine (C110), increases the temperature at which liquid-liquid phase separation occurs (LLPS), lowers protein solubility and suggests an important role for this amino acid in maintaining its long-term resistance to condensation. Here we demonstrate that mutation of the cysteine does not alter the structure or solubility (liquidus) line for the protein, but dramatically increases the protein crystal nucleation rate following LLPS, suggesting that the free cysteine has a vital role in suppressing crystallization in the human lens.