Proteins at interfaces play important roles in cell biology, immunology, bioengineering, and biomimetic material design. Many biological processes are based on interfacial protein action, ranging ...from cellular communication to immune responses and the protein-driven mineralization of bone. Despite the importance of interfacial proteins, comparatively little is known about their structure. The standard methods for studying crystalline or solution-phase proteins (X-ray diffraction and NMR spectroscopy) are not well-suited for studying proteins at interfaces, and for these proteins we still lack a corresponding technique that can provide the same level of structural resolution. This is not surprising in view of the challenges involved in probing the structure of proteins within monomolecular films assembled at a very thin interface in situ. Vibrational sum-frequency generation (SFG) spectroscopy has the potential to overcome this challenge and investigate the structure and dynamics of proteins at interfaces at the molecular level with subpicosecond time resolution. While SFG studies were initially limited to simple model peptides, the past decade has seen a dramatic advancement of experimental techniques and data analysis methods that has made it possible to also study interfacial proteins and their folding, binding, orientation, hydration, and dynamics. In this review, we first explain the principles of SFG spectroscopy and the experimental and theoretical methods to measure and analyze protein SFG spectra. Then we give an extensive overview of the interfacial proteins studied to date with SFG. We highlight representative examples to demonstrate recent advances in probing the structure of proteins at the interfaces of liquids, membranes, minerals, and synthetic materials.
Understanding the assembly of proteins at the air-water interface (AWI) informs the formation of protein films, emulsion properties, and protein aggregation. Determination of protein conformation and ...orientation at an interface is difficult to resolve with a single experimental or simulation technique alone. To date, the interfacial structure of even one of the most widely studied proteins, lysozyme, at the AWI remains unresolved. In this study, molecular dynamics (MD) simulations are used to determine if the protein adopts a side-on, head-on, or axial orientation at the AWI with two different forcefields, GROMOS-53a6 + SPC/E and a99SB-disp + TIP4P-D. Vibrational sum frequency generation (SFG) spectroscopy experiments and spectral SFG calculations validate consistency between the structure determined from MD and experiments. Overall, we show with strong agreement that lysozyme adopts an axial conformation at pH 7. Further, we provide molecular-level insight as to how pH influences the binding domains of lysozyme resulting in side-on adsorption near the isoelectric point of the lysozyme.
The aggregation of the intrinsically disordered protein alpha-synuclein (αS) into amyloid fibrils is thought to play a central role in the pathology of Parkinson's disease. Using a combination of ...techniques (AFM, UV-CD, XRD, and amide-I 1D- and 2D-IR spectroscopy) we show that the structure of αS fibrils varies as a function of ionic strength: fibrils aggregated in low ionic-strength buffers (NaCl ≤ 25 mM) have a significantly different structure than fibrils grown in higher ionic-strength buffers. The observations for fibrils aggregated in low-salt buffers are consistent with an extended conformation of αS molecules, forming hydrogen-bonded intermolecular β-sheets that are loosely packed in a parallel fashion. For fibrils aggregated in high-salt buffers (including those prepared in buffers with a physiological salt concentration) the measurements are consistent with αS molecules in a more tightly-packed, antiparallel intramolecular conformation, and suggest a structure characterized by two twisting stacks of approximately five hydrogen-bonded intermolecular β-sheets each. We find evidence that the high-frequency peak in the amide-I spectrum of αS fibrils involves a normal mode that differs fundamentally from the canonical high-frequency antiparallel β-sheet mode. The high sensitivity of the fibril structure to the ionic strength might form the basis of differences in αS-related pathologies.
The conversion of biomass into green fuels and chemicals is of great societal interest. Engineers have been designing new cellulase enzymes for the breakdown of otherwise insoluble cellulose ...materials. A barrier to the rational design of new enzymes has been our lack of a molecular picture of how cellulase binding occurs. A critical factor is the attachment via the enzyme’s carbohydrate binding module (CBM). To elucidate the structural and mechanistic details of cellulase adsorption, we have combined experimental data from sum frequency generation spectroscopy with molecular dynamics simulations to probe the equilibrium structure and surface alignment of a 14-residue peptide mimicking the CBM. The data show that binding is driven by hydrogen bonding and that tyrosine side chains within the CBM align the cellulase with the registry of the cellulose surface. Such an alignment is favorable for the translocation and effective cellulose breakdown and is therefore likely an important parameter for the design of novel enzymes.
The vibrational coupling between protein backbone modes and the role of water interactions are important topics in biomolecular spectroscopy. Our work reports the first study of the coupling between ...amide I and amide A modes within peptides and proteins with secondary structure and water contacts. We use two-color two-dimensional infrared (2D IR) spectroscopy and observe cross peaks between amide I and amide A modes. In experiments with peptides with different secondary structures and side chains, we observe that the spectra are sensitive to secondary structure. Water interactions affect the cross peaks, which may be useful as probes for the accessibility of protein sites to hydration water. Moving to two-color 2D IR spectra of proteins, the data demonstrate that the cross peaks integrate the sensitivities of both amide I and amide A spectra and that a two-color detection scheme may be a promising tool for probing secondary structures in proteins.
Ice-nucleation active (INA) bacteria can promote the growth of ice more effectively than any other known material. Using specialized ice-nucleating proteins (INPs), they obtain nutrients from plants ...by inducing frost damage and, when airborne in the atmosphere, they drive ice nucleation within clouds, which may affect global precipitation patterns. Despite their evident environmental importance, the molecular mechanisms behind INP-induced freezing have remained largely elusive. We investigate the structural basis for the interactions between water and the ice-nucleating protein InaZ from the INA bacterium Pseudomonas syringae. Using vibrational sum-frequency generation (SFG) and two-dimensional infrared spectroscopy, we demonstrate that the ice-active repeats of InaZ adopt a β-helical structure in solution and at water surfaces. In this configuration, interaction between INPs and water molecules imposes structural ordering on the adjacent water network. The observed order of water increases as the interface is cooled to temperatures close to the melting point of water. Experimental SFG data combined with molecular-dynamics simulations and spectral calculations show that InaZ reorients at lower temperatures. This reorientation can enhance water interactions, and thereby the effectiveness of ice nucleation.
The development of methods that allow a structural interpretation of linear and non-linear vibrational spectra is of great importance, both for spectroscopy and for optimizing force-field quality. ...The experimentally measured signals are ensemble averages over all accessible configurations, which complicates spectral calculations. To account for this, we present a recipe for calculating vibrational amide-I spectra of proteins based on metadynamics molecular dynamics simulations. For each frame, a one-exciton Hamiltonian is set up for the backbone amide groups, in which the couplings are estimated with the transition-charge coupling model for non-nearest neighbors, and with a parametrized map of ab initio calculations that give the coupling as a function of the dihedral angles for nearest neighbors. The local-mode frequency variations due to environmental factors such as hydrogen bonds are modelled by exploiting the linear relationship between the amide C-O bond length and the amide-I frequency. The spectra are subsequently calculated while taking into account the equilibrium statistical weights of the frames that are determined using a previously-published reweighting procedure. By implementing all these steps in an efficient Fortran code, the spectra can be averaged over very large amounts of structures, thereby extensively covering the phase space of proteins. Using this recipe, the spectral responses of 2.5 million frames of a metadynamics simulation of the miniprotein Trp-cage are averaged to reproduce the experimental temperature-dependent IR spectra very well. The spectral calculations provide new insight into the origin of the various spectral signatures (which are typically challenging to disentangle in the congested amide-I region), and allow for a direct structural interpretation of the experimental spectra and for validation of the molecular dynamics simulations of ensembles.The development of methods that allow a structural interpretation of linear and non-linear vibrational spectra is of great importance, both for spectroscopy and for optimizing force-field quality. The experimentally measured signals are ensemble averages over all accessible configurations, which complicates spectral calculations. To account for this, we present a recipe for calculating vibrational amide-I spectra of proteins based on metadynamics molecular dynamics simulations. For each frame, a one-exciton Hamiltonian is set up for the backbone amide groups, in which the couplings are estimated with the transition-charge coupling model for non-nearest neighbors, and with a parametrized map of ab initio calculations that give the coupling as a function of the dihedral angles for nearest neighbors. The local-mode frequency variations due to environmental factors such as hydrogen bonds are modelled by exploiting the linear relationship between the amide C-O bond length and the amide-I frequency. The spectra are subsequently calculated while taking into account the equilibrium statistical weights of the frames that are determined using a previously-published reweighting procedure. By implementing all these steps in an efficient Fortran code, the spectra can be averaged over very large amounts of structures, thereby extensively covering the phase space of proteins. Using this recipe, the spectral responses of 2.5 million frames of a metadynamics simulation of the miniprotein Trp-cage are averaged to reproduce the experimental temperature-dependent IR spectra very well. The spectral calculations provide new insight into the origin of the various spectral signatures (which are typically challenging to disentangle in the congested amide-I region), and allow for a direct structural interpretation of the experimental spectra and for validation of the molecular dynamics simulations of ensembles.
We propose a computational investigation on the interaction mechanisms between SARS-CoV-2 spike protein and possible human cell receptors. In particular, we make use of our newly developed numerical ...method able to determine efficiently and effectively the relationship of complementarity between portions of protein surfaces. This innovative and general procedure, based on the representation of the molecular isoelectronic density surface in terms of 2D Zernike polynomials, allows the rapid and quantitative assessment of the geometrical shape complementarity between interacting proteins, which was unfeasible with previous methods. Our results indicate that SARS-CoV-2 uses a dual strategy: in addition to the known interaction with angiotensin-converting enzyme 2, the viral spike protein can also interact with sialic-acid receptors of the cells in the upper airways.
Silaffin peptide R5 is key for the biogenesis of silica cell walls of diatoms. Biosilification by the R5 peptide has potential in biotechnology, drug development, and materials science due to its ...ability to precipitate stable, high fidelity silica sheets and particles. A true barrier for the design of novel peptide-based architectures for wider applications has been the limited understanding of the interfacial structure of R5 when precipitating silica nanoparticles. While R5–silica interactions have been studied in detail at flat surfaces, the structure within nanophase particles is still being debated. We herein elucidate the conformation of R5 in its active form within silica particles by combining interface-specific vibrational spectroscopy data with solid-state NMR torsion angles using theoretical spectra. Our calculations show that R5 is structured and undergoes a conformational transition from a strand-type motif in solution to a more curved, contracted structure when interacting with silica precursors.
Abstract
The amyloid aggregation of
α
-synuclein (
α
S), related to Parkinson’s disease, can be catalyzed by lipid membranes. Despite the importance of lipid surfaces, the 3D-structure and ...orientation of lipid-bound
α
S is still not known in detail. Here, we report interface-specific vibrational sum-frequency generation (VSFG) experiments that reveal how monomeric
α
S binds to an anionic lipid interface over a large range of
α
S-lipid ratios. To interpret the experimental data, we present a frame-selection method ("
ViscaSelect
”) in which out-of-equilibrium molecular dynamics simulations are used to generate structural hypotheses that are compared to experimental amide-I spectra via excitonic spectral calculations. At low and physiological
α
S concentrations, we derive flat-lying helical structures as previously reported. However, at elevated and potentially disease-related concentrations, a transition to interface-protruding
α
S structures occurs. Such an upright conformation promotes lateral interactions between
α
S monomers and may explain how lipid membranes catalyze the formation of
α
S amyloids at elevated protein concentrations.