Characterization of Low-Frequency Modes in Aqueous Peptides Using Far-Infrared Spectroscopy and Molecular Dynamics Simulation

Real-time measure of solvation free energy changes upon liquid-liquid phase separation of α-elastin

Biological condensates are known to retain a large fraction of water to remain in a liquid and reversible state. Local solvation contributions from water hydrating hydrophilic and hydrophobic protein surfaces were proposed to play a prominent role for the formation of condensates through liquid-liquid phase separation (LLPS). However, although the total free energy is accessible by calorimetry, the partial solvent contributions to the free energy changes upon LLPS remained experimentally inaccessible so far. Here, we show that the recently developed THz calorimetry approach allows to quantify local hydration enthalpy and entropy changes upon LLPS of α-elastin in real time, directly from experimental THz spectroscopy data. We find that hydrophobic solvation dominates the entropic solvation term, whereas hydrophilic solvation mainly contributes to the enthalpy. Both terms are in the order of hundreds of kJ/mol, which is more than one order of magnitude larger than the total free energy changes at play during LLPS. However, since we show that entropy/enthalpy mostly compensates, a small entropy/enthalpy imbalance is sufficient to tune LLPS. Theoretically, a balance was proposed before. Here we present experimental evidence based on our spectroscopic approach. We finally show that LLPS can be steered by inducing small changes of solvation entropy/enthalpy compensation via concentration or temperature in α-elastin.

Nanoantenna enhanced terahertz interaction of biomolecules

Terahertz time-domain spectroscopy (THz-TDS) is a non-invasive, non-contact and label-free technique for biological and chemical sensing as THz-spectra are less energetic and lie in the characteristic vibration frequency regime of proteins and DNA molecules. However, THz-TDS is less sensitive for the detection of micro-organisms of size equal to or less than λ/100 (where, λ is the wavelength of the incident THz wave), and molecules in extremely low concentration solutions (like, a few femtomolar). After successful high-throughput fabrication of nanostructures, nanoantennas were found to be indispensable in enhancing the sensitivity of conventional THz-TDS. These nanostructures lead to strong THz field enhancement when in resonance with the absorption spectrum of absorptive molecules, causing significant changes in the magnitude of the transmission spectrum, therefore, enhancing the sensitivity and allowing the detection of molecules and biomaterials in extremely low concentration solutions. Herein, we review the recent developments in ultra-sensitive and selective nanogap biosensors. We have also provided an in-depth review of various high-throughput nanofabrication techniques. We also discussed the physics behind the field enhancements in the sub-skin depth as well as sub-nanometer sized nanogaps. We introduce finite-difference time-domain (FDTD) and molecular dynamics (MD) simulation tools to study THz biomolecular interactions. Finally, we provide a comprehensive account of nanoantenna enhanced sensing of viruses (like, H1N1) and biomolecules such as artificial sweeteners which are addictive and carcinogenic.

Application of far-infrared spectroscopy to the structural identification of protein materials

Far-IR spectroscopy was applied to monitor the structure of two types of silk fibroins and the results indicate that they both show several sharp characteristic peaks, which are totally different from those of globular proteins.

Peptide-Protein Binding Investigated by Far-IR Spectroscopy and Molecular Dynamics Simulations

Molecular dynamics (MD) simulations and far-infrared (far-IR) spectroscopy were combined to study peptide binding by the second PDZ domain (PDZ1) of MAGI1, which has been identified as an important target for the Human Papilloma Virus. PDZ1 recognizes and binds to the C-terminal end of the E6 protein from high-risk Human Papilloma Virus. The far-IR spectra of two forms of the protein, an unbound APO form and a HOLO form (where the PDZ1 is bound to an 11-residue peptide derived from the C terminus of HPV16 E6), were obtained. MD simulations were used to determine the most representative structure of each form and these were used to compute their respective IR spectra by normal mode analysis. Far-UV circular dichroism spectroscopy was used to confirm the secondary structure content and the stability through temperature-dependent studies. Both the experimental and calculated far-IR spectra showed a red shift of the low-frequency peaks upon peptide binding. The calculations show that this is coincident with an increased number of hydrogen bonds formed as the peptide augments the protein β-sheet. We further identified the contribution of surface-bound water molecules to bands in the far-IR and, through the calculations, identified potential pathways for allosteric communication. Together, these results demonstrate the utility of combining far-IR experiments and MD studies to study peptide binding by proteins.

Silk: Optical Properties over 12.6 Octaves THz-IR-Visible-UV Range

Domestic (Bombyx mori) and wild (Antheraea pernyi) silk fibers were characterised over a wide spectral range from THz 8 cm -1 ( λ = 1.25 mm, f = 0.24 THz) to deep-UV 50 × 10 3 cm - 1 ( λ = 200 nm, f = 1500 THz) wavelengths or over a 12.6 octave frequency range. Spectral features at β-sheet, α-coil and amorphous fibroin were analysed at different spectral ranges. Single fiber cross sections at mid-IR were used to determine spatial distribution of different silk constituents and revealed an α-coil rich core and more broadly spread β-sheets in natural silk fibers obtained from wild Antheraea pernyi moths. Low energy T-ray bands at 243 and 229 cm -1 were observed in crystalline fibers of domestic and wild silk fibers, respectively, and showed no spectral shift down to 78 K temperature. A distinct 20±4 cm-1 band was observed in the crystalline Antheraea pernyi silk fibers. Systematic analysis and assignment of the observed spectral bands is presented. Water solubility and biodegradability of silk, required for bio-medical and sensor applications, are directly inferred from specific spectral bands.

Far-infrared amide IV-VI spectroscopy of isolated 2- and 4-Methylacetanilide

Delocalized molecular vibrations in the far-infrared and THz ranges are highly sensitive to the molecular structure, as well as to intra- and inter-molecular interactions. Thus, spectroscopic studies of biomolecular structures can greatly benefit from an extension of the conventional mid-infrared to the far-infrared wavelength range. In this work, the conformer-specific gas-phase far-infrared spectra of two aromatic molecules containing the peptide -CO-NH- link, namely, 2- and 4-Methylacetanilide, are investigated. The planar conformations with trans configuration of the peptide link have only been observed in the supersonic-jet expansion. The corresponding far-infrared signatures associated with the vibrations of the peptide -CO-NH- moiety, the so-called amide IV-VI bands, have been assigned and compared with the results of density functional theory frequency calculations based on the anharmonic vibrational second-order perturbation theory approach. The analysis of the experimental and theoretical data shows that the amide IV-VI bands are highly diagnostic for the geometry of the peptide moiety and the molecular backbone. They are also strongly blue-shifted upon formation of the NH⋯O-C hydrogen bonding, which is, for example, responsible for the formation of secondary protein structures. Furthermore, the amide IV-VI bands are also diagnostic for the cis configuration of the peptide link, which can be present in cyclic peptides. The experimental gas-phase data presented in this work can assist the vibrational assignment of similar biologically important systems, either isolated or in natural environments.

THz frequency spectrum of protein-solvent interaction energy using a recurrence plot-based Wiener-Khinchin method

The dynamics of a protein and the water surrounding it are coupled via nonbonded energy interactions. This coupling can exhibit a complex, nonlinear, and nonstationary nature. The THz frequency spectrum for this interaction energy characterizes both the vibration spectrum of the water hydrogen bond network, and the frequency range of large amplitude modes of proteins. We use a Recurrence Plot based Wiener-Khinchin method RPWK to calculate this spectrum, and the results are compared to those determined using the classical auto-covariance-based Wiener-Khinchin method WK. The frequency spectra for the total nonbonded interaction energy extracted from molecular dynamics simulations between the β-Lactamase Inhibitory Protein BLIP, and water molecules within a 10 Å distance from the protein surface, are calculated at 150, 200, 250, and 310 K, respectively. Similar calculations are also performed for the nonbonded interaction energy between the residues 49ASP, 53TYR, and 142PHE in BLIP, with water molecules within 10 Å from each residue respectively at 150, 200, 250, and 310 K. A comparison of the results shows that RPWK performs better than WK, and is able to detect some frequency data points that WK fails to detect. This points to the importance of using methods capable of taking the complex nature of the protein-solvent energy landscape into consideration, and not to rely on standard linear methods. In general, RPWK can be a valuable addition to the analysis tools for protein molecular dynamics simulations. Proteins 2016; 84:1549-1557. © 2016 Wiley Periodicals, Inc.

Formulations for microprojection/microneedle vaccine delivery: Structure, strength and release profiles

To develop novel methods for vaccine delivery, the skin is viewed as a high potential target, due to the abundance of immune cells that reside therein. One method, the use of dissolving microneedle technologies, has the potential to achieve this, with a range of formulations now being employed. Within this paper we assemble a range of methods (including FT-FIR using synchrotron radiation, nanoindentation and skin delivery assays) to systematically examine the effect of key bulking agents/excipients - sugars/polyols - on the material form, structure, strength, failure properties, diffusion and dissolution for dissolving microdevices. We investigated concentrations of mannitol, sucrose, trehalose and sorbitol from 1:1 to 30:1 with carboxymethylcellulose (CMC), although mannitol did not form our micro-structures so was discounted early in the study. The other formulations showed a variety of crystalline (sorbitol) and amorphous (sucrose, trehalose) structures, when investigated using Fourier transform far infra-red (FT-FIR) with synchrotron radiation. The crystalline structures had a higher elastic modulus than the amorphous formulations (8-12GPa compared to 0.05-11GPa), with sorbitol formulations showing a bimodal distribution of results including both amorphous and crystalline behaviour. In skin, diffusion properties were similar among all formulations with dissolution occurring within 5s for our small projection array structures (~100μm in length). Overall, slight variations in formulation can significantly change the ability of our projections to perform their required function, making the choice of bulking/vaccine stabilising agents of great importance for these devices.

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.

Far-infrared spectroscopy of salt penetration into a collagen fiber scaffold

We employed far-infrared spectroscopy to observe the amount of salt that penetrates into collagen fiber masses. The absorption properties of collagen sheets prepared from tilapia skin, bovine skin, rat tail, and sea cucumber dermis were measured using a transmission Fourier transform spectrometer in a band from approximately 100 to 700 cm(-1). We confirmed that the absorbance spectra of the four types of dried collagen sheet show good agreement, even though the amino acid compositions differed. The absorbance peaks observed in the band corresponded to collective vibrations of plural functional groups such as methylene and imino groups in collagen. When salt solution was added to the collagen sheets and then dried, the spectral shapes of the sheets at approximately 166 cm(-1) were clearly different from those of the plain collagen sheets. The differential absorbance between wavenumbers 166 cm(-1) and 250 cm(-1) sensitively reflected the difference between higher-order structures, and the salt diffusion (crystallization) depended on the collagen fiber condition. From these results, we consider that spectral changes can be used for the numerical evaluation of salt penetration into a collagen fiber scaffold.

The hydration of serine: multipole moments versus point charges

Next-generation force fields must incorporate improved electrostatic potentials in order to increase the reliability of their predictions. A crucial decision toward this goal is to abandon point charges in favour of multipole moments centered on nuclear sites. Here we compare the geometries generated by quantum topological multipole moments with those generated by four popular point charge models (TAFF, OPLS-AA, MMFF94x and PFROSST) for a hydrated serine. A main feature of this study is the dual comparison made, both at static level (geometry optimisation via energy minimisation) and at dynamic level (via molecular dynamics and radial/spatial distribution function analysis). At static level, multipolar electrostatics best reproduces the ab initio reference geometry. At dynamic level, multipolar electrostatics produces more structure than point charge electrostatics does, over the whole range. From our previous work on liquid water [Int. J. Quantum. Chem., 2004, 99, 685], where agreement with experiment only occurs when using multipole moments, we deduce that our predictions for hydrated serine will also be closer to experiment when using multipolar electrostatics. The spatial distribution function shows that only multipolar electrostatics shows pronounced structure at long range. Even at short range there are many regions where waters appear in the system governed by multipolar electrostatics but not in that governed by point charges.

The glassy state of crambin and the THz time scale protein-solvent fluctuations possibly related to protein function

BACKGROUND

THz experiments have been used to characterize the picosecond time scale fluctuations taking place in the model, globular protein crambin.

RESULTS

Using both hydration and temperature as an experimental parameter, we have identified collective fluctuations (<= 200 cm(-1)) in the protein. Observation of the protein dynamics in the THz spectrum from both below and above the glass transition temperature (Tg) has provided unique insight into the microscopic interactions and modes that permit the solvent to effectively couple to the protein thermal fluctuations.

CONCLUSIONS

Our findings suggest that the solvent dynamics on the picosecond time scale not only contribute to protein flexibility but may also delineate the types of fluctuations that are able to form within the protein structure.

Using THz time-scale infrared spectroscopy to examine the role of collective, thermal fluctuations in the formation of myoglobin allosteric communication pathways and ligand specificity

In this investigation we use THz time-scale spectroscopy to conduct an initial set of studies on myoglobin with the aim of providing further insight into the global, collective thermal fluctuations in the protein that have been hypothesized to play a prominent role in the dynamic formation of transient ligand channels as well as in shaping the molecular level basis for ligand discrimination. Using the two ligands O2 and CO, we have determined that the perturbation from the heme-ligand complex has a strong influence on the characteristics of the myoglobin collective dynamics that are excited upon binding. Further, the differences detected in the collective protein motions in Mb-O2 compared with those in Mb-CO appear to be intimately tied with the pathways of long-range allosteric communication in the protein, which ultimately determine the trajectories selected by the respective ligands on the path to and from the heme-binding cavity.