51
|
Suvlu D, Thirumalai D, Rasaiah JC. Water-Mediated Interactions Determine Helix Formation of Peptides in Open Nanotubes. J Phys Chem B 2021; 125:817-824. [PMID: 33464101 DOI: 10.1021/acs.jpcb.0c10178] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Water-mediated interactions (WMIs) play diverse roles in molecular biology. They are particularly relevant in geometrically confined spaces such as the interior of the chaperonin, at the interface between ligands and their binding partners, and in the ribosome tunnel. Inspired in part by the geometry of the ribosome tunnel, we consider confinement effects on the stability of peptides. We describe results from replica exchange molecular dynamics simulations of a system containing a 23-alanine or 23-serine polypeptide confined to nonpolar and polar nanotubes in the gas phase and when open to a water reservoir. We quantify the effect of water in determining the preferred conformational states of these polypeptides by calculating the difference in the solvation free energy for the helix and coil states in the open nanotube in the two phases. Our simulations reveal several possibilities. We find that nanoscopic confinement preferentially stabilizes the helical state of polypeptides with hydrophobic side chains, which is explained by the entropic stabilization mechanism proposed on the basis of polymer physics. Polypeptide chains with hydrophilic side chains can adopt helical structures within nanotubes, but helix formation is sensitive to the nature of the nanotube due to WMIs. We elaborate on the potential implications of our findings to the stability of peptides in the ribosome tunnel.
Collapse
|
52
|
Mugnai ML, Thirumalai D. Step-Wise Hydration of Magnesium by Four Water Molecules Precedes Phosphate Release in a Myosin Motor. J Phys Chem B 2021; 125:1107-1117. [PMID: 33481593 DOI: 10.1021/acs.jpcb.0c10004] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular motors, such as myosin, kinesin, and dynein, convert the energy released by the hydrolysis of ATP into mechanical work, thus allowing them to undergo directional motion on cytoskeletal tracks. A pivotal step in the chemomechanical transduction in myosin motors occurs after they bind to the actin filament, which triggers the release of phosphate (Pi, product of ATP hydrolysis) and the rotation of the lever arm. Here, we investigate the mechanism of phosphate release in myosin VI using extensive molecular dynamics simulations involving multiple trajectories of several μs. Because the escape of phosphate is expected to occur on time-scales on the order of milliseconds or more in myosin VI, we observed Pi release only if the trajectories were initiated with a rotated phosphate inside the nucleotide binding pocket. We discovered that although Pi populates the traditional "back door" route, phosphate exits through various other gateways, thus establishing the heterogeneity in the escape routes. Remarkably, we observed that the release of phosphate is preceded by a stepwise hydration of the ADP-bound magnesium ion. The release of the anion occurred only after four water molecules hydrated the cation (Mg2+). By performing comparative structural analyses, we show that hydration of magnesium is the key step in the phosphate release in a number of ATPases and GTPases. Nature may have evolved hydration of Mg2+ as a general molecular switch for Pi release, which is a universal step in the catalytic cycle of many machines that share little sequence or structural similarity.
Collapse
|
53
|
Koculi E, Thirumalai D. Retardation of Folding Rates of Substrate Proteins in the Nanocage of GroEL. Biochemistry 2021; 60:460-464. [PMID: 33464880 DOI: 10.1021/acs.biochem.0c00903] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Escherichia coli ATP-consuming chaperonin machinery, a complex between GroEL and GroES, has evolved to facilitate folding of substrate proteins (SPs) that cannot do so spontaneously. A series of kinetic experiments show that the SPs are encapsulated in the GroEL/ES nanocage for a short duration. If confinement of the SPs is the mechanism by which GroEL/ES facilitates folding, it follows that the assisted folding rate, relative to the bulk value, should always be enhanced. Here, we show that this is not the case for the folding of rhodanese in the presence of the full machinery of GroEL/ES and ATP. The assisted folding rate of rhodanese decreases. On the basis of our finding and those reported in other studies, we suggest that the ATP-consuming chaperonin machinery has evolved to optimize the product of the folding rate and the yield of the folded SPs on the biological time scale. Neither the rate nor the yield is separately maximized.
Collapse
|
54
|
Gajalakshmi S, Asharani IV, Thirumalai D. Synthesis and Antioxidant Activity of Functionalized Pyridinyl-Methylthiosemicarbazide Derivatives. Polycycl Aromat Compd 2020. [DOI: 10.1080/10406638.2020.1866033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
55
|
Houser JR, Hayden CC, Thirumalai D, Stachowiak JC. A Förster Resonance Energy Transfer-Based Sensor of Steric Pressure on Membrane Surfaces. J Am Chem Soc 2020; 142:20796-20805. [PMID: 33237768 PMCID: PMC8274331 DOI: 10.1021/jacs.0c09802] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cellular membranes are densely covered by proteins. Steric pressure generated by protein collisions plays a significant role in shaping and curving biological membranes. However, no method currently exists for measuring steric pressure at membrane surfaces. Here, we developed a sensor based on Förster resonance energy transfer (FRET), which uses the principles of polymer physics to precisely detect changes in steric pressure. The sensor consists of a polyethylene glycol chain tethered to the membrane surface. The polymer has a donor fluorophore at its free end, such that FRET with acceptor fluorophores in the membrane provides a real-time readout of polymer extension. As a demonstration of the sensor, we measured the steric pressure generated by a model protein involved in membrane bending, the N-terminal homology domain (ENTH) of Epsin1. As the membrane becomes crowded by ENTH proteins, the polymer chain extends, increasing the fluorescence lifetime of the donor. Drawing on polymer theory, we use this change in lifetime to calculate steric pressure as a function of membrane coverage by ENTH, validating theoretical equations of state. Further, we find that ENTH's ability to break up larger vesicles into smaller ones correlates with steric pressure rather than the chemistry used to attach ENTH to the membrane surface. This result addresses a long-standing question about the molecular mechanisms of membrane remodeling. More broadly, this sensor makes it possible to measure steric pressure in situ during diverse biochemical events that occur on membrane surfaces, such as membrane remodeling, ligand-receptor binding, assembly of protein complexes, and changes in membrane organization.
Collapse
|
56
|
Mugnai ML, Templeton C, Elber R, Thirumalai D. Role of Long-range Allosteric Communication in Determining the Stability and Disassembly of SARS-COV-2 in Complex with ACE2. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2020:2020.11.30.405340. [PMID: 33299995 PMCID: PMC7724663 DOI: 10.1101/2020.11.30.405340] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Severe acute respiratory syndrome (SARS) and novel coronavirus disease (COVID-19) are caused by two closely related beta-coronaviruses, SARS-CoV and SARS-CoV-2, respectively. The envelopes surrounding these viruses are decorated with spike proteins, whose receptor binding domains (RBDs) initiate invasion by binding to the human angiotensin-converting enzyme 2 (ACE2). Subtle changes at the interface with ACE2 seem to be responsible for the enhanced affinity for the receptor of the SARS-CoV-2 RBD compared to SARS-CoV RBD. Here, we use Elastic Network Models (ENMs) to study the response of the viral RBDs and ACE2 upon dissassembly of the complexes. We identify a dominant detachment mode, in which the RBD rotates away from the surface of ACE2, while the receptor undergoes a conformational transition which stretches the active-site cleft. Using the Structural Perturbation Method, we determine the network of residues, referred to as the Allostery Wiring Diagram (AWD), which drives the large-scale motion activated by the detachment of the complex. The AWD for SARS-CoV and SARS-CoV-2 are remarkably similar, showing a network that spans the interface of the complex and reaches the active site of ACE2, thus establishing an allosteric connection between RBD binding and receptor catalytic function. Informed in part by the AWD, we used Molecular Dynamics simulations to probe the effect of interfacial mutations in which SARS-CoV-2 residues are replaced by their SARS-CoV counterparts. We focused on a conserved glycine (G502 in SARS-CoV-2, G488 in SARS-CoV) because it belongs to a region that initiates the dissociation of the complex along the dominant detachment mode, and is prominent in the AWD. Molecular Dynamics simulations of SARS-CoV-2 wild-type and G502P mutant show that the affinity for the human receptor of the mutant is drastically diminished. Our results suggest that in addition to residues that are in direct contact with the interface those involved in long range allosteric communication are also a determinant of the stability of the RBD-ACE2 complex.
Collapse
|
57
|
Kirmizialtin S, Pitici F, Cardenas AE, Elber R, Thirumalai D. Dramatic Shape Changes Occur as Cytochrome c Folds. J Phys Chem B 2020; 124:8240-8248. [PMID: 32840372 DOI: 10.1021/acs.jpcb.0c05802] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Extensive experimental studies on the folding of cytochrome c (Cyt c) make this small protein an ideal target for atomic detailed simulations for the purposes of quantitatively characterizing the structural transitions and the associated time scales for folding to the native state from an ensemble of unfolded states. We use previously generated atomically detailed folding trajectories by the stochastic difference equation in length to calculate the time-dependent changes in the small-angle X-ray scattering (SAXS) profiles. Excellent agreement is obtained between experiments and simulations for the time-dependent SAXS spectra, allowing us to identify the structures of the folding intermediates, which shows that Cyt c reaches the native state by a sequential folding mechanism. Using the ensembles of structures along the folding pathways, we show that compaction and the sphericity of Cyt c change dramatically from the prolate ellipsoid shape in the unfolded state to the spherical native state. Our data, which are in unprecedented quantitative agreement with all aspects of time-resolved SAXS experiments, show that hydrophobic collapse and amide group protection coincide on the 100 microseconds time scale, which is in accordance with ultrafast hydrogen/deuterium exchange studies. Based on these results, we propose that compaction of polypeptide chains, accompanied by dramatic shape changes, is a universal characteristic of globular proteins, regardless of the underlying folding mechanism.
Collapse
|
58
|
Sinha S, Malmi-Kakkada AN, Li X, Samanta HS, Thirumalai D. Correction: Spatially heterogeneous dynamics of cells in a growing tumor spheroid: comparison between theory and experiments. SOFT MATTER 2020; 16:6442. [PMID: 32602498 DOI: 10.1039/d0sm90123g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Correction for 'Spatially heterogeneous dynamics of cells in a growing tumor spheroid: comparison between theory and experiments' by Sumit Sinha et al., Soft Matter, 2020, 16, 5294-5304, DOI: .
Collapse
|
59
|
Sinha S, Malmi-Kakkada AN, Li X, Samanta HS, Thirumalai D. Spatially heterogeneous dynamics of cells in a growing tumor spheroid: comparison between theory and experiments. SOFT MATTER 2020; 16:5294-5304. [PMID: 32462163 DOI: 10.1039/c9sm02277e] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Collective cell movement, characterized by multiple cells that are in contact for substantial periods of time and undergo correlated motion, plays a central role in cancer and embryogenesis. Recent imaging experiments have provided time-dependent traces of individual cells, thus providing an unprecedented picture of tumor spheroid growth. By using simulations of a minimal cell model, we analyze the experimental data that map the movement of cells in a fibrosarcoma tumor spheroid embedded in a collagen matrix. Both simulations and experiments show that cells in the core of the spheroid exhibit subdiffusive glassy dynamics (mean square displacement, Δ(t) ≈ tα with α < 1), whereas cells in the periphery exhibit superdiffusive motion, Δ(t) ≈ tα with α > 1. The motion of most of the cells near the periphery is highly persistent and correlated directional motion due to cell doubling and apoptosis rates, thus explaining the observed superdiffusive behavior. The α values for cells in the core and periphery, extracted from simulations and experiments, are in near quantitative agreement with each other, which is surprising given that no parameter in the model was used to fit the measurements. The qualitatively different dynamics of cells in the core and periphery is captured by the fourth order susceptibility, introduced to characterize metastable states in glass forming systems. Analyses of the velocity autocorrelation of individual cells show remarkable spatial heterogeneity with no two cells exhibiting similar behavior. The prediction that α should depend on the location of the cells in the tumor is amenable to experimental testing. The highly heterogeneous dynamics of cells in the tumor spheroid provides a plausible mechanism for the origin of intratumor heterogeneity.
Collapse
|
60
|
Nguyen HT, Thirumalai D. Charge Density of Cation Determines Inner versus Outer Shell Coordination to Phosphate in RNA. J Phys Chem B 2020; 124:4114-4122. [PMID: 32342689 DOI: 10.1021/acs.jpcb.0c02371] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Divalent cations are often required to fold RNA, which is a highly charged polyanion. Condensation of ions, such as Mg2+ or Ca2+, in the vicinity of RNA renormalizes the effective charges on the phosphate groups, thus minimizing the intra RNA electrostatic repulsion. The prevailing view is that divalent ions bind diffusively in a nonspecific manner. In sharp contrast, we arrive at the exact opposite conclusion using a theory for the interaction of ions with the phosphate groups using RISM theory in conjunction with simulations based on an accurate three-interaction-site RNA model. The divalent ions bind in a nucleotide-specific manner using either the inner (partially dehydrated) or outer (fully hydrated) shell coordination. The high charge density Mg2+ ion has a preference to bind to the outer shell, whereas the opposite is the case for Ca2+. Surprisingly, we find that bridging interactions, involving ions that are coordinated to two or more phosphate groups, play a crucial role in maintaining the integrity of the folded state. Their importance could become increasingly prominent as the size of the RNA increases. Because the modes of interaction of divalent ions with DNA are likely to be similar, we propose that specific inner and outer shell coordination could play a role in DNA condensation, and perhaps genome organization as well.
Collapse
|
61
|
Cho HW, Mugnai ML, Kirkpatrick TR, Thirumalai D. Fragile-to-strong crossover, growing length scales, and dynamic heterogeneity in Wigner glasses. Phys Rev E 2020; 101:032605. [PMID: 32290023 DOI: 10.1103/physreve.101.032605] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Accepted: 02/26/2020] [Indexed: 06/11/2023]
Abstract
Colloidal particles, which are ubiquitous, have become ideal testing grounds for the structural glass transition theories. In these systems glassy behavior arises as the density of the particles is increased. Thus, soft colloidal particles with varying degree of softness capture diverse glass-forming properties, observed normally in molecular glasses. Brownian dynamics simulations for a binary mixture of micron-sized charged colloidal suspensions show that tuning the softness of the interaction potential, achievable by changing the monovalent salt concentration results in a continuous transition from fragile to strong behavior. Remarkably, this is found in a system where the well characterized interaction potential between the colloidal particles is isotropic. We also show that the predictions of the random first-order transition (RFOT) theory quantitatively describes the universal features such as the growing correlation length, ξ∼(ϕ_{K}/ϕ-1)^{-ν} with ν=2/3 where ϕ_{K}, the analog of the Kauzmann temperature, depends on the salt concentration. As anticipated by the RFOT predictions, we establish a causal relationship between the growing correlation length and a steep increase in the relaxation time and dynamic heterogeneity as the system is compressed. The broad range of fragility observed in Wigner glasses is used to draw analogies with molecular and polymer glasses. The large variations in the fragility are normally found only when the temperature dependence of the viscosity is examined for a large class of diverse glass-forming materials. In sharp contrast, this is vividly illustrated in a single system that can be experimentally probed. Our work also shows that the RFOT predictions are accurate in describing the dynamics over the entire density range, regardless of the fragility of the glasses.
Collapse
|
62
|
Li X, Thirumalai D. Share, but unequally: a plausible mechanism for emergence and maintenance of intratumour heterogeneity. J R Soc Interface 2020; 16:20180820. [PMID: 30958159 DOI: 10.1098/rsif.2018.0820] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Intratumour heterogeneity (ITH), referring to the coexistence of different cell subpopulations in a single tumour, has been a major puzzle in cancer research for almost half a century. The lack of understanding of the underlying mechanism of ITH hinders progress in developing effective therapies for cancers. Based on the findings in a recent quantitative experiment on pancreatic cancer, we developed a general evolutionary model for one type of cancer, accounting for interactions between different cell populations through paracrine or juxtacrine factors. We show that the emergence of a stable heterogeneous state in a tumour requires an unequal allocation of paracrine growth factors (public goods) between cells that produce them and those that merely consume them. Our model provides a quantitative explanation of recent in vitro experimental studies in pancreatic cancer in which insulin-like growth factor II (IGF-II) plays the role of public goods. The calculated phase diagrams as a function of exogenous resources and fraction of growth factor producing cells show ITH persists only in a narrow range of concentration of exogenous IGF-II. Remarkably, maintenance of ITH requires cooperation among tumour cell subpopulations in harsh conditions, specified by lack of exogenous IGF-II, whereas surplus exogenous IGF-II elicits competition. Our theory also quantitatively accounts for measured in vivo tumour growth in glioblastoma multiforme (GBM). The predictions for GBM tumour growth as a function of the fraction of tumour cells are amenable to experimental tests. The mechanism for ITH also provides hints for devising efficacious therapies.
Collapse
|
63
|
Hathcock D, Tehver R, Hinczewski M, Thirumalai D. Myosin V executes steps of variable length via structurally constrained diffusion. eLife 2020; 9:51569. [PMID: 31939739 PMCID: PMC7054003 DOI: 10.7554/elife.51569] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2019] [Accepted: 01/14/2020] [Indexed: 11/16/2022] Open
Abstract
The molecular motor myosin V transports cargo by stepping on actin filaments, executing a random diffusive search for actin binding sites at each step. A recent experiment suggests that the joint between the myosin lever arms may not rotate freely, as assumed in earlier studies, but instead has a preferred angle giving rise to structurally constrained diffusion. We address this controversy through comprehensive analytical and numerical modeling of myosin V diffusion and stepping. When the joint is constrained, our model reproduces the experimentally observed diffusion, allowing us to estimate bounds on the constraint energy. We also test the consistency between the constrained diffusion model and previous measurements of step size distributions and the load dependence of various observable quantities. The theory lets us address the biological significance of the constrained joint and provides testable predictions of new myosin behaviors, including the stomp distribution and the run length under off-axis force.
Collapse
|
64
|
Thirumalai D, Lorimer GH, Hyeon C. Iterative annealing mechanism explains the functions of the GroEL and RNA chaperones. Protein Sci 2019; 29:360-377. [PMID: 31800116 DOI: 10.1002/pro.3795] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2019] [Revised: 11/21/2019] [Accepted: 11/22/2019] [Indexed: 12/16/2022]
Abstract
Molecular chaperones are ATP-consuming machines, which facilitate the folding of proteins and RNA molecules that are kinetically trapped in misfolded states. Unassisted folding occurs by the kinetic partitioning mechanism according to which folding to the native state, with low probability as well as misfolding to one of the many metastable states, with high probability, occur rapidly. GroEL is an all-purpose stochastic machine that assists misfolded substrate proteins to fold. The RNA chaperones such as CYT-19, which are ATP-consuming enzymes, help the folding of ribozymes that get trapped in metastable states for long times. GroEL does not interact with the folded proteins but CYT-19 disrupts both the folded and misfolded ribozymes. The structures of GroEL and RNA chaperones are strikingly different. Despite these differences, the iterative annealing mechanism (IAM) quantitatively explains all the available experimental data for assisted folding of proteins and ribozymes. Driven by ATP binding and hydrolysis and GroES binding, GroEL undergoes a catalytic cycle during which it samples three allosteric states, T (apo), R (ATP bound), and R″ (ADP bound). Analyses of the experimental data show that the efficiency of the GroEL-GroES machinery and mutants is determined by the resetting rate k R ″ → T , which is largest for the wild-type (WT) GroEL. Generalized IAM accurately predicts the folding kinetics of Tetrahymena ribozyme and its variants. Chaperones maximize the product of the folding rate and the steady-state native state fold by driving the substrates out of equilibrium. Neither the absolute yield nor the folding rate is optimized.
Collapse
|
65
|
Chakrabarti S, Jarzynski C, Thirumalai D. Processivity, Velocity, and Universal Characteristics of Nucleic Acid Unwinding by Helicases. Biophys J 2019; 117:867-879. [PMID: 31400912 PMCID: PMC6731385 DOI: 10.1016/j.bpj.2019.07.021] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2019] [Revised: 06/20/2019] [Accepted: 07/12/2019] [Indexed: 12/20/2022] Open
Abstract
Helicases are components of the cellular replisome that are essential for unwinding double-strand nucleic acids during the process of replication. Intriguingly, most helicases are inefficient and require either oligomerization or assistance from other partner proteins to increase the processivity of unwinding in the presence of the replication fork, which acts as a barrier to progress. Single-molecule force spectroscopy has emerged as a promising experimental technique to probe how relieving this barrier on the helicase can allow for increased efficiency of unwinding. However, there exists no comprehensive theoretical framework to provide unique interpretations of the underlying helicase kinetics from the force spectroscopy data. This remains a major confounding issue in the field. Here, we develop a mathematical framework and derive analytic expressions for the velocity and run length of a general model of finitely processive helicases, the two most commonly measured experimental quantities. We show that in contrast to the unwinding velocity, the processivity exhibits a universal increase in response to external force, irrespective of the underlying architecture and unwinding kinetics of the helicase. Our work provides the first, to our knowledge, explanation to a wide array of experiments and suggests that helicases may have evolved to maximize processivity rather than speed. To demonstrate the use of our theory on experimental data, we analyze velocity and processivity data on the T7 helicase and provide unique inferences on the kinetics of the helicase. Our results show that T7 is a weakly active helicase that destabilizes the fork ahead by less than 1 kBT and back steps very frequently while unwinding DNA. Our work generates fundamental insights into the force response of helicases and provides a widely applicable method for inferring the underlying helicase kinetics from force spectroscopy data.
Collapse
|
66
|
Shi G, Thirumalai D. Conformational heterogeneity in human interphase chromosome organization reconciles the FISH and Hi-C paradox. Nat Commun 2019; 10:3894. [PMID: 31467267 PMCID: PMC6715811 DOI: 10.1038/s41467-019-11897-0] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2018] [Accepted: 08/06/2019] [Indexed: 11/19/2022] Open
Abstract
Hi-C experiments are used to infer the contact probabilities between loci separated by varying genome lengths. Contact probability should decrease as the spatial distance between two loci increases. However, studies comparing Hi-C and FISH data show that in some cases the distance between one pair of loci, with larger Hi-C readout, is paradoxically larger compared to another pair with a smaller value of the contact probability. Here, we show that the FISH-Hi-C paradox can be resolved using a theory based on a Generalized Rouse Model for Chromosomes (GRMC). The FISH-Hi-C paradox arises because the cell population is highly heterogeneous, which means that a given contact is present in only a fraction of cells. Insights from the GRMC is used to construct a theory, without any adjustable parameters, to extract the distribution of subpopulations from the FISH data, which quantitatively reproduces the Hi-C data. Our results show that heterogeneity is pervasive in genome organization at all length scales, reflecting large cell-to-cell variations. Studies comparing Hi-C and FISH data show that in some cases the distance between one pair of loci is paradoxically larger compared to another pair with a smaller value of the contact probability. Here the authors use a theory based on a Generalized Rouse Model for Chromosomes to resolve this paradox.
Collapse
|
67
|
Thirumalai D, Samanta HS, Maity H, Reddy G. Universal Nature of Collapsibility in the Context of Protein Folding and Evolution. Trends Biochem Sci 2019; 44:675-687. [DOI: 10.1016/j.tibs.2019.04.003] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2018] [Revised: 04/03/2019] [Accepted: 04/04/2019] [Indexed: 12/21/2022]
|
68
|
Samanta HS, Mugnai ML, Kirkpatrick TR, Thirumalai D. Giant Casimir Nonequilibrium Forces Drive Coil to Globule Transition in Polymers. J Phys Chem Lett 2019; 10:2788-2793. [PMID: 31066561 DOI: 10.1021/acs.jpclett.9b00695] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We develop a theory to probe the effect of nonequilibrium fluctuation-induced forces on the size of a polymer confined between two horizontal, thermally conductive plates subject to a constant temperature gradient, ∇ T. We assume that (a) the solvent is good and (b) the distance between the plates is large so that in the absence of a thermal gradient the polymer is a coil, whose size scales with the number of monomers as Nν, with ν ≈ 0.6. We find that above a critical temperature gradient, ∇ Tc ≈ N-5/4, a favorable attractive monomer-monomer interaction due to the giant Casimir force (GCF) overcomes the chain conformational entropy, resulting in a coil-globule transition. Our predictions can be verified using light-scattering experiments with polymers, such as polystyrene or polyisoprene in organic solvents in which the GCF is attractive.
Collapse
|
69
|
Hori N, Denesyuk NA, Thirumalai D. Ion Condensation onto Ribozyme Is Site Specific and Fold Dependent. Biophys J 2019; 116:2400-2410. [PMID: 31130233 DOI: 10.1016/j.bpj.2019.04.037] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2019] [Revised: 04/12/2019] [Accepted: 04/23/2019] [Indexed: 01/20/2023] Open
Abstract
The highly charged RNA molecules, with each phosphate carrying a single negative charge, cannot fold into well-defined architectures with tertiary interactions in the absence of ions. For ribozymes, divalent cations are known to be more efficient than monovalent ions in driving them to a compact state, although Mg2+ ions are needed for catalytic activities. Therefore, how ions interact with RNA is relevant in understanding RNA folding. It is often thought that most of the ions are territorially and nonspecifically bound to the RNA, as predicted by the counterion condensation theory. Here, we show using simulations of Azoarcus ribozyme, based on an accurate coarse-grained three-site interaction model with explicit divalent and monovalent cations, that ion condensation is highly specific and depends on the nucleotide position. The regions with high coordination between the phosphate groups and the divalent cations are discernible even at very low Mg2+ concentrations when the ribozyme does not form tertiary interactions. Surprisingly, these regions also contain the secondary structural elements that nucleate subsequently in the self-assembly of RNA, implying that ion condensation is determined by the architecture of the folded state. These results are in sharp contrast to interactions of ions (monovalent and divalent) with rigid charged rods, in which ion condensation is uniform and position independent. The differences are explained in terms of the dramatic nonmonotonic shape fluctuations in the ribozyme as it folds with increasing Mg2+ or Ca2+ concentration.
Collapse
|
70
|
Baul U, Chakraborty D, Mugnai ML, Straub JE, Thirumalai D. Sequence Effects on Size, Shape, and Structural Heterogeneity in Intrinsically Disordered Proteins. J Phys Chem B 2019; 123:3462-3474. [PMID: 30913885 PMCID: PMC6920032 DOI: 10.1021/acs.jpcb.9b02575] [Citation(s) in RCA: 100] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Intrinsically disordered proteins (IDPs) lack well-defined three-dimensional structures, thus challenging the archetypal notion of structure-function relationships. Determining the ensemble of conformations that IDPs explore under physiological conditions is the first step toward understanding their diverse cellular functions. Here, we quantitatively characterize the structural features of IDPs as a function of sequence and length using coarse-grained simulations. For diverse IDP sequences, with the number of residues ( NT) ranging from 20 to 441, our simulations not only reproduce the radii of gyration ( Rg) obtained from experiments, but also predict the full scattering intensity profiles in excellent agreement with small-angle X-ray scattering experiments. The Rg values are well-described by the standard Flory scaling law, Rg = Rg0 NTν, with ν ≈ 0.588, making it tempting to assert that IDPs behave as polymers in a good solvent. However, clustering analysis reveals that the menagerie of structures explored by IDPs is diverse, with the extent of heterogeneity being highly sequence-dependent, even though ensemble-averaged properties, such as the dependence of Rg on chain length, may suggest synthetic polymer-like behavior in a good solvent. For example, we show that for the highly charged Prothymosin-α, a substantial fraction of conformations is highly compact. Even if the sequence compositions are similar, as is the case for α-Synuclein and a truncated construct from the Tau protein, there are substantial differences in the conformational heterogeneity. Taken together, these observations imply that metrics based on net charge or related quantities alone cannot be used to anticipate the phases of IDPs, either in isolation or in complex with partner IDPs or RNA. Our work sets the stage for probing the interactions of IDPs with each other, with folded protein domains, or with partner RNAs, which are critical for describing the structures of stress granules and biomolecular condensates with important cellular functions.
Collapse
|
71
|
Thirumalai D, Hyeon C, Zhuravlev PI, Lorimer GH. Symmetry, Rigidity, and Allosteric Signaling: From Monomeric Proteins to Molecular Machines. Chem Rev 2019; 119:6788-6821. [DOI: 10.1021/acs.chemrev.8b00760] [Citation(s) in RCA: 51] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
72
|
Thirumalai D, Hyeon C. Signalling networks and dynamics of allosteric transitions in bacterial chaperonin GroEL: implications for iterative annealing of misfolded proteins. Philos Trans R Soc Lond B Biol Sci 2019; 373:rstb.2017.0182. [PMID: 29735736 DOI: 10.1098/rstb.2017.0182] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/01/2018] [Indexed: 12/14/2022] Open
Abstract
Signal transmission at the molecular level in many biological complexes occurs through allosteric transitions. Allostery describes the responses of a complex to binding of ligands at sites that are spatially well separated from the binding region. We describe the structural perturbation method, based on phonon propagation in solids, which can be used to determine the signal-transmitting allostery wiring diagram (AWD) in large but finite-sized biological complexes. Application to the bacterial chaperonin GroEL-GroES complex shows that the AWD determined from structures also drives the allosteric transitions dynamically. From both a structural and dynamical perspective these transitions are largely determined by formation and rupture of salt-bridges. The molecular description of allostery in GroEL provides insights into its function, which is quantitatively described by the iterative annealing mechanism. Remarkably, in this complex molecular machine, a deep connection is established between the structures, reaction cycle during which GroEL undergoes a sequence of allosteric transitions, and function, in a self-consistent manner.This article is part of a discussion meeting issue 'Allostery and molecular machines'.
Collapse
|
73
|
Samanta HS, Thirumalai D. Origin of superdiffusive behavior in a class of nonequilibrium systems. Phys Rev E 2019; 99:032401. [PMID: 30999548 DOI: 10.1103/physreve.99.032401] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Indexed: 06/09/2023]
Abstract
Experiments and simulations have established that dynamics in a class of living and abiotic systems that are far from equilibrium exhibit superdiffusive behavior at long times, which in some cases (for example, an evolving tumor) is preceded by slow glass-like dynamics. By using the evolution of a collection of tumor cells, driven by mechanical forces and subject to cell birth and apoptosis, as a case study we show theoretically that on short timescales the mean-square displacement is subdiffusive due to jamming, whereas at long times it is superdiffusive. The results obtained by using a stochastic quantization method, which is needed because of the absence of the fluctuation-dissipation theorem, show that the superdiffusive behavior is universal and impervious to the nature of cell-cell interactions. Surprisingly, the theory also quantitatively accounts for the nontrivial dynamics observed in simulations of a model soap foam characterized by creation and destruction of spherical bubbles, which suggests that the two nonequilibrium systems belong to the same universality class. The theoretical prediction for the superdiffusion exponent is in excellent agreement with simulations for collective motion of tumor cells and dynamics associated with soap bubbles.
Collapse
|
74
|
Mondal T, Dutta S, De S, Thirumalai D, Koley D. Donor Stabilized Diatomic Gr.14 E 2 (E = C-Pb) Molecule D-E 2-D (D = NHC, aNHC, NNHC, NHSi, NHGe, cAAC, cAASi, cAAGe): A Theoretical Insight. J Phys Chem A 2019; 123:565-581. [PMID: 30562460 DOI: 10.1021/acs.jpca.8b11005] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Quantum chemical calculations have been carried out to explore the detailed electronic structure and bonding scenario in various bis-donor stabilized E2 compounds (E = C-Pb). Our computational findings reveal that the thermodynamic stabilities of the E2 core gradually decrease as we move down the group. A linear D-E-E'-D framework is observed for C2 systems, while the heavier group 14 analogues possess trans-bent geometries. Consideration of few compounds as viable targets for synthesis is suggested by their corresponding calculated formation energies. In addition, the thermodynamic stabilities of C2 systems notably increase with the saturation of the donor ring framework and are even more pronounced for boron-substituted saturated NHD ligand. QTAIM calculations affirmed that the covalent nature of E-E' bonds shifts toward the donor-acceptor region as one traverses from top to bottom along group 14. The E-D and E'-D bonds in the C2 systems have covalent nature, whereas those in Si2-Pb2 systems are characterized by donor-acceptor bonds. In addition, we have computed proton affinities and vertical ionization potentials (VIPs) of these compounds. An excellent correlation was obtained between calculated VIPs and orbital energies of HOMOs. Furthermore, in the present study, we also explored the effect of bis-donors in the stabilization of heterodiatomic SiC compounds. Our calculations indicate that a typical bonding description of the SiC(D)2 compounds should be represented by a combination of a classical double bond between C-D with significant donor-acceptor interaction in Si-D, i.e., D → Si═C═D. The SiC(D)2 systems are found to be less stable than the corresponding dicarbon compounds C2(D)2, but they show significant stabilization compared to the corresponding disilicon systems Si2(D)2.
Collapse
|
75
|
Denesyuk NA, Hori N, Thirumalai D. Molecular Simulations of Ion Effects on the Thermodynamics of RNA Folding. J Phys Chem B 2018; 122:11860-11867. [PMID: 30468380 DOI: 10.1021/acs.jpcb.8b08142] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
How ions affect RNA folding thermodynamics and kinetics is an important but a vexing problem that remains unsolved. Experiments have shown that the free-energy change, Δ G( c), of RNA upon folding varies with the salt concentration ( c) as, Δ G( c) = k c ln c + const, where the coefficient k c is proportional to the difference in the ion preferential coefficient, ΔΓ. We performed simulations of a coarse-grained model, by modeling electrostatic interactions implicitly and with explicit representation of ions, to elucidate the molecular underpinnings of the relationship between Δ G and ΔΓ. The simulations quantitatively reproduce the heat capacity for a pseudoknot, thus validating the model. We show that Δ G( c), calculated directly from ΔΓ, varies linearly with ln c ( c < 0.2 M), for a hairpin and the pseudoknot, demonstrating a molecular link between the two quantities. Explicit ion simulations also show the linear dependence of Δ G( c) on ln c at all c with k c = 2 kB T, except that Δ G( c) values are shifted by ∼2 kcal/mol higher than experiments. The discrepancy is due to an underestimation of Γ for both the folded and unfolded states while giving accurate values for ΔΓ. The predictions for the salt dependence of ΔΓ are amenable to test using single-molecule pulling experiments. The framework provided here can be used to obtain accurate thermodynamics for other RNA molecules as well.
Collapse
|
76
|
Toan NM, Thirumalai D. Forced-rupture of cell-adhesion complexes reveals abrupt switch between two brittle states. J Chem Phys 2018; 148:123332. [PMID: 29604893 DOI: 10.1063/1.5011056] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Cell adhesion complexes (CACs), which are activated by ligand binding, play key roles in many cellular functions ranging from cell cycle regulation to mediation of cell extracellular matrix adhesion. Inspired by single molecule pulling experiments using atomic force spectroscopy on leukocyte function-associated antigen-1 (LFA-1), expressed in T-cells, bound to intercellular adhesion molecules (ICAM), we performed constant loading rate (rf) and constant force (F) simulations using the self-organized polymer model to describe the mechanism of ligand rupture from CACs. The simulations reproduce the major experimental finding on the kinetics of the rupture process, namely, the dependence of the most probable rupture forces (f*s) on ln rf (rf is the loading rate) exhibits two distinct linear regimes. The first, at low rf, has a shallow slope, whereas the slope at high rf is much larger, especially for a LFA-1/ICAM-1 complex with the transition between the two occurring over a narrow rf range. Locations of the two transition states (TSs) extracted from the simulations show an abrupt change from a high value at low rf or constant force, F, to a low value at high rf or F. This unusual behavior in which the CACs switch from one brittle (TS position is a constant over a range of forces) state to another brittle state is not found in forced-rupture in other protein complexes. We explain this novel behavior by constructing the free energy profiles, F(Λ)s, as a function of a collective reaction coordinate (Λ), involving many key charged residues and a critical metal ion (Mg2+). The TS positions in F(Λ), which quantitatively agree with the parameters extracted using the Bell-Evans model, change abruptly at a critical force, demonstrating that it, rather than the molecular extension, is a good reaction coordinate. Our combined analyses using simulations performed in both the pulling modes (constant rf and F) reveal a new mechanism for the two loading regimes observed in the rupture kinetics in CACs.
Collapse
|
77
|
Samanta HS, Chakraborty D, Thirumalai D. Charge fluctuation effects on the shape of flexible polyampholytes with applications to intrinsically disordered proteins. J Chem Phys 2018; 149:163323. [PMID: 30384718 DOI: 10.1063/1.5035428] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Random polyampholytes (PAs) contain positively and negatively charged monomers that are distributed randomly along the polymer chain. The interaction between charges is assumed to be given by the Debye-Huckel potential. We show that the size of the PA is determined by an interplay between electrostatic interactions, giving rise to the polyelectrolyte effect due to net charge per monomer (σ) and an effective attractive PA interaction due to charge fluctuations, δσ. The interplay between these terms gives rise to non-monotonic dependence of the radius of gyration, R g , on the inverse Debye length, κ, when PA effects are important ( δ σ σ > 1 ). In the opposite limit, R g decreases monotonically with increasing κ. Simulations of PA chains, using a charged bead-spring model, further corroborate our theoretical predictions. The simulations unambiguously show that conformational heterogeneity manifests itself among sequences that have identical PA parameters. A clear implication is that the phases of PA sequences, and by inference intrinsically disordered proteins (IDPs), cannot be determined using only the bare PA parameters (σ and δσ). The theory is used to calculate the changes in R g on N, the number of residues for a set of IDPs. For a certain class of IDPs, with N between 24 and 441, the size grows as R g ∼ N 0.6, which agrees with data from small angle X-ray scattering experiments.
Collapse
|
78
|
Goldtzvik Y, Mugnai ML, Thirumalai D. Dynamics of Allosteric Transitions in Dynein. Structure 2018; 26:1664-1677.e5. [PMID: 30270176 DOI: 10.1016/j.str.2018.08.005] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2018] [Revised: 07/19/2018] [Accepted: 08/08/2018] [Indexed: 12/20/2022]
Abstract
Cytoplasmic dynein, whose motor domain belongs to the AAA+ family, walks on microtubules toward the minus end. Using the available structures in different nucleotide states, we performed simulations of a coarse-grained model to elucidate the dynamics of allosteric transitions. Binding of ATP closes the cleft between the AAA1 and AAA2 domains, triggering conformational changes in the rest of the motor domain, thus forming the pre-power stroke state. Interactions with the microtubule, modeled implicitly, enhance ADP release rate, and the formation of the post-power stroke state. The dynamics of the linker (LN), which reversibly changes from a straight to a bent state, is heterogeneous. Persistent interactions between the LN and the insert loops in the AAA2 domain prevent the formation of pre-power stroke state when ATP is bound to AAA3, thus locking dynein in a repressed non-functional state. Application of mechanical force to the LN restores motility in the repressed state.
Collapse
|
79
|
Hori N, Denesyuk NA, Thirumalai D. Frictional Effects on RNA Folding: Speed Limit and Kramers Turnover. J Phys Chem B 2018; 122:11279-11288. [PMID: 30179471 DOI: 10.1021/acs.jpcb.8b07129] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
We investigated frictional effects on the folding rates of a human telomerase hairpin (hTR HP) and H-type pseudoknot from the Beet Western Yellow Virus (BWYV PK) using simulations of the Three Interaction Site (TIS) model for RNA. The heat capacity from TIS model simulations, calculated using temperature replica exchange simulations, reproduces nearly quantitatively the available experimental data for the hTR HP. The corresponding results for BWYV PK serve as predictions. We calculated the folding rates ( kF) from more than 100 folding trajectories for each value of the solvent viscosity (η) at a fixed salt concentration of 200 mM. By using the theoretical estimate (∝ √N where N is the number of nucleotides) for folding free energy barrier, kF data for both the RNAs are quantitatively fit using one-dimensional Kramers's theory with two parameters specifying the curvatures in the unfolded basin and the barrier top. In the high-friction regime (η ≳ 10-5 Pa·s), for both HP and PK, kF values decrease as 1/η, whereas in the low friction regime, kF values increase as η increases, leading to a maximum folding rate at a moderate viscosity (∼10-6 Pa·s), which is the Kramers turnover. From the fits, we find that the speed limit to RNA folding at water viscosity is between 1 and 4 μs, which is in accord with our previous theoretical prediction as well as results from several single molecule experiments. Both the RNA constructs fold by parallel pathways. Surprisingly, we find that the flux through the pathways could be altered by changing solvent viscosity, a prediction that is more easily testable in RNA than in proteins.
Collapse
|
80
|
Chen J, Thirumalai D. Interface Residues That Drive Allosteric Transitions Also Control the Assembly of l-Lactate Dehydrogenase. J Phys Chem B 2018; 122:11195-11205. [PMID: 30102042 DOI: 10.1021/acs.jpcb.8b06430] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The allosteric enzyme, l-lactate dehydrogenase (LDH), is activated by fructose 1,6-metaphosphate (FBP) to reduce pyruvate to lactate. The molecular details of the FBP-driven transition from the low affinity T state to the high affinity R state in LDH, a tetramer composed of identical subunits, are not known. The dynamics of the T → R allosteric transition, investigated using Brownian dynamics (BD) simulations of the self-organized polymer (SOP) model, revealed that coordinated rotations of the subunits drive the T → R transition. We used the structural perturbation method (SPM), which requires only the static structure, to identify the allostery wiring diagram (AWD), a network of residues that transmits signals across the tetramer, as LDH undergoes the T → R transition. Interestingly, the residues that play a major role in the dynamics, which are predominantly localized at the interfaces, coincide with the AWD identified using the SPM. Although the allosteric pathways are highly heterogeneous, on the basis of our simulations, we surmise that predominantly the conformational changes in the T → R transition start from the region near the active site, comprised of helix αC, helix α1/2G, helix α3G, and helix α2F, and proceed to other structural units, thus completing the global motion. Brownian dynamics simulations of the tetramer assembly, triggered by a temperature quench from the fully disrupted conformations, show that the bottleneck for assembly is the formation of the correct orientational registry between the subunits, requiring contacts between the interface residues. Surprisingly, these residues are part of the AWD, which was identified using the SPM. Taken together, our results show that LDH, and perhaps other multidomain proteins, may have evolved to stabilize distinct states of allosteric enzymes using precisely the same AWD that also controls the functionally relevant allosteric transitions.
Collapse
|
81
|
Shi G, Liu L, Hyeon C, Thirumalai D. Interphase human chromosome exhibits out of equilibrium glassy dynamics. Nat Commun 2018; 9:3161. [PMID: 30089831 PMCID: PMC6082855 DOI: 10.1038/s41467-018-05606-6] [Citation(s) in RCA: 81] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2017] [Accepted: 07/05/2018] [Indexed: 01/01/2023] Open
Abstract
Fingerprints of the three-dimensional organization of genomes have emerged using advances in Hi-C and imaging techniques. However, genome dynamics is poorly understood. Here, we create the chromosome copolymer model (CCM) by representing chromosomes as a copolymer with two epigenetic loci types corresponding to euchromatin and heterochromatin. Using novel clustering techniques, we establish quantitatively that the simulated contact maps and topologically associating domains (TADs) for chromosomes 5 and 10 and those inferred from Hi-C experiments are in good agreement. Chromatin exhibits glassy dynamics with coherent motion on micron scale. The broad distribution of the diffusion exponents of the individual loci, which quantitatively agrees with experiments, is suggestive of highly heterogeneous dynamics. This is reflected in the cell-to-cell variations in the contact maps. Chromosome organization is hierarchical, involving the formation of chromosome droplets (CDs) on genomic scale, coinciding with the TAD size, followed by coalescence of the CDs, reminiscent of Ostwald ripening.
Collapse
|
82
|
Chakraborty D, Hori N, Thirumalai D. Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA. J Chem Theory Comput 2018; 14:3763-3779. [PMID: 29870236 PMCID: PMC6423546 DOI: 10.1021/acs.jctc.8b00091] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.
Collapse
|
83
|
Pantelopulos GA, Straub JE, Thirumalai D, Sugita Y. Structure of APP-C99 1-99 and implications for role of extra-membrane domains in function and oligomerization. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2018; 1860:1698-1708. [PMID: 29702072 DOI: 10.1016/j.bbamem.2018.04.002] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2018] [Revised: 04/07/2018] [Accepted: 04/09/2018] [Indexed: 01/30/2023]
Abstract
The 99 amino acid C-terminal fragment of Amyloid Precursor Protein APP-C99 (C99) is cleaved by γ-secretase to form Aβ peptide, which plays a critical role in the etiology of Alzheimer's Disease (AD). The structure of C99 consists of a single transmembrane domain flanked by intra and intercellular domains. While the structure of the transmembrane domain has been well characterized, little is known about the structure of the flanking domains and their role in C99 processing by γ-secretase. To gain insight into the structure of full-length C99, REMD simulations were performed for monomeric C99 in model membranes of varying thickness. We find equilibrium ensembles of C99 from simulation agree with experimentally-inferred residue insertion depths and protein backbone chemical shifts. In thin membranes, the transmembrane domain structure is correlated with extra-membrane structural states and the extra-membrane domain structural states become less correlated to each other. Mean and variance of the transmembrane and G37G38 hinge angles are found to increase with thinning membrane. The N-terminus of C99 forms β-strands that may seed aggregation of Aβ on the membrane surface, promoting amyloid formation. In thicker membranes the N-terminus forms α-helices that interact with the nicastrin domain of γ-secretase. The C-terminus of C99 becomes more α-helical as the membrane thickens, forming structures that may be suitable for binding by cytoplasmic proteins, while C-terminal residues essential to cytotoxic function become α-helical as the membrane thins. The heterogeneous but discrete extra-membrane domain states analyzed here open the path to new investigations of the role of C99 structure and membrane in amyloidogenesis. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.
Collapse
|
84
|
Dutta S, Maity B, Thirumalai D, Koley D. Computational Investigation of Carbene–Phosphinidenes: Correlation between 31P Chemical Shifts and Bonding Features to Estimate the π-Backdonation of Carbenes. Inorg Chem 2018. [DOI: 10.1021/acs.inorgchem.8b00174] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
85
|
Liu Z, Thirumalai D. Denaturants Alter the Flux through Multiple Pathways in the Folding of PDZ Domain. J Phys Chem B 2018; 122:1408-1416. [PMID: 29303586 DOI: 10.1021/acs.jpcb.7b11408] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Although we understand many aspects of how small proteins (number of residues less than about hundred) fold, it is a major challenge to quantitatively describe how large proteins self-assemble. To partially overcome this challenge, we performed simulations using the self-organized polymer model with side chains (SOP-SC) in guanidinium chloride (GdmCl), using the molecular transfer model (MTM), to describe the folding of the 110-residue PDZ3 domain. The simulations reproduce the folding thermodynamics accurately including the melting temperature (Tm), the stability of the folded state with respect to the unfolded state. We show that the calculated dependence of ln kobs (kobs is the relaxation rate) has the characteristic chevron shape. The slopes of the chevron plots are in good agreement with experiments. We show that PDZ3 folds by four major pathways populating two metastable intermediates, in accord with the kinetic partitioning mechanism. The structure of one of the intermediates, populated after polypeptide chain collapse, is structurally similar to an equilibrium intermediate. Surprisingly, the connectivities between the intermediates and hence, the fluxes through the pathways depend on the concentration of GdmCl. The results are used to predict possible outcomes for unfolding of PDZ domain subject to mechanical forces. Our study demonstrates that, irrespective of the size or topology, simulations based on MTM and SOP-SC offer a theoretical framework for describing the folding of proteins, mimicking precisely the conditions used in experiments.
Collapse
|
86
|
Asharani IV, Thirumalai D, Sivakumar A. Dendrimer encapsulated Silver nanoparticles as novel catalysts for reduction of aromatic nitro compounds. ACTA ACUST UNITED AC 2017. [DOI: 10.1088/1757-899x/263/2/022010] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
87
|
Lee Y, Thirumalai D, Hyeon C. Ultrasensitivity of Water Exchange Kinetics to the Size of Metal Ion. J Am Chem Soc 2017; 139:12334-12337. [PMID: 28853881 DOI: 10.1021/jacs.7b04198] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Metal ions play a vital role in many biological processes. An important factor in these processes is the dynamics of exchange between ion bound-water molecules and the bulk. Although structural and dynamical properties of labile waters bound to metal ions, such as Na+ and Ca2+, can be elucidated using molecular dynamics simulations, direct evaluation of rates of exchange of waters rigidly bound to high charge density Mg2+, has been elusive. Here, we report a universal relationship, allowing us to determine the water exchange time on metal ions as a function of valence and hydration radius. The proposed relationship, which covers times spanning 14 orders of magnitude, highlights the ultrasensitivity of water lifetime to the ion size, as exemplified by divalent ions, Ca2+ (∼100 ps) and Mg2+ (∼1.5 μs). We show that even when structures, characterized by radial distributions are similar, a small difference in hydration radius leads to a qualitatively different (associative or dissociative) mechanism of water exchange. Our work provides a theoretical basis for determination of hydration radius, which is critical for accurately modeling the water dynamics around multivalent ions, and hence in describing all electrostatically driven events such as ribozyme folding and catalysis.
Collapse
|
88
|
Samanta HS, Hinczewski M, Thirumalai D. Optimal information transfer in enzymatic networks: A field theoretic formulation. Phys Rev E 2017; 96:012406. [PMID: 29347079 DOI: 10.1103/physreve.96.012406] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2017] [Indexed: 06/07/2023]
Abstract
Signaling in enzymatic networks is typically triggered by environmental fluctuations, resulting in a series of stochastic chemical reactions, leading to corruption of the signal by noise. For example, information flow is initiated by binding of extracellular ligands to receptors, which is transmitted through a cascade involving kinase-phosphatase stochastic chemical reactions. For a class of such networks, we develop a general field-theoretic approach to calculate the error in signal transmission as a function of an appropriate control variable. Application of the theory to a simple push-pull network, a module in the kinase-phosphatase cascade, recovers the exact results for error in signal transmission previously obtained using umbral calculus [Hinczewski and Thirumalai, Phys. Rev. X 4, 041017 (2014)2160-330810.1103/PhysRevX.4.041017]. We illustrate the generality of the theory by studying the minimal errors in noise reduction in a reaction cascade with two connected push-pull modules. Such a cascade behaves as an effective three-species network with a pseudointermediate. In this case, optimal information transfer, resulting in the smallest square of the error between the input and output, occurs with a time delay, which is given by the inverse of the decay rate of the pseudointermediate. Surprisingly, in these examples the minimum error computed using simulations that take nonlinearities and discrete nature of molecules into account coincides with the predictions of a linear theory. In contrast, there are substantial deviations between simulations and predictions of the linear theory in error in signal propagation in an enzymatic push-pull network for a certain range of parameters. Inclusion of second-order perturbative corrections shows that differences between simulations and theoretical predictions are minimized. Our study establishes that a field theoretic formulation of stochastic biological signaling offers a systematic way to understand error propagation in networks of arbitrary complexity.
Collapse
|
89
|
Abstract
Myosin VI (MVI) is the only known member of the myosin superfamily that, upon dimerization, walks processively toward the pointed end of the actin filament. The leading head of the dimer directs the trailing head forward with a power stroke, a conformational change of the motor domain exaggerated by the lever arm. Using a unique coarse-grained model for the power stroke of a single MVI, we provide the molecular basis for its motility. We show that the power stroke occurs in two major steps. First, the motor domain attains the poststroke conformation without directing the lever arm forward; and second, the lever arm reaches the poststroke orientation by undergoing a rotational diffusion. From the analysis of the trajectories, we discover that the potential that directs the rotating lever arm toward the poststroke conformation is almost flat, implying that the lever arm rotation is mostly uncoupled from the motor domain. Because a backward load comparable to the largest interhead tension in a MVI dimer prevents the rotation of the lever arm, our model suggests that the leading-head lever arm of a MVI dimer is uncoupled, in accord with the inference drawn from polarized total internal reflection fluorescence (polTIRF) experiments. Without any adjustable parameter, our simulations lead to quantitative agreement with polTIRF experiments, which validates the structural insights. Finally, in addition to making testable predictions, we also discuss the implications of our model in explaining the broad step-size distribution of the MVI stepping pattern.
Collapse
|
90
|
Samanta HS, Zhuravlev PI, Hinczewski M, Hori N, Chakrabarti S, Thirumalai D. Protein collapse is encoded in the folded state architecture. SOFT MATTER 2017; 13:3622-3638. [PMID: 28447708 DOI: 10.1039/c7sm00074j] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Folded states of single domain globular proteins are compact with high packing density. The radius of gyration, Rg, of both the folded and unfolded states increase as Nν where N is the number of amino acids in the protein. The values of the Flory exponent ν are, respectively, ≈⅓ and ≈0.6 in the folded and unfolded states, coinciding with those for homopolymers. However, the extent of compaction of the unfolded state of a protein under low denaturant concentration (collapsibility), conditions favoring the formation of the folded state, is unknown. We develop a theory that uses the contact map of proteins as input to quantitatively assess collapsibility of proteins. Although collapsibility is universal, the propensity to be compact depends on the protein architecture. Application of the theory to over two thousand proteins shows that collapsibility depends not only on N but also on the contact map reflecting the native structure. A major prediction of the theory is that β-sheet proteins are far more collapsible than structures dominated by α-helices. The theory and the accompanying simulations, validating the theoretical predictions, provide insights into the differing conclusions reached using different experimental probes assessing the extent of compaction of proteins. By calculating the criterion for collapsibility as a function of protein length we provide quantitative insights into the reasons why single domain proteins are small and the physical reasons for the origin of multi-domain proteins. Collapsibility of non-coding RNA molecules is similar β-sheet proteins structures adding support to "Compactness Selection Hypothesis".
Collapse
|
91
|
Thirumalai D. Theory and computations in biology: Kamal's legacy. Phys Biol 2017; 14:010401. [PMID: 28177933 DOI: 10.1088/1478-3975/aa5564] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
92
|
Thirumalai D, Shi G. Chromatin Is Stretched but Intact When the Nucleus Is Squeezed through Constrictions. Biophys J 2017; 112:411-412. [PMID: 27989542 DOI: 10.1016/j.bpj.2016.11.902] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2016] [Revised: 11/19/2016] [Accepted: 11/21/2016] [Indexed: 11/18/2022] Open
|
93
|
Hori N, Denesyuk NA, Thirumalai D. Multistep Folding Kinetics of Group I Intron RNA Studied by Mg 2+ -Concentration Jump Simulations. Biophys J 2017. [DOI: 10.1016/j.bpj.2016.11.2006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022] Open
|
94
|
Reddy G, Thirumalai D. Collapse Precedes Folding in Denaturant-Dependent Assembly of Ubiquitin. J Phys Chem B 2017; 121:995-1009. [DOI: 10.1021/acs.jpcb.6b13100] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
95
|
Suvlu D, Samaratunga S, Thirumalai D, Rasaiah JC. Thermodynamics of Helix-Coil Transitions of Polyalanine in Open Carbon Nanotubes. J Phys Chem Lett 2017; 8:494-499. [PMID: 28060517 DOI: 10.1021/acs.jpclett.6b02620] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Understanding structure formation in polypeptide chains and synthetic polymers encapsulated in pores is important in biology and nanotechnology. We present replica exchange molecular dynamics studies of the phase diagram for α-helix formation of capped polyalanine in nanotubes (NT) open to a water reservoir as a function of the NT diameter and hydrophobicity. A helix forms only in a narrow range of diameters, which surprisingly is comparable to the width of the ribosome tunnel. Increasing the hydrophobicity enhances helicity in the NT. Helix formation in polyalanine is driven by a small negative enthalpy and a positive entropy change at ≈300 K, in contrast to the large negative entropy change that destabilizes the helix and favors the coiled state in bulk water. There is an anticorrelation between water density inside the nanotube and structure formation. Confinement-induced helix formation depends on amino acid sequence. There is complete absence of helix in polyglutamine and polyserine confined to a open carbon nanotube.
Collapse
|
96
|
Chakrabarti S, Hinczewski M, Thirumalai D. Phenomenological and microscopic theories for catch bonds. J Struct Biol 2017; 197:50-56. [PMID: 27046010 PMCID: PMC5580263 DOI: 10.1016/j.jsb.2016.03.022] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2016] [Revised: 03/18/2016] [Accepted: 03/30/2016] [Indexed: 12/15/2022]
Abstract
Lifetimes of bound states of protein complexes or biomolecule folded states typically decrease when subject to mechanical force. However, a plethora of biological systems exhibit the counter-intuitive phenomenon of catch bonding, where non-covalent bonds become stronger under externally applied forces. The quest to understand the origin of catch-bond behavior has led to the development of phenomenological and microscopic theories that can quantitatively recapitulate experimental data. Here, we assess the successes and limitations of such theories in explaining experimental data. The most widely applied approach is a phenomenological two-state model, which fits all of the available data on a variety of complexes: actomyosin, kinetochore-microtubule, selectin-ligand, and cadherin-catenin binding to filamentous actin. With a primary focus on the selectin family of cell-adhesion complexes, we discuss the positives and negatives of phenomenological models and the importance of evaluating the physical relevance of fitting parameters. We describe a microscopic theory for selectins, which provides a structural basis for catch bonds and predicts a crucial allosteric role for residues Asn82-Glu88. We emphasize the need for new theories and simulations that can mimic experimental conditions, given the complex response of cell adhesion complexes to force and their potential role in a variety of biological contexts.
Collapse
|
97
|
|
98
|
Vu HT, Chakrabarti S, Hinczewski M, Thirumalai D. Discrete Step Sizes of Molecular Motors Lead to Bimodal Non-Gaussian Velocity Distributions under Force. PHYSICAL REVIEW LETTERS 2016; 117:078101. [PMID: 27564000 DOI: 10.1103/physrevlett.117.078101] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/17/2015] [Indexed: 06/06/2023]
Abstract
Fluctuations in the physical properties of biological machines are inextricably linked to their functions. Distributions of run lengths and velocities of processive molecular motors, like kinesin-1, are accessible through single-molecule techniques, but rigorous theoretical models for these probabilities are lacking. Here, we derive exact analytic results for a kinetic model to predict the resistive force (F)-dependent velocity [P(v)] and run length [P(n)] distribution functions of generic finitely processive molecular motors. Our theory quantitatively explains the zero force kinesin-1 data for both P(n) and P(v) using the detachment rate as the only parameter. In addition, we predict the F dependence of these quantities. At nonzero F, P(v) is non-Gaussian and is bimodal with peaks at positive and negative values of v, which is due to the discrete step size of kinesin-1. Although the predictions are based on analyses of kinesin-1 data, our results are general and should hold for any processive motor, which walks on a track by taking discrete steps.
Collapse
|
99
|
Hori N, Denesyuk NA, Thirumalai D. Salt Effects on the Thermodynamics of a Frameshifting RNA Pseudoknot under Tension. J Mol Biol 2016; 428:2847-59. [PMID: 27315694 DOI: 10.1016/j.jmb.2016.06.002] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2016] [Revised: 06/06/2016] [Accepted: 06/07/2016] [Indexed: 12/24/2022]
Abstract
Because of the potential link between -1 programmed ribosomal frameshifting and response of a pseudoknot (PK) RNA to force, a number of single-molecule pulling experiments have been performed on PKs to decipher the mechanism of programmed ribosomal frameshifting. Motivated in part by these experiments, we performed simulations using a coarse-grained model of RNA to describe the response of a PK over a range of mechanical forces (fs) and monovalent salt concentrations (Cs). The coarse-grained simulations quantitatively reproduce the multistep thermal melting observed in experiments, thus validating our model. The free energy changes obtained in simulations are in excellent agreement with experiments. By varying f and C, we calculated the phase diagram that shows a sequence of structural transitions, populating distinct intermediate states. As f and C are changed, the stem-loop tertiary interactions rupture first, followed by unfolding of the 3'-end hairpin (I⇌F). Finally, the 5'-end hairpin unravels, producing an extended state (E⇌I). A theoretical analysis of the phase boundaries shows that the critical force for rupture scales as (logCm)(α) with α=1(0.5) for E⇌I (I⇌F) transition. This relation is used to obtain the preferential ion-RNA interaction coefficient, which can be quantitatively measured in single-molecule experiments, as done previously for DNA hairpins. A by-product of our work is the suggestion that the frameshift efficiency is likely determined by the stability of the 5'-end hairpin that the ribosome first encounters during translation.
Collapse
|
100
|
Hinczewski M, Thirumalai D. Noise Control in Gene Regulatory Networks with Negative Feedback. J Phys Chem B 2016; 120:6166-77. [DOI: 10.1021/acs.jpcb.6b02093] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|