Effects of Mutating Trp42 Residue on γD-Crystallin Stability

Effect of mutations on the folding and stability of γD-crystallin protein

Interprotein interactions between the partially unfolded states of γD-crystallin (γD-crys) protein are known to cause cataracts. Therefore, understanding the unfolding pathways of native γD-crys is extremely crucial to delineate their aggregation mechanism. In this study, we have performed extensive all-atom Molecular Dynamics simulations with explicit solvent to understand the role of the critical residues that drive the stability of the motifs and domains of γD-crys in its wild type and mutant forms. Our findings show that while the individual motifs of wild type are not stable in the native form, the individual domains remain structurally stable at 425K. This enhanced stability of the domain was attributed to the hydrophobic interactions between the motifs. Single and double point mutations of the domains with negatively charged aspartic and glutamic acid amino acid residues (I3E, W42D, W42E, I3D/W42D, I3E/W42E, and L92D/W157D) decreases the structural stability, leading to unfolding of individual domains of γD-crys. We believe that our study sheds light on the weakest links of γD-crys, along with the role of interactions stabilizing the domains. Further, this study bolsters and provides a better understanding of the domain swapping mechanism of aggregation of γD-crys.Communicated by Ramaswamy H. Sarma.

Conformational stability of the deamidated and mutated human βB2-crystallin

Previous studies propose that genetic mutations and post-translational modifications in protein crystallins promote protein aggregation and are considered significant risk factors for cataract formation. The βB2-crystallin (HβB2C) forms a high proportion of proteins in the human eye lens. Different congenital mutations and post-translational deamidations in βB2-crystallin have been reported and linked to cataract formation. In this work, we employed extensive all-atom molecular dynamics simulations to evaluate the conformational stability of deamidated and mutated HβB2C. Our results show critical changes in the protein surface and its native contacts due to a modification in the conformational equilibrium of these proteins. The double deamidated (Q70E/Q162E) and single deamidated (Q70E) impact the well compact conformation of the HβB2C. These post-translational modifications allow the exposure of the protein hydrophobic interface, which lead to the exposure of electronegative residues. On the other hand, our mutational studies showed that the S143F mutation modifies the hydrogen-bond network of an antiparallel β-sheet, unfolding the C-terminal domain. Interestingly, the chain termination mutation (Q155X) does not unfold the N-terminal domain. However, the resultant conformation is more compact and avoids the exposure of the hydrophobic interface. Our results provide valuable information about the first steps of HβB2C unfolding in the presence of deamidated amino acids that have been reported to appear during aging. The findings reported in this work are essential for the general knowledge of the initial steps in the cataract formation mechanism, which may be helpful for the further development of molecules with pharmacological potential against cataract disease.

Coordinate-Dependent Drift-Diffusion Reveals the Kinetic Intermediate Traps of Top7-Based Proteins

The computer-designed Top7 served as a scaffold to produce immunoreactive proteins by grafting of the 2F5 HIV-1 antibody epitope (Top7-2F5) followed by biotinylation (Top7-2F5-biotin). The resulting nonimmunoglobulin affinity proteins were effective in inducing and detecting the HIV-1 antibody. However, the grafted Top7-2F5 design led to protein aggregation, as opposed to the soluble biotinylated Top7-2F5-biotin. The structure-based model predicted that the thermodynamic cooperativity of Top7 increases after grafting and biotin-labeling, reducing their intermediate state populations. In this work, the folding kinetic traps that might contribute to the aggregation propensity are investigated by the diffusion theory. Since the engineered proteins have similar sequence and structural homology, they served as protein models to study the kinetic intermediate traps that were uncovered by characterizing the position-dependent drift-velocity (v(Q)) and the diffusion (D(Q)) coefficients. These coordinate-dependent coefficients were taken into account to obtain the folding and transition path times over the free energy transition states containing the intermediate kinetic traps. This analysis may be useful to predict the aggregated kinetic traps of scaffold-epitope proteins that might compose novel diagnostic and therapeutic platforms.

A novel cataract-causing mutation Ile82Met of γA crystallin trends to aggregate with unfolding intermediate

Cataract is the most common pathogenic ophthalmic disease leading to blindness in children worldwide. Genetic disorder is the leading cause of congenital cataract, among which crystallin mutations have a high incidence. There are few reports on γA-crystallin, one critical member of crystallin superfamilies. In this study, we identified a novel pathogenic mutation (Ile82Met) in γA-crystallin from a three-generation Chinese family with cataract, and investigated the potential molecular mechanism in detail. To elucidate the pathogenic mechanism of I82M mutant, spectroscopic and solubility experiments were performed to determine the difference between the purified γA-crystallin wild type (WT) and I82M mutant under both physiological conditions and environmental stresses (UV irradiation, thermal denaturation or chemical denaturation). The I82M mutant did not affect the secondary/tertiary structure of monomeric γA-crystallin under physiological status, but decreased protein stability and increased aggregatory potency under the stressful treatment. Surprisingly, the chemical denaturation caused I82M to switch from the two-state unfolding of γA-crystallin to three-state unfolding involving an unfolding intermediate. This study expands the genetic variation map of cataract, and provides novel insights into the pathomechanism, in particular, filling in a gap in the understanding of γA-crystallin mutants causing cataract.

Biotin-painted proteins have thermodynamic stability switched by kinetic folding routes

Biotin-labeled proteins are widely used as tools to study protein-protein interactions and proximity in living cells. Proteomic methods broadly employ proximity-labeling technologies based on protein biotinylation in order to investigate the transient encounters of biomolecules in subcellular compartments. Biotinylation is a post-translation modification in which the biotin molecule is attached to lysine or tyrosine residues. So far, biotin-based technologies proved to be effective instruments as affinity and proximity tags. However, the influence of biotinylation on aspects such as folding, binding, mobility, thermodynamic stability, and kinetics needs to be investigated. Here, we selected two proteins [biotin carboxyl carrier protein (BCCP) and FKBP3] to test the influence of biotinylation on thermodynamic and kinetic properties. Apo (without biotin) and holo (biotinylated) protein structures were used separately to generate all-atom structure-based model simulations in a wide range of temperatures. Holo BCCP contains one biotinylation site, and FKBP3 was modeled with up to 23 biotinylated lysines. The two proteins had their estimated thermodynamic stability changed by altering their energy landscape. In all cases, after comparison between the apo and holo simulations, differences were observed on the free-energy profiles and folding routes. Energetic barriers were altered with the density of states clearly showing changes in the transition state. This study suggests that analysis of large-scale datasets of biotinylation-based proximity experiments might consider possible alterations in thermostability and folding mechanisms imposed by the attached biotins.

Insights into the binding of morin to human γD-crystallin

Crystallin aggregation in the eye lens is one of the leading causes of cataract formation. The increase in the human γD-crystallin (HγDC) aggregation propensity has been associated with the oligomerization of its partially folded and fully unfolded structure. A recent study demonstrated that the binding of flavonoid morin (MOR) to HγDC inhibits the fibrillation of this protein. In this work, we carry out an exhaustive search for the possible binding site of MOR on HγDC by combining an ensemble docking approach with the Wrap 'N' Shake protocol. In agreement with previous results, we found a potential MOR-binding site in the cleft formed between the N-terminal and C-terminal domains of HγDC. MOR preference for the cleft residues was observed even with the interface-opened intermediate conformers of HγDC. In addition, metadynamics simulations were carried out to corroborate the stabilizing activity of MOR on HγDC structure and to identify the structural regions implicated during the unfolding inhibition. Overall, this study provides relevant insights into the identification of new HγDC aggregation inhibitors.

Replica exchange molecular dynamics simulations reveal self-association sites in M-crystallin caused by mutations provide insights of cataract

Crystallins are ubiquitous, however, prevalence is seen in eye lens. Eye lens crystallins are long-lived and structural intactness is required for maintaining lens transparency and protein solubility. Mutations in crystallins often lead to cataract. In this study, we performed mutations at specific sites of M-crystallin, a close homologue of eye lens crystallin and studied by using replica exchange molecular dynamics simulation with generalized Born implicit solvent model. Mutations were made on the Ca²⁺ binding residues (K34D and S77D) and in the hydrophobic core (W45R) which is known to cause congenital cataract in homologous γD-crystallin. The chosen mutations caused large motion of the N-terminal Greek key, concomitantly broke the interlocking Greek keys interactions and perturbed the compact core resulting in several folded and partially unfolded states. Partially unfolded states exposed large hydrophobic patches that could act as precursors for self-aggregation. Accumulation of such aggregates is the potential cause of cataract in homologous eye lens crystallins.

Exploring the folding process of human βB2-crystallin using multiscale molecular dynamics and the Markov state model

Adequate knowledge of protein conformations is crucial for understanding their function and their association properties with other proteins. The cataract disease is correlated with conformational changes in key proteins called crystallins. These changes are due to mutations or post-translational modifications that may lead to protein unfolding, and thus the formation of aggregate states. Human βB2-crystallin (HβB2C) is found in high proportion in the eye lens, and its mutations are related to some cataracts. HβB2C also associates into dimers, tetramers, and other higher-order supramolecular complexes. However, it is the only protein of the βγ-crystallin family that has been found in an extended conformation. Therefore, we hypothesize that the extended conformation is not energetically favourable and that HβB2C may adopt a closed (completely folded) conformation, similar to the other members of the βγ-crystallin family. To corroborate this hypothesis, we performed extensive molecular dynamics simulations of HβB2C in its monomeric and dimeric conformations, using all-atom and coarse-grained scales. We employed Markov state model (MSM) analysis to characterize the conformational and kinetically relevant states in the folding process of monomeric HβB2C. The MSM analysis clearly shows that HβB2C adopts a completely folded structure, and this conformation is the most kinetically and energetically favourable one. In contrast, the extended conformations are kinetically unstable and energetically unfavourable. Our MSM analysis also reveals a key metastable state, which is particularly interesting because it is from this state that the folded state is reached. The folded state is stabilized by the formation of two salt bridges between the residue-pairs E74-R187 and R97-E166 and the two hydrophobic residue-pairs V59-L164 and V72-V151. Furthermore, free energy surface (FES) analysis revealed that the HβB2C dimer with both monomers in a closed conformation (face-en-face dimer) is energetically more stable than the domain-swapped dimer (crystallographic structure). The results presented in this report shed light on the molecular details of the folding mechanism of HβB2C in an aqueous environment and may contribute to interpreting different experimental findings. Finally, a detailed knowledge of HβB2C folding may be key to the rational design of potential molecules to treat cataract disease.