The mechanisms leading to aggregation of the crystallin proteins of the eye lens remain largely unknown. We use atomistic multiscale molecular simulations to model the solution-state conformational dynamics of γD-crystallin and its cataract-related W42R variant at both infinite dilution and physiologically relevant concentrations. We find that the W42R variant assumes a distinct conformation in solution that leaves the Greek key domains of the native fold largely unaltered but lacks the hydrophobic interdomain interface that is key to the stability of wild-type γD-crystallin. At physiologically relevant concentrations, exposed hydrophobic regions in this alternative conformation become primary sites for enhanced interprotein interactions leading to large-scale aggregation.
Background The objective of this study was to determine whether the cataract-related G18V variant of human γS-crystallin has increased exposure of hydrophobic residues that could explain its aggregation propensity and/or recognition by αB-crystallin. Methods We used an ANS fluorescence assay and NMR chemical shift perturbation to experimentally probe exposed hydrophobic surfaces. These results were compared to flexible docking simulations of ANS molecules to the proteins, starting with the solution-state NMR structures of γS-WT and γS-G18V. Results γS-G18V exhibits increased ANS fluorescence, suggesting increased exposed hydrophobic surface area. The specific residues involved in ANS binding were mapped by NMR chemical shift perturbation assays, revealing ANS binding sites in γS-G18V that are not present in γS-WT. Molecular docking predicts three binding sites that are specific to γS-G18V corresponding to the exposure of a hydrophobic cavity located at the interdomain interface, as well as two hydrophobic patches near a disordered loop containing solvent-exposed cysteines, all but one of which is buried in γS-WT. Conclusions Although both proteins display non-specific binding, more residues are involved in ANS binding to γS-G18V, and the affected residues are localized in the N-terminal domain and the nearby interdomain interface, proximal to the mutation site. General Significance Characterization of changes in exposed hydrophobic surface area between wild-type and variant proteins can help elucidate the mechanisms of aggregation propensity and chaperone recognition, presented here in the context of cataract formation. Experimental data and simulations provide complementary views of the interactions between proteins and the small molecule probes commonly used to study aggregation.
Efficient computational modeling of biological systems characterized by high concentrations of macromolecules often relies on rigid-body Brownian Dynamics or Monte Carlo (MC) simulations. However, the accuracy of rigid-body models is limited by the fixed conformation of the simulated biomolecules. Multi-conformation Monte Carlo (mcMC) simulations of protein solutions incorporate conformational flexibility via a conformational swap trial move within a predetermined library of discrete protein structures, thereby alleviating artifacts arising from the use of a single protein conformation. Here, we investigate the impact of the number of distinct protein structures in the conformational library and the extent of conformational sampling used in its generation on structural observables computed from simulations of hen egg white lysozyme (HEWL), human γD-Crystallin, and bovine γB-Crystallin solutions. We find that the importance of specific protocols for the construction of the protein structure library is strongly dependent on the nature of the simulated system.
We reveal that the axial stiffness of amyloid fibrils is inversely correlated with their cross-sectional area. Because amyloid fibrils’ stiffness is determined by hydrogen bond (H-bond) density with a linear correlation, our finding implies that amyloid fibrils with larger radial sizes are generally softer and have lower density H-bond networks. In silico calculations show that the stiffness–size relationship of amyloid fibrils is, indeed, driven by the packing densities of residues and H-bonds. Our results suggest that polypeptide chains which form amyloid fibrils with narrow cross sections can optimize packing densities in the fibrillar core structure, in contrast to those forming wide amyloid fibrils. Consequently, the density of residues and H-bonds that contribute to mechanical stability is higher in amyloid fibrils with narrow cross sections. This size dependence of nanomechanics appears to be a global property of amyloid fibrils, just like the well-known cross-β sheet topology.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.