BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Changbong Hyeon

Changbong Hyeon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softBiomoleculesBiological Physicscond-mat.stat-mechMolecular Networksquant-phphysics.atom-phSubcellular Processes

Catalog footprint

What is connected

31works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2023arXiv

Fundamental trade-off between the speed of light and the Fano factor of photon current in three-level lambda systems

Electromagnetically induced slow-light medium is a promising system for quantum memory devices, but controlling its noise level remains a major challenge to overcome. This work considers the simplest model for such medium, comprised of three-level $Λ$-systems interacting with bosonic bath, and provides a new fundamental trade-off relation in light-matter interaction between the group velocity of light and the Fano factor of photon current due to radiative transitions. Considering the steady state limits of a newly derived Lindblad-type equation, we find that the Fano factor of the photon current maximizes to 3 at the minimal group velocity of light, which holds true universally regardless of detailed values of parameters characterizing the medium.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Polymer brush-induced depletion interactions and clustering of membrane proteins

We investigate the effect of mobile polymer brushes on proteins embedded in biological membranes by employing both Asakura-Oosawa type of theoretical model and coarse-grained molecular dynamics simulations. The brush polymer-induced depletion attraction between proteins changes non-monotonically with the size of brush. The depletion interaction, which is determined by the ratio of protein size to the grafting distance between brush polymers, increases linearly with brush size as long as the polymer brush height is shorter than the protein size. When the brush height exceeds the protein size, however, the depletion attraction among proteins is slightly reduced. We also explore the possibility of brush polymer-induced assembly of a large protein cluster, which can be related to one of many molecular mechanisms underlying recent experimental observations of integrin nanocluster formation and signaling.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The origin of loose bound of the thermodynamic uncertainty relation in a dissipative two-level quantum system

The thermodynamic uncertainty relations (TURs), originally discovered for classical systems, dictate the trade-off between dissipation and fluctuations of irreversible current, specifying a minimal bound that constrains the two quantities. In a series of effort to extend the relation to the one under more generalized conditions, it has been noticed that the bound is less tight in open quantum processes. To study the origin of the loose bounds, we consider an external field-driven transition dynamics of two-level quantum system weakly coupled to the bosonic bath as a model of open quantum system. The model makes it explicit that the imaginary part of quantum coherence, which contributes to dissipation to the environment, is responsible for loosening the TUR bound by suppressing the relative fluctuations in the irreversible current of transitions, whereas the real part of the coherence tightens it. Our study offers a better understanding of how quantum nature affects the TUR bound.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Thermodynamic optimality of glycolytic oscillations

Temporal order in living matters reflects the self-organizing nature of dynamical processes driven out of thermodynamic equilibrium. Because of functional reason, the period of a biochemical oscillation must be tuned to a specific value with precision; however, according to the thermodynamic uncertainty relation (TUR), the precision of oscillatory period is constrained by the thermodynamic cost of generating it. After reviewing the basics of chemical oscillations using Brusselator as a model system, we study the glycolytic oscillation generated by octameric phosphofructokinase (PFK), which is known to display a period of several minutes. By exploring the phase space of glycolytic oscillations, we find that the glycolytic oscillation under the cellular condition is realized in a cost effective manner. Specifically, over the biologically relevant range of parameter values of glycolysis and octameric PFK, the entropy production from the glycolytic oscillation is minimal when the oscillation period is (5 - 10) minutes. Further, the glycolytic oscillation is found at work near the phase boundary of limit cycles, suggesting that a moderate increase of glucose injection rate leads to the loss of oscillatory dynamics, which is reminiscent of the loss of pulsatile insulin release resulting from elevated blood glucose level.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Thermodynamic uncertainty relation to assess biological processes

We review the trade-offs between speed, fluctuations, and thermodynamic cost involved with biological processes in nonequilibrium states, and discuss how optimal these processes are in light of the universal bound set by the thermodynamic uncertainty relation (TUR). The values of the uncertainty product $\mathcal{Q}$ of TUR, which can be used as a measure of the precision of enzymatic processes realized for a given thermodynamic cost, are suboptimal when the substrate concentration $[S]$ is at the Michaelis constant ($K_\text{M}$), and some of the key biological processes are found to work around this condition. We illustrate the utility of $\mathcal{Q}$ in assessing how close the molecular motors and biomass producing machineries are to the TUR bound, and for the cases of biomass production (or biological copying processes) we discuss how their optimality quantified in terms of $\mathcal{Q}$ is balanced with the error rate in the information transfer process. We also touch upon the trade-offs in other error-minimizing processes in biology, such as gene regulation and chaperone-assisted protein folding. A spectrum of $\mathcal{Q}$ recapitulating the biological processes surveyed here provides glimpses into how biological systems are evolved to optimize and balance the conflicting functional requirements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermodynamic cost, speed, fluctuations, and error reduction of biological copy machines

Due to large fluctuations in cellular environments, transfer of information in biological processes without regulation is inherently error-prone. The mechanistic details of error-reducing mechanisms in biological copying processes have been a subject of active research; however, how error reduction of a process is balanced with its thermodynamic cost and dynamical properties remain largely unexplored. Here, we study the error reducing strategies in light of the recently discovered thermodynamic uncertainty relation (TUR) that sets a physical bound to the cost-precision trade-off relevant in general dissipative processes. We found that the two representative copying processes, DNA replication by the exonuclease-deficient T7 DNA polymerase and mRNA translation by the \textit{E. coli} ribosome, reduce the error rates to biologically acceptable levels while also optimizing the processes close to the physical limit dictated by TUR.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Compressing $Θ$-chain in slit geometry

When compressed in a slit of width $D$, a $Θ$-chain that displays the scaling of size $R_0$ (diameter) with respect to the number of monomers $N$, $R_0\sim aN^{1/2}$, expands in the lateral direction as $R_{\parallel}\sim a N^ν(a/D)^{2ν-1}$. Provided that the $Θ$ condition is strictly maintained throughout the compression, the well-known scaling exponent of $Θ$-chain in 2 dimensions, $ν=4/7$, is anticipated in a perfect confinement. However, numerics shows that upon increasing compression from $R_0/D<1$ to $R_0/D\gg 1$, $ν$ gradually deviates from $ν=1/2$ and plateaus at $ν=3/4$, the exponent associated with the self-avoiding walk in two dimensions. Using both theoretical considerations and numerics, we argue that it is highly nontrivial to maintain the $Θ$ condition under confinement because of two major effects. First, as the dimension is reduced from 3 to 2 dimensions, the contributions of higher order virial terms, which can be ignored in 3 dimensions at large $N$, become significant. Second and more importantly, the geometrical confinement, which is regarded as an applied external field, alters the second virial coefficient ($B_2$) changes from $B_2=0$ ($Θ$ condition) in free space to $B_2>0$ (good-solvent condition) in confinement. Our study provides practical insight into how confinement affects the conformation of a single polymer chain.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Theoretical Perspectives on Biological Machines

Many biological functions are executed by molecular machines, which consume energy and convert it into mechanical work. Biological machines have evolved to transport cargo, facilitate folding of proteins and RNA, remodel chromatin and replicate DNA. A common aspect of these machines is that their functions are driven out of equilibrium. It is a challenge to provide a general framework for understanding the functions of biological machines, such as molecular motors, molecular chaperones, and helicases. Using these machines as prototypical examples, we describe a few general theoretical methods providing insights into their functions. Although the theories rely on coarse-graining of these complex systems they have proven useful in not only accounting for many in vitro experiments but also addressing questions such as how the trade-off between precision, energetic costs and optimal performances are balanced. However, many complexities associated with biological machines will require one to go beyond current theoretical methods. Simple point mutations in the enzyme could drastically alter functions, making the motors bi-directional or result in unexpected diseases or dramatically restrict the capacity of molecular chaperones to help proteins fold. These examples are reminders that while the search for principles of generality in biology is intellectually stimulating, one also ought to keep in mind that molecular details must be accounted for to develop a deeper understanding of processes driven by biological machines. Going beyond generic descriptions of in vitro behavior to making genuine understanding of in vivo functions will likely remain a major challenge for some time to come. The combination of careful experiments and the use of physical principles will be useful in elucidating the rules governing the workings of biological machines.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Symmetry, Rigidity, and Allosteric Signaling: From Monomeric Proteins to Molecular Machines

Allosteric signaling in biological molecules, which may be viewed as specific action at a distance due to localized perturbation upon binding of ligands or changes in environmental cues, is pervasive in biology. Phenomenological MWC and KNF models galvanized research for over five decades, and these models continue to be the basis for describing the allostery in an array of systems. However, understanding allosteric signaling and the associated dynamics between distinct allosteric states is challenging. In this review, we first describe symmetry and rigidity as essential requirements for allosteric proteins. The general features, with MWC and KNF as two extreme scenarios, emerge when allosteric signaling is viewed from an energy landscape perspective. To go beyond the general theories, we describe computational tools to predict the network of residues that carry allosteric signals. Methods to obtain molecular insights into the dynamics of allosteric transitions are briefly mentioned. The utility of the methods is illustrated by applications to systems ranging from monomeric proteins to multisubunit proteins. Finally, the role allostery plays in the functions of ATP-consuming molecular machines, bacterial chaperonin GroEL and molecular motors, is described. Although universal molecular principles governing allosteric signaling do not exist, we can draw the following general conclusions. (1) Multiple pathways connecting allosteric states are highly heterogeneous. (2) Allosteric signaling is exquisitely sensitive to the architecture of the system, which implies that the capacity for allostery is encoded in the structure itself. (3) The mechanical modes that connect distinct allosteric states are robust to sequence variations. (4) Extensive investigations of allostery in Hemoglobin and more recently GroEL, show that a network of salt-bridge rearrangements serves as allosteric switches.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Contact statistics highlight distinct organizing principles of proteins and RNA

Although both RNA and proteins have densely packed native structures, chain organizations of these two biopolymers are fundamentally different. Motivated by the recent discoveries in chromatin folding that interphase chromosomes have territorial organization with signatures pointing to metastability, we analyzed the biomolecular structures deposited in the Protein Data Bank and found that the intrachain contact probabilities, $P(s)$ as a function of the arc length $s$, decay in power-law $\sim s^{-γ}$ over the intermediate range of $s$, $10\lesssim s\lesssim 110$. We found that the contact probability scaling exponent is $γ\approx 1.11$ for large RNA ($N>110$), $γ\approx 1.41$ for small sized RNA ($N<110$), and $γ\approx 1.65$ for proteins. Given that Gaussian statistics is expected for a fully equilibrated chain in polymer melts, the deviation of $γ$ value from $γ=1.5$ for the subchains of large RNA in the native state suggests that the chain configuration of RNA is not fully equilibrated. It is visually clear that folded structures of large sized RNA ($N\gtrsim 110$) adopt crumpled structures, partitioned into modular multi-domains assembled by proximal sequences along the chain, whereas the polypeptide chain of folded proteins looks better mixed with the rest of the structure. Our finding of $γ\approx 1$ for large RNA might be an ineluctable consequence of the hierarchical ordering of the secondary to tertiary elements in the folding process.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ultraslow water-mediated transmembrane interactions regulate the activation of A$_{\text{2A}}$ adenosine receptor

Water molecules inside G-protein coupled receptor have recently been spotlighted in a series of crystal structures. To decipher the dynamics and functional roles of internal waters in GPCR activity, we studied A$_{\text{2A}}$ adenosine receptor using $μ$sec-molecular dynamics simulations. Our study finds that the amount of water flux across the transmembrane (TM) domain varies depending on the receptor state, and that the water molecules of the TM channel in the active state flow three times slower than those in the inactive state. Depending on the location in solvent-protein interface as well as the receptor state, the average residence time of water in each residue varies from $\sim\mathcal{O}(10^2)$ psec to $\sim\mathcal{O}(10^2)$ nsec. Especially, water molecules, exhibiting ultraslow relaxation ($\sim\mathcal{O}(10^2)$ nsec) in the active state, are found around the microswitch residues that are considered activity hotspots for GPCR function. A continuous allosteric network spanning the TM domain, arising from water-mediated contacts, is unique in the active state, underscoring the importance of slow waters in the GPCR activation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Confinement-induced glassy dynamics in a model for chromosome organization

Recent experiments showing scaling of the intrachromosomal contact probability, $P(s)\sim s^{-1}$ with the genomic distance $s$, are interpreted to mean a self-similar fractal-like chromosome organization. However, scaling of $P(s)$ varies across organisms, requiring an explanation. We illustrate dynamical arrest in a highly confined space as a discriminating marker for genome organization, by modeling chromosome inside a nucleus as a homopolymer confined to a sphere of varying sizes. Brownian dynamics simulations show that the chain dynamics slows down as the polymer volume fraction ($ϕ$) inside the confinement approaches a critical value $ϕ_c$. The universal value of $ϕ_c^{\infty}\approx 0.44$ for a sufficiently long polymer ($N\gg 1$) allows us to discuss genome dynamics using $ϕ$ as a single parameter. Our study shows that the onset of glassy dynamics is the reason for the segregated chromosome organization in human ($N\approx 3\times 10^9$, $ϕ\gtrsimϕ_c^{\infty}$), whereas chromosomes of budding yeast ($N\approx 10^8$, $ϕ<ϕ_c^{\infty}$) are equilibrated with no clear signature of such organization.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of macromolecular crowding on the collapse of biopolymers

Experiments show that macromolecular crowding modestly reduces the size of intrinsically disordered proteins (IDPs) even at volume fraction ($ϕ$) similar to that in the cytosol whereas DNA undergoes a coil-to-globule transition at very small $ϕ$. We show using a combination of scaling arguments and simulations that the polymer size $\overline{R}_g(ϕ)$ depends on $x = \overline{R}_g(0)/D$ where $D$ is the $ϕ$-dependent distance between the crowders. If $x\lesssim \mathcal{O}(1)$, there is only a small decrease in $\overline{R}_g(ϕ)$ as $ϕ$ increases. When $x\gg \mathcal{O}(1)$, a cooperative coil-to-globule transition is induced. Our theory quantitatively explains a number of experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Extracting folding landscape characteristics of biomolecules using mechanical forces

In recent years single molecule force spectroscopy has opened a new avenue to provide profiles of the complex energy landscape of biomolecules. In this field, quantitative analyses of the data employing sound theoretical models, have played a major role in interpreting data and anticipating outcomes of experiments. Here, we explain how by using temperature as a variable in mechanical unfolding of biomolecules in force spectroscopy, the roughness of the energy landscape can be measured without making any assumptions about the underlying reaction coordinate. Estimates of other aspects of the energy landscape such as free energy barriers or the transition state (TS) locations could depend on the precise model used to analyze the experimental data. We illustrate the inherent difficulties in obtaining the transition state location from loading rate or force-dependent unfolding rates. Because the transition state moves as the force or the loading rate is varied, which is reminiscent of the Hammond effect, it is in general difficult to invert the experimental data. The independence of the TS location on force holds good only for brittle or hard biomolecules whereas the TS location changes considerably if the molecule is soft or plastic. Finally, we discuss the goodness of the end-to-end distance (or pulling) coordinate of the molecule as a surrogate reaction coordinate in a situation like force-induced ligand-unbinding from a gated molecular receptor as well as force-quench refolding of an RNA hairpin.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Time traces of individual kinesin motors suggest functional heterogeneity

Conventional analysis of in vitro assays of motor proteins rests on the assumption that all proteins with the same chemical composition function identically; however molecule-to-molecule variation is often seen even in well-controlled experiments. In an effort to obtain a statistically meaningful set of time traces that simultaneously avoid any experimental artifacts, we performed quantum-dot labeled kinesin experiments on both surface and levitated microtubules. Similar to glassy systems, we found that mean velocities of individual kinesin motors vary widely from one motor to another, the variation of which is greater than that expected from the stochastic variation of stepping times. In the presence of heterogeneity, an ensemble-averaged quantity such as diffusion constant or randomness parameter is ill-defined. We propose to analyze heterogeneous data from single molecule measurements by decomposing them into homogeneous subensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Unexpected swelling of stiff DNA in a polydisperse crowded environment

We investigate the conformations of DNA-like stiff chains, characterized by contour length ($L$) and persistence length ($l_p$), in a variety of crowded environments containing mono disperse soft spherical (SS) and spherocylindrical (SC) particles, mixture of SS and SC, and a milieu mimicking the composition of proteins in $E. coli.$ cytoplasm. The stiff chain, whose size modestly increases in SS crowders up to $ϕ\approx 0.1$, is considerably more compact at low volume fractions ($ϕ\leq 0.2$) in monodisperse SC particles than in a medium containing SS particles. A 1:1 mixture of SS and SC crowders induces greater chain compaction than the pure SS or SC crowders at the same $ϕ$ with the effect being highly non-additive. We also discover a counter-intuitive result that polydisperse crowding environment, mimicking the composition of a cell lysate, swells the DNA-like polymer, which is in stark contrast to the size reduction of flexible polymer in the same milieu. Trapping of the stiff chain in a fluctuating tube-like environment created by large-sized crowders explains the dramatic increase in size and persistence length of the stiff chain. In the polydisperse medium, mimicking the cellular environment, the size of the DNA (or related RNA) is determined by $L/l_p$. At low $L/l_p$ the size of the polymer is unaffected whereas there is a dramatic swelling at intermediate value of $L/l_p$. We use these results to provide insights into recent experiments on crowding effects on RNA, and also make testable predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Communication over the network of binary switches regulates the activation of A$_{2A}$ adenosine receptor

Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately regulated by the type of ligands that bind to the orthosteric or allosteric binding sites. To glean the structural and dynamical origin of ligand-dependent modulation of GPCR activity, we performed total $\sim$ 5 $μ$sec molecular dynamics simulations of A$_{2A}$ adenosine receptor (A$_{2A}$AR) in its apo, antagonist-bound, and agonist-bound forms in an explicit water and membrane environment, and examined the corresponding dynamics and correlation between the 10 key structural motifs that serve as the allosteric hotspots in intramolecular signaling network. We dubbed these 10 structural motifs "binary switches" as they display molecular interactions that switch between two distinct states. By projecting the receptor dynamics on these binary switches that yield $2^{10}$ microstates, we show that (i) the receptors in apo, antagonist-bound, and agonist-bound states explore vastly different conformational space; (ii) among the three receptor states the apo state explores the broadest range of microstates; (iii) in the presence of the agonist, the active conformation is maintained through coherent couplings among the binary switches; and (iv) to be most specific, our analysis shows that W246, located deep inside the binding cleft, can serve as both an agonist sensor and actuator of ensuing intramolecular signaling for the receptor activation.Finally, our analysis of multiple trajectories generated by inserting an agonist to the apo state underscores that the transition of the receptor from inactive to active form requires the disruption of ionic-lock in the DRY motif.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Development and applications of Coarse Grained models for RNA

In contrast to proteins much less attention has been focused on development of computational models for describing RNA molecules, which are being recognized as playing key roles in many cellular functions. Current atomically detailed force fields are not accurate enough to capture the properties of even simple nucleic acid constructs. In this article, we review our efforts to develop coarse-grained (CG) models that capture the underlying physics for the particular length scale of interest. Two models are discussed. One of them is the Three Interaction Site (TIS) model in which each nucleotide is represented by three beads corresponding to sugar, phosphate, and base. The other is the Self-Organized Polymer (SOP) model in which each nucleotide is represented as a single interaction center. Applications of the TIS model to study the complexity of hairpin formation and the effects of crowding in shifting equilibrium between two conformations in human telomerase pseudoknot are described. The work on crowding illustrates a direct link to the activity of telomerase. We use the SOP model to describe the response of {\it Tetrahymena} ribozyme to force. The simulated unfolding pathways agree well with single molecule pulling experiments. We also review predictions for the unfolding pathways for {\it Azoarcus} ribozyme. The success of the CG applications to describe dynamics in RNA gives hope that more complex processes involving RNA-protein interactions can be tackled using variants of the proposed models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamical Transition and Heterogeneous Hydration Dynamics in RNA

Enhanced dynamical fluctuations of RNAs, facilitated by a network of water molecules with strong interactions with RNA, are suspected to be critical in their ability to respond to a variety of cellular signals. Using atomically detailed molecular dynamics simulations at various temperatures of purine (adenine)- and preQ$_1$ sensing riboswitch aptamers, we show that water molecules in the vicinity of RNAs undergo complex dynamics depending on the local structures of the RNAs. The overall lifetimes of hydrogen bonds (HBs) of surface bound waters are more than at least 1-2 orders of magnitude longer than bulk water. Slow hydration dynamics, revealed in non-Arrhenius behavior of the relaxation time, arises from high activation barriers to break water hydrogen bonds with a nucleotide and by reduced diffusion of water. The relaxation kinetics at specific locations in the two RNAs show a broad spectrum of time scales reminiscent of glass-like behavior, suggesting that the hydration dynamics is highly heterogeneous. Both RNAs undergo dynamic transition at $T = T_D \gtrsim 200$ K as assessed by the mean square fluctuation of hydrogen atoms $\langle x^2\rangle$, which undergoes an abrupt harmonic-to-anharmonic transition at $T_D$. The near universal value of $T_D$ found for these RNAs and previously for tRNA is strongly correlated with changes in hydration dynamics as $T$ is altered. Hierarchical dynamics of waters associated with the RNA surface, revealed in the motions of distinct classes of water with well-separated time scales, reflects the heterogeneous local environment on the molecular surface of RNA. At low temperatures slow water dynamics predominates over structural transitions. Our study demonstrates that the complex interplay of dynamics between water and local environment in the RNA structures could be a key determinant of the functional activities of RNA.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mapping the intramolecular signal transduction of G-protein coupled receptors

G-protein coupled receptors (GPCRs), a major gatekeeper of extracellular signals on plasma membrane, are unarguably one of the most important therapeutic targets. Given the recent discoveries of allosteric modulations, an allosteric wiring diagram of intramolecular signal transductions would be of great use to glean the mechanism of receptor regulation. Here, by evaluating betweenness centrality ($C_B$) of each residue, we calculate maps of information flow in GPCRs and identify key residues for signal transductions and their pathways. Compared with preexisting approaches, the allosteric hotspots that our $C_B$-based analysis detects for A$_{2A}$ adenosine receptor (A$_{2A}$AR) and bovine rhodopsin are better correlated with biochemical data. In particular, our analysis outperforms other methods in locating the rotameric microswitches, which are generally deemed critical for mediating orthosteric signaling in class A GPCRs. For A$_{2A}$AR, the inter-residue cross-correlation map, calculated using equilibrium structural ensemble from molecular dynamics simulations, reveals that strong signals of long-range transmembrane communications exist only in the agonist-bound state. A seemingly subtle variation in structure, found in different GPCR subtypes or imparted by agonist bindings or a point mutation at an allosteric site, can lead to a drastic difference in the map of signaling pathways and protein activity. The signaling map of GPCRs provides valuable insights into allosteric modulations as well as reliable identifications of orthosteric signaling pathways.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chain Length Determines the Folding Rates of RNA

We show that the folding rates (k_F) of RNA are determined by N, the number of nucleotides. By assuming that the distribution of free energy barriers separating the folded and the unfolded states is Gaussian, which follows from central limit theorem arguments and polymer physics concepts, we show that k_F ~ k_0 exp(-alpha N^0.5). Remarkably, the theory fits the experimental rates spanning over seven orders of magnitude with k_0 ~ 1.0 (microsec)^{-1}. An immediate consequence of our finding is that the speed limit of RNA folding is about one microsecond just as it is in the folding of globular proteins.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hidden Complexity in the Isomerization Dynamics of Holliday Junctions

A plausible consequence of rugged energy landscapes of biomolecules is that functionally competent folded states may not be unique, as is generally assumed. Indeed, molecule-to-molecule variations in the dynamics of enzymes and ribozymes under folding conditions have recently been identified in single molecule experiments. However, systematic quantification and the structural origin of the observed complex behavior remain elusive. Even for a relatively simple case of isomerization dynamics in Holliday Junctions (HJs), molecular heterogeneities persist over a long observation time (Tobs ~ 40 sec). Here, using concepts in glass physics and complementary clustering analysis, we provide a quantitative method to analyze the smFRET data probing the isomerization in HJ dynamics. We show that ergodicity of HJ dynamics is effectively broken; as a result, the conformational space of HJs is partitioned into a folding network of kinetically disconnected clusters. While isomerization dynamics in each cluster occurs rapidly as if the associated conformational space is fully sampled, distinct patterns of time series belonging to different clusters do not interconvert on Tobs. Theory suggests that persistent heterogeneity of HJ dynamics is a consequence of internal multiloops with varying sizes and flexibilities frozen by Mg2+ ions. An annealing experiment using Mg2+ pulse that changes the Mg2+ cocentration from high to low to high values lends support to this idea by explicitly showing that interconversions can be driven among trajectories with different patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multiple barriers in forced rupture of protein complexes

Curvatures in the most probable rupture force ($f^*$) versus log-loading rate ($\log{r_f}$) observed in dynamic force spectroscopy (DFS) on biomolecular complexes are interpreted using a one-dimensional free energy profile with multiple barriers or a single barrier with force-dependent transition state. Here, we provide a criterion to select one scenario over another. If the rupture dynamics occurs by crossing a single barrier in a physical free energy profile describing unbinding, the exponent $ν$, from $(1- f^*/f_c)^{1/ν}\sim(\log r_f)$ with $f_c$ being a critical force in the absence of force, is restricted to $0.5 \leq ν\leq 1$. For biotin-ligand complexes and leukocyte-associated antigen-1 bound to intercellular adhesion molecules, which display large curvature in the DFS data, fits to experimental data yield $ν<0.5$, suggesting that ligand unbinding is associated with multiple-barrier crossing.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

RNA under Tension: Folding Landscapes, Kinetic Partitioning Mechanism, and Molecular Tensegrity

Non-coding RNA sequences play a great role in controlling a number of cellular functions, thus raising the need to understand their complex conformational dynamics in quantitative detail. In this perspective, we first show that single molecule pulling experiments when combined with with theory and simulations can be used to quantitatively explore the folding landscape of nucleic acid hairpins, and riboswitches with tertiary interactions. Applications to riboswitches, which are non-coding RNA elements that control gene expression by undergoing dynamical conformational changes in response to binding of metabolites, lead to an organization principle that assembly of RNA is determined by the stability of isolated helices. We also point out the limitations of single molecule pulling experiments, with molecular extension as the only accessible parameter, in extracting key parameters of the folding landscapes of RNA molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Weak temporal signals can synchronize and accelerate the transition dynamics of biopolymers under tension

In addition to thermal noise, which is essential to promote conformational transitions in biopolymers, cellular environment is replete with a spectrum of athermal fluctuations that are produced from a plethora of active processes. To understand the effect of athermal noise on biological processes, we studied how a small oscillatory force affects the thermally induced folding and unfolding transition of an RNA hairpin, whose response to constant tension had been investigated extensively in both theory and experiments. Strikingly, our molecular simulations performed under overdamped condition show that even at a high (low) tension that renders the hairpin (un)folding improbable, a weak external oscillatory force at a certain frequency can synchronously enhance the transition dynamics of RNA hairpin and increase the mean transition rate. Furthermore, the RNA dynamics can still discriminate a signal with resonance frequency even when the signal is mixed among other signals with nonresonant frequencies. In fact, our computational demonstration of thermally induced resonance in RNA hairpin dynamics is a direct realization of the phenomena called stochastic resonance (SR) and resonant activation (RA). Our study, amenable to experimental tests using optical tweezers, is of great significance to the folding of biopolymers in vivo that are subject to the broad spectrum of cellular noises.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A structural perspective on the dynamics of kinesin motors

Despite significant fluctuation under thermal noise, biological machines in cells perform their tasks with exquisite precision. Using molecular simulation of a coarse-grained model and theoretical arguments we envisaged how kinesin, a prototype of biological machines, generates force and regulates its dynamics to sustain persistent motor action. A structure based model, which can be versatile in adapting its structure to external stresses while maintaining its native fold, was employed to account for several features of kinesin dynamics along the biochemical cycle. This analysis complements our current understandings of kinesin dynamics and connections to experiments. We propose a thermodynamic cycle for kinesin that emphasizes the mechanical and regulatory role of the neck-linker and clarify issues related the motor directionality, and the difference between the external stalling force and the internal tension responsible for the head-head coordination. The comparison between the thermodynamic cycle of kinesin and macroscopic heat engines highlights the importance of structural change as the source of work production in biomolecular machines.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Capturing the essence of folding and functions of biomolecules using Coarse-Grained Models

The distances over which biological molecules and their complexes can function range from a few nanometres, in the case of folded structures, to millimetres, for example during chromosome organization. Describing phenomena that cover such diverse length, and also time scales, requires models that capture the underlying physics for the particular length scale of interest. Theoretical ideas, in particular, concepts from polymer physics, have guided the development of coarse-grained models to study folding of DNA, RNA, and proteins. More recently, such models and their variants have been applied to the functions of biological nanomachines. Simulations using coarse-grained models are now poised to address a wide range of problems in biology.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Compaction and tensile forces determine the accuracy of folding landscape parameters from single molecule pulling experiments

We establish a framework for assessing whether the transition state location of a biopolymer, which can be inferred from single molecule pulling experiments, corresponds to the ensemble of structures that have equal probability of reaching either the folded or unfolded states (Pfold = 0.5). Using results for the forced-unfolding of a RNA hairpin, an exactly soluble model and an analytic theory, we show that Pfold is solely determined by s, an experimentally measurable molecular tensegrity parameter, which is a ratio of the tensile force and a compaction force that stabilizes the folded state. Applications to folding landscapes of DNA hairpins and leucine zipper with two barriers provide a structural interpretation of single molecule experimental data. Our theory can be used to assess whether molecular extension is a good reaction coordinate using measured free energy profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring the energy landscape of biopolymers using single molecule force spectroscopy and molecular simulations

In recent years, single molecule force techniques have opened a new avenue to decipher the folding landscapes of biopolymers by allowing us to watch and manipulate the dynamics of individual proteins and nucleic acids. In single molecule force experiments, quantitative analyses of measurements employing sound theoretical models and molecular simulations play central role more than any other field. With a brief description of basic theories for force mechanics and molecular simulation technique using self-organized polymer (SOP) model, this chapter will discuss various issues in single molecule force spectroscopy (SMFS) experiments, which include pulling speed dependent unfolding pathway, measurement of energy landscape roughness, the in uence of molecular handles in optical tweezers on measurement and molecular motion, and folding dynamics of biopolymers under force quench condition.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Theoretical Perspectives on Protein Folding

Understanding how monomeric proteins fold under in vitro conditions is crucial to describing their functions in the cellular context. Significant advances both in theory and experiments have resulted in a conceptual framework for describing the folding mechanisms of globular proteins. The experimental data and theoretical methods have revealed the multifaceted character of proteins. Proteins exhibit universal features that can be determined using only the number of amino acid residues (N) and polymer concepts. The sizes of proteins in the denatured and folded states, cooperativity of the folding transition, dispersions in the melting temperatures at the residue level, and time scales of folding are to a large extent determined by N. The consequences of finite N especially on how individual residues order upon folding depends on the topology of the folded states. Such intricate details can be predicted using the Molecular Transfer Model that combines simulations with measured transfer free energies of protein building blocks from water to the desired concentration of the denaturant. By watching one molecule fold at a time, using single molecule methods, the validity of the theoretically anticipated heterogeneity in the folding routes, and the N-dependent time scales for the three stages in the approach to the native state have been established. Despite the successes of theory, of which only a few examples are documented here, we conclude that much remains to be done to solve the "protein folding problem" in the broadest sense.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Refolding dynamics of stretched biopolymers upon force quench

Single molecule force spectroscopy methods can be used to generate folding trajectories of biopolymers from arbitrary regions of the folding landscape. We illustrate the complexity of the folding kinetics and generic aspects of the collapse of RNA and proteins upon force quench, using simulations of an RNA hairpin and theory based on the de Gennes model for homopolymer collapse. The folding time, $τ_F$, depends asymmetrically on $δf_S = f_S - f_m$ and $δf_Q = f_m - f_Q$ where $f_S$ ($f_Q$) is the stretch (quench) force, and $f_m$ is the transition mid-force of the RNA hairpin. In accord with experiments, the relaxation kinetics of the molecular extension, $R(t)$, occurs in three stages: a rapid initial decrease in the extension is followed by a plateau, and finally an abrupt reduction in $R(t)$ that occurs as the native state is approached. The duration of the plateau increases as $λ=τ_Q/τ_F$ decreases (where $τ_Q$ is the time in which the force is reduced from $f_S$ to $f_Q$). Variations in the mechanisms of force quench relaxation as $λ$ is altered are reflected in the experimentally measurable time-dependent entropy, which is computed directly from the folding trajectories. An analytical solution of the de Gennes model under tension reproduces the multistage stage kinetics in $R(t)$. The prediction that the initial stages of collapse should also be a generic feature of polymers is validated by simulation of the kinetics of toroid (globule) formation in semiflexible (flexible) homopolymers in poor solvents upon quenching the force from a fully stretched state. Our findings give a unified explanation for multiple disparate experimental observations of protein folding.

0 citations0 reviewsNo code linkedOpen access