1
|
Jayapurna I, Ruan Z, Eres M, Jalagam P, Jenkins S, Xu T. Sequence Design of Random Heteropolymers as Protein Mimics. Biomacromolecules 2023; 24:652-660. [PMID: 36638823 PMCID: PMC9930114 DOI: 10.1021/acs.biomac.2c01036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Random heteropolymers (RHPs) have been computationally designed and experimentally shown to recapitulate protein-like phase behavior and function. However, unlike proteins, RHP sequences are only statistically defined and cannot be sequenced. Recent developments in reversible-deactivation radical polymerization allowed simulated polymer sequences based on the well-established Mayo-Lewis equation to more accurately reflect ground-truth sequences that are experimentally synthesized. This led to opportunities to perform bioinformatics-inspired analysis on simulated sequences to guide the design, synthesis, and interpretation of RHPs. We compared batches on the order of 10000 simulated RHP sequences that vary by synthetically controllable and measurable RHP characteristics such as chemical heterogeneity and average degree of polymerization. Our analysis spans across 3 levels: segments along a single chain, sequences within a batch, and batch-averaged statistics. We discuss simulator fidelity and highlight the importance of robust segment definition. Examples are presented that demonstrate the use of simulated sequence analysis for in-silico iterative design to mimic protein hydrophobic/hydrophilic segment distributions in RHPs and compare RHP and protein sequence segments to explain experimental results of RHPs that mimic protein function. To facilitate the community use of this workflow, the simulator and analysis modules have been made available through an open source toolkit, the RHPapp.
Collapse
Affiliation(s)
- Ivan Jayapurna
- Department of Materials Science and Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| | - Zhiyuan Ruan
- Department of Materials Science and Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| | - Marco Eres
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Prajna Jalagam
- Department of Materials Science and Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| | - Spencer Jenkins
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Ting Xu
- Department of Materials Science and Engineering, University of California, Berkeley, Berkeley, California 94720, United States.,Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States.,Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Palácio G, Pulcinelli SH, Santilli CV. Fingerprint of semi-crystalline structure memory in the thermal and ionic conduction properties of amorphous ureasil-polyether hybrid solid electrolytes. RSC Adv 2022; 12:5225-5235. [PMID: 35425554 PMCID: PMC8981479 DOI: 10.1039/d1ra09138g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Accepted: 02/04/2022] [Indexed: 11/21/2022] Open
Abstract
Correlations among the structure, thermal properties, and ionic conductivity of solid polymer electrolytes (SPEs) were studied using a ureasil-polyethylene oxide (U-PEO) organic-inorganic hybrid prepared according to a simple sol-gel route, employing a low molecular weight PEO macromer (M w = 1900 g mol-1). The behavior of an amorphous sample loaded with lithium triflate (LiTFSI) at an optimum ratio between ether oxygen and lithium (EO/Li+ = 15) was compared with that of a semicrystalline sample prepared without salt loading. The temperature range investigated by differential scanning calorimetry (DSC), Raman spectroscopy, small angle X-ray scattering (SAXS), and complex impedance spectroscopy covered both the glass transition and the melting temperature of the U-PEO. The gauche to trans conformational transformation of the (O-C-C-O)Li+ sequence showed similarity between the temperature evolution of the semi-crystalline U-PEO and amorphous U-PEO:Li+ samples, providing an indication of the local structural memory of crystalline state in the amorphous SPE. The linear thermal expansion of the average correlation distance between the siloxane crosslink nodes and the long-distance period of the lamellar semi-crystalline edifice were determined by SAXS. Comparison of the expansion curves suggested that although the siloxane nodes were excluded from the PEO crystalline edifice, the sharp expansion of the amorphous region between the lamellae during melting permitted modulation of the free volume of the hybrid network. In addition, the temperature-induced Li+-EO decomplexation observed by Raman spectroscopy explained the change of the average activation energy of the conduction process revealed by the different Arrhenius regimes. These results evidence the key role of the ionic conductivity decoupling from the segmental motion of chain pair channels on the improvement of ion mobility through the free volume between chains. This concept may inspire materials chemistry researchers to design optimized structures of polymer electrolytes with minimized structural memory of crystaline building blocks and improved ionic conductivity.
Collapse
Affiliation(s)
- Gustavo Palácio
- Chemistry Institute of the São Paulo State University, UNESP 14800-060 Araraquara São Paulo Brazil
| | - Sandra H Pulcinelli
- Chemistry Institute of the São Paulo State University, UNESP 14800-060 Araraquara São Paulo Brazil
| | - Celso V Santilli
- Chemistry Institute of the São Paulo State University, UNESP 14800-060 Araraquara São Paulo Brazil
| |
Collapse
|
3
|
Hilburg SL, Ruan Z, Xu T, Alexander-Katz A. Behavior of Protein-Inspired Synthetic Random Heteropolymers. Macromolecules 2020. [DOI: 10.1021/acs.macromol.0c01886] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Shayna L. Hilburg
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Zhiyuan Ruan
- Department of Materials Science & Engineering, University of California Berkeley, Berkeley, California 94720, United States
| | - Ting Xu
- Department of Materials Science & Engineering, University of California Berkeley, Berkeley, California 94720, United States
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, United States
- Tsinghua−Berkeley Shenzhen Institute, University of California Berkeley, Berkeley, California 94720, United States
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Alfredo Alexander-Katz
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
4
|
Weber JK, Pande VS. Percolation-like phase transitions in network models of protein dynamics. J Chem Phys 2016; 142:215105. [PMID: 26049529 DOI: 10.1063/1.4921989] [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/28/2023] Open
Abstract
In broad terms, percolation theory describes the conditions under which clusters of nodes are fully connected in a random network. A percolation phase transition occurs when, as edges are added to a network, its largest connected cluster abruptly jumps from insignificance to complete dominance. In this article, we apply percolation theory to meticulously constructed networks of protein folding dynamics called Markov state models. As rare fluctuations are systematically repressed (or reintroduced), we observe percolation-like phase transitions in protein folding networks: whole sets of conformational states switch from nearly complete isolation to complete connectivity in a rapid fashion. We analyze the general and critical properties of these phase transitions in seven protein systems and discuss how closely dynamics on protein folding landscapes relate to percolation on random lattices.
Collapse
Affiliation(s)
- Jeffrey K Weber
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Vijay S Pande
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
5
|
Nazarov LI, Tamm MV, Avetisov VA, Nechaev SK. A statistical model of intra-chromosome contact maps. SOFT MATTER 2015; 11:1019-1025. [PMID: 25521815 DOI: 10.1039/c4sm02519a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
A statistical model describing a fine structure of the intra-chromosome maps obtained by a genome-wide chromosome conformation capture method (Hi-C) is proposed. The model combines hierarchical chain folding with a quenched heteropolymer structure of primary chromatin sequences. It is conjectured that the observed Hi-C maps are statistical averages over many different ways of hierarchical genome folding. It is shown that the existence of a quenched primary structure coupled with hierarchical folding induces a full range of features observed in experimental Hi-C maps: hierarchical elements, chess-board intermittency and large-scale compartmentalization.
Collapse
Affiliation(s)
- Leonid I Nazarov
- Physics Department, M. V. Lomonosov Moscow State University, 119992 Moscow, Russia
| | | | | | | |
Collapse
|
6
|
Tiana G, Sutto L. Equilibrium properties of realistic random heteropolymers and their relevance for globular and naturally unfolded proteins. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:061910. [PMID: 22304119 DOI: 10.1103/physreve.84.061910] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/05/2011] [Indexed: 05/31/2023]
Abstract
Random heteropolymers do not display the typical equilibrium properties of globular proteins, but are the starting point to understand the physics of proteins and, in particular, to describe their non-native states. So far, they have been studied with mean-field models in the thermodynamic limit, or with computer simulations of very small chains on lattice. After describing a self-adjusting parallel-tempering technique to sample efficiently the low-energy states of frustrated systems without the need of tuning the system-dependent parameters of the algorithm, we apply it to random heteropolymers moving in continuous space. We show that if the mean interaction between monomers is negative, the usual description through the random-energy model is nearly correct, provided that it is extended to account for noncompact conformations. If the mean interaction is positive, such a simple description breaks out and the system behaves in a way more similar to Ising spin glasses. The former case is a model for the denatured state of globular proteins, the latter of naturally unfolded proteins, whose equilibrium properties thus result as qualitatively different.
Collapse
Affiliation(s)
- G Tiana
- Department of Physics, Università degli Studi di Milano and Istituto Nazionale di Fisica Nucleare, via Celoria 16, I-20133 Milano, Italy
| | | |
Collapse
|
7
|
A probabilistic approach to the effect of water hydrogen bonds on the kinetics of protein folding and protein denaturation. Adv Colloid Interface Sci 2010; 154:77-90. [PMID: 20163782 DOI: 10.1016/j.cis.2010.01.009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2009] [Accepted: 01/20/2010] [Indexed: 11/23/2022]
Abstract
Previously, we presented a review of our kinetic models for the nucleation mechanism of protein folding and for the protein thermal denaturation in a barrierless way. A protein was treated as a random heteropolymer consisting of hydrophobic, hydrophilic, and neutral beads. As a crucial idea of the model, an overall potential around the cluster of native residues wherein a residue performs a chaotic motion was considered as the combination of the average dihedral, effective pairwise, and confining potentials. The overall potential as a function of the distance from the cluster has a double well shape. This allowed one to develop kinetic models for the nucleation mechanism of protein folding (NMPF) and barrierless protein denaturation (BPD) by using the mean first passage time analysis. In the original models, however, hydrogen bonding effects were taken into account only indirectly which affected the accuracy of the models because hydrogen bonding does play a crucial role in the folding, stability, and denaturation of proteins. To improve the NMPF and BPD models and explicitly take into account the hydrogen bonding "water-water" and "water-protein residue", we have developed a probabilistic hydrogen bond (PHB) model for the effect of hydrogen bond networks of water molecules around two solute particles (immersed in water) on their interaction, and have then combined the PHB model with the NMPF and BPD models. In this paper, that can be regarded as sequel of our previous review, we analyze the modified NMPF and BPD models that explicitly take into account the effect of water-water hydrogen bonding on these processes. As expected, the application of the modified models to the folding/unfolding of two model proteins (one short, consisting of 124 residues and the other large, consisting of 2500 residues) demonstrate that the hydrogen bond networks play a very important role in the protein folding/unfolding phenomena.
Collapse
|
8
|
Djikaev YS, Ruckenstein E. Effect of hydrogen bond networks on the nucleation mechanism of protein folding. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:061918. [PMID: 20365201 DOI: 10.1103/physreve.80.061918] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2009] [Revised: 10/23/2009] [Indexed: 05/29/2023]
Abstract
We have recently developed a kinetic model for the nucleation mechanism of protein folding (NMPF) in terms of ternary nucleation by using the first passage time analysis. A protein was considered as a random heteropolymer consisting of hydrophobic, hydrophilic (some of which are negatively or positively ionizable), and neutral beads. The main idea of the NMPF model consisted of averaging the dihedral potential in which a selected residue is involved over all possible configurations of all neighboring residues along the protein chain. The combination of the average dihedral, effective pairwise (due to Lennard-Jones-type and electrostatic interactions), and confining (due to the polymer connectivity constraint) potentials gives rise to an overall potential around the cluster that, as a function of the distance from the cluster center, has a double-well shape. This allows one to evaluate the protein folding time. In the original NMPF model hydrogen bonding was not taken into account explicitly. To improve the NMPF model and make it more realistic, in this paper we modify our (previously developed) probabilistic hydrogen bond model and combine it with the former. Thus, a contribution due to the disruption of hydrogen bond networks around the interacting particles (cluster of native residues and residue in the protein unfolded part) appears in the overall potential field around a cluster. The modified model is applied to the folding of the same model proteins that were examined in the original model: a short protein consisting of 124 residues (roughly mimicking bovine pancreatic ribonuclease) and a long one consisting of 2500 residues (as a representative of large proteins with superlong polypeptide chains), at pH=8.3 , 7.3, and 6.3. The hydrogen bond contribution now plays a dominant role in the total potential field around the cluster (except for very short distances thereto where the repulsive energy tends to infinity). It is by an order of magnitude stronger for hydrophobic residues than for hydrophilic ones. The range of "residue-cluster" distances, at which the hydrogen bond effect exists, is twice as long for hydrophobic residues as for hydrophilic ones.
Collapse
Affiliation(s)
- Y S Djikaev
- Department of Chemical and Biological Engineering, SUNY at Buffalo, Buffalo, New York 14260, USA.
| | | |
Collapse
|
9
|
Djikaev YS, Ruckenstein E. The role of hydrogen bond networks in the barrierless thermal denaturation of a native protein. J Chem Phys 2009; 131:045105. [DOI: 10.1063/1.3179688] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Djikaev YS, Ruckenstein E. Thermal denaturation of a native protein via spinodal decomposition in the framework of first-passage-time analysis. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 78:011909. [PMID: 18763984 DOI: 10.1103/physreve.78.011909] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/10/2008] [Indexed: 05/26/2023]
Abstract
In this paper we present a kinetic model for the thermal unfolding of a native protein. Due to a sufficiently large temperature increase or decrease, the rate with which a cluster of native residues within a protein emits residues becomes larger than the rate of absorption of residues from the unfolded part of the protein in the whole range of cluster sizes up to the size of the whole protein. This leads to the unfolding of the protein in a barrierless way, i.e., as spinodal decomposition. Using the formalism of the first passage time analysis [previously applied also to the problem of protein folding via nucleation by the authors, J. Phys. Chem. B 111, 886 (2007); J. Chem. Phys. 126, 175103 (2007)], one can determine the temperature dependence of the rates of emission and absorption of residues by the cluster. Knowing these rates as functions of temperature and cluster size, one can find the threshold temperatures of cold and hot barrierless denaturation as well as the unfolding times at temperatures lower and higher, respectively, than those threshold values. For a numerical illustration, the method is applied to the thermal unfolding of a model protein consisting of 2500 residues.
Collapse
Affiliation(s)
- Y S Djikaev
- Department of Chemical and Biological Engineering, SUNY at Buffalo, Buffalo, New York 14260, USA.
| | | |
Collapse
|
11
|
Suzuki Y, Noel JK, Onuchic JN. An analytical study of the interplay between geometrical and energetic effects in protein folding. J Chem Phys 2008; 128:025101. [PMID: 18205476 DOI: 10.1063/1.2812956] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Analytical studies have several advantages for an understanding of the mechanisms of protein folding such as the interplay between geometrical and energetic effects. In this paper, we introduce a Gaussian filament with a C(alpha) structure-based (Go) potential as a new theoretical scheme based on a Hamiltonian approach. This model takes into account geometrical information in a realistic fashion without the need of phenomenological descriptions. In order to make this model more appropriate for comparison with protein folding simulations and experiments, we introduce a many-body interaction into the potential term to enhance cooperativity. We apply our new analytical model to a beta-hairpin-type peptide and compare our results with a molecular dynamics simulation of a structure-based model.
Collapse
Affiliation(s)
- Yoko Suzuki
- Department of Physics, School of Sciences and Engineering, Meisei University, 2-1-1 Hodokubo, Hino-shi, Tokyo 191-8506, Japan.
| | | | | |
Collapse
|
12
|
Djikaev YS, Ruckenstein E. Effect of ionized protein residues on the nucleation pathway of protein folding. J Chem Phys 2008; 128:025103. [DOI: 10.1063/1.2820771] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Djikaev YS, Ruckenstein E. Temperature effects on the nucleation mechanism of protein folding and on the barrierless thermal denaturation of a native protein. Phys Chem Chem Phys 2008; 10:6281-300. [DOI: 10.1039/b807399f] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
14
|
Djikaev YS, Ruckenstein E. A ternary nucleation model for the nucleation pathway of protein folding. J Chem Phys 2007; 126:175103. [PMID: 17492888 DOI: 10.1063/1.2727469] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recently [Y. S. Djikaev and E. Ruckenstein, J. Phys. Chem. B 111, 886 (2007)], the authors proposed a kinetic model for the nucleation mechanism of protein folding where a protein was modeled as a heteropolymer consisting of hydrophobic and hydrophilic beads and the composition of the growing cluster of protein residues was assumed to be constant and equal to the overall protein composition. Here, they further develop the model by considering a protein as a three-component heteropolymer and by allowing the composition of the growing cluster of protein residues to vary independently of the overall one. All the bonds in the heteropolymer (now consisting of hydrophobic, hydrophilic, and neutral beads) have the same constant length, and all the bond angles are equal and fixed. As a crucial idea of the model, an overall potential around the cluster wherein a residue performs a chaotic motion is considered to be a combination of the average dihedral and average pairwise potentials assigned to the bead. The overall potential as a function of the distance from the cluster center has a double well shape which allows one to determine its emission and absorption rates by using a first passage time analysis. Knowing these rates as functions of three independent variables of a ternary cluster, one can develop a self-consistent kinetic theory for the nucleation mechanism of folding of a protein using a ternary nucleation formalism and evaluate the size and composition of the nucleus and the protein folding time. As an illustration, the model is applied to the folding of bovine pancreatic ribonuclease consisting of 124 amino acids whereof 40 are hydrophobic, 81 hydrophilic, and 3 neutral. With a reasonable choice of diffusion coefficients of the residues in the native state and potential parameters, the model predicts folding times in the range of 1-100 s.
Collapse
Affiliation(s)
- Y S Djikaev
- Department of Chemical and Biological Engineering, SUNY at Buffalo, Buffalo, New York 14260, USA.
| | | |
Collapse
|
15
|
Abstract
A nucleation-like pathway of protein folding involves the formation of a cluster containing native residues that grows by including residues from the unfolded part of the protein. This pathway is examined by using a heteropolymer as a protein model. The model heteropolymer consists of hydrophobic and hydrophilic beads with fixed bond lengths and bond angles. The total energy of the heteropolymer is determined by the pairwise repulsive/attractive interactions between nonlinked beads and by the contribution from the dihedral angles involved. The parameters of these interactions can be rigorously defined, unlike the ill-defined surface tension of a cluster of protein residues that constitutes the basis of a previous nucleation model. The main idea underlying the new model consists of averaging the dihedral potential of a selected residue over all possible configurations of all neighboring residues along the protein chain. The resulting average dihedral potential depends on the distance between the selected residue and the cluster center. Its combination with the average pairwise potential of the selected residue and with a confining potential caused by the bonds between the residues leads to an overall potential around the cluster that has a double-well shape. Residues in the inner (closer to the cluster) well are considered as belonging to the folded cluster, whereas those in the outer well are treated as belonging to the unfolded part of the protein. Transitions of residues from the inner well into the outer one and vice versa are considered as elementary emission and absorption events, respectively. The double-well character of the potential well around the cluster allows one to determine the rates of both emission and absorption of residues by the cluster using a first passage time analysis. Once these rates are found as functions of the cluster size, one can develop a self-consistent kinetic theory for the nucleation mechanism of folding of a protein. The model allows one to evaluate the size of the nucleus and the protein folding time. The latter is evaluated as the sum of the times necessary for the first nucleation event to occur and for the nucleus to grow to the maximum size (of the folded protein). Depending on the diffusion coefficients of the native residues in the range from 10(-6) to 10(-8) cm2/s, numerical calculations for a protein of 2500 residues suggest that the folding time ranges from several seconds to several hundreds of seconds.
Collapse
Affiliation(s)
- Y S Djikaev
- Department of Chemical and Biological Engineering, State University of New York at Buffalo, Buffalo, New York 14260, USA.
| | | |
Collapse
|
16
|
Goldschmidt YY. Large-time dynamics and aging of a polymer chain in a random potential. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:021804. [PMID: 17025463 DOI: 10.1103/physreve.74.021804] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2006] [Indexed: 05/12/2023]
Abstract
We study the out-of-equilibrium large-time dynamics of a Gaussian polymer chain in a quenched random potential. The dynamics studied is a simple Langevin dynamics commonly referred to as the Rouse model. The equations for the two-time correlation and response functions are derived within the Gaussian variational approximation. In order to implement this approximation faithfully, we employ the supersymmetric representation of the Martin-Siggia-Rose dynamical action. For a short-ranged correlated random potential the equations are solved analytically in the limit of large times using certain assumptions concerning the asymptotic behavior. Two possible dynamical behaviors are identified depending upon the time separation: a stationary regime and an aging regime. In the stationary regime time translation invariance holds and so does the fluctuation dissipation theorem. The aging regime which occurs for large time separations of the two-time correlation functions is characterized by a history dependence and the breakdown of certain equilibrium relations. The large-time limit of the equations yields equations among the order parameters that are similar to the equations obtained in statics using replicas. In particular the aging solution corresponds to the broken replica solution. But there is a difference in one equation that leads to important consequences for the solution. The stationary regime corresponds to the motion of the polymer inside a local minimum of the random potential, whereas in the aging regime the polymer hops between different minima. As a by-product we also solve exactly the dynamics of a chain in a random potential with quadratic correlations.
Collapse
Affiliation(s)
- Yadin Y Goldschmidt
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| |
Collapse
|
17
|
Shakhnovich E. Protein folding thermodynamics and dynamics: where physics, chemistry, and biology meet. Chem Rev 2006; 106:1559-88. [PMID: 16683745 PMCID: PMC2735084 DOI: 10.1021/cr040425u] [Citation(s) in RCA: 253] [Impact Index Per Article: 14.1] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Eugene Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford Street, Cambridge, Massachusetts 02138, USA.
| |
Collapse
|
18
|
Konkoli Z, Hertz J. Embedding a native state into a random heteropolymer model: the dynamic approach. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 67:051915. [PMID: 12786186 DOI: 10.1103/physreve.67.051915] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2002] [Indexed: 05/24/2023]
Abstract
We study a random heteropolymer model with Langevin dynamics, in the supersymmetric formulation. Employing a procedure similar to one that has been used in static calculations, we construct an ensemble in which the affinity of the system for a native state is controlled by a "selection temperature" T0. In the limit of high T0, the model reduces to a random heteropolymer, while for T0-->0 the system is forced into the native state. Within the Gaussian variational approach that we employed previously for the random heteropolymer, we explore the phases of the system for high and low T0. For high T0, the system exhibits a (dynamical) spin-glass phase, like that found for the random heteropolymer, below a temperature T(g). For low T0, we find an ordered phase, characterized by a nonzero overlap with the native state, below a temperature T(n) proportional to 1/T(0)>T(g). However, the random-globule phase remains locally stable below T(n), down to the dynamical glass transition at T(g). Thus, in this model, folding is rapid for temperatures between T(g) and T(n), but below T(g) the system can get trapped in conformations uncorrelated with the native state. At a lower temperature, the ordered phase can also undergo a dynamical glass transition, splitting into substates separated by large barriers.
Collapse
Affiliation(s)
- Z Konkoli
- Department of Applied Physics, Chalmers University of Technology and Göteborg University, SE 412 96 Göteborg, Sweden
| | | |
Collapse
|
19
|
Konkoli Z, Hertz J, Franz S. Random heteropolymer dynamics. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2001; 64:051910. [PMID: 11735971 DOI: 10.1103/physreve.64.051910] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/25/2001] [Indexed: 05/23/2023]
Abstract
We study the Langevin dynamics of the standard random heteropolymer model by mapping the problem to a supersymmetric field theory using the Martin-Siggia-Rose formalism. The resulting model is solved nonperturbatively employing a Gaussian variational approach. In constructing the solution, we assume that the chain is very long and impose the translational invariance which is expected to be present in the bulk of the globule by averaging over the center of mass coordinate. In this way we derive equations of motion for the correlation and response functions C(t,t') and R(t,t'). The order parameters are extracted from the asymptotic behavior of these functions. We find a dynamical phase diagram with frozen (glassy) and melted (ergodic) phases. In the glassy phase the system fails to reach equilibrium and exhibits aging of the type found in p-spin glasses. Within the approximations used in this study, the random heteropolymer model can be mapped to the problem of a manifold in a random potential with power law correlations.
Collapse
Affiliation(s)
- Z Konkoli
- NORDITA, Blegdamsvej 17, DK 2100 København, Denmark
| | | | | |
Collapse
|
20
|
|
21
|
Shiferaw Y, Goldschmidt YY. Localization of a polymer in random media: relation to the localization of a quantum particle. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2001; 63:051803. [PMID: 11414926 DOI: 10.1103/physreve.63.051803] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2000] [Indexed: 05/23/2023]
Abstract
In this paper we consider in detail the connection between the problem of a polymer in a random medium and that of a quantum particle in a random potential. We are interested in a system of finite volume where the polymer is known to be localized inside a low minimum of the potential. We show how the end-to-end distance of a polymer that is free to move can be obtained from the density of states of the quantum particle using extreme value statistics. We give a physical interpretation to the recently discovered one-step replica-symmetry-breaking solution for the polymer [Phys. Rev. E 61, 1729 (2000)] in terms of the statistics of localized tail states. Numerical solutions of the variational equations for chains of different length are performed and compared with quenched averages computed directly by using the eigenfunctions and eigenenergies of the Schrödinger equation for a particle in a one-dimensional random potential. The quantities investigated are the radius of gyration of a free Gaussian chain, its mean square distance from the origin and the end-to-end distance of a tethered chain. The probability distribution for the position of the chain is also investigated. The glassiness of the system is explained and is estimated from the variance of the measured quantities.
Collapse
Affiliation(s)
- Y Shiferaw
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennyslvania 15260, USA
| | | |
Collapse
|
22
|
Pitard E, Shakhnovich EI. Mode-coupling theory for heteropolymers. PHYSICAL REVIEW E 2001; 63:041501. [PMID: 11308842 DOI: 10.1103/physreve.63.041501] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/1999] [Revised: 03/31/2000] [Indexed: 11/07/2022]
Abstract
We study the Langevin dynamics of a heteropolymer by means of a mode-coupling approximation scheme, giving rise to a set of coupled integro-differential equations relating the response and correlation functions. The analysis shows that there is a regime at low temperature characterized by out-of-equilibrium dynamics, with violation of time-translational invariance and of the fluctuation-dissipation theorem. The onset of aging dynamics at low temperatures gives insight into the nature of the slow dynamics of a disordered polymer. We also introduce a renormalization-group treatment of our mode-coupling equations, which supports our analysis, and might be applicable to other systems.
Collapse
Affiliation(s)
- E Pitard
- Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford Street, Cambridge, MA 02138, USA
| | | |
Collapse
|
23
|
Schulz M. Thermodynamics and mechanical properties of random copolymer networks. J Chem Phys 2000. [DOI: 10.1063/1.1324712] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Enoki T, Tanaka K, Watanabe T, Oya T, Sakiyama T, Takeoka Y, Ito K, Wang G, Annaka M, Hara K, Du R, Chuang J, Wasserman K, Masamune S, Tanaka T. Frustrations in polymer conformation in gels and their minimization through molecular imprinting. PHYSICAL REVIEW LETTERS 2000; 85:5000-5003. [PMID: 11102171 DOI: 10.1103/physrevlett.85.5000] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2000] [Indexed: 05/23/2023]
Abstract
We report an experimental realization of a gel system in which frustrations exist and can be minimized, thus meeting two crucial criteria predicted to enable memory of conformations in polymers. The gels consist of a thermosensitive major monomer component and two minor components. One minor component is positively charged and will form complexes around negatively charged target molecules placed in solution. The complexes can be imprinted into the gel by then cross-linking the second minor component, which will form cross-links additional to those in the major polymer matrix. The complexes are destroyed and reformed upon swelling and reshrinking of the gels, showing that memorization has been achieved.
Collapse
Affiliation(s)
- T Enoki
- Department of Physics and Center for Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
25
|
Wilder J, Shakhnovich EI. Proteins with selected sequences: a heteropolymeric study. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 2000; 62:7100-10. [PMID: 11102067 DOI: 10.1103/physreve.62.7100] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/05/2000] [Indexed: 11/07/2022]
Abstract
Protein sequences are expected not to be random but selected in order to form a stable native structure that is kinetically accessible. Therefore our model contains a selective temperature in sequence space (see [S. Ramanathan and E. Shakhnovich, Phys. Rev. E 50, 1303 (1994)] ) to optimize the sequence for the target conformation statistically. Replica calculations, which go beyond quadratic approximations in the field-theoretical Hamiltonian, are presented. A phase diagram indicating the temperatures and selective temperatures at which transitions to a frozen globule, i.e., the native state, occur is obtained. It is shown that going beyond the quadratic approximation in the field Hamiltonian is very important, since it results in a significant change of the phase diagram. Moreover, we suggest that a one-step replica permutation symmetry scheme is sufficient to solve the model. In addition to this we present a result for the sequence correlation function along the chain in the case of a short-ranged potential between the monomers. A correlation function between monomers that form a contact in the native state is given depending on the temperature and the interaction parameter.
Collapse
Affiliation(s)
- J Wilder
- Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford Street, Cambridge, Massachusetts 02138, USA.
| | | |
Collapse
|
26
|
Alvarez-Lorenzo C, Guney O, Oya T, Sakai Y, Kobayashi M, Enoki T, Takeoka Y, Ishibashi T, Kuroda K, Tanaka K, Wang G, Grosberg AY, Masamune S, Tanaka T. Polymer Gels That Memorize Elements of Molecular Conformation. Macromolecules 2000. [DOI: 10.1021/ma000603v] [Citation(s) in RCA: 103] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Carmen Alvarez-Lorenzo
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Orhan Guney
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Taro Oya
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Yasuzo Sakai
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Masatoshi Kobayashi
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Takashi Enoki
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Yukikazu Takeoka
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Toru Ishibashi
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Kenichi Kuroda
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Kazunori Tanaka
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Guoqiang Wang
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Alexander Yu. Grosberg
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Satoru Masamune
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| | - Toyoichi Tanaka
- Department of Physics and Center for Materials Science and Engineering and Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139
| |
Collapse
|
27
|
Rostiashvili VG, Vilgis TA. Dynamics of structural models with a long-range interaction: glassy versus nonglassy behavior. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 2000; 62:1560-1576. [PMID: 11088618 DOI: 10.1103/physreve.62.1560] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/06/2000] [Indexed: 05/23/2023]
Abstract
By making use of the Langevin dynamics and its generating functional (GF) formulation, the influence of the long-range nature of the interaction on the tendency of the glass formation is systematically investigated. In doing so, two types of models are considered: (i) the nondisordered model with a pure repulsive type of interaction, and (ii) the model with a randomly distributed strength of interaction (a quenched disordered model). The long-ranged potential of interaction is scaled with a number of particles N in such a way as to enable for the GF the saddle-point treatment as well as the systematic 1/N expansion around it. We show that the nondisordered model has no glass transition, which is in line with the mean-field limit of the mode-coupling theory (MCT) predictions. On the other hand, the model with a long-range interaction that has a quenched disorder leads to MC equations which are generic for the p-spin glass model and polymeric manifold in a random media.
Collapse
Affiliation(s)
- VG Rostiashvili
- Max-Planck-Institut fur Polymerforschung, Postfach 3148, D-55021 Mainz, Germany and Institute of Chemical Physics, Russian Academy of Science, 142432, Chernogolovka, Moscow region, Russia
| | | |
Collapse
|
28
|
Timoshenko EG, Kuznetsov YA. Analysis of stability of macromolecular clusters in dilute heteropolymer solutions. J Chem Phys 2000. [DOI: 10.1063/1.481417] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
29
|
Du R, Grosberg AY, Tanaka T, Rubinstein M. Unexpected scenario of glass transition in polymer globules: an exactly enumerable model. PHYSICAL REVIEW LETTERS 2000; 84:2417-2420. [PMID: 11018899 DOI: 10.1103/physrevlett.84.2417] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/1999] [Indexed: 05/23/2023]
Abstract
We introduce a lattice model of glass transition in polymer globules. This model exhibits ergodicity breaking in which the disjoint regions of phase space do not arise uniformly, but as small chambers whose number increases exponentially with polymer density. Chamber sizes obey power law distribution, making phase space similar to a fractal foam. This clearly demonstrates the importance of the phase space geometry and topology in describing any glass-forming system, such as semicompact polymers during protein folding.
Collapse
Affiliation(s)
- R Du
- Department of Physics and Center for Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | | | | | | |
Collapse
|
30
|
Buchler NEG, Goldstein RA. Surveying determinants of protein structure designability across different energy models and amino-acid alphabets: A consensus. J Chem Phys 2000. [DOI: 10.1063/1.480893] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
Buchler NEG, Goldstein RA. Universal correlation between energy gap and foldability for the random energy model and lattice proteins. J Chem Phys 1999. [DOI: 10.1063/1.479951] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Angerman HJ, Shakhnovich E. Freezing in polyampholytes globules: Influence of the long-range nature of the interaction. J Chem Phys 1999. [DOI: 10.1063/1.479356] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Iori G, Marinari E, Parisi G. Random self-interacting chains: a mechanism for protein folding. ACTA ACUST UNITED AC 1999. [DOI: 10.1088/0305-4470/24/22/019] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
34
|
Gutin AM, Grosberg AY, Shakhnovich EI. Globular state of branched random heteropolymers. ACTA ACUST UNITED AC 1999. [DOI: 10.1088/0305-4470/26/5/023] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
35
|
Struglia MV. Conformational properties of random heteropolymers in the folded phase. ACTA ACUST UNITED AC 1999. [DOI: 10.1088/0305-4470/28/6/005] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
36
|
|
37
|
|
38
|
Bouchaud JP, Mezard M, Parisi G, Yedidia JS. Polymers with long-range self-repulsion: a variational approach. ACTA ACUST UNITED AC 1999. [DOI: 10.1088/0305-4470/24/17/011] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
39
|
|
40
|
Derrida B, Higgs PG. Low-temperature properties of directed walks with random self interactions. ACTA ACUST UNITED AC 1999. [DOI: 10.1088/0305-4470/27/16/014] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
41
|
Goldbart PM, Zippelius A. Issues of replica-symmetry breaking for the amorphous solid state of vulcanized macromolecules. ACTA ACUST UNITED AC 1999. [DOI: 10.1088/0305-4470/27/19/011] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
42
|
|
43
|
Simmons ER, Chakraborty AK. Adhesion promotion at a homopolymer–solid interface using random heteropolymers. J Chem Phys 1998. [DOI: 10.1063/1.477534] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Srebnik S, Chakraborty AK, Bratko D. Random heteropolymer adsorption on disordered multifunctional surfaces: Effect of specific intersegment interactions. J Chem Phys 1998. [DOI: 10.1063/1.477285] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Gutin A, Sali A, Abkevich V, Karplus M, Shakhnovich EI. Temperature dependence of the folding rate in a simple protein model: Search for a “glass” transition. J Chem Phys 1998. [DOI: 10.1063/1.476053] [Citation(s) in RCA: 69] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
46
|
Bratko D, Chakraborty AK, Shakhnovich E. Adsorption of random copolymers on disordered surfaces. ACTA ACUST UNITED AC 1998. [DOI: 10.1016/s1089-3156(98)00023-3] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
47
|
Takada S, Wolynes PG. Microscopic theory of critical folding nuclei and reconfiguration activation barriers in folding proteins. J Chem Phys 1997. [DOI: 10.1063/1.475256] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
48
|
Affiliation(s)
- Andrey V. Dobrynin
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290
| |
Collapse
|
49
|
Bratko D, Chakraborty AK, Shakhnovich EI. The structure of a random heteropolymer in a disordered medium: Ensemble growth simulation. J Chem Phys 1997. [DOI: 10.1063/1.473223] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Roan JR, Shakhnovich EI. Dynamics of heteropolymers in dilute solution: Effective equation of motion and relaxation spectrum. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1996; 54:5340-5357. [PMID: 9965719 DOI: 10.1103/physreve.54.5340] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|