Publications

2018
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. Inorganic Chemistry 2018, 57, 3993–4008. Publisher's VersionAbstract

Abstract Image

Detailed investigations of the electronic structure and bonding scenario in different carbene–phosphinidenes have been presented using state-of-the-art computational methods (BP86/def2-TZVPP//BP86/def2-SVP). We have endeavored to find the correlation of the calculated 31P chemical shifts with different bonding parameters of compounds to access the relative π-acceptor strengths of the carbenes. 31P chemical shifts exhibit a weak correlation with σ-polarizations of Ccarb–P bonds toward phosphorus; however excellent correlations are obtained in the case of π-polarizations of Ccarb–P bonds toward the carbene carbon (Ccarb) and NPA charges on phosphorus atoms. 31P chemical shifts also show excellent correlations with the electron densities and energy densities of Ccarb–P bonds at BCPs, as suggested by QTAIM calculations. Moreover, EDA-NOCV analysis is implemented to gain brief insight into the bonding scenario in this class of compounds. Good correlation exists between the interaction energies between the carbene and PPh fragments and 31P chemical shifts. Additionally, we have investigated the correlations of calculated 31P chemical shifts with different bonding parameters of the corresponding free carbenes. The bonding scenario in different carbene-substituted phosphinidenes is also explored to see how the bonding situation depends on various substituents on phosphinidenes. The other substituted carbene–phosphinidenes show correlations similar to those of carbene–phenylphosphinidenes.

computational_investigation_of_carbenephosphinidenes.pdf
Thirumalai, D. ; Toan, N. M. Forced-rupture of cell-adhesion complexes reveals abrupt switch between two brittle states. Journal of Chemical Physics 2018. Publisher's Version forced-rupture_of_cell-adhesion_complexes_reveals_abrupt_switch_between_two_brittle_states.pdf
Thirumalai, D. ; Mugnai, M. L. ; Goldtzvik, Y. Dynamics of Allosteric Transitions in Dynein. bioRxiv 2018. Publisher's Version dynamics_of_allosteric_transitions_in_dynei.pdf
Nguyen, H. T. ; Thirumalai, D. Calculation of Ion-Dependent RNA Folding Free Energy using Coarse-Grained Simulation. Biophysical Journal 2018. Publisher's Version calculation_of_ion-dependent_rna_folding_free_energy_using_coarse-grained_simulation.pdf
Thirumalai, D. ; Li, X. ; Samanta, H. S. ; Sinha, S. ; Malmi-Kakkada, A. N. Cell Growth Rate Dictates the Onset of Glass to Fluidlike Transition and Long Time Superdiffusion in an Evolving Cell Colony. Physical Review 2018. Publisher's Version cell_growth_rate_dictates.pdf
Thirumalai, D. ; Hyeon, C. Signalling networks and dynamics of allosteric transitions in bacterial chaperonin GroEL: implications for iterative annealing of misfolded proteins. The Royal Society Publishing 2018. Publisher's VersionAbstract

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’.

signalling_networks_and_dynamics_of.pdf
Thirumalai, D. ; Mugnai, M. L. Hydration of Magnesium is Required for Myosin VI Phosphate Release. Biophysical Journal 2018, 318A. Publisher's Version hydration_of_magnesium_is_required_for_myosin_vi_phosphate_release.pdf
Thirumalai, D. ; Liu, Z. Denaturants alter the flux through multiple pathways in the folding of PDZ domain. The Journal of Physical Chemistry B, 2018 - ACS Publications 2018. Publisher's Version Denaturants Alter the flux through multiple pathways in the folding of pdz domain.pdf
2017
Thirumalai, D. ; Lee, Y. ; Hyeon, C. Ultrasensitivity of water exchange kinetics to the size of metal ion. Journal of the American Chemical Society 2017, 139, 12334-12337. Publisher's Version ultrasensitivity-of-water-exchange-kinetics-to-the-size-of-metal-ion.pdf
Thirumalai, D. ; Suvlu, D. ; Samaratunga, S. ; Rasaiah, J. Thermodynamics of Helix–Coil Transitions of Polyalanine in Open Carbon Nanotubes. The Journal of Physical Chemistry Letters 2017, 8 494-499. Publisher's Version thermodynamics-of-helix-coil-transitions-of-polyalanine-in-open-carbon-nanotubes.pdf
Thirumalai, D. ; Hyeon, C. Ripping RNA by force using Gaussian Network Models. The Journal of Physical Chemistry B 2017, 121, 3515-3522. Publisher's Version ripping-rna-by-force-using-gaussian-network-models.pdf
Thirumalai, D. ; Reddy, G. Collapse precedes folding in denaturant-dependent assembly of ubiquitin. The Journal of Physical Chemistry B 2017, 121, 995-2009. Publisher's Version collapse-precedes-folding-in-denaturant-assembly-of-ubiquitin.pdf
Thirumalai, D. ; Shi, G. Chromatin is stretched but intact when the nucleus is squeezed through constrictions. Biophysical Journal 2017, 112, 411-412. Publisher's Version chromatin-is-stretched-but-intact-when-the-nucleus-is-squeezed-through-constrictions.pdf
Thirumalai, D. ; Chakraborti, S. ; Hinczewski, M. Phenomenological and microscopic theories for catch bonds. Journal of Structural Biology 2017, 197, 50-56. Publisher's Version phenomenological-and-microscopic-theories-for-catch-bonds.pdf
Mugnai, M. L. ; Thirumalai, D. Kinematics of the lever arm swing in myosin VI. Proc Natl Acad Sci U S A 2017, 114, E4389–E4398. DOI:10.1073/pnas.1615708114Abstract

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.

kinematics-of-the-lever-arm-swing-in-myosin_vi.pdf
Samanta, H. S. ; Hinczewski, M. ; Thirumalai, D. Optimal information transfer in enzymatic networks: A field theoretic formulation. Phys Rev E 2017, 96, 012406. DOI:10.1103/PhysRevE.96.012406Abstract

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)]. 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.

optimal-information-transfer-in-enzymatic-networks-a-field-theoretic-formulation.pdf
Samanta, H. S. ; Zhuravlev, P. I. ; Hinczewski, M. ; Hori, N. ; Chakrabarti, S. ; Thirumalai, D. Protein collapse is encoded in the folded state architecture. Soft Matter 2017, 13, 3622-3638. DOI:10.1039/C7SM00074JAbstract

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”.

protein-collapse-is-encoded-in-the-folded-state-architecture.pdf
2016
Vu, H. T. ; Chakrabarti, S. ; Hinczewski, M. ; Thirumalai, D. Discrete Step Sizes of Molecular Motors Lead to Bimodal Non-Gaussian Velocity Distributions under Force. Phys. Rev. Lett. 2016, 117, 078101. DOI:10.1103/PhysRevLett.117.078101Abstract

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 FP(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.

discrete-step-sizes-of-molecular-motors-lead-to-bimodal-non-gaussian-velocity-distributions-under-force.pdf
Hinczewski, M. ; Hyeon, C. ; Thirumalai, D. Directly measuring single-molecule heterogeneity using force spectroscopy. Proc Natl Acad Sci U S A 2016, 113, E3852-61.Abstract
One of the most intriguing results of single-molecule experiments on proteins and nucleic acids is the discovery of functional heterogeneity: the observation that complex cellular machines exhibit multiple, biologically active conformations. The structural differences between these conformations may be subtle, but each distinct state can be remarkably long-lived, with interconversions between states occurring only at macroscopic timescales, fractions of a second or longer. Although we now have proof of functional heterogeneity in a handful of systems-enzymes, motors, adhesion complexes-identifying and measuring it remains a formidable challenge. Here, we show that evidence of this phenomenon is more widespread than previously known, encoded in data collected from some of the most well-established single-molecule techniques: atomic force microscopy or optical tweezer pulling experiments. We present a theoretical procedure for analyzing distributions of rupture/unfolding forces recorded at different pulling speeds. This results in a single parameter, quantifying the degree of heterogeneity, and also leads to bounds on the equilibration and conformational interconversion timescales. Surveying 10 published datasets, we find heterogeneity in 5 of them, all with interconversion rates slower than 10 s(-1) Moreover, we identify two systems where additional data at realizable pulling velocities is likely to find a theoretically predicted, but so far unobserved crossover regime between heterogeneous and nonheterogeneous behavior. The significance of this regime is that it will allow far more precise estimates of the slow conformational switching times, one of the least understood aspects of functional heterogeneity.
directly-measuring-single-molecule-heterogeneity-using-force-spectroscopy.pdf
Liu, Z. ; Reddy, G. ; Thirumalai, D. Folding PDZ2 Domain Using the Molecular Transfer Model. J Phys Chem B 2016. folding-pdz2-domain-using-the-molecular-transfer-model.pdf

Pages