A network of conformational transitions in an unfolding process of HP-35 revealed by high-temperature MD simulation and a Markov state model

An understanding of protein folding/unfolding processes has important implications for all biological processes, including protein degradation, protein translocation, aging, and diseases. All-atom molecular dynamics(MD) simulations are uniquely suitable for it because of their atomic level resolution and accuracy. However, limited by computational capabilities, nowadays even for small and fast-folding proteins, all-atom MD simulations of protein folding still presents a great challenge. An alternative way is to study unfolding process using MD simulations at high temperature. High temperature provides more energy to overcome energetic barriers to unfolding, and information obtained from studying unfolding can shed light on the mechanism of folding. In the present study, a 1000-ns MD simulation at high temperature(500 K)was performed to investigate the unfolding process of a small protein, chicken villin headpiece(HP-35). To infer the folding mechanism, a Markov state model was also built from our simulation, which maps out six macrostates during the folding/unfolding process as well as critical transitions between them, revealing the folding mechanism unambiguously.     

Folding and unfolding of an elastinlike oligopeptide: "inverse temperature transition," reentrance, and hydrogen-bond dynamics

The temperature-dependent behavior of a solvated oligopeptide, GVG(VPGVG), is investigated. Spectroscopic measurements, thermodynamic measurements, and molecular dynamics simulations find that this elastinlike octapeptide behaves as a two-state system that undergoes an "inverse temperature" folding transition and reentrant unfolding close to the boiling point of water. A molecular picture of these processes is presented, emphasizing changes in the dynamics of hydrogen bonding at the protein/water interface and peptide backbone librational entropy.     

基于TorchMD的粗粒化分子动力模拟研究

粗粒化模型通过简化原子性质以及原子间的相互作用实现生物大分子长时间尺度的分子动力学模拟. 深度学习通过模拟人类的认知过程实现海量数据的准确分类和回归过程. 本论文将这两种技术进行融合,利用基于深度学习的粗粒化分子动力学模拟技术研究分子在不同状态之间的变化过程,并提出基于TorchMD的分子动力学模拟的分析框架. 在本工作中,MFDP聚类算法被用于在三维的CV变量空间中进行聚类,并确定分子的若干主要状态,在完成聚类的同时,给出各类中的代表分子构象,并给出类之间的分子构象. 这为后续利用String算法分析分子在不同状态间的转换路径打下基础. 通过String算法,迭代搜索得到分子在不同状态之间的变化路径以及对应的势能变化曲线. 通过与已有文献的结果进行对比,验证了基于TorchMD的粗粒化分子动力学模拟的理论框架可以在相对较短的时间尺度里研究分子的变化过程.     

Nanosecond-time-resolved infrared spectroscopic study of fast relaxation kinetics of protein folding by means of laser-induced temperature-jump

Elucidating the initial kinetics of folding pathways is critical to the understanding of the protein folding mechanism. Transient infrared spectroscopy has proved a powerful tool to probe the folding kinetics. Herein we report the construction of a nanosecond laser-induced temperature-jump （T-jump） technique coupled to a nanosecond timeresolved transient mid-infrared （mid-IR） spectrometer system capable of investigating the protein folding kinetics with a temporal resolution of 50 ns after deconvolution of the instrumental response function. The mid-IR source is a liquid N2 cooled CO laser covering a spectral range of 5.0μm （2000 cm^-1）-6.5μm （1540 cm^-1）. The heating pulse was generated by a high pressure H2 Raman shifter at wavelength of 1.9μm. The maximum temperature-jump could reach as high as 26±1℃. The fast folding/unfolding dynamics of cytochrome C was investigated by the constructed system, providing an example.     

Protein folding kinetics and thermodynamics from atomistic simulations

Determining protein folding kinetics and thermodynamics from all-atom molecular dynamics (MD) simulations without using experimental data represents a formidable scientific challenge because simulations can easily get trapped in local minima on rough free energy landscapes. This necessitates the computation of multiple simulation trajectories, which can be independent from each other or coupled in some manner, as, for example, in the replica exchange MD method. Here we present results obtained with a new analysis tool that allows the deduction of faithful kinetics data from a heterogeneous ensemble of simulation trajectories. The method is demonstrated on the decapeptide Chignolin for which we predict folding and unfolding time constants of 1.0 +/- 0.3 and 2.6 +/- 0.4 micros, respectively. We also derive the energetics of folding, and calculate a realistic melting curve for Chignolin.     

蛋白质空间结构的实验技术和理论方法

文章主要介绍几种蛋白质空间结构的实验测定方法，在现代生物学研究中，最常用的方法包括X射线晶体学、二维核磁共振（2D-NMR）和低温冷冻电镜，近几年发展起来的单分子技术在生物大分子动态结构的研究中应用越来越多，这些方法都有它们特定的时间和空间分辨率，所测定的结构及其动力学受环境热运动涨落的影响也非常不同，文章对这些问题作了较详细的分析，在蛋白质结构的理论方法方面，介绍了一个新的折叠理论及其与现有折叠模型的关系．讨论了模拟计算在研究蛋白质构象变化和动力学方面的应用，同时强调了分子动力学和蒙特卡罗方法．指出粗粒化模型是研究的热点之一，对生物学中经常遇到的多长度多时问尺度问题提供了一个可行的解决方案。     

A challenge to chaotic itinerancy from brain dynamics

Brain hermeneutics and chaotic itinerancy proposed by Tsuda are attractive characterizations of perceptual dynamics in the mammalian olfactory system. This theory proposes that perception occurs at the interface between itinerant neural representation and interaction with the environment. Quantifiable application of these dynamics has been hampered by the lack of definable history and action processes which characterize the changes induced by behavioral state, attention, and learning. Local field potentials measured from several brain areas were used to characterize dynamic activity patterns for their use as representations of history and action processes. The signals were recorded from olfactory areas (olfactory bulb, OB, and pyriform cortex) and hippocampal areas (entorhinal cortex and dentate gyrus, DG) in the brains of rats. During odor-guided behavior the system shows dynamics at three temporal scales. Short time-scale changes are system-wide and can occur in the space of a single sniff. They are predictable, associated with learned shifts in behavioral state and occur periodically on the scale of the intertrial interval. These changes occupy the theta (2-12 Hz), beta (15-30 Hz), and gamma (40-100 Hz) frequency bands within and between all areas. Medium time-scale changes occur relatively unpredictably, manifesting in these data as alterations in connection strength between the OB and DG. These changes are strongly correlated with performance in associated trial blocks (5-10 min) and may be due to fluctuations in attention, mood, or amount of reward received. Long time-scale changes are likely related to learning or decline due to aging or disease. These may be modeled as slow monotonic processes that occur within or across days or even weeks or years. The folding of different time scales is proposed as a mechanism for chaotic itinerancy, represented by dynamic processes instead of static connection strengths. Thus, the individual maintains continuity of experience within the stability of fast periodic and slow monotonic processes, while medium scale events alter experience and performance dramatically but temporarily. These processes together with as yet to be determined action effects from motor system feedback are proposed as an instantiation of brain hermeneutics and chaotic itinerancy.     

Folding of lattice protein model chains with fixed ends

Using Monte Carlo simulations, we have studied the folding dynamics and thermodynamics of geometricall yconstrained lattice protein model chains. The constraints are realized by fixing one or both terminals of the chains. By comparing the results with that of the free-end chains, we find that the folding behaviours of the end-constralned chains are not completely similar to that of the free-end chains. Both kinds of constraints on the chain ends affect the folding dynamics of the chains: i.e., the folding rate, but not the thermodynamics. The thermodynamic behaviour of the one-end-fixed chains shows less difference from that of the free-end chains, while the thermodynamic behaviour of the two-end-fixed chains has obvious difference from that of the free-end chains. The origin of these differences comes from the differences of the ergodicitv of the chains in the conformational space.     

飞秒物理、飞秒化学和飞秒生物学

飞秒激光技术因其极高的时间分辨特性而被广泛应用于研究多种材料的超快过程,文章从几个侧面就飞秒技术在物理学,化学及生物学等方面的应用作了介绍,在飞秒物理方面,介绍了飞秒技术在研究半导体量子阱材料,纳米材料的性质及高次谐波产生等方面的研究进展,飞秒化学则主要介绍了飞秒技术在研究光化学反应,光解离过程、键的断裂及结合以及相关的动力学过程的应用;在生物方面,则介绍利用飞秒技术研究光合作用中的能量传递过程,视觉系统中的光致异构化过程以及DNA中的电荷传递及质子传递等过程的研究现状。     

Noisy optical pulses enhance the temporal resolution of pump-probe spectroscopy

Time-resolved measurements of quantum dynamics are based on the availability of controlled events that are shorter than the typical evolution time scale of the processes to be observed. Here we introduce the concept of noise-enhanced pump-probe spectroscopy, allowing the measurement of dynamics significantly shorter than the average pulse duration by exploiting randomly varying, partially coherent light fields consisting of bunched colored noise. These fields are shown to be superior by more than a factor of 10 to frequency-stabilized fields, with important implications for time-resolved experiments at x-ray free-electron lasers and, in general, for measurements at the frontiers of temporal resolution (e.g., attosecond spectroscopy). As an example application, the concept is used to explain the recent experimental observation of vibrational wave-packet motion in D(2)(+) on time scales shorter than the average pulse duration.     

Dominant pathways in protein folding

We present a method to investigate the kinetics of protein folding and the dynamics underlying the formation of secondary and tertiary structures during the entire reaction. By writing the solution of the Fokker-Planck equation in terms of a path integral, we derive a Hamilton-Jacobi variational principle from which we are able to compute the most probable pathway of folding. The method is applied to the folding of the Villin headpiece subdomain simulated using a Go model. An initial collapsing phase driven by the initial configuration is followed by a rearrangement phase, in which secondary structures are formed and all computed paths display strong similarities. This completely general method does not require the prior knowledge of any reaction coordinate and is an efficient tool to perform simulations of the entire folding process with available computers.     

Hydrophobic and Ionic-Interactions in Bulk and Confined Water with Implications for Collapse and Folding of Proteins

Water and water-mediated interactions determine the thermodynamics and kinetics of protein folding, protein aggregation and self-assembly in confined spaces. To obtain insights into the role of water in the context of folding problems, we describe computer simulations of a few related model systems. The dynamics of collapse of eicosane shows that upon expulsion of water the linear hydrocarbon chain adopts an ordered helical hairpin structure with 1.5 turns. The structure of dimer of eicosane molecules has two well ordered helical hairpins that are stacked perpendicular to each other. As a prelude to studying folding in confined spaces we used simulations to understand changes in hydrophobic and ionic interactions in nano-sized water droplets. Solvation of hydrophobic and charged species change drastically in nano-scale water droplets. Hydrophobic species are localized at the boundary. The tendency of ions to be at the boundary where water density is low increases as the charge density decreases. The interactions between hydrophobic, polar, and charged residue are also profoundly altered in confined spaces. Using the results of computer simulations and accounting for loss of chain entropy upon confinement we argue and then demonstrate, using simulations in explicit water, that ordered states of generic amphiphilic peptide sequences should be stabilized in cylindrical nanopores.     

Phase contrast MRI with improved temporal resolution by view sharing: k-space related velocity mapping properties.

Phase contrast techniques in combination with k-space segmented CINE imaging are widely used for the quantitative assessment of blood flow or tissue motion. The temporal resolution of the corresponding pulse sequences plays an important role concerning the potential of the method to fully detect time resolved flow or motion patterns. A further improvement of temporal or spatial resolution in phase contrast CINE MRI can be achieved by the application of view sharing. Based on simulations with point-spread-functions resulting from different cyclic flow or motion patterns an analysis of view sharing techniques in combination with phase contrast MRI is presented. Velocity mapping properties and the role of different k-space regions concerning the resulting values in the phase images and thus encoded velocities were investigated. It could be shown that the velocity induced phase shifts in phase contrast techniques are mainly encoded in the central sections of k-space which makes view sharing also suitable for velocity mapping. As a result the use of appropriate sampling and data acquisition schemes permits the assessment of flow or motion patterns with significantly improved temporal resolution without loss of functional information. In addition phantom measurements with an oscillation phantom were performed in order to validate the simulation results and to demonstrate the potential of view sharing techniques to accelerate phase contrast imaging and improve the detection of the underlying flow or motion dynamics.     

Simple theory of protein folding kinetics

We present a simple model of protein folding dynamics that captures key qualitative elements recently seen in all-atom simulations. The goals of this theory are to serve as a simple formalism for gaining deeper insight into the physical properties seen in detailed simulations as well as to serve as a model to easily compare why these simulations suggest a different kinetic mechanism than previous simple models. Specifically, we find that non-native contacts play a key role in determining the mechanism, which can shift dramatically as the energetic strength of non-native interactions is changed. For proteinlike non-native interactions, our model finds that the native state is a kinetic hub, connecting the strength of relevant interactions directly to the nature of folding kinetics.     

Kinetic characterization of strongly coupled systems

We propose a simple method to determine the local coupling strength Gamma experimentally, by linking the individual particle dynamics with the local density and crystal structure of a 2D plasma crystal. By measuring particle trajectories with high spatial and temporal resolution we obtain the first maps of Gamma and temperature at individual particle resolution. We employ numerical simulations to test this new method, and discuss the implications to characterize strongly coupled systems.     

生物大分子多尺度理论和计算方法

分子模拟是研究生物大分子的重要手段. 过去二十年来, 人们将分子模拟与实验研究相结合, 揭示出生物大分子结构和动力学方面的诸多重要性质. 传统分子模拟主要采用全原子分子模型或各种粗粒化的分子模型. 在实际应用中, 传统分子模拟方法通常存在精度或效率瓶颈, 一定程度上限制了其应用范围. 近年来, 多尺度分子模型越来越受到人们的关注. 多尺度分子模型基于统计力学原理, 将全原子模型和粗粒化模型相耦合, 有望克服传统分子模拟方法中的精度/效率瓶颈, 进而拓展分子模拟在生物大分子研究中的应用范围. 根据模型之间的耦合方式, 近年来发展起来的多尺度分子模拟方法可归纳为如下四种类型: 混合分辨多尺度模型、并行耦合多尺度模型、单向耦合多尺度模型、以及自学习多尺度模型. 本文将对上述四类多尺度模型做简要介绍, 并讨论其主要优缺点、应用范围以及进一步发展方向.     

Microsecond resolution of quasiparticle tunneling in the single-Cooper-pair transistor

We present radio-frequency measurements on a single-Cooper-pair-transistor in which individual quasiparticle poisoning events were observed with microsecond temporal resolution. Thermal activation of the quasiparticle dynamics is investigated, and consequently, we are able to determine energetics of the poisoning and unpoisoning processes. In particular, we are able to assign an effective quasiparticle temperature to parametrize the poisoning rate.     

A method of improving overall resolution in ultrasonic array imaging using spatio-temporal deconvolution

In this paper a beamforming method for ultrasonic array imaging is presented that performs both spatial and temporal deconvolution based on a minimum mean square error (MMSE) criteria. The presented MMSE receive mode beamformer performs a regularized inversion of the propagation operator for the ultrasonic array system at hand. The MMSE beamformer accounts for the transmit and receive processes, defined in terms of finite array element sizes, transmit focusing laws and electrical transducer characteristics. The MMSE beamformer is compared to the traditional delay-and-sum (DAS) beamformer with respect to both resolution and signal-to-noise ratio. The two algorithms are compared using both simulated and measured data. The simulated data was obtained using ultrasonic field simulations and the measured data was acquired using a linear phased array imaging wire targets in water. The results show that the MMSE beamformer has superior temporal and lateral resolution compared to DAS. It is also shown that the MMSE beamformer can be expressed as a filter bank, which enables parallel processing at high frame rates.