Energy Transport Pathways in Proteins: A Non-equilibrium Molecular Dynamics Simulation Study

Temperature Dependence of Thermal Conductivity of Proteins: Contributions of Thermal Expansion and Grüneisen Parameter

The thermal conductivity of many materials depends on temperature due to several factors, including variation of heat capacity with temperature, changes in vibrational dynamics with temperature, and change in volume with temperature. For proteins some, but not all, of these influences on the variation of thermal conductivity with temperature have been investigated in the past. In this study, we examine the influence of change in volume, and corresponding changes in vibrational dynamics, on the temperature dependence of the thermal conductivity. Using a measured value for the coefficient of thermal expansion and recently computed values for the Grüneisen parameter of proteins we find that the thermal conductivity increases with increasing temperature due to change in volume with temperature. We compare the impact of thermal expansion on the variation of the thermal conductivity with temperature found in this study with contributions of heat capacity and anharmonic coupling examined previously. Using values of thermal transport coefficients computed for proteins we also model heating of water in a protein solution following photoexcitation.

Allosteric Communication Mediated by Protein Contact Clusters: A Dynamical Model

Describing the puzzling phenomenon of long-range communication between distant protein sites, allostery is of paramount importance in biomolecular regulation and signal transduction. It is commonly assumed to arise from a conformational rearrangement of the protein, although the underlying dynamical process has remained largely elusive. This study introduces a dynamical model of allosteric communication based on "contact clusters"─localized groups of highly correlated contacts that facilitate interactions between secondary structures. The model shows that allostery involves a multistep process with cooperative contact changes within clusters and communication between distant clusters mediated by rigid secondary structures. Considering time-dependent experiments on a photoswitchable PDZ3 domain, extensive (in total ∼500 μs) molecular dynamics simulations are conducted that directly monitor the photoinduced allosteric transition. The structural reorganization is illustrated by the time evolution of the contact clusters and the ligand, which effects the nonlocal coupling between distant clusters. A time scale analysis reveals dynamics from nano- to microseconds, which are in excellent agreement with the experimentally measured time scales. While the simulation of larger systems may require enhanced sampling techniques, it is expected that the general picture of allostery mediated by communicating contact clusters will still be applicable.

Intermolecular Protein-Water Coupling Impedes the Coupling Between the Amide A and Amide I Mode in Interfacial Proteins

The coupling between different vibrational modes in proteins is essential for chemical dynamics and biological functions and is linked to the propagation of conformational changes and pathways of allosteric communication. However, little is known about the influence of intermolecular protein-H2O coupling on the vibrational coupling between amide A (NH) and amide I (C═O) bands. Here, we investigate the NH/CO coupling strength in various peptides with different secondary structures at the lipid cell membrane/H2O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy (SFG-VS) in which a femtosecond infrared pump is used to excite the amide A band, and SFG-VS is used to probe transient spectral evolution in the amide A and amide I bands. Our results reveal that the NH/CO coupling strength strongly depends on the bandwidth of the amide I mode and the coupling of proteins with water molecules. A large extent of protein-water coupling significantly reduces the delocalization of the amide I mode along the peptide chain and impedes the NH/CO coupling strength. A large NH/CO coupling strength is found to show a strong correlation with the high energy transfer rate found in the light-harvesting proteins of green sulfur bacteria, which may understand the mechanism of energy transfer through a molecular system and assist in controlling vibrational energy transfer by engineering the molecular structures to achieve high energy transfer efficiency.

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.

Energy Transport in Class B GPCRs: Role of Protein-Water Dynamics and Activation

We compute energy exchange networks (EENs) through glucagon-like peptide-1 receptor (GLP-1R), a class B G-protein-coupled receptor (GPCR), in inactive and two active states, one activated by a peptide ligand and the other by a small molecule agonist, from results of molecular dynamics simulations. The reorganized network upon activation contains contributions from structural as well as from dynamic changes and corresponding entropic contributions to the free energy of activation, which are estimated in terms of the change in rates of energy transfer across non-covalent contacts. The role of water in the EENs and in activation of GLP-1R is also investigated. The dynamics of water in contact with the central polar network of the transmembrane region is found to be significantly slower for both activated states compared to the inactive state. This result is consistent with the contribution of water molecules to activation of GLP-1R previously suggested and resembles water dynamics in parts of the transmembrane region found in earlier studies of rhodopsin-like GPCRs.

Computational Study on the Thermal Conductivity of a Protein

Protein molecules are thermally fluctuating and tightly packed amino acid residues strongly interact with each other. Such interactions are characterized in terms of heat current at the atomic level. We calculated the thermal conductivity of a small globular protein, villin headpiece subdomain, based on the linear response theory using equilibrium molecular dynamics simulation. The value of its thermal conductivity was 0.3 ± 0.01 [W m^-1 K^-1], which is in good agreement with experimental and computational studies on the other proteins in the literature. Heat current along the main chain was dominated by local vibrations in the polypeptide bonds, with amide I, II, III, and A bands on the Fourier transform of the heat current autocorrelation function.

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.

Kanamycin-Mediated Conformational Dynamics of Escherichia coli Outer Membrane Protein TolC

TolC is a member of the outer membrane efflux proteins (OEPs) family and acts as an exit duct to export proteins, antibiotics, and substrate molecules across the Escherichia coli cell membrane. Export of these molecules is evidenced to be brought about through the reversible interactions and binding of substrate-specific drug molecules or antibiotics with TolC and by being open for transport, which afterward leads to cross-resistance. Hence, the binding of kanamycin with TolC was monitored through molecular docking (MD), the structural fluctuations and conformational changes to the atomic level. The results were further supported from the steady-state fluorescence binding and isothermal titration calorimetry (ITC) studies. Binding of kanamycin with TolC resulted in a concentration dependent fluorescence intensity quenching with 7 nm blue shift. ITC binding data maintains a single binding site endothermic energetic curve with binding parameters indicating an entropy driven binding process. The confirmational changes resulting from this binding were monitored by a circular dichroism (CD) study, and the results showed insignificant changes in the α-helix and β-sheets secondary structure contents, but the tertiary structure shows inclusive changes in the presence of kanamycin. The experimental data substaintially correlates the RMSD, R _g, and RMSF results. The resulting conformational changes of the TolC-kanamycin complexation was stabilized through H-bonding and other interactions.

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.

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.

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.