Ranjith Padinhateeri

Elucidating the role of disulfide bond on amyloid formation and fibril reversibility of somatostatin-14: Relevance to its storage and secretion

Background: Peptide/protein hormones are stored as amyloids within endocrine secretory granules. Results: Disulfide bond cleavage enhances conformational dynamics and aggregation kinetics in somatostatin-14, resulting in amyloid fibrils with increased resistance to denaturing conditions and decreased reversibility. Conclusion: Disulfide bond could be a key modulating factor in somatostatin-14 amyloid formation associated with secretory granule biogenesis. Significance: Defective disulfide bonding might cause dysregulation of hormone storage/secretion. The storage of protein/peptide hormones within subcellular compartments and subsequent release are crucial for their native function, and hence these processes are intricately regulated in mammalian systems. Several peptide hormones were recently suggested to be stored as amyloids within endocrine secretory granules. This leads to an apparent paradox where storage requires formation of aggregates, and their function requires a supply of non-aggregated peptides on demand. The precise mechanism behind amyloid formation by these hormones and their subsequent release remain an open question. To address this, we examined aggregation and fibril reversibility of a cyclic peptide hormone somatostatin (SST)-14 using various techniques. After proving that SST gets stored as amyloid in vivo, we investigated the role of native structure in modulating its conformational dynamics and self-association by disrupting the disulfide bridge (Cys3–Cys14) in SST. Using two-dimensional NMR, we resolved the initial structure of somatostatin-14 leading to aggregation and further probed its conformational dynamics in silico. The perturbation in native structure (S-S cleavage) led to a significant increase in conformational flexibility and resulted in rapid amyloid formation. The fibrils formed by disulfide-reduced noncyclic SST possess greater resistance to denaturing conditions with decreased monomer releasing potency. MD simulations reveal marked differences in the intermolecular interactions in SST and noncyclic SST providing plausible explanation for differential aggregation and fibril reversibility observed experimentally in these structural variants. Our findings thus emphasize that subtle changes in the native structure of peptide hormone(s) could alter its conformational dynamics and amyloid formation, which might have significant implications on their reversible storage and secretion.

Random hydrolysis controls the dynamic instability of microtubules.

Uncovering mechanisms that control the dynamics of microtubules is fundamental for our understanding of multiple cellular processes such as chromosome separation and cell motility. Building on previous theoretical work on the dynamic instability of microtubules, we propose here a stochastic model that includes all relevant biochemical processes that affect the dynamics of microtubule plus-end, namely, the binding of GTP-bound monomers, unbinding of GTP- and GDP-bound monomers, and hydrolysis of GTP monomers. The inclusion of dissociation processes, present in our approach but absent from many previous studies, is essential to guarantee the thermodynamic consistency of the model. Our theoretical method allows us to compute all dynamic properties of microtubules explicitly. Using experimentally determined rates, it is found that the cap size is ∼3.6 layers, an estimate that is compatible with several experimental observations. In the end, our model provides a comprehensive description of the dynamic instability of microtubules that includes not only the statistics of catastrophes but also the statistics of rescues.

Characterization of Amyloid Formation by Glucagon-Like Peptides: Role of Basic Residues in Heparin-Mediated Aggregation

Glycosaminoglycans (GAGs) have been reported to play a significant role in amyloid formation of a wide range of proteins/peptides either associated with diseases or native biological functions. The exact mechanism by which GAGs influence amyloid formation is not clearly understood. Here, we studied two closely related peptides, glucagon-like peptide 1 (GLP1) and glucagon-like peptide 2 (GLP2), for their amyloid formation in the presence and absence of the representative GAG heparin using various biophysical and computational approaches. We show that the aggregation and amyloid formation by these peptides follow distinct mechanisms: GLP1 follows nucleation-dependent aggregation, whereas GLP2 forms amyloids without any significant lag time. Investigating the role of heparin, we also found that heparin interacts with GLP1, accelerates its aggregation, and gets incorporated within its amyloid fibrils. In contrast, heparin neither affects the aggregation kinetics of GLP2 nor gets embedded within its fibrils. Furthermore, we found that heparin preferentially influences the stability of the GLP1 fibrils over GLP2 fibrils. To understand the specific nature of the interaction of heparin with GLP1 and GLP2, we performed all-atom MD simulations. Our in silico results show that the basic-nonbasic-basic (B-X-B) motif of GLP1 (K28-G29-R30) facilitates the interaction between heparin and peptide monomers. However, the absence of such a motif in GLP2 could be the reason for a significantly lower strength of interaction between GLP2 and heparin. Our study not only helps to understand the role of heparin in inducing protein aggregation but also provides insight into the nature of heparin-protein interaction.

Nucleosome positioning in a model of active chromatin remodeling enzymes

Accounting for enzyme-mediated active sliding, disassembly, and sequence-dependent positioning of nucleosomes, we simulate nucleosome occupancy over cell-cycle-scale times using a stochastic kinetic model. We show that ATP-dependent active nucleosome sliding and nucleosome removal processes are essential to obtain in vivo-like nucleosome positioning. While active sliding leads to dense nucleosome filling, sliding events alone cannot ensure sequence-dependent nucleosome positioning: Active nucleosome removal is the crucial remodeling event that drives positioning. We also show that remodeling activity changes nucleosome dynamics from glassy to liquid-like, and that remodeling dramatically influences exposure dynamics of promoter regions.

Nucleosome positioning and kinetics near transcription-start-site barriers are controlled by interplay between active remodeling and DNA sequence

We investigate how DNA sequence, ATP-dependent chromatin remodeling and nucleosome-depleted ‘barriers’ co-operate to determine the kinetics of nucleosome organization, in a stochastic model of nucleosome positioning and dynamics. We find that ‘statistical’ positioning of nucleosomes against ‘barriers’, hypothesized to control chromatin structure near transcription start sites, requires active remodeling and therefore cannot be described using equilibrium statistical mechanics. We show that, unlike steady-state occupancy, DNA site exposure kinetics near a barrier is dominated by DNA sequence rather than by proximity to the barrier itself. The timescale for formation of positioning patterns near barriers is proportional to the timescale for active nucleosome eviction. We also show that there are strong gene-to-gene variations in nucleosome positioning near barriers, which are eliminated by averaging over many genes. Our results suggest that measurement of nucleosome kinetics can reveal information about sequence-dependent regulation that is not apparent in steady-state nucleosome occupancy.

Stretching and Bending Fluctuations of Short DNA Molecules

Recent measurements of the distribution of end-to-end distance in short DNA molecules infer cooperative stretching fluctuations. The assumptions underlying the analysis can be questioned if transient, thermally induced defects producing a localized decrease in bending stiffness are present in thermal equilibrium, such as regions in which DNA melts locally (bubbles), sustains large-angle bends (kinks), or can locally transform into an alternative (S-DNA) state. We study a generalized discrete worm-like chain model for DNA, capable of describing these experiments, showing that the model yields accurate fits to available experimental data. Our results indicate that DNA bending arising from such localized defects, rather than solely stretching, can be an equal contributor to end-to-end distance fluctuations for 35-bp DNA and contributes nontrivially to such fluctuations at all scales below the persistence length. The analysis suggests that such fluctuations should exhibit a scale-dependent cooperativity, specifically relevant in determining the behavior of short chains, but which saturates rapidly to a length-independent value for longer DNA, to ensure a consistent physical description of DNA across multiple scales. Our approach provides a minimal, yet accurate, coarse-grained description of DNA at the subpersistence length scales of current experimental interest.

Investigating the intrinsic aggregation potential of evolutionarily conserved segments in p53.

Protein aggregation and amyloid formation are known to play a role both in diseases and in biological functions. Transcription factor p53 plays a major role in tumor suppression by maintaining genomic stability. Recent studies have suggested that amyloid formation of p53 could lead to its loss of physiological function as a tumor suppressor. Here, we investigated the intrinsic amyloidogenic nature of wild-type p53 using sequence analysis. We used bioinformatics and aggregation prediction algorithms to establish the evolutionarily conserved nature of aggregation-prone sequences in wild-type p53. Further, we analyzed the amyloid forming capacity of conserved and aggregation-prone p53-derived peptides PILTIITL and YFTLQI in vitro using various biophysical techniques, including all atom molecular dynamics simulation. Finally, we probed the seeding ability of the PILTIITL peptide on p53 aggregation in vitro and in cells. Our data demonstrate the intrinsic amyloid forming ability of a sequence stretch of the p53 DNA binding domain (DBD) and its aggregation templating behavior on full-length and p53 core domain. Therefore, p53 aggregation, instigated through an amyloidogenic segment in its DBD, could be a putative driving force for p53 aggregation in vivo.

Collective force generated by multiple biofilaments can exceed the sum of forces due to individual ones

Collective dynamics and force generation by cytoskeletal filaments are crucial in many cellular processes. Investigating growth dynamics of a bundle of N independent cytoskeletal filaments pushing against a wall, we show that chemical switching (ATP/GTP hydrolysis) leads to a collective phenomenon that is currently unknown. Obtaining force-velocity relations for different models that capture chemical switching, we show, analytically and numerically, that the collective stall force of N filaments is greater than N times the stall force of a single filament. Employing an exactly solvable toy model, we analytically prove the above result for N = 2. We, further, numerically show the existence of this collective phenomenon, for , in realistic models (with random and sequential hydrolysis) that simulate actin and microtubule bundle growth. We make quantitative predictions for the excess forces, and argue that this collective effect is related to the non-equilibrium nature of chemical switching.

Molecular interpretation of ACTH-β-endorphin coaggregation: relevance to secretory granule biogenesis.

Peptide/protein hormones could be stored as non-toxic amyloid-like structures in pituitary secretory granules. ACTH and β-endorphin are two of the important peptide hormones that get co-stored in the pituitary secretory granules. Here, we study molecular interactions between ACTH and β-endorphin and their colocalization in the form of amyloid aggregates. Although ACTH is known to be a part of ACTH-β-endorphin aggregate, ACTH alone cannot aggregate into amyloid under various plausible conditions. Using all atom molecular dynamics simulation we investigate the early molecular interaction events in the ACTH-β-endorphin system, β-endorphin-only system and ACTH-only system. We find that β-endorphin and ACTH formed an interacting unit, whereas negligible interactions were observed between ACTH molecules in ACTH-only system. Our data suggest that ACTH is not only involved in interaction with β-endorphin but also enhances the stability of mixed oligomers of the entire system.

Theoretical estimates of exposure timescales of protein binding sites on DNA regulated by nucleosome kinetics

It is being increasingly realized that nucleosome organization on DNA crucially regulates DNA–protein interactions and the resulting gene expression. While the spatial character of the nucleosome positioning on DNA has been experimentally and theoretically studied extensively, the temporal character is poorly understood. Accounting for ATPase activity and DNA-sequence effects on nucleosome kinetics, we develop a theoretical method to estimate the time of continuous exposure of binding sites of non-histone proteins (e.g. transcription factors and TATA binding proteins) along any genome. Applying the method to Saccharomyces cerevisiae, we show that the exposure timescales are determined by cooperative dynamics of multiple nucleosomes, and their behavior is often different from expectations based on static nucleosome occupancy. Examining exposure times in the promoters of GAL1 and PHO5, we show that our theoretical predictions are consistent with known experiments. We apply our method genome-wide and discover huge gene-to-gene variability of mean exposure times of TATA boxes and patches adjacent to TSS (+1 nucleosome region); the resulting timescale distributions have non-exponential tails.