1
|
Pérez-Ramírez HA, Moncho-Jordá A, Odriozola G. Phenol release from pNIPAM hydrogels: scaling molecular dynamics simulations with dynamical density functional theory. SOFT MATTER 2022; 18:8271-8284. [PMID: 36278506 DOI: 10.1039/d2sm01083f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
We employed molecular dynamic simulations (MD) and the Bennett's acceptance ratio method to compute the free energy of transfer, ΔGtrans, of phenol, methane, and 5-fluorouracil (5-FU), between bulk water and water-pNIPAM mixtures of different polymer volume fractions, ϕp. For this purpose, we first calculate the solvation free energies in both media to obtain ΔGtrans. Phenol and 5-FU (a medication used to treat cancer) attach to the pNIPAM surface so that they show negative values of ΔGtrans irrespective of temperature (above or below the lower critical solution temperature of pNIPAM, Tc). Conversely, methane switches the ΔGtrans sign when considering temperatures below (positive) and above (negative) Tc. In all cases, and contrasting with some theoretical predictions, ΔGtrans maintains a linear behavior with the pNIPAM concentration up to large polymer densities. We have also employed MD to compute the diffusion coefficient, D, of phenol in water-pNIPAM mixtures as a function of ϕp in the diluted limit. Both ΔGtrans and D as a function of ϕp are required inputs to obtain the release halftime of hollow pNIPAM microgels through Dynamic Density Functional Theory (DDFT). Our scaling strategy captures the experimental value of 2200 s for 50 μm radius microgels with no cavity, for ϕp ≃ 0.83 at 315 K.
Collapse
Affiliation(s)
- H A Pérez-Ramírez
- Área de Física de Procesos Irreversibles, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana-Azcapotzalco, Avenida San Pablo 180, 02200 Ciudad de México, Mexico.
| | - A Moncho-Jordá
- Departamento de Física Aplicada, Universidad de Granada, Campus Fuentenueva S/N, 18071 Granada, Spain
| | - G Odriozola
- Área de Física de Procesos Irreversibles, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana-Azcapotzalco, Avenida San Pablo 180, 02200 Ciudad de México, Mexico.
| |
Collapse
|
2
|
Winogradoff D, Chou HY, Maffeo C, Aksimentiev A. Percolation transition prescribes protein size-specific barrier to passive transport through the nuclear pore complex. Nat Commun 2022; 13:5138. [PMID: 36050301 PMCID: PMC9437005 DOI: 10.1038/s41467-022-32857-1] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Accepted: 08/18/2022] [Indexed: 11/19/2022] Open
Abstract
Nuclear pore complexes (NPCs) control biomolecular transport in and out of the nucleus. Disordered nucleoporins in the complex's pore form a permeation barrier, preventing unassisted transport of large biomolecules. Here, we combine coarse-grained simulations of experimentally derived NPC structures with a theoretical model to determine the microscopic mechanism of passive transport. Brute-force simulations of protein transport reveal telegraph-like behavior, where prolonged diffusion on one side of the NPC is interrupted by rapid crossings to the other. We rationalize this behavior using a theoretical model that reproduces the energetics and kinetics of permeation solely from statistics of transient voids within the disordered mesh. As the protein size increases, the mesh transforms from a soft to a hard barrier, enabling orders-of-magnitude reduction in permeation rate for proteins beyond the percolation size threshold. Our model enables exploration of alternative NPC architectures and sets the stage for uncovering molecular mechanisms of facilitated nuclear transport.
Collapse
Affiliation(s)
- David Winogradoff
- grid.35403.310000 0004 1936 9991Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA ,grid.35403.310000 0004 1936 9991Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA
| | - Han-Yi Chou
- grid.35403.310000 0004 1936 9991Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA
| | - Christopher Maffeo
- grid.35403.310000 0004 1936 9991Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA ,grid.35403.310000 0004 1936 9991Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA
| | - Aleksei Aksimentiev
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA. .,Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA. .,Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
| |
Collapse
|
3
|
Bezemek Z, Spiliopoulos K. Large deviations for interacting multiscale particle systems. Stoch Process Their Appl 2022. [DOI: 10.1016/j.spa.2022.09.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
4
|
Moncho-Jordá A, Jódar-Reyes AB, Kanduč M, Germán-Bellod A, López-Romero JM, Contreras-Cáceres R, Sarabia F, García-Castro M, Pérez-Ramírez HA, Odriozola G. Scaling Laws in the Diffusive Release of Neutral Cargo from Hollow Hydrogel Nanoparticles: Paclitaxel-Loaded Poly(4-vinylpyridine). ACS NANO 2020; 14:15227-15240. [PMID: 33174725 DOI: 10.1021/acsnano.0c05480] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We study the nonequilibrium diffusive release of electroneutral molecular cargo encapsulated inside hollow hydrogel nanoparticles. We propose a theoretical model that includes osmotic, steric, and short-range polymer-cargo attractions to determine the effective cargo-hydrogel interaction, ueff*, and the effective diffusion coefficient of the cargo inside the polymer network, Deff*. Using dynamical density functional theory (DDFT), we investigate the scaling of the characteristic release time, τ1/2, with the key parameters involved in the process, namely, ueff*, Deff*, and the swelling ratio. This effort represents a full study of the problem, covering a broad range of cargo sizes and providing predictions for repulsive and attractive polymer shells. Our calculations show that the release time through repulsive polymer networks scales with q2eβueff*/Deff* for βueff* ≫ 1. In this case, the cargo molecules are excluded from the shell of the hydrogel. For attractive shells, the polymer retains the cargo molecules on its internal surface and its interior, and the release time grows exponentially with the attraction strength. The DDFT calculations are compared to an analytical model for the mean first passage time, which provides an excellent quantitative description of the kinetics for both repulsive and attractive shells without fitting parameters. Finally, we apply the method to reproduce experimental results on the release of paclitaxel from hollow poly(4-vinylpyridine) nanoparticles and find that the slow release of the drug can be explained in terms of the strong binding attraction between the drug and the polymer.
Collapse
Affiliation(s)
- Arturo Moncho-Jordá
- Departamento de Física Aplicada, Universidad de Granada, Campus Fuentenueva S/N, 18071 Granada, Spain
- Instituto Carlos I de Física Teórica y Computacional, Facultad de Ciencias, Universidad de Granada, Campus Fuentenueva S/N, 18071 Granada, Spain
| | - Ana B Jódar-Reyes
- Departamento de Física Aplicada, Universidad de Granada, Campus Fuentenueva S/N, 18071 Granada, Spain
- Excellence Research Unit "Modeling Nature" (MNat), Universidad de Granada, Campus Fuentenueva S/N, 18071 Granada, Spain
| | - Matej Kanduč
- Jožef Stefan Institute, Jamova 39, SI-1000 Ljubljana, Slovenia
| | - Alicia Germán-Bellod
- Departamento de Física Aplicada, Universidad de Granada, Campus Fuentenueva S/N, 18071 Granada, Spain
| | - Juan M López-Romero
- Departamento de Química Orgánica, Facultad de Ciencias, Universidad de Málaga, 29071 Málaga, Spain
| | - Rafael Contreras-Cáceres
- Departamento de Química en Ciencias Farmacéuticas, Facultad de Farmacia, Universidad Complutense de Madrid, Plaza Ramón y Cajal, 28040 Madrid, Spain
| | - Francisco Sarabia
- Departamento de Química Orgánica, Facultad de Ciencias, Universidad de Málaga, 29071 Málaga, Spain
| | - Miguel García-Castro
- Departamento de Química Orgánica, Facultad de Ciencias, Universidad de Málaga, 29071 Málaga, Spain
| | - Héctor A Pérez-Ramírez
- Física de Procesos Irreversibles, Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana-Azcapotzalco, Avenida San Pablo 180, 02200 Ciudad de México, Mexico
| | - Gerardo Odriozola
- Física de Procesos Irreversibles, Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana-Azcapotzalco, Avenida San Pablo 180, 02200 Ciudad de México, Mexico
| |
Collapse
|
5
|
Sigg D, Voelz VA, Carnevale V. Microcanonical coarse-graining of the kinetic Ising model. J Chem Phys 2020; 152:084104. [PMID: 32113343 PMCID: PMC7042020 DOI: 10.1063/1.5139228] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2019] [Accepted: 02/06/2020] [Indexed: 11/15/2022] Open
Abstract
We propose a scheme for coarse-graining the dynamics of the 2-D kinetic Ising model onto the microcanonical ensemble. At subcritical temperatures, 2-D and higher-dimensional Ising lattices possess two basins of attraction separated by a free energy barrier. Projecting onto the microcanonical ensemble has the advantage that the dependence of the crossing rate constant on environmental conditions can be obtained from a single Monte Carlo trajectory. Using various numerical methods, we computed the forward rate constants of coarse-grained representations of the Ising model and compared them with the true value obtained from brute force simulation. While coarse-graining preserves detailed balance, the computed rate constants for barrier heights between 5 kT and 9 kT were consistently 50% larger than the true value. Markovianity testing revealed loss of dynamical memory, which we propose accounts for coarse-graining error. Committor analysis did not support the alternative hypothesis that microcanonical projection is incompatible with an optimal reaction coordinate. The correct crossing rate constant was obtained by spectrally decomposing the diffusion coefficient near the free energy barrier and selecting the slowest (reactive) component. The spectral method also yielded the correct rate constant in the 3-D Ising lattice, where coarse-graining error was 6% and memory effects were diminished. We conclude that microcanonical coarse-graining supplemented by spectral analysis of short-term barrier fluctuations provides a comprehensive kinetic description of barrier crossing in a non-inertial continuous-time jump process.
Collapse
Affiliation(s)
| | - Vincent A. Voelz
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Vincenzo Carnevale
- Institute for Computational Molecular Science, College of Science and Technology, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
6
|
Li Y, Xu Y, Kurths J. First-passage-time distribution in a moving parabolic potential with spatial roughness. Phys Rev E 2019; 99:052203. [PMID: 31212431 DOI: 10.1103/physreve.99.052203] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2019] [Indexed: 06/09/2023]
Abstract
In this paper, we investigate the first-passage-time distribution (FPTD) within a time-dependent parabolic potential in the presence of roughness with two methods: the Kramers theory and a nonsingular integral equation. By spatially averaging, the rough potential is equivalent to the combination of an effective smooth potential and an effective diffusion coefficient. Based on the Kramers theory, we first obtain Kramers approximations (KAs) of FPTD for both smooth and rough potentials. As expected, KA is valid only for high barriers and small external forces, and generally applicable for high barriers in rough potentials. To overcome the shortcoming of KA, a probability asymptotic approximation (PAA) based on an integral equation is proposed, which uses the transient probability density function (PDF) of the natural boundary conditions instead of the absorbing boundary conditions. We find that PAA fits very well even for large external forces. This enables us to analytically solve the FPTD for large external forces and low barriers as a strong extension to KA. In addition, we show that in the presence of a rough potential, the PAA of FPTD is in good agreement with numerical simulations for low barrier potentials. The PAA makes it possible to investigate the first-passage problem with ultrafast varying potentials and short exiting time. Thus, KA and PAA are complementary in determining the FPTD both for various barriers and external forces. Finally, the mean first-passage time (MFPT) is studied, which illustrates that the PAA of MFPT is effective almost in the whole range of external forces, while the KA of MFPT is valid only for small external forces.
Collapse
Affiliation(s)
- Yongge Li
- Center for Mathematical Sciences & School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Yong Xu
- Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072, China
- MIIT Key Laboratory of Dynamics and Control of Complex Systems, Northwestern Polytechnical University, Xi'an 710072, China
| | - Jürgen Kurths
- Potsdam Institute for Climate Impact Research, Potsdam 14412, Germany
- Human and Animal Physiology Department, Saratov State University, Saratov 410000, Russia
| |
Collapse
|
7
|
Knips A, Zacharias M. Both DNA global deformation and repair enzyme contacts mediate flipping of thymine dimer damage. Sci Rep 2017; 7:41324. [PMID: 28128222 PMCID: PMC5269681 DOI: 10.1038/srep41324] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2016] [Accepted: 12/16/2016] [Indexed: 01/31/2023] Open
Abstract
The photo-induced cis-syn-cyclobutane pyrimidine (CPD) dimer is a frequent DNA lesion. In bacteria photolyases efficiently repair dimers employing a light-driven reaction after flipping out the CPD damage to the active site. How the repair enzyme identifies a damaged site and how the damage is flipped out without external energy is still unclear. Employing molecular dynamics free energy calculations, the CPD flipping process was systematically compared to flipping undamaged nucleotides in various DNA global states and bound to photolyase enzyme. The global DNA deformation alone (without protein) significantly reduces the flipping penalty and induces a partially looped out state of the damage but not undamaged nucleotides. Bound enzyme further lowers the penalty for CPD damage flipping with a lower free energy of the flipped nucleotides in the active site compared to intra-helical state (not for undamaged DNA). Both the reduced penalty and partial looping by global DNA deformation contribute to a significantly shorter mean first passage time for CPD flipping compared to regular nucleotides which increases the repair likelihood upon short time encounter between repair enzyme and DNA.
Collapse
Affiliation(s)
- Alexander Knips
- Physik-Department T38, Technische Universität München, James-Franck-Str. 1, D-85748 Garching, Germany
| | - Martin Zacharias
- Physik-Department T38, Technische Universität München, James-Franck-Str. 1, D-85748 Garching, Germany
| |
Collapse
|
8
|
Li Y, Xu Y, Kurths J, Yue X. Lévy-noise-induced transport in a rough triple-well potential. Phys Rev E 2016; 94:042222. [PMID: 27841518 DOI: 10.1103/physreve.94.042222] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2016] [Indexed: 06/06/2023]
Abstract
Rough energy landscape and noisy environment are two common features in many subjects, such as protein folding. Due to the wide findings of bursting or spiking phenomenon in biology science, small diffusions mixing large jumps are adopted to model the noisy environment that can be properly described by Lévy noise. We combine the Lévy noise with the rough energy landscape, modeled by a potential function superimposed by a fast oscillating function, and study the transport of a particle in a rough triple-well potential excited by Lévy noise, rather than only small perturbations. The probabilities of a particle staying in the middle well are considered under different amplitudes of roughness to find out how roughness affects the steady-state probability density function. Variations in the mean first passage time from the middle well to the right well have been investigated with respect to Lévy parameters and amplitudes of the roughness. In addition, we have examined the influences of roughness on the splitting probabilities of the first escape from the middle well. We uncover that the roughness can enhance significantly the first escape of a particle from the middle well, especially for different skewness parameters, but weak differences are found for stability index and noise intensity on the probabilities a particle staying in the middle well and splitting probability to the right.
Collapse
Affiliation(s)
- Yongge Li
- Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072, China
- Potsdam Institute for Climate Impact Research, Potsdam 14412, Germany
| | - Yong Xu
- Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072, China
| | - Jürgen Kurths
- Potsdam Institute for Climate Impact Research, Potsdam 14412, Germany
- Institute of Applied Physics of the Russian Academy of Sciences, 603950 Nizhny Novgorod, Russia
| | - Xiaole Yue
- Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072, China
| |
Collapse
|
9
|
Roy M, Grazioli G, Andricioaei I. Rate turnover in mechano-catalytic coupling: A model and its microscopic origin. J Chem Phys 2016; 143:045105. [PMID: 26233168 DOI: 10.1063/1.4926664] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
A novel aspect in the area of mechano-chemistry concerns the effect of external forces on enzyme activity, i.e., the existence of mechano-catalytic coupling. Recent experiments on enzyme-catalyzed disulphide bond reduction in proteins under the effect of a force applied on the termini of the protein substrate reveal an unexpected biphasic force dependence for the bond cleavage rate. Here, using atomistic molecular dynamics simulations combined with Smoluchowski theory, we propose a model for this behavior. For a broad range of forces and systems, the model reproduces the experimentally observed rates by solving a reaction-diffusion equation for a "protein coordinate" diffusing in a force-dependent effective potential. The atomistic simulations are used to compute, from first principles, the parameters of the model via a quasiharmonic analysis. Additionally, the simulations are also used to provide details about the microscopic degrees of freedom that are important for the underlying mechano-catalysis.
Collapse
Affiliation(s)
- Mahua Roy
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - Gianmarc Grazioli
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - Ioan Andricioaei
- Department of Chemistry, University of California, Irvine, California 92697, USA
| |
Collapse
|
10
|
A phase transition in the first passage of a Brownian process through a fluctuating boundary with implications for neural coding. Proc Natl Acad Sci U S A 2013; 110:E1438-43. [PMID: 23536302 DOI: 10.1073/pnas.1212479110] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Finding the first time a fluctuating quantity reaches a given boundary is a deceptively simple-looking problem of vast practical importance in physics, biology, chemistry, neuroscience, economics, and industrial engineering. Problems in which the bound to be traversed is itself a fluctuating function of time include widely studied problems in neural coding, such as neuronal integrators with irregular inputs and internal noise. We show that the probability p(t) that a Gauss-Markov process will first exceed the boundary at time t suffers a phase transition as a function of the roughness of the boundary, as measured by its Hölder exponent H. The critical value occurs when the roughness of the boundary equals the roughness of the process, so for diffusive processes the critical value is Hc = 1/2. For smoother boundaries, H > 1/2, the probability density is a continuous function of time. For rougher boundaries, H < 1/2, the probability is concentrated on a Cantor-like set of zero measure: the probability density becomes divergent, almost everywhere either zero or infinity. The critical point Hc = 1/2 corresponds to a widely studied case in the theory of neural coding, in which the external input integrated by a model neuron is a white-noise process, as in the case of uncorrelated but precisely balanced excitatory and inhibitory inputs. We argue that this transition corresponds to a sharp boundary between rate codes, in which the neural firing probability varies smoothly, and temporal codes, in which the neuron fires at sharply defined times regardless of the intensity of internal noise.
Collapse
|
11
|
Dupuis P, Spiliopoulos K. Large deviations for multiscale diffusion via weak convergence methods. Stoch Process Their Appl 2012. [DOI: 10.1016/j.spa.2011.12.006] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
12
|
Sheu SY, Yang DY. Determination of Protein Surface Hydration Shell Free Energy of Water Motion: Theoretical Study and Molecular Dynamics Simulation. J Phys Chem B 2010; 114:16558-66. [DOI: 10.1021/jp105164t] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Sheh-Yi Sheu
- Department of Life Sciences and Institute of Genome Sciences and Institute of Biomedical Informatics, National Yang-Ming University, Taipei 112, Taiwan, Institute of Atomic and Molecular Science, Academia Sinica, Taipei 106, Taiwan, and Division of Biomolecular Sensing, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| | - Dah-Yen Yang
- Department of Life Sciences and Institute of Genome Sciences and Institute of Biomedical Informatics, National Yang-Ming University, Taipei 112, Taiwan, Institute of Atomic and Molecular Science, Academia Sinica, Taipei 106, Taiwan, and Division of Biomolecular Sensing, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| |
Collapse
|
13
|
A multiscale model linking ion-channel molecular dynamics and electrostatics to the cardiac action potential. Proc Natl Acad Sci U S A 2009; 106:11102-6. [PMID: 19549851 DOI: 10.1073/pnas.0904505106] [Citation(s) in RCA: 105] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Ion-channel function is determined by its gating movement. Yet, molecular dynamics and electrophysiological simulations were never combined to link molecular structure to function. We performed multiscale molecular dynamics and continuum electrostatics calculations to simulate a cardiac K(+) channel (I(Ks)) gating and its alteration by mutations that cause arrhythmias and sudden death. An all-atom model of the I(Ks) alpha-subunit KCNQ1, based on the recent Kv1.2 structure, is used to calculate electrostatic energies during gating. Simulations are compared with experiments where varying degrees of positive charge-added via point mutation-progressively reduce current. Whole-cell simulations show that mutations cause action potential and ECG QT interval prolongation, consistent with clinical phenotypes. This framework allows integration of multiscale observations to study the molecular basis of excitation and its alteration by disease.
Collapse
|
14
|
Mondal D, Ghosh PK, Ray DS. Noise-induced transport in a rough ratchet potential. J Chem Phys 2009; 130:074703. [DOI: 10.1063/1.3076934] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Stepanova M. Dynamics of essential collective motions in proteins: theory. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 76:051918. [PMID: 18233698 DOI: 10.1103/physreve.76.051918] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/20/2006] [Revised: 07/23/2007] [Indexed: 05/25/2023]
Abstract
A general theoretical background is introduced for characterization of conformational motions in protein molecules, and for building reduced coarse-grained models of proteins, based on the statistical analysis of their phase trajectories. Using the projection operator technique, a system of coupled generalized Langevin equations is derived for essential collective coordinates, which are generated by principal component analysis of molecular dynamic trajectories. The number of essential degrees of freedom is not limited in the theory. An explicit analytic relation is established between the generalized Langevin equation for essential collective coordinates and that for the all-atom phase trajectory projected onto the subspace of essential collective degrees of freedom. The theory introduced is applied to identify correlated dynamic domains in a macromolecule and to construct coarse-grained models representing the conformational motions in a protein through a few interacting domains embedded in a dissipative medium. A rigorous theoretical background is provided for identification of dynamic correlated domains in a macromolecule. Examples of domain identification in protein G are given and employed to interpret NMR experiments. Challenges and potential outcomes of the theory are discussed.
Collapse
Affiliation(s)
- Maria Stepanova
- National Institute for Nanotechnology, National Research Council of Canada, Department of Electrical and Computer Engineering, University of Alberta, 11421 Saskatchewan Drive, Edmonton, Alberta, Canada T6G 2M9
| |
Collapse
|
16
|
Banushkina P, Meuwly M. Diffusive dynamics on multidimensional rough free energy surfaces. J Chem Phys 2007; 127:135101. [DOI: 10.1063/1.2775444] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
17
|
Danielsson J, Banushkina P, Nutt DR, Meuwly M. Computer simulations of structures, energetics and dynamics of myoglobin ··· ligand complexes. INT REV PHYS CHEM 2006. [DOI: 10.1080/01442350600798253] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
|
18
|
Banushkina P, Meuwly M. Hierarchical Numerical Solution of Smoluchowski Equations with Rough Potentials. J Chem Theory Comput 2005; 1:208-14. [DOI: 10.1021/ct0499480] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
19
|
Affiliation(s)
- Yiqing Shen
- Department of Physics (M/C 273) and Department of Bioengineering (M/C 063), University of Illinois at Chicago, 845 West Taylor Street, Chicago, Illinois 60607
| | - Serguei V. Kuznetsov
- Department of Physics (M/C 273) and Department of Bioengineering (M/C 063), University of Illinois at Chicago, 845 West Taylor Street, Chicago, Illinois 60607
| | - Anjum Ansari
- Department of Physics (M/C 273) and Department of Bioengineering (M/C 063), University of Illinois at Chicago, 845 West Taylor Street, Chicago, Illinois 60607
| |
Collapse
|
20
|
Ansari A, Kuznetsov SV, Shen Y. Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins. Proc Natl Acad Sci U S A 2001; 98:7771-6. [PMID: 11438730 PMCID: PMC35417 DOI: 10.1073/pnas.131477798] [Citation(s) in RCA: 165] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Elucidating the mechanism of folding of polynucleotides depends on accurate estimates of free energy surfaces and a quantitative description of the kinetics of structure formation. Here, the kinetics of hairpin formation in single-stranded DNA are measured after a laser temperature jump. The kinetics are modeled as configurational diffusion on a free energy surface obtained from a statistical mechanical description of equilibrium melting profiles. The effective diffusion coefficient is found to be strongly temperature-dependent in the nucleation step as a result of formation of misfolded loops that do not lead to subsequent zipping. This simple system exhibits many of the features predicted from theoretical studies of protein folding, including a funnel-like energy surface with many folding pathways, trapping in misfolded conformations, and non-Arrhenius folding rates.
Collapse
Affiliation(s)
- A Ansari
- Department of Physics (M/C 273), University of Illinois, 845 West Taylor Street, Chicago, IL 60607, USA.
| | | | | |
Collapse
|