Structural mechanism of HP1⍺-dependent transcriptional repression and chromatin compaction
Not Available
Engineering topics
Publications and source records attributed to Singharoy, Abhishek.
Not Available
Small, diffusible redox proteins play an essential role in electron transfer (ET) in respiration and photosynthesis, sustaining life on Earth by shuttling electrons between membrane-bound complexes via finely tuned and reversible interactions. Ensemble kinetic studies show transient ET complexes form in two distinct stages: an “encounter” complex largely mediated by electrostatic interactions, which subsequently, through subtle reorganization of the binding interface, forms a “productive” ET complex stabilized by additional hydrophobic interactions around the redox-active cofactors. Here, using single-molecule force spectroscopy (SMFS) we dissected the transient ET complexes formed between the photosynthetic reaction center-light harvesting complex 1 (RC-LH1) of Rhodobacter sphaeroides and its native electron donor cytochrome c 2 (cyt c 2 ). Importantly, SMFS resolves the distribution of interaction forces into low (~150 pN) and high (~330 pN) components, with the former more susceptible to salt concentration and to alteration of key charged residues on the RC. Thus, the low force component is suggested to reflect the contribution of electrostatic interactions in forming the initial encounter complex, whereas the high force component reflects the additional stabilization provided by hydrophobic interactions to the productive ET complex. Employing molecular dynamics simulations, we resolve five intermediate states that comprise the encounter, productive ET and leaving complexes, predicting a weak interaction between cyt c 2 and the LH1 ring near the RC-L subunit that could lie along the exit path for oxidized cyt c 2 . The multimodal nature of the interactions of ET complexes captured here may have wider implications for ET in all domains of life.
Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, in this work, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.
Abstract Mutations in the TP53 tumor suppressor gene occur in >80% of the triple-negative or basal-like breast cancer. To test whether neomorphic functions of specific TP53 missense mutations contribute to phenotypic heterogeneity, we characterized phenotypes of non-transformed MCF10A-derived cell lines expressing the ten most common missense mutant p53 proteins and observed a wide spectrum of phenotypic changes in cell survival, resistance to apoptosis and anoikis, cell migration, invasion and 3D mammosphere architecture. The p53 mutants R248W, R273C, R248Q, and Y220C are the most aggressive while G245S and Y234C are the least, which correlates with survival rates of basal-like breast cancer patients. Interestingly, a crucial amino acid difference at one position—R273C vs. R273H—has drastic changes on cellular phenotype. RNA-Seq and ChIP-Seq analyses show distinct DNA binding properties of different p53 mutants, yielding heterogeneous transcriptomics profiles, and MD simulation provided structural basis of differential DNA binding of different p53 mutants. Integrative statistical and machine-learning-based pathway analysis on gene expression profiles with phenotype vectors across the mutant cell lines identifies quantitative association of multiple pathways including the Hippo/YAP/TAZ pathway with phenotypic aggressiveness. Further, comparative analyses of large transcriptomics datasets on breast cancer cell lines and tumors suggest that dysregulation of the Hippo/YAP/TAZ pathway plays a key role in driving the cellular phenotypes towards basal-like in the presence of more aggressive p53 mutants. Overall, our study describes distinct gain-of-function impacts on protein functions, transcriptional profiles, and cellular behaviors of different p53 missense mutants, which contribute to clinical phenotypic heterogeneity of triple-negative breast tumors.
Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge or artificial intelligence driven, and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary- structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models, which contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly, while refining them against the density maps. Here, we introduce such an adaptive decision making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, and the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two protein systems, Adenylate Kinase and Carbon Monoxide Dehydrogenase. For the test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2 - 3 Å, map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that is inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from Chimpanzee adenovirus. The overhead for decision making remaining low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install the R-MDFF on different computing environments, from local Linux based workstations to High Performance Computing (HPC) environments.
The mechanism of rotatory catalysis in ATPhydrolyzing molecular motors remains an unresolved puzzle in biological energy transfer. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, knowledge on how the coupling between the chemical and mechanical steps within motors enforces directional rotatory movements remains fragmentary. Even more contentious is to pinpoint the rate-limiting step of a multistep rotation process. Here, using vacuolar or V 1 -type hexameric ATPase as an exemplary rotational motor, we present a model of the complete 4-step conformational cycle involved in rotatory catalysis. First, using Xray crystallography, a new intermediate or “dwell” is identified, which enables the release of an inorganic phosphate (or P i ) after ATP hydrolysis. Using molecular dynamics simulations, this new dwell is placed in a sequence with three other crystal structures to derive a putative cyclic rotation path. Free-energy simulations are employed to estimate the rate of the hexameric protein transformations and delineate allosteric effects that allow new reactant ATP entry only after hydrolysis product exit. An analysis of transfer entropy brings to light how the side-chain-level interactions transcend into larger-scale reorganizations, highlighting the role of the ubiquitous arginine-finger residues in coupling chemical and mechanical information. An inspection of all known rates encompassing the 4-step rotation mechanism implicates the overcoming of the ADP interactions with V 1 -ATPase to be the rate-limiting step of motor action.
Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.
Vaccines derived from chimpanzee adenovirus Y25 (ChAdOx1), human adenovirus type 26 (HAdV-D26), and human adenovirus type 5 (HAdV-C5) are critical in combatting the severe acute respiratory coronavirus 2 (SARS-CoV-2) pandemic. As part of the largest vaccination campaign in history, ultrarare side effects not seen in phase 3 trials, including thrombosis with thrombocytopenia syndrome (TTS), a rare condition resembling heparin-induced thrombocytopenia (HIT), have been observed. This study demonstrates that all three adenoviruses deployed as vaccination vectors versus SARS-CoV-2 bind to platelet factor 4 (PF4), a protein implicated in the pathogenesis of HIT. We have determined the structure of the ChAdOx1 viral vector and used it in state-of-the-art computational simulations to demonstrate an electrostatic interaction mechanism with PF4, which was confirmed experimentally by surface plasmon resonance. These data confirm that PF4 is capable of forming stable complexes with clinically relevant adenoviruses, an important step in unraveling the mechanisms underlying TTS.
Cryoelectron microscopy requires molecular modeling for refinement of structures. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins. Here, we introduce CryoFold, a pipeline of molecular dynamics simulations that determines ensembles of protein structures by integrating density data of varying sparsity at 3–5 Å resolution with sequence information and coarse-grained topological knowledge of the protein folds. We present six examples, folding proteins between 72 and 2,000 residues, including large membrane and multi-domain systems, and results from two Electron Microscopy Data Bank (EMDB) competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python graphical user interface (GUI).