Nonequilibrium molecular dynamics simulation of the energy transport through a peptide helix

Thermal Energy Transport through Nonbonded Native Contacts in Protein

Within the protein interior, where we observe various types of interactions, nonuniform flow of thermal energy occurs along the polypeptide chain and through nonbonded native contacts, leading to inhomogeneous transport efficiencies from one site to another. The folded native protein serves not merely as thermal transfer medium but, more importantly, as sophisticated molecular nanomachines in cells. Therefore, we are particularly interested in what sort of "communication" is mediated through native contacts in the folded proteins and how such features are quantitatively depicted in terms of local transport coefficients of heat currents. To address the issue, we introduced a concept of inter-residue thermal conductivity and characterized the nonuniform thermal transport properties of a small globular protein, HP36, using equilibrium molecular dynamics simulation and the Green-Kubo formula. We observed that the thermal transport of the protein was dominated by that along the polypeptide chain, while the local thermal conductivity of nonbonded native contacts decreased in the order of H-bonding > π-stacking > electrostatic > hydrophobic contacts. Furthermore, we applied machine learning techniques to analyze the molecular mechanism of protein thermal transport. As a result, the contact distance, variance in contact distance, and H-bonding occurrence probability during MD simulations are found to be the top three important determinants for predicting local thermal transport coefficients.

Study on Thermophysical Properties and Phase Change Regulation Mechanism of Optically-Controlled Phase Change Materials: Synthesis, Crystal Structure and Molecular Dynamics

Optically-controlled phase change materials, which are prepared by introducing molecular photoswitches into traditional phase change materials (PCMs), can convert and store solar energy into photochemical enthalpy and phase change enthalpy. However, the thermophysical properties of optically controlled PCMs, which are crucial in the practical, are rarely paid attention to. 4-(phenyldiazenyl)phenyl decanoate (Azo-A-10) is experimentally prepared as an optically-controlled PCMs, whose energy storage density is 210.0 kJ·kg-1, and the trans single crystal structure is obtained. The density, phase transition temperature, thermal conductivity, and other parameters in trans state are measured experimentally. Furthermore, a microscopic model of Azo-A-10 is established, and the thermophysical properties are analyzed based on molecular dynamics. The results show that the microstructure parameter (order parameters) and thermophysical properties (density, radial distribution function, self-diffusion coefficient, phase change temperature, and thermal conductivity) of partially or completely isomerized Azo-A-10, which are challenging to observe in experiments, can be predicted by molecular dynamics simulation. The optically-controlled phase change mechanism can be clarified according to the differences in microstructure. The optically-controlled switchability of thermophysical properties of an optically-controlled PCM is analyzed. This study provides ideas for the improvement, development, and application of optically-controlled PCMs in the future.

Vibrational Energy Landscapes and Energy Flow in GPCRs

We construct and analyze disconnectivity graphs to provide the first graphical representation of the vibrational energy landscape of a protein, in this study β2AR, a G-protein coupled receptor (GPCR), in active and inactive states. The graphs, which indicate the relative free energy of each residue and the minimum free energy barriers for energy transfer between them, reveal important composition, structural and dynamic properties that mediate the flow of energy. Prolines and glycines, which contribute to GPCR plasticity and function, are identified as bottlenecks to energy transport along the backbone from which alternative pathways for energy transport via nearby noncovalent contacts emerge, seen also in the analysis of first passage time (FPT) distributions presented here. Striking differences between the disconnectivity graphs and FPT distributions for the inactive and active states of β2AR are found where structural and dynamic changes occur upon activation, contributing to allosteric regulation.

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose-Einstein statistics into an "effective temperature" in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

Inhibition of vibrational energy flow within an aromatic scaffold via heavy atom effect

The regulation of intramolecular vibrational energy redistribution (IVR) to influence energy flow within molecular scaffolds provides a way to steer fundamental processes of chemistry, such as chemical reactivity in proteins and design of molecular diodes. Using two-dimensional infrared (2D IR) spectroscopy, changes in the intensity of vibrational cross-peaks are often used to evaluate different energy transfer pathways present in small molecules. Previous 2D IR studies of para-azidobenzonitrile (PAB) demonstrated that several possible energy pathways from the N3 to the cyano-vibrational reporters were modulated by Fermi resonance, followed by energy relaxation into the solvent [Schmitz et al., J. Phys. Chem. A 123, 10571 (2019)]. In this work, the mechanisms of IVR were hindered via the introduction of a heavy atom, selenium, into the molecular scaffold. This effectively eliminated the energy transfer pathway and resulted in the dissipation of the energy into the bath and direct dipole-dipole coupling between the two vibrational reporters. Several structural variations of the aforementioned molecular scaffold were employed to assess how each interrupted the energy transfer pathways, and the evolution of 2D IR cross-peaks was measured to assess the changes in the energy flow. By eliminating the energy transfer pathways through isolation of specific vibrational transitions, through-space vibrational coupling between an azido (N3) and a selenocyanato (SeCN) probe is facilitated and observed for the first time. Thus, the rectification of this molecular circuitry is accomplished through the inhibition of energy flow using heavy atoms to suppress the anharmonic coupling and, instead, favor a vibrational coupling pathway.

Locating dynamic contributions to allostery via determining rates of vibrational energy transfer

Determining rates of energy transfer across non-covalent contacts for different states of a protein can provide information about dynamic and associated entropy changes during transitions between states. We investigate the relationship between rates of energy transfer across polar and nonpolar contacts and contact dynamics for the β₂-adrenergic receptor, a rhodopsin-like G-protein coupled receptor, in an antagonist-bound inactive state and agonist-bound active state. From structures sampled during molecular dynamics (MD) simulations, we find the active state to have, on average, a lower packing density, corresponding to generally more flexibility and greater entropy than the inactive state. Energy exchange networks (EENs) are computed for the inactive and active states from the results of the MD simulations. From the EENs, changes in the rates of energy transfer across polar and nonpolar contacts are found for contacts that remain largely intact during activation. Change in dynamics of the contact, and entropy associated with the dynamics, can be estimated from the change in rates of energy transfer across the contacts. Measurement of change in the rates of energy transfer before and after the transition between states thereby provides information about dynamic contributions to activation and allostery.

Energy Transport and Its Function in Heptahelical Transmembrane Proteins

Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β₂ adrenergic receptor (β₂AR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [J. Chem. Phys.2020, 152, 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure. Despite their overall structural similarity, the three 7TM proteins reveal quite different strategies to redistribute excess energy. While bR quickly removes the energy using the TM7 helix as a "lightning rod", Rho exhibits a rather poor energy dissipation, which might eventually require the hydrolysis of the Schiff base between the protein and the retinal chromophore to prevent overheating. Heating the ligand adrenaline of β₂AR, the resulting energy transport network of the protein is found to change significantly upon switching from the active state to the inactive state. While the energy flow may highlight aspects of the inter-residue couplings of β₂AR, it seems not particularly suited to explain allosteric phenomena.

Dependence of Vibrational Energy Transfer on Distance in a Four-Helix Bundle Protein: Equidistant Increments with the Periodicity of α Helices

Vibrational energy exchanges between various degrees of freedom are critical to barrier-crossing processes in proteins. Heme proteins are highly suitable for studies of the vibrational energy exchanges in proteins. The migration of excess energy released by heme in a protein moiety can be observed using time-resolved anti-Stokes ultraviolet resonance Raman spectroscopy. The anti-Stokes resonance Raman intensity of a tryptophan residue is an excellent probe for the excess energy and the spatial resolution of a single amino acid residue can be achieved. Here, we studied dependence of vibrational energy transfer on the distance in cytochrome b₅₆₂, which is a heme-containing, four-helix bundle protein. The vibrational energy transfer from the heme group to a single tryptophan residue introduced by site-directed mutagenesis was examined for different heme-tryptophan distances by a quasi-constant length with the periodicity of α helices. Taken together with structural data obtained by molecular dynamics simulations, the energy transfer could be well described by the model of classical thermal diffusion, which suggests that continuum media provide a good approximation of the protein interior, of which the atomic packing density is very high.

Learning the Hydrophobic, Hydrophilic, and Aromatic Character of Amino Acids from Thermal Relaxation and Interfacial Thermal Conductance

In this study, the thermal relaxation of the 20 naturally occurring amino acids in water and in the protein lysozyme is investigated using transient nonequilibrium molecular dynamics simulations. By modeling the thermal relaxation process, the relaxation times of the amino acids in water occurs over a time scale covering 2-5 ps. For the hydrophobic amino acids, the relaxation time is controlled by the size of the hydrocarbon side chain, while for hydrophilic amino acids, the number of hydrogen bonds does not significantly affect the time scales of the heat dissipation. Our results show that the interfacial thermal conductance at the amino acid-water interface is in the range of 40-80 MW m^-2 K^-1. Hydrophobic and aromatic amino acids tend to have a lower interfacial thermal conductance. Notably, we show that amino acids can be correlated with their thermal relaxation times and molar masses, into simply connected phases with the same hydrophilicity, hydrophobicity, and aromaticity. The thermal relaxation slows down by a factor of up to five in the protein relative to that in water. In the case of the hydrophobic amino acids in the protein lysozyme, the slow down in the thermal relaxation relative to that in water appears to be controlled primarily by the size of the side chain.

Through bonds or contacts? Mapping protein vibrational energy transfer using non-canonical amino acids

Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the β-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway.

Asymmetric dynamics of dimeric SARS-CoV-2 and SARS-CoV main proteases in an apo form: Molecular dynamics study on fluctuations of active site, catalytic dyad, and hydration water

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been widely spread around the world. It is necessary to examine the viral proteins that play a notorious role in the invasion of our body. The main protease (3CL^pro) facilitates the maturation of the coronavirus. It is thought that the dimerization of 3CL^pro leads to its catalytic activity; the detailed mechanism has, however, not been suggested. Furthermore, the structural differences between the predecessor SARS-CoV 3CL^pro and SARS-CoV-2 3CL^pro have not been fully understood. Here, we show the structural and dynamical differences between the two main proteases, and demonstrate the relationship between the dimerization and the activity via atomistic molecular dynamics simulations. Simulating monomeric and dimeric 3CL^pro systems for each protease, we show that (i) global dynamics between the two different proteases are not conserved, (ii) the dimerization stabilizes the catalytic dyad and hydration water molecules behind the dyad, and (iii) the substrate-binding site (active site) and hydration water molecules in each protomer fluctuate asymmetrically. We then speculate the roles of hydration water molecules in their catalytic activity.

Nonequilibrium molecular dynamics simulations of infrared laser-induced dissociation of a tetrameric Aβ42 β-barrel in a neuronal membrane model

Experimental studies have reported that the amyloid-β proteins can form pores in cell membranes, and this could be one possible source of toxicity in Alzheimer's disease. Dissociation of these pores could therefore be a potential therapeutic approach. It is known that high photon density free-electron laser experiments and laser-induced nonequilibrium molecular dynamics simulations (NEMD) can dissociate amyloid fibrils at specific frequencies in vitro. Our question is whether NEMD simulations can dissociate amyloid pores in a bilayer mimicking a neuronal membrane, and as an example, we select a tetrameric Aβ42 β-barrel. Our simulations shows that the resonance between the laser field and the amide I vibrational mode of the barrel destabilises all intramolecular and intermolecular hydrogen bonds of Aβ42 and converts the β-barrel to a random/coil disordered oligomer. Starting from this disordered oligomer, extensive standard MD simulations shows sampling of disordered Aβ42 states without any increase of β-sheet and reports that the orientational order of lipids is minimally disturbed. Interestingly, the frequency to be employed to dissociate this beta-barrel is specific to the amino acid sequence. Taken together with our previous simulation results, this study indicates that infrared laser irradiation can dissociate amyloid fibrils and oligomers in bulk solution and in a membrane environment without affecting the surrounding molecules, offering therefore a promising way to retard the progression of AD.

MSMPathfinder: Identification of Pathways in Markov State Models

Markov state models represent a popular means to interpret biomolecular processes in terms of memoryless transitions between metastable conformational states. To gain insight into the underlying mechanism, it is instructive to determine all relevant pathways between initial and final states of the process. Currently available methods, such as Markov chain Monte Carlo and transition path theory, are convenient for identifying the most frequented pathways. They are less suited to account for the typically huge amount of pathways with low probability which, though, may dominate the cumulative flux of the reaction. On the basis of a systematic construction of all possible pathways, the here proposed method MSMPathfinder is able to characterize the multitude of unique pathways (say, up to 10¹⁰) in a complex system and to quantitatively calculate their correct weights and associated waiting times with predefined accuracy. Adopting the chiral transitions of a peptide helix and the folding of the villin headpiece as model problems, mechanisms and associated waiting times of these processes are discussed using a kinetic network representation. The analysis reveals that the waiting time distribution may yield only little insight into the diversity of pathways, because the measured folding times do typically not reflect the most probable path lengths but rather the cumulative effect of many different pathways.

Energy Transfer across Nonpolar and Polar Contacts in Proteins: Role of Contact Fluctuations

Molecular dynamics simulations of the villin headpiece subdomain HP36 have been carried out to examine relations between rates of vibrational energy transfer across non-covalently bonded contacts and equilibrium structural fluctuations, with focus on van der Waals contacts. Rates of energy transfer across van der Waals contacts vary inversely with the variance of the contact length, with the same constant of proportionality for all nonpolar contacts of HP36. A similar relation is observed for hydrogen bonds, but the proportionality depends on contact pairs, with hydrogen bonds stabilizing the α-helices all exhibiting the same constant of proportionality, one that is distinct from those computed for other polar contacts. Rates of energy transfer across van der Waals contacts are found to be up to 2 orders of magnitude smaller than rates of energy transfer across polar contacts.

Stochastic simulation of nonequilibrium heat conduction in extended molecular junctions

Understanding phononic heat transport processes in molecular junctions is a central issue in the developing field of nanoscale heat conduction. Here, we present a Langevin dynamics simulation framework to investigate heat transport processes in molecular junctions at and beyond the linear response regime and apply it to saturated and unsaturated linear hydrocarbon chains connecting two gold substrates. Thermal boundary conditions represented by Markovian noise and damping are filtered through several (up to four) gold layers to provide a realistic and controllable bath spectral density. Classical simulations using the full universal force field are compared with quantum calculations that use only the harmonic part of this field. The close agreement found at about room temperature between these very different calculations suggests that heat transport at such temperatures is dominated by lower frequency vibrations whose dynamics is described well by classical mechanics. The results obtained for alkanedithiol molecules connecting gold substrates agree with previous quantum calculations based on the Landauer formula and match recent experimental measurements [e.g., thermal conductance around 20 pW/K for alkanedithiols in single-molecule junctions (SMJs)]. Heat conductance simulations on polyynes of different lengths illuminate the effects of molecular conjugation on thermal transport. The difference between alkanes and polyynes is not large but correlates with the larger rigidity and stronger mode localization that characterize the polyyne structure. This computational approach has been recently used [R. Chen, I. Sharony, and A. Nitzan, J. Phys. Chem. Lett. 11, 4261-4268 (2020)] to unveil local atomic heat currents and phononic interference effect in aromatic-ring based SMJs.

Energy flow and intersubunit signalling in GSAM: A non-equilibrium molecular dynamics study

Non-equilibrium molecular dynamics simulations of vibrational energy flow induced by the imposition of a thermal gradient have been performed on the μ2-dimeric enzyme glutamate-1-semialdehyde aminomutase (GSAM), the key enzyme in the biosynthesis of chlorophyll, in order to identify energy transport pathways and to elucidate their role as potential allosteric communication networks for coordinating functional dynamics, specifically the negative cooperativity observed in the motion of the two active site gating loops. Fully atomistic MD simulations of thermal diffusion were executed with a GROMACS simulation package on a fully solvated GSAM enzyme by heating various active site target ligands (initially, catalytic intermediates and cofactors) to 300K while holding the remainder of the protein and the solvent bath at 10K and monitoring the temperature T(t) of all the enzyme residues as a function of time over a 1ns observation window. Energy is observed to be deposited in a relatively small number of discrete chains of residues most of which contribute to specific structural or biochemical functionality. Thermal linkages between all thermally active chains were established by isolating a specific pair of chains and performing a thermal diffusion simulation on the pair, one held at 300K and the other at 10K, with the rest of the protein frozen in its initial atomic configuration and thus thermally unresponsive. Proceeding in this way, it was possible to map out multiple pathways of vibrational energy flow leading from one of the active sites through a network of contiguous residues, many of which were evolutionarily conserved and linked by hydrogen bonds, into the other active site and ultimately to the other gating loop.

Infrared Laser-Induced Amyloid Fibril Dissociation: A Joint Experimental/Theoretical Study on the GNNQQNY Peptide

Neurodegenerative diseases are usually characterized by plaques made of well-ordered aggregates of distinct amyloid proteins. Dissociating these very stable amyloid plaques is a critical clinical issue. In this study, we present a joint mid-infrared free electron laser experiment/nonequilibrium molecular dynamics simulation to understand the dissociation process of a representative example GNNQQNY fibril. By tuning the laser frequency to the amide I band of the fibril, the resonance takes place and dissociation is occurred. With the calculated and observed wide-angle X-ray scattering profiles and secondary structures before and after laser irradiation being identical, we can propose a dissociation mechanism with high confidence from our simulations. We find that dissociation starts in the core of the fibrils by fragmenting the intermolecular hydrogen bonds and separating the peptides and then propagates to the fibril extremities leading to the formation of unstructured expanded oligomers. We suggest that this should be a generic mechanism of the laser-induced dissociation of amyloid fibrils.

Tuning Molecular Vibrational Energy Flow within an Aromatic Scaffold via Anharmonic Coupling

From guiding chemical reactivity in synthesis or protein folding to the design of energy diodes, intramolecular vibrational energy redistribution harnesses the power to influence the underlying fundamental principles of chemistry. To evaluate the ability to steer these processes, the mechanism and time scales of intramolecular vibrational energy redistribution through aromatic molecular scaffolds have been assessed by utilizing two-dimensional infrared (2D IR) spectroscopy. 2D IR cross peaks reveal energy relaxation through an aromatic scaffold from the azido- to the cyano-vibrational reporters in para-azidobenzonitrile (PAB) and para-(azidomethyl)benzonitrile (PAMB) prior to energy relaxation into the solvent. The rates of energy transfer are modulated by Fermi resonances, which are apparent by the coupling cross peaks identified within the 2D IR spectrum. Theoretical vibrational mode analysis allowed the determination of the origins of the energy flow, the transfer pathway, and a direct comparison of the associated transfer rates, which were in good agreement with the experimental results. Large variations in energy-transfer rates, approximately 1.9 ps for PAB and 23 ps for PAMB, illustrate the importance of strong anharmonic coupling, i.e., Fermi resonance, on the transfer pathways. In particular, vibrational energy rectification is altered by Fermi resonances of the cyano- and azido-modes allowing control of the propensity for energy flow.

Topology, landscapes, and biomolecular energy transport

While ubiquitous, energy redistribution remains a poorly understood facet of the nonequilibrium thermodynamics of biomolecules. At the molecular level, finite-size effects, pronounced nonlinearities, and ballistic processes produce behavior that diverges from the macroscale. Here, we show that transient thermal transport reflects macromolecular energy landscape architecture through the topological characteristics of molecular contacts and the nonlinear processes that mediate dynamics. While the former determines transport pathways via pairwise interactions, the latter reflects frustration within the landscape for local conformational rearrangements. Unlike transport through small-molecule systems, such as alkanes, nonlinearity dominates over coherent processes at even quite short time- and length-scales. Our exhaustive all-atom simulations and novel local-in-time and space analysis, applicable to both theory and experiment, permit dissection of energy migration in biomolecules. The approach demonstrates that vibrational energy transport can probe otherwise inaccessible aspects of macromolecular dynamics and interactions that underly biological function.

Understanding vibrational energy transfer in macromolecules has been challenging to both theory and experiment. Here the authors use non-equilibrium molecular dynamics to reveal the relationship between heat transport in a model peptide, emergent nonlinearity, and the underlying free energy landscape.

A spin-1 representation for dual-funnel energy landscapes

The interconversion between the left- and right-handed helical folds of a polypeptide defines a dual-funneled free energy landscape. In this context, the funnel minima are connected through a continuum of unfolded conformations, evocative of the classical helix-coil transition. Physical intuition and recent conjectures suggest that this landscape can be mapped by assigning a left- or right-handed helical state to each residue. We explore this possibility using all-atom replica exchange molecular dynamics and an Ising-like model, demonstrating that the energy landscape architecture is at odds with a two-state picture. A three-state model-left, right, and unstructured-can account for most key intermediates during chiral interconversion. Competing folds and excited conformational states still impose limitations on the scope of this approach. However, the improvement is stark: Moving from a two-state to a three-state model decreases the fit error from 1.6 kBT to 0.3 kBT along the left-to-right interconversion pathway.

Energy transport pathway in proteins: Insights from non-equilibrium molecular dynamics with elastic network model

Intra-molecular energy transport between distant functional sites plays important roles in allosterically regulating the biochemical activity of proteins. How to identify the specific intra-molecular signaling pathway from protein tertiary structure remains a challenging problem. In the present work, a non-equilibrium dynamics method based on the elastic network model (ENM) was proposed to simulate the energy propagation process and identify the specific signaling pathways within proteins. In this method, a given residue was perturbed and the propagation of energy was simulated by non-equilibrium dynamics in the normal modes space of ENM. After that, the simulation results were transformed from the normal modes space to the Cartesian coordinate space to identify the intra-protein energy transduction pathways. The proposed method was applied to myosin and the third PDZ domain (PDZ3) of PSD-95 as case studies. For myosin, two signaling pathways were identified, which mediate the energy transductions form the nucleotide binding site to the 50 kDa cleft and the converter subdomain, respectively. For PDZ3, one specific signaling pathway was identified, through which the intra-protein energy was transduced from ligand binding site to the distant opposite side of the protein. It is also found that comparing with the commonly used cross-correlation analysis method, the proposed method can identify the anisotropic energy transduction pathways more effectively.

MELD-Path Efficiently Computes Conformational Transitions, Including Multiple and Diverse Paths

The molecular actions of proteins occur along reaction coordinates. Current computer methods have limited ability to explore them. We describe a fast protocol called MELD-path that (1) efficiently samples relevant conformational states via MELD, an accelerator of Molecular Dynamics (MD), (2) seeds multiple short MD trajectories from MELD states, and then (3) constructs Markov State Models (MSM) that give the routes and kinetics. We tested the method against extensive (multi μs) MD simulations of the right-handed- to left-handed-helix transition of a 9-mer peptide of AIB, the symmetry of which allows us to establish convergence. MELD-path finds all the metastable states, their correct relative populations, and the full ensemble of routes, not just a single assumed route. For this transition, we find a very broad route structure. MELD-path is highly parallelizable and efficient, yielding the full route map in a few days of computation. We believe MELD-path could be a general and rapid way to explore mechanistic processes in biomolecules on the computer.

Probing the communication of deoxythymidine triphosphate in HIV-1 reverse transcriptase by communication maps and interaction energy studies

We calculate communication maps for HIV-1 Reverse Transcriptase (RT) to elucidate energy transfer pathways between deoxythymidine triphosphate (dTTP) and other parts of the protein.

Molecular Dynamics Study on the Photothermal Actuation of a Glassy Photoresponsive Polymer Reinforced with Gold Nanoparticles with Size Effect

We investigated the optical and thermal actuation behavior of densely cross-linked photoresponsive polymer (PRP) and polymer nanocomposites containing gold nanoparticles (PRP/Au) using all-atom molecular dynamics (MD) simulations. The modeled molecular structures contain a large number of photoreactive mesogens with linear orientation. Flexible side chains are interconnected through covalent bonds under periodic boundary conditions. A switchable dihedral potential was applied on a diazene moiety to describe the photochemical trans-to-cis isomerization. To quantify the photoinduced molecular reorientation and its effect on the macroscopic actuation of the neat PRP and PRP/Au materials, we characterized the photostrain and other material properties including elastic stiffness and thermal stability according to the photoisomerization ratio of the reactive groups. We particularly examined the effect of nanoparticle size on the photothermal actuation by varying the diameter of the nanofiller (10-20 Å) under the same volume fraction of 1.62%. The results indicated that the insertion of the gold nanoparticles enlarges the photostrain of the material while enhancing its mechanical stiffness and thermal stability. When the diameter of the nanoparticle reaches a size similar to or smaller than the length of the mesogen, the interfacial energy between the nanofiller and the surrounding polymer matrix does not significantly affect the alignment of the mesogens, but rather the adsorption energy at the interface generates a stable interphase layer. Hence, these improvements were more effective as the size of the gold nanoparticle decreased. The present findings suggest a wider analysis of the nanofiller-reinforced PRP composites and could be a guide for the mechanical design of the PRP actuator system.

Picosecond dissociation of amyloid fibrils with infrared laser: A nonequilibrium simulation study

Recently, mid-infrared free-electron laser technology has been developed to dissociate amyloid fibrils. Here, we present a theoretical framework for this type of experiment based on laser-induced nonequilibrium all-atom molecular dynamics simulations. We show that the fibril is destroyed due to the strong resonance between its amide I vibrational modes and the laser field. The effects of laser irradiation are determined by a balance between fibril formation and dissociation. While the overall rearrangements of the fibril finish over short time scales, the interaction between the peptides and the solvent continues over much longer times indicating that the waters play an important role in the dissociation process. Our results thus provide new insights into amyloid fibril dissociation by laser techniques and open up new venues to investigate the complex phenomena associated with amyloidogenesis.

Hierarchical Biomolecular Dynamics: Picosecond Hydrogen Bonding Regulates Microsecond Conformational Transitions

Biomolecules exhibit structural dynamics on a number of time scales, including picosecond (ps) motions of a few atoms, nanosecond (ns) local conformational transitions, and microsecond (μs) global conformational rearrangements. Despite this substantial separation of time scales, fast and slow degrees of freedom appear to be coupled in a nonlinear manner; for example, there is theoretical and experimental evidence that fast structural fluctuations are required for slow functional motion to happen. To elucidate a microscopic mechanism of this multiscale behavior, Aib peptide is adopted as a simple model system. Combining extensive molecular dynamics simulations with principal component analysis techniques, a hierarchy of (at least) three tiers of the molecule's free energy landscape is discovered. They correspond to chiral left- to right-handed transitions of the entire peptide that happen on a μs time scale, conformational transitions of individual residues that take about 1 ns, and the opening and closing of structure-stabilizing hydrogen bonds that occur within tens of ps and are triggered by sub-ps structural fluctuations. Providing a simple mechanism of hierarchical dynamics, fast hydrogen bond dynamics is found to be a prerequisite for the ns local conformational transitions, which in turn are a prerequisite for the slow global conformational rearrangement of the peptide. As a consequence of the hierarchical coupling, the various processes exhibit a similar temperature behavior which may be interpreted as a dynamic transition.

Scaling Rules for Vibrational Energy Transport in Globular Proteins

Computational studies of vibrational energy flow in biomolecules have to date mapped out transport pathways on a case-by-case basis. To provide a more general approach, we derive scaling rules for vibrational energy transport in a globular protein, which are identified from extensive nonequilibrium molecular dynamics simulations of vibrational energy flow in the villin headpiece subdomain HP36. We parametrize a master equation based on inter-residue, residue-solvent, and heater-residue energy-transfer rates, which closely reproduces the results of the all-atom simulations. From that fit, two scaling rules emerge, one for energy transport along the protein backbone which relies on a diffusion model and another for energy transport between tertiary contacts, which is based on a harmonic model. Requiring only the calculation of mean and variance of relatively few atomic distances, the approach holds the potential to predict the pathways and time scales of vibrational energy flow in large proteins.

2D-IR spectroscopy of hydrogen-bond-mediated vibrational excitation transfer

Inter-molecular vibrational energy transfer in the hydrogen-bonded complexes of methyl acetate and 4-cyanophenol is studied by dual-frequency 2D-IR spectroscopy.

Picosecond infrared laser-induced all-atom nonequilibrium molecular dynamics simulation of dissociation of viruses

Laser-induced all-atom nonequilibrium molecular dynamics simulation of virus dissociation.

Picosecond melting of peptide nanotubes using an infrared laser: a nonequilibrium simulation study

Resonance between carboxylate bond vibrations and laser frequency results in melting of nanotube.

Communication maps of vibrational energy transport through Photoactive Yellow Protein

We calculate communication maps for Photoactive Yellow Protein (PYP) from the purple phototropic eubacterium Halorhodospira halophile and use them to elucidate energy transfer pathways from the chromophore through the rest of the protein in the ground and excited state. The calculations reveal that in PYP excess energy from the chromophore flows mainly to regions of the surrounding residues that hydrogen bond to the chromophore. In addition, quantum mechanics/molecular mechanics and molecular dynamics (MD) simulations of the dielectric response of the protein and solvent environment due to charge rearrangement on the chromophore following photoexcitation are also presented, with both approaches yielding similar time constants for the response. Results of MD simulations indicate that the residues hydrogen bonding to the chromophore make the largest contribution to the response.

The promoting vibration in human heart lactate dehydrogenase is a preferred vibrational channel

We examine whether the rate-promoting vibration of lactate dehydrogenase is a preferred axis of thermal energy transfer. While it seems plausible that such a mechanistically important motion is also a favored direction of energy transfer, none of the previous studies of rate-promoting vibrations in enzymatic catalysis have addressed this question. It is equally likely that the promoting vibration, though catalytically important, has no different properties than any other axis in the protein. Resolution of this issue is important for two reasons: First, if energy is transferred along this axis in a preferred fashion, it shows that the protein is engineered in a way that transfers thermal energy into a motion that is coupled to the chemical step. Second, the discovery of a preferred direction of thermal transfer provides a potential route to experimental verification of the promoting vibration concept. Our computational experiments are specifically designed to mimic potential laser experiment with the deposition of thermal energy in an active-site chromophore with subsequent measurement of temperature at various points in the protein. Our results indicate that the promoting vibration is indeed a preferred channel of energy transfer. In addition, we study the vibrational structure of the protein via the dynamical structure factor to show preferred vibrational motion along the promoting vibration axis is an inherent property of the protein structure via thermal fluctuations.

Temperature dependence of the heat diffusivity of proteins

In a combined experimental-theoretical study, we investigated the transport of vibrational energy from the surrounding solvent into the interior of a heme protein, the sperm whale myoglobin double mutant L29W-S108L, and its dependence on temperature from 20 to 70 K. The hindered libration of a CO molecule that is not covalently bound to any part of the protein but is trapped in one of its binding pockets (the Xe4 pocket) was used as the local thermometer. Energy was deposited into the solvent by IR excitation. Experimentally, the energy transfer rate increased from (30 ps)(-1) at 20 K to (8 ps)(-1) at 70 K. This temperature trend is opposite to what is expected, assuming that the mechanism of heat transport is similar to that in glasses. In order to elucidate the mechanism and its temperature dependence, nonequilibrium molecular dynamics (MD) simulations were performed, which, however, predicted an essentially temperature-independent rate of vibrational energy flow. We tentatively conclude that the MD potentials overestimate the coupling between the protein and the CO molecule, which appears to be the rate-limiting step in the real system at low temperatures. Assuming that this coupling is anharmonic in nature, the observed temperature trend can readily be explained.