Response of villin headpiece-capped gold nanoparticles to ultrafast laser heating

Plasmonic-Driven Regulation of Biomolecular Activity In Situ

Selective and remote manipulation of activity for biomolecules, including protein, DNA, and lipids, is crucial to elucidate their molecular function and to develop biomedical applications. While advances in tool development, such as optogenetics, have significantly impacted these directions, the requirement for genetic modification significantly limits their therapeutic applications. Plasmonic nanoparticle heating has brought new opportunities to the field, as hot nanoparticles are unique point heat sources at the nanoscale. In this review, we summarize fundamental engineering problems such as plasmonic heating and the resulting biomolecular responses. We highlight the biological responses and applications of manipulating biomolecules and provide perspectives for future directions in the field.

Temporal Resolution of Activity-Related Solvation Dynamics in the TIM Barrel Enzyme Murine Adenosine Deaminase

Murine adenosine deaminase (mADA) is a prototypic system for studying the thermal activation of active site chemistry within the TIM barrel family of enzyme reactions. Previous temperature-dependent hydrogen deuterium exchange studies under various conditions have identified interconnected thermal networks for heat transfer from opposing protein-solvent interfaces to active site residues in mADA. One of these interfaces contains a solvent exposed helix-loop-helix moiety that presents the hydrophobic face of its long α-helix to the backside of bound substrate. Herein we pursue the time and temperature dependence of solvation dynamics at the surface of mADA, for comparison to established kinetic parameters that represent active site chemistry. We first created a modified protein devoid of native tryptophans with close to native kinetic behavior. Single site-specific tryptophan mutants were back inserted into each of the four positions where native tryptophans reside. Measurements of nanosecond fluorescence relaxation lifetimes and Stokes shift decays, that reflect time dependent environmental reoroganization around the photo-excited state of Trp*, display minimal temperature dependences. These regions serve as controls for the behavior of a new single tryptophan inserted into a solvent exposed region near the helix-loop-helix moiety located behind the bound substrate, Lys54Trp. This installed Trp displays a significantly elevated value for Ea ( k Stokes shift ) ; further, when Phe61 within the long helix positioned behind bound substrate is replaced by a series of aliphatic hydrophobic side chains, the trends in Ea ( k Stokes shift ) mirror the earlier reported impact of the same series of function-altering hydrophobic side chains on the activation energy of catalysis, Ea ( k cat ) .The reported experimental findings implicate a solvent initiated and rapid (>ns) protein restructuring that contributes to the enthalpic activation barrier to catalysis in mADA.

Glutathione Self-Assembles into a Shell of Hydrogen-Bonded Intermolecular Aggregates on "Naked" Silver Nanoparticles

A detailed understanding of the molecular structure in nanoparticle ligand capping layers is crucial for their efficient incorporation into modern scientific and technological applications. Peptide ligands render the nanoparticles as biocompatible materials. Glutathione, a γ-ECG tripeptide, self-assembles into aggregates on the surface of ligand-free silver nanoparticles through intermolecular hydrogen bonding and forms a few nanometer-thick shells. Two-dimensional nonlinear infrared (2DIR) spectroscopy suggests that aggregates adopt a conformation resembling the β-sheet secondary structure. The shell thickness was evaluated with localized surface plasmon resonance spectroscopy and X-ray photoelectron spectroscopy. The amount of glutathione on the surface was obtained with spectrophotometry of a thiol-reactive probe. Our results suggest that the shell consists of ∼15 stacked molecular layers. These values correspond to the inter-sheet distances, which are significantly shorter than those in amyloid fibrils with relatively bulky side chains, but are comparable to glycine-rich silk fibrils, where the side chains are compact. The tight packing of the glutathione layers can be facilitated by hydrogen-bonded carboxylic acid dimers of glycine and the intermolecular salt bridges between the zwitterionic γ-glutamyl groups. The structure of the glutathione aggregates was studied by 2DIR spectroscopy of the amide-I vibrational modes using ¹³C isotope labeling of the cysteine carbonyl. Isotope dilution experiments revealed the coupling of modes forming vibrational excitons along the cysteine chain. The coupling along the γ-glutamyl exciton chain was estimated from these values. The obtained coupling strengths are slightly lower than those of native β-sheets, yet they appear large enough to point onto an ordered conformation of the peptides within the aggregate. Analysis of the excitons' anharmonicities and the strength of the transition dipole moments generally is in agreement with these observations.

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.

Non-Arrhenius Reaction-Diffusion Kinetics for Protein Inactivation over a Large Temperature Range

Understanding protein folding and unfolding has been a long-standing fundamental question and has important applications in manipulating protein activity in biological systems. Experimental investigations of protein unfolding have been predominately conducted by small temperature perturbations (e.g., temperature jump), while molecular simulations are limited to small time scales (microseconds) and high temperatures to observe unfolding. Thus, it remains unclear how fast a protein unfolds irreversibly and loses function (i.e., inactivation) across a large temperature range. In this work, using nanosecond pulsed heating of individual plasmonic nanoparticles to create precise localized heating, we examine the protein inactivation kinetics at extremely high temperatures. Connecting this with protein inactivation measurements at low temperatures, we observe that the kinetics of protein unfolding is less sensitive to temperature change at the higher temperatures, which significantly departs from the Arrhenius behavior extrapolated from low temperatures. To account for this effect, we propose a reaction-diffusion model that modifies the temperature-dependence of protein inactivation by introducing a diffusion limit. Analysis of the reaction-diffusion model provides general guidelines in the behavior of protein inactivation (reaction-limited, transition, diffusion-limited) across a large temperature range from physiological temperature to extremely high temperatures. We further demonstrate that the reaction-diffusion model is particularly useful for designing optimal operating conditions for protein photoinactivation. The experimentally validated reaction-diffusion kinetics of protein unfolding is an important step toward understanding protein-inactivation kinetics over a large temperature range. It has important applications including molecular hyperthermia and calls for future studies to examine this model for other protein molecules.

Intermolecular energy flows between surface molecules on metal nanoparticles

Three model systems are designed to investigate energy transport between molecules on metal nanoparticle surfaces. Energy is rapidly transferred from one carbon monoxide (CO) molecule to another CO molecule or an organic molecule on adjacent surface sites of 2 nm Pt particles within a few picoseconds. On the contrary, energy flow from a surface organic molecule to an adjacent CO molecule is significantly slower and, in fact, within experimental sensitivity and uncertainty the transfer is not observed. The energy transport on particle surfaces (about 2 km s-1) is almost ten times faster than inside a molecule (200 m s-1). The seemingly perplexing observations can be well explained by the combination of electron/vibration and vibration/vibration coupling mechanisms, which mediate molecular energy dynamics on metal nanoparticle surfaces: the strong electron/vibration coupling rapidly converts CO vibrational energy into heat that can be immediately sensed by nearby molecules; but the vibration/vibration coupling dissipates the vibrational excitation in the organic molecule as low-frequency intramolecular vibrations that may or may not couple to surface electronic motions.

Characterizing the Surface Coverage of Protein-Gold Nanoparticle Bioconjugates

Functional enzyme-nanoparticle bioconjugates are increasingly important in biomedical and biotechnology applications such as drug delivery and biosensing. Optimization of the function of such bioconjugates requires careful control and characterization of their structures and activity, but current methods are inadequate for this purpose. A key shortcoming of existing approaches is the lack of an accurate method for quantitating protein content of bioconjugates for low (monolayer) surface coverages. In this study, an integrated characterization methodology for protein-gold nanoparticle (AuNP) bioconjugates is developed, with a focus on site-specific attachment and surface coverage of protein on AuNPs. Single-cysteine-containing mutants of dihydrofolate reductase are covalently attached to AuNPs with diameters of 5, 15, and 30 nm, providing a range of surface curvature. Site-specific attachment to different regions of the protein surface is investigated, including attachment to a flexible loop versus a rigid α helix. Characterization methods include SDS-PAGE, UV-vis spectrophotometry, dynamic light scattering, and a novel fluorescence-based method for accurate determination of low protein concentration on AuNPs. An accurate determination of both protein and AuNP concentration in conjugate samples allows for the calculation of the surface coverage. We find that surface coverage is related to the surface curvature of the AuNP, with a higher surface coverage observed for higher surface curvature. The combination of these characterization methods is important for understanding the functionality of protein-AuNP bioconjugates, particularly enzyme activity.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Folding stability modulation of the villin headpiece helical subdomain by 4-fluorophenylalanine and 4-methylphenylalanine

HP36, the helical subdomain of villin headpiece, contains a hydrophobic core composed of three phenylalanine residues (Phe47, Phe51, and Phe58). Hydrophobic effects and electrostatic interactions were shown to be the critical factors in stabilizing this core and the global structure. To assess the interactions among Phe47, Phe51, and Phe58 residues and investigate how they affect the folding stability, we implanted 4-fluorophenylalanine (Z) and 4-methylphenylalanine (X) into the hydrophobic core of HP36. We chemically synthesized HP36 and its seven variants including four single mutants whose Phe51 or Phe58 was replaced with Z or X, and three double mutants whose Phe51 and Phe58 were both substituted. Circular dichroism and nuclear magnetic resonance measurements show that the variants exhibit a native HP36 like fold, of which F51Z and three double mutants are more stable than the wild type. Molecular modeling provided detailed interaction energy within the phenylalanine residues, revealing that electrostatic interactions dominate the stability modulation upon the introduction of 4-fluorophenylalanine and 4-methylphenylalanine. Our results show that these two non-natural amino acids can successfully tune the interactions in a relatively compact hydrophobic core and the folding stability without inducing dramatic steric effects. Such an approach may be applied to other folded motifs or proteins.

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.

Protein dynamics in the solid state from 2H NMR line shape analysis: a consistent perspective

Deuterium line shape analysis of CD3 groups has emerged as a particularly useful tool for studying microsecond-millisecond protein motions in the solid state. The models devised so far consist of several independently conceived simple jump-type motions. They are comprised of physical quantities encoded in their simplest form; improvements are only possible by adding yet another simple motion, thereby changing the model. The various treatments developed are case-specific; hence comparison among the different systems is not possible. Here we develop a new methodology for (2)H NMR line shape analysis free of these limitations. It is based on the microscopic-order-macroscopic-disorder (MOMD) approach. In MOMD motions are described by diffusion tensors, spatial restrictions by potentials/ordering tensors, and geometric features by relative tensor orientations. Jump-type motions are recovered in the limit of large orientational potentials. Model improvement is accomplished by monitoring the magnitude, symmetry, and orientation of the various tensors. The generality of MOMD makes possible comparison among different scenarios. CD3 line shapes from the Chicken Villin Headpiece Subdomain and the Streptomyces Subtilisin Inhibitor are used as experimental examples. All of these spectra are reproduced by using rhombic local potentials constrained for simplicity to be given by the L = 2 spherical harmonics, and by axial diffusion tensors. Potential strength and rhombicity are found to be ca. 2-3 k(B)T. The diffusion tensor is tilted at 120° from the C-CD3 axis. The perpendicular (parallel) correlation times for local motion are 0.1-1.0 ms (3.3-30 μs). Activation energies in the 1.1-8.0 kcal/mol range are estimated. Future prospects include extension to the (2)H relaxation limit, application to the (15)N and (13)C NMR nuclei, and accounting for collective motions and anisotropic media.