Anisotropic energy flow and allosteric ligand binding in albumin

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.

Ballistic Brownian Motion of Nanoconfined DNA

Theoretical treatments of polymer dynamics in liquid generally start with the basic assumption that motion at the smallest scale is heavily overdamped; therefore, inertia can be neglected. We report on the Brownian motion of tethered DNA under nanoconfinement, which was analyzed by molecular dynamics simulation and nanoelectrochemistry-based single-electron shuttle experiments. Our results show a transition into the ballistic Brownian motion regime for short DNA in sub-5 nm gaps, with quality coefficients as high as 2 for double-stranded DNA, an effect mainly attributed to a drastic increase in stiffness. The possibility for DNA to enter the underdamped regime could have profound implications on our understanding of the energetics of biomolecular engines such as the replication machinery, which operates in nanocavities that are a few nanometers wide.

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen-deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein-solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation (Ea) for the Stokes shifts decay rates and the millisecond C-H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

Time-resolved spectroscopic mapping of vibrational energy flow in proteins: Understanding thermal diffusion at the nanoscale

Vibrational energy exchange between various degrees of freedom is critical to barrier-crossing processes in proteins. Hemeproteins are well suited for studying vibrational energy exchange in proteins because the heme group is an efficient photothermal converter. The released energy by heme following photoexcitation shows migration in a protein moiety on a picosecond timescale, which is observed using time-resolved ultraviolet resonance Raman spectroscopy. The anti-Stokes ultraviolet resonance Raman intensity of a tryptophan residue is an excellent probe for the vibrational energy in proteins, allowing the mapping of energy flow with the spatial resolution of a single amino acid residue. This Perspective provides an overview of studies on vibrational energy flow in proteins, including future perspectives for both methodologies and applications.

Time-resolved infra-red studies of photo-excited porphyrins in the presence of nucleic acids and in HeLa tumour cells: insights into binding site and electron transfer dynamics

Cationic porphyrins based on the 5,10,15,20-meso-(tetrakis-4-N-methylpyridyl) core (TMPyP4) have been studied extensively over many years due to their strong interactions with a variety of nucleic acid structures, and their potential use as photodynamic therapeutic agents and telomerase inhibitors. In this paper, the interactions of metal-free TMPyP4 and Pt(II)TMPyP4 with guanine-containing nucleic acids are studied for the first time using time-resolved infrared spectroscopy (TRIR). In D2O solution (where the metal-free form exists as D2TMPyP4) both compounds yielded similar TRIR spectra (between 1450-1750 cm-1) following pulsed laser excitation in their Soret B-absorption bands. Density functional theory calculations reveal that vibrations centred on the methylpyridinium groups are responsible for the dominant feature at ca. 1640 cm-1. TRIR spectra of D2TMPyP4 or PtTMPyP4 in the presence of guanosine 5'-monophosphate (GMP), double-stranded {d(GC)5}2 or {d(CGCAAATTTGCG)}2 contain negative-going signals, 'bleaches', indicative of binding close to guanine. TRIR signals for D2TMPyP4 or PtTMPyP bound to the quadruplex-forming cMYC sequence {d(TAGGGAGGG)}2T indicate that binding occurs on the stacked guanines. For D2TMPyP4 bound to guanine-containing systems, the TRIR signal at ca. 1640 cm-1 decays on the picosecond timescale, consistent with electron transfer from guanine to the singlet excited state of D2TMPyP4, although IR marker bands for the reduced porphyrin/oxidised guanine were not observed. When PtTMPyP is incorporated into HeLa tumour cells, TRIR studies show protein binding with time-dependent ps/ns changes in the amide absorptions demonstrating TRIR's potential for studying light-activated molecular processes not only with nucleic acids in solution but also in biological cells.

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.

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.

Allosteric Type and Pathways Are Governed by the Forces of Protein-Ligand Binding

Allostery is central to many cellular processes, by up- or down-regulating target function. However, what determines the allosteric type remains elusive and currently it is impossible to predict whether the allosteric compounds would activate or inhibit target function before experimental studies. We demonstrated that the allosteric type and allosteric pathways are governed by the forces imposed by ligand binding to target protein using the anisotropic network model and developed an allosteric type prediction method (AlloType). AlloType correctly predicted 13 of the 16 allosteric systems in the data set with experimentally determined protein and complex structures as well as verified allosteric types, which was also used to identify allosteric pathways. When applied to glutathione peroxidase 4, a protein with no complex structure information, AlloType could still be able to predict the allosteric type of the recently reported allosteric activators, demonstrating its potential application in designing specific allosteric drugs and uncovering allosteric mechanisms.

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.

Diffusivelike Motions in a Solvent-Free Protein-Polymer Hybrid

The interaction between proteins and hydration water stabilizes protein structure and promotes functional dynamics, with water translational motions enabling protein flexibility. Engineered solvent-free protein-polymer hybrids have been shown to preserve protein structure, function, and dynamics. Here, we used neutron scattering, protein and polymer perdeuteration, and molecular dynamics simulations to explore how a polymer dynamically replaces water. Even though relaxation rates and vibrational properties are strongly modified in polymer coated compared to hydrated proteins, liquidlike polymer dynamics appear to plasticize the conjugated protein in a qualitatively similar way as do hydration-water translational motions.

Acceleration of catalysis in dihydrofolate reductase by transient, site-specific photothermal excitation

We have studied the role of protein dynamics in chemical catalysis in the enzyme dihydrofolate reductase (DHFR), using a pump-probe method that employs pulsed-laser photothermal heating of a gold nanoparticle (AuNP) to directly excite a local region of the protein structure and transient absorbance to probe the effect on enzyme activity. Enzyme activity is accelerated by pulsed-laser excitation when the AuNP is attached close to a network of coupled motions in DHFR (on the FG loop, containing residues 116-132, or on a nearby alpha helix). No rate acceleration is observed when the AuNP is attached away from the network (distal mutant and His-tagged mutant) with pulsed excitation, or for any attachment site with continuous wave excitation. We interpret these results within an energy landscape model in which transient, site-specific addition of energy to the enzyme speeds up the search for reactive conformations by activating motions that facilitate this search.

The hydrogen-bond configuration modulates the energy transfer efficiency in helical protein nanotubes

Energy transport in proteins is critical to a variety of physical, chemical, and biological processes in living organisms. While strenuous efforts have been made to study vibrational energy transport in proteins, thermal transport processes across the most fundamental building blocks of proteins, i.e. helices, are not well understood. This work studies energy transport in a group of "isomer" helices. The π-helix is shown to have the highest thermal conductivity, 110% higher than that of the α-helix and 207% higher than that of the 310-helix. The H-bond connectivity is found to govern thermal transport mechanisms including the phonon spectral energy density, dispersion, mode-specific transport, group velocity, and relaxation time. The energy transport is strongly correlated with the H-bond strength which is also modulated by the H-bond connectivity. These fundamental insights provide a novel perspective for understanding energy transfer in proteins and guiding a rational molecule-level design of novel materials with configurable H-bonds.

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

To facilitate the observation of biomolecular energy transport in real time and with single-residue resolution, recent experiments by Baumann et al. ( Angew. Chem. Int. Ed. 2019 , 58 , 2899 , DOI: 10.1002/anie.201812995 ) have used unnatural amino acids β-(1-azulenyl)alanine (Azu) and azidohomoalanine (Aha) to site-specifically inject and probe vibrational energy in proteins. To aid the interpretation of such experiments, non-equilibrium molecular dynamics simulations of the anisotropic energy flow in proteins TrpZip2 and PDZ3 domains are presented. On this account, an efficient simulation protocol is established that accurately mimics the excitation and probing steps of Azu and Aha. The simulations quantitatively reproduce the experimentally found cooling times of the solvated proteins at room temperature and predict that the cooling slows by a factor 2 below the glass temperature of water. In PDZ3, vibrational energy is shown to travel from the initially excited peptide ligand via a complex network of inter-residue contacts and backbone transport to distal regions of the protein. The supposed connection of these energy transport pathways with pathways of allosteric communication is discussed.

Computational simulations determining disulfonic stilbene derivative bioavailability within human serum albumin

Disulfonic stilbene (DS) derivatives are a member of the large family of compounds widely employed in medicine and biology as modulators for membrane transporters or inhibitors of a protein involved in DNA repair. They constitute interesting compounds that have not yet been investigated within the bioavailability framework. No crystallographic structures exist involving such compounds embedded in the most common drug carrier, human serum albumin (HSA). The present work studies, for the first time, the physico-chemical features driving the inclusion of three DS derivatives (amino, nitro and acetamido, named DADS, DNDS and DATDS, respectively) within the four common HSA binding sites using combined molecular docking and molecular dynamics simulations. A careful analysis of each ligand within each of the studied binding sites is carried out, highlighting specific interactions and key residues playing a role in stabilizing the ligand within each pocket. The comparison between DADS, DNDS and DATDS reveals that depending on the binding site, the conclusions are rather different. For instance, the IB binding site shows a specificity to DADS compounds while IIIA is the most favorable site for DNDS and DATDS.

Site-Resolved Observation of Vibrational Energy Transfer Using a Genetically Encoded Ultrafast Heater

Allosteric information transfer in proteins has been linked to distinct vibrational energy transfer (VET) pathways in a number of theoretical studies. Experimental evidence for such pathways, however, is sparse because site-selective injection of vibrational energy into a protein, that is, localized heating, is required for their investigation. Here, we solved this problem by the site-specific incorporation of the non-canonical amino acid β-(1-azulenyl)-l-alanine (AzAla) through genetic code expansion. As an exception to Kasha's rule, AzAla undergoes ultrafast internal conversion and heating after S₁ excitation while upon S₂ excitation, it serves as a fluorescent label. We equipped PDZ3, a protein interaction domain of postsynaptic density protein 95, with this ultrafast heater at two distinct positions. We indeed observed VET from the incorporation sites in the protein to a bound peptide ligand on the picosecond timescale by ultrafast IR spectroscopy. This approach based on genetically encoded AzAla paves the way for detailed studies of VET and its role in a wide range of proteins.

Vibrational Energy in Proteins Correlates with Topology

The exchange of vibrational energy in proteins is crucial for their function. Here, we establish a connection between quantities related to it with geometry-based properties such as the proteins' residues coordination number. This relation is proven by molecular simulation in a neuro-pharmacologically relevant transmembrane receptor. The connection demonstrated here paves the way to studies of protein allostery and conformational changes based solely on protein structure.

Allostery and cooperativity in multimeric proteins: bond-to-bond propensities in ATCase

Aspartate carbamoyltransferase (ATCase) is a large dodecameric enzyme with six active sites that exhibits allostery: its catalytic rate is modulated by the binding of various substrates at distal points from the active sites. A recently developed method, bond-to-bond propensity analysis, has proven capable of predicting allosteric sites in a wide range of proteins using an energy-weighted atomistic graph obtained from the protein structure and given knowledge only of the location of the active site. Bond-to-bond propensity establishes if energy fluctuations at given bonds have significant effects on any other bond in the protein, by considering their propagation through the protein graph. In this work, we use bond-to-bond propensity analysis to study different aspects of ATCase activity using three different protein structures and sources of fluctuations. First, we predict key residues and bonds involved in the transition between inactive (T) and active (R) states of ATCase by analysing allosteric substrate binding as a source of energy perturbations in the protein graph. Our computational results also indicate that the effect of multiple allosteric binding is non linear: a switching effect is observed after a particular number and arrangement of substrates is bound suggesting a form of long range communication between the distantly arranged allosteric sites. Second, cooperativity is explored by considering a bisubstrate analogue as the source of energy fluctuations at the active site, also leading to the identification of highly significant residues to the T ↔ R transition that enhance cooperativity across active sites. Finally, the inactive (T) structure is shown to exhibit a strong, non linear communication between the allosteric sites and the interface between catalytic subunits, rather than the active site. Bond-to-bond propensity thus offers an alternative route to explain allosteric and cooperative effects in terms of detailed atomistic changes to individual bonds within the protein, rather than through phenomenological, global thermodynamic arguments.

Photoprogramming Allostery in Human Serum Albumin

Developing strategies to interfere with allosteric interactions in proteins not only promises to deepen our understanding of vital cellular processes but also allows their regulation using external triggers. Light is particularly attractive as a trigger being spatiotemporally selective and compatible with the physiological environment. Here, we engineered a hybrid protein in which irradiation with light opens a new allosteric communication route that is not inherent to the natural system. We select human serum albumin, a promiscuous protein responsible for transporting a variety of ligands in plasma, and show that by covalently incorporating a synthetic photoswitch to subdomain IA we achieve optical control of the ligand binding in subdomain IB. Molecular dynamics simulations confirm the allosteric nature of the interactions between IA and IB in the engineered protein. Specifically, upon illumination, photoconversion of the switch is found to correlate with a less-coordinated motion of the two subdomains and an increased flexibility of the binding pocket in subdomain IB, whose fluctuations are cooperatively enhanced by the presence of ligands, ultimately facilitating their release. Our combined experimental and computational work demonstrates how harnessing artificial molecular switches enables photoprogramming the allosteric regulation of binding activities in such a prominent protein.

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.

Two Photon Spectroscopy Can Serve as a Marker of Protein Denaturation Pathway

Rhodamine group of molecules are widely used dyes for imaging of biological molecules. Application of these dyes however includes a limitation that these molecules absorb in the visible range of the spectrum, which does not fall in the 'biologically transparent window' (BTW). Two photon absorption (TPA) process could come up with an alternate solution to this as these dyes could be excited in the near infrared (NIR) window to extract similar information. To validate this we have investigated TPA cross section (TPACS, σ₂) of two rhodamine dyes, namely Rhodamine 6G (R6G), Rhodamine B (RhB), site selectively bound with a model protein, bovine serum albumin (BSA), by exciting at 800 nm. Two photon spectroscopy and imaging confirms the binding of the dye to the protein. The decreases in TPACS with increasing temperature at a fixed BSA concentration excellently follows the temperature induced structural transition of BSA as the protein transforms from a molten globule to unfolded conformation beyond 60 °C, which has previously been established through circular dichroism (CD) measurements. The thus established resemblance in TPACS and CD measurement trends thus strongly affirms the suitability of TPA process in protein imaging and as an alternative marker to tracking its conformational transformations using NIR radiation.

Spiers Memorial Lecture. Introductory lecture: the impact of structure on photoinduced processes in nucleic acids and proteins

Light is an important environmental variable and most organisms have evolved means to sense, exploit or avoid it and to repair detrimental effects on their genome. In general, light absorption is the task of specific chromophores, however other biomolecules such as oligonucleotides also do so which can result in undesired outcomes such as mutations and cancer. Given the biological importance of light-induced processes and applications for imaging, optogenetics, photodynamic therapy or photovoltaics, there is a great interest in understanding the detailed molecular mechanisms of photoinduced processes in proteins and nucleic acids. The processes are typically characterized by time-resolved spectroscopic approaches or computation, inferring structural information on transient species from stable ground state structures. Recently, however, structure determination of excited states or other short-lived species has become possible with the advent of X-ray free-electron lasers. This review gives an overview of the impact of structure on the understanding of photoinduced processes in macromolecules, focusing on systems presented at this Faraday Discussion meeting.

The effect of surface charge and pH on the physiological behaviour of cobalt, copper, manganese, antimony, zinc and titanium oxide nanoparticles in vitro

The precise knowledge on various interactions of metal nanoparticles (NP) in a living organism is scarce. It is expected that metals can bind to nucleic acids, peptides and proteins (e.g. enzymes), and modify the functioning of vital cellular compartments after entering the organism. The predictive factors for quantitative nanostructure-activity relationship (QNAR) analysis could enhance efficient and harmless usage of nanoparticles (NPs) in the industry as well in the medicine. The studies value the composition of the NP corona determined by time, temperature and source of protein which has been found to implicate the physiological behaviour of NPs. One has largely been ignored: the NPs specific isoelectric point (IEP) and pH at the state of measurement. Herein, this study investigates the effect of pH and surface charge of six metal oxide (MeOx) NPs on time dependency of cytotoxicity. Several aspects of the characterization of ultrafine particles in the actual test system which is the most relevant for the interpretation of the toxicological data are referred: (i) the difference of pH in the room temperature and in the incubation conditions (ii) the difference of dispersions in MilliQ and complete cell media; (iii) the need to exemplify also the pH and isoelectric point when the hydrodynamic size is measured; (iv) the importance of time due to the time-dependent equilibration and changes of NPs corona. The surface charge determines the formation of corona and could be modified by pH. MeOx NPs without fully charge equilibrated corona might play the main role of MeOx NPs entering into the cell and consequently the time dependent manifestation of the cellular effect.

The Wisdom of Networks: A General Adaptation and Learning Mechanism of Complex Systems: The Network Core Triggers Fast Responses to Known Stimuli; Innovations Require the Slow Network Periphery and Are Encoded by Core-Remodeling

I hypothesize that re-occurring prior experience of complex systems mobilizes a fast response, whose attractor is encoded by their strongly connected network core. In contrast, responses to novel stimuli are often slow and require the weakly connected network periphery. Upon repeated stimulus, peripheral network nodes remodel the network core that encodes the attractor of the new response. This "core-periphery learning" theory reviews and generalizes the heretofore fragmented knowledge on attractor formation by neural networks, periphery-driven innovation, and a number of recent reports on the adaptation of protein, neuronal, and social networks. The core-periphery learning theory may increase our understanding of signaling, memory formation, information encoding and decision-making processes. Moreover, the power of network periphery-related "wisdom of crowds" inventing creative, novel responses indicates that deliberative democracy is a slow yet efficient learning strategy developed as the success of a billion-year evolution. Also see the video abstract here: https://youtu.be/IIjP7zWGjVE.

A Universal Pattern in the Percolation and Dissipation of Protein Structural Perturbations

Understanding the extent to which information is transmitted through the intramolecular interaction network of proteins upon a perturbation, that is, an allosteric effect, has long remained an unsolved problem. Through an analysis of high-resolution NMR data from the literature on 28 different proteins and 49 structural perturbations, we show that the extent of induced structural changes through mutations and molecular events including protein-protein, protein-peptide, protein-ligand binding, and post-translational modifications exhibit a near-universal exponential functional form. The extent of percolation into the protein structures can be up to 20-25 Å despite no apparent change in the 3D structures. These observations are also consistent with theoretical expectations, elementary graph theoretic analysis of protein structures, detailed molecular dynamics simulations, and experimental double-mutant cycles. Our analysis highlights that most molecular events would contribute to allosteric effects independent of protein structure, topology, or identity and provides a simple avenue to test and potentially model their effects.

Ultrafast Vibrational Dynamics of Membrane-Bound Peptides at the Lipid Bilayer/Water Interface

Vibrational energy transfer (VET) of proteins at cell membrane plays critical roles in controlling the protein functionalities, but its detection is very challenging. By using a surface-sensitive femtosecond time-resolved sum-frequency generation vibrational spectroscopy with infrared pump, the detection of the ultrafast VET in proteins at cell membrane has finally become possible. The vibrational relaxation time of the N-H groups is determined to be 1.70(±0.05) ps for the α-helix located in the hydrophobic core of the lipid bilayer and 0.9(±0.05) ps for the membrane-bound β-sheet structure. The N-H groups with strong hydrogen bonding gain faster relaxation time. By pumping the amide A band and probing amide I band, the vibrational relaxation from N-H mode to C=O mode through two pathways (direct coupling and through intermediate states) is revealed. The ratio of the pathways depends on the NH⋅⋅⋅O=C hydrogen-bonding strength. Strong hydrogen bonding favors the coupling through intermediate states.

What Mutagenesis Can and Cannot Reveal About Allostery

Allosteric regulation of protein function is recognized to be widespread throughout biology; however, knowledge of allosteric mechanisms, the molecular changes within a protein that couple one binding site to another, is limited. Although mutagenesis is often used to probe allosteric mechanisms, we consider herein what the outcome of a mutagenesis study truly reveals about an allosteric mechanism. Arguably, the best way to evaluate the effects of a mutation on allostery is to monitor the allosteric coupling constant (Qax), a ratio of the substrate binding constants in the absence versus presence of an allosteric effector. A range of substitutions at a given residue position in a protein can reveal when a particular substitution causes gain-of-function, which addresses a key challenge in interpreting mutation-dependent changes in the magnitude of Qax. Thus, whole-protein mutagenesis studies offer an acceptable means of identifying residues that contribute to an allosteric mechanism. With this focus on monitoring Qax, and keeping in mind the equilibrium nature of allostery, we consider alternative possibilities for what an allosteric mechanism might be. We conclude that different possible mechanisms (rotation-of-solid-domains, movement of secondary structure, side-chain repacking, changes in dynamics, etc.) will result in different findings in whole-protein mutagenesis studies.

Comparative Study of the Collective Dynamics of Proteins and Inorganic Nanoparticles

Molecular dynamics simulations of ubiquitin in water/glycerol solutions are used to test the suggestion by Karplus and coworkers that proteins in their biologically active state should exhibit a dynamics similar to 'surface-melted' inorganic nanoparticles (NPs). Motivated by recent studies indicating that surface-melted inorganic NPs are in a 'glassy' state that is an intermediate dynamical state between a solid and liquid, we probe the validity and significance of this proposed analogy. In particular, atomistic simulations of ubiquitin in solution based on CHARMM36 force field and pre-melted Ni NPs (Voter-Chen Embedded Atom Method potential) indicate a common dynamic heterogeneity, along with other features of glass-forming (GF) liquids such as collective atomic motion in the form of string-like atomic displacements, potential energy fluctuations and particle displacements with long range correlations ('colored' or 'pink' noise), and particle displacement events having a power law scaling in magnitude, as found in earthquakes. On the other hand, we find the dynamics of ubiquitin to be even more like a polycrystalline material in which the α-helix and β-sheet regions of the protein are similar to crystal grains so that the string-like collective atomic motion is concentrated in regions between the α-helix and β-sheet domains.

Prediction of allosteric sites and mediating interactions through bond-to-bond propensities

Allostery is a fundamental mechanism of biological regulation, in which binding of a molecule at a distant location affects the active site of a protein. Allosteric sites provide targets to fine-tune protein activity, yet we lack computational methodologies to predict them. Here we present an efficient graph-theoretical framework to reveal allosteric interactions (atoms and communication pathways strongly coupled to the active site) without a priori information of their location. Using an atomistic graph with energy-weighted covalent and weak bonds, we define a bond-to-bond propensity quantifying the non-local effect of instantaneous bond fluctuations propagating through the protein. Significant interactions are then identified using quantile regression. We exemplify our method with three biologically important proteins: caspase-1, CheY, and h-Ras, correctly predicting key allosteric interactions, whose significance is additionally confirmed against a reference set of 100 proteins. The almost-linear scaling of our method renders it suitable for high-throughput searches for candidate allosteric sites.

Human plasma lipocalins and serum albumin: Plasma alternative carriers?

Lipocalins are an evolutionarily conserved family of proteins that bind and transport a variety of exogenous and endogenous ligands. Lipocalins share a conserved eight anti-parallel β-sheet structure. Among the different lipocalins identified in humans, α-1-acid glycoprotein (AGP), apolipoprotein D (apoD), apolipoprotein M (apoM), α1-microglobulin (α1-m) and retinol-binding protein (RBP) are plasma proteins. In particular, AGP is the most important transporter for basic and neutral drugs, apoD, apoM, and RBP mainly bind endogenous molecules such as progesterone, pregnenolone, bilirubin, sphingosine-1-phosphate, and retinol, while α1-m binds the heme. Human serum albumin (HSA) is a monomeric all-α protein that binds endogenous and exogenous molecules like fatty acids, heme, and acidic drugs. Changes in the plasmatic levels of lipocalins and HSA are responsible for the onset of pathological conditions associated with an altered drug transport and delivery. This, however, does not necessary result in potential adverse effects in patients because many drugs can bind both HSA and lipocalins, and therefore mutual compensatory binding mechanisms can be hypothesized. Here, molecular and clinical aspects of ligand transport by plasma lipocalins and HSA are reviewed, with special attention to their role as alterative carriers in health and disease.

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.

Using THz Spectroscopy, Evolutionary Network Analysis Methods, and MD Simulation to Map the Evolution of Allosteric Communication Pathways in c-Type Lysozymes

It is now widely accepted that protein function is intimately tied with the navigation of energy landscapes. In this framework, a protein sequence is not described by a distinct structure but rather by an ensemble of conformations. And it is through this ensemble that evolution is able to modify a protein's function by altering its landscape. Hence, the evolution of protein functions involves selective pressures that adjust the sampling of the conformational states. In this work, we focus on elucidating the evolutionary pathway that shaped the function of individual proteins that make-up the mammalian c-type lysozyme subfamily. Using both experimental and computational methods, we map out specific intermolecular interactions that direct the sampling of conformational states and accordingly, also underlie shifts in the landscape that are directly connected with the formation of novel protein functions. By contrasting three representative proteins in the family we identify molecular mechanisms that are associated with the selectivity of enhanced antimicrobial properties and consequently, divergent protein function. Namely, we link the extent of localized fluctuations involving the loop separating helices A and B with shifts in the equilibrium of the ensemble of conformational states that mediate interdomain coupling and concurrently moderate substrate binding affinity. This work reveals unique insights into the molecular level mechanisms that promote the progression of interactions that connect the immune response to infection with the nutritional properties of lactation, while also providing a deeper understanding about how evolving energy landscapes may define present-day protein function.

Direct observation of ultrafast collective motions in CO myoglobin upon ligand dissociation

The hemoprotein myoglobin is a model system for the study of protein dynamics. We used time-resolved serial femtosecond crystallography at an x-ray free-electron laser to resolve the ultrafast structural changes in the carbonmonoxy myoglobin complex upon photolysis of the Fe-CO bond. Structural changes appear throughout the protein within 500 femtoseconds, with the C, F, and H helices moving away from the heme cofactor and the E and A helices moving toward it. These collective movements are predicted by hybrid quantum mechanics/molecular mechanics simulations. Together with the observed oscillations of residues contacting the heme, our calculations support the prediction that an immediate collective response of the protein occurs upon ligand dissociation, as a result of heme vibrational modes coupling to global modes of the protein.

Allosteric communication pathways and thermal rectification in PDZ-2 protein: a computational study

Allosteric communication in proteins is a fundamental and yet unresolved problem of structural biochemistry. Previous findings, from computational biology ( Ota, N.; Agard, D. A. J. Mol. Biol. 2005 , 351 , 345 - 354 ), have proposed that heat diffuses in a protein through cognate protein allosteric pathways. This work studied heat diffusion in the well-known PDZ-2 protein, and confirmed that this protein has two cognate allosteric pathways and that heat flows preferentially through these. Also, a new property was also observed for protein structures: heat diffuses asymmetrically through the structures. The underling structure of this asymmetrical heat flow was a normal length hydrogen bond (∼2.85 Å) that acted as a thermal rectifier. In contrast, thermal rectification was compromised in short hydrogen bonds (∼2.60 Å), giving rise to symmetrical thermal diffusion. Asymmetrical heat diffusion was due, on a higher scale, to the local, structural organization of residues that, in turn, was also mediated by hydrogen bonds. This asymmetrical/symmetrical energy flow may be relevant for allosteric signal communication directionality in proteins and for the control of heat flow in materials science.

Fast infrared spectroscopy of protein dynamics: advancing sensitivity and selectivity

2D-IR spectroscopy has matured to a powerful technique to study the structure and dynamics of peptides, but its extension to larger proteins is still in its infancy, the major limitations being sensitivity and selectivity. Site-selective information requires measuring single vibrational probes at sub-millimolar concentrations where most proteins are still stable, which is a severe challenge for conventional (FT)IR spectroscopy. Besides its ultrafast time-resolution, a so far largely underappreciated potential of 2D-IR spectroscopy lies in its sensitivity gain. The present paper sets the goals and outlines strategies how to use that sensitivity gain together with properly designed vibrational labels to make IR spectroscopy a versatile tool to study a wide class of proteins.

Ultrafast myoglobin structural dynamics observed with an X-ray free-electron laser

Light absorption can trigger biologically relevant protein conformational changes. The light-induced structural rearrangement at the level of a photoexcited chromophore is known to occur in the femtosecond timescale and is expected to propagate through the protein as a quake-like intramolecular motion. Here we report direct experimental evidence of such 'proteinquake' observed in myoglobin through femtosecond X-ray solution scattering measurements performed at the Linac Coherent Light Source X-ray free-electron laser. An ultrafast increase of myoglobin radius of gyration occurs within 1 picosecond and is followed by a delayed protein expansion. As the system approaches equilibrium it undergoes damped oscillations with a ~3.6-picosecond time period. Our results unambiguously show how initially localized chemical changes can propagate at the level of the global protein conformation in the picosecond timescale.

Observing Vibrational Energy Flow in a Protein with the Spatial Resolution of a Single Amino Acid Residue

One of the challenges in physical chemistry has been understanding how energy flows in a condensed phase from the microscopic viewpoint. To address this, space-resolved information at the molecular scale is required but has been lacking due to experimental difficulties. We succeeded in the real-time mapping of the vibrational energy flow in a protein with the spatial resolution of a single amino acid residue by combining time-resolved resonance Raman spectroscopy and site-directed single-Trp mutagenesis. Anti-Stokes Raman intensities of the Trp residues at different sites exhibited different temporal evolutions, reflecting propagation of the energy released by the heme group. A classical heat transport model was not able to reproduce the entire experimental data set, showing that we need a molecular-level description to explain the energy flow in a protein. The systematic application of our general methodology to proteins with different structural motifs may provide a greatly increased understanding of the energy flow in proteins.