Coherent Neutron Scattering and Collective Dynamics in the Protein,GFP

Collective dynamics are considered to be one of the major properties of soft materials, including biological macromolecules. We present coherent neutron scattering studies of the low-frequency vibrations, the so-called boson peak, in fully deuterated green fluorescent protein (GFP). Our analysis revealed unexpectedly low coherence of the atomic motions in GFP. This result implies a low amount of in-phase collective motion of the secondary structural units contributing to the boson peak vibrations and fast conformational fluctuations on the picosecond timescale. These observations are in contrast to earlier studies of polymers and glass-forming systems, and suggest that random or out-of-phase motions of the β-strands contribute greater than two-thirds of the intensity to the low-frequency vibrational spectra of GFP.     

Rigidity,Secondary Structure,and the Universality of the Boson Peak in Proteins

Complementary neutron- and light-scattering results on nine proteins and amino acids reveal the role of rigidity and secondary structure in determining the time- and lengthscales of low-frequency collective vibrational dynamics in proteins. These dynamics manifest in a spectral feature, known as the boson peak (BP), which is common to all disordered materials. We demonstrate that BP position scales systematically with structural motifs, reflecting local rigidity: disordered proteins appear softer than α-helical proteins; which are softer than β-sheet proteins. Our analysis also reveals a universal spectral shape of the BP in proteins and amino acid mixtures; superimposable on the shape observed in typical glasses. Uniformity in the underlying physical mechanism, independent of the specific chemical composition, connects the BP vibrations to nanometer-scale heterogeneities, providing an experimental benchmark for coarse-grained simulations, structure/rigidity relationships, and engineering of proteins for novel applications.     

Rigidity,Secondary Structure,and the Universality of the Boson Peak in Proteins

Complementary neutron- and light-scattering results on nine proteins and amino acids reveal the role of rigidity and secondary structure in determining the time- and lengthscales of low-frequency collective vibrational dynamics in proteins. These dynamics manifest in a spectral feature, known as the boson peak (BP), which is common to all disordered materials. We demonstrate that BP position scales systematically with structural motifs, reflecting local rigidity: disordered proteins appear softer than α-helical proteins; which are softer than β-sheet proteins. Our analysis also reveals a universal spectral shape of the BP in proteins and amino acid mixtures; superimposable on the shape observed in typical glasses. Uniformity in the underlying physical mechanism, independent of the specific chemical composition, connects the BP vibrations to nanometer-scale heterogeneities, providing an experimental benchmark for coarse-grained simulations, structure/rigidity relationships, and engineering of proteins for novel applications.     

Hydration effect on low-frequency protein dynamics observed in simulated neutron scattering spectra

Hydration effects on protein dynamics were investigated by comparing the frequency dependence of the calculated neutron scattering spectra between full and minimal hydration states at temperatures between 100 and 300 K. The protein boson peak is observed in the frequency range 1-4 meV at 100 K in both states. The peak frequency in the minimal hydration state shifts to lower than that in the full hydration state. Protein motions with a frequency higher than 4 meV were shown to undergo almost harmonic motion in both states at all temperatures simulated, whereas those with a frequency lower than 1 meV dominate the total fluctuations above 220 K and contribute to the origin of the glass-like transition. At 300 K, the boson peak becomes buried in the quasielastic contributions in the full hydration state but is still observed in the minimal hydration state. The boson peak is observed when protein dynamics are trapped within a local minimum of its energy surface. Protein motions, which contribute to the boson peak, are distributed throughout the whole protein. The fine structure of the dynamics structure factor is expected to be detected by the experiment if a high resolution instrument (<∼20 μeV) is developed in the near future.     

Biological functions of low-frequency vibrations (phonons). III. Helical structures and microenvironment.

Low-frequency vibrations in biomacromolecules possess significant biological functions. In this paper, the alpha-helix element is compared with a mass-distributed spring. Based on this, a set of intuitive and easily handled equations are derived for predicting the fundamental frequencies of helical structures in protein molecules. As shown in the equations, the fundamental frequency depends not only on the constituents of a helix itself but also on its microenvironment. The calculated results agree with the observations. The calculations also demonstrate that the low-frequency vibrations with wave number of approximately 30 cm-1 do not necessarily arise from motions that involve either all or very large portions of the protein molecule as previously thought; a piece of helix containing more than 10 residues and surrounded by a proper microenvironment can also generate such low-frequency motions. Furthermore , we illustrate that the low-frequency motions are closely related to the native state of a protein molecule. Upon denaturation, which is accompanied by a radical change of the relevant microenvironment, the original fundamental frequency also disappears. Consequently, this kind of special frequency termed activating low frequency can serve as a dynamic criterion in identifying whether a biomacromolecule is in its native state. The energy of a phonon excited by this kind of low-frequency vibration is of the same order of magnitude as the average enthalpy value per residue measured during conformational change in some protein molecules. Therefore, there must be some intrinsic relation between the allosteric transitions of protein molecules and their low-frequency motions.     

Low-temperature molecular dynamics simulations of horse heart cytochrome c and comparison with inelastic neutron scattering data

Molecular dynamics (MD) simulation combined with inelastic neutron scattering can provide information about the thermal dynamics of proteins, especially the low-frequency vibrational modes responsible for large movement of some parts of protein molecules. We performed several 30-ns MD simulations of cytochrome c (Cyt c) in a water box for temperatures ranging from 110 to 300 K and compared the results with those from experimental inelastic neutron scattering. The low-frequency vibrational modes were obtained via dynamic structure factors, S(Q, ω), obtained both from inelastic neutron scattering experiments and calculated from MD simulations for Cyt c in the same range of temperatures. The well known thermal transition in structural movements of Cyt c is clearly seen in MD simulations; it is, however, confined to unstructured fragments of loops Ω₁ and Ω₂; movement of structured loop Ω₃ and both helical ends of the protein is resistant to thermal disturbance. Calculated and experimental S(Q, ω) plots are in qualitative agreement for low temperatures whereas above 200 K a boson peak vanishes from the calculated plots. This may be a result of loss of crystal structure by the protein–water system compared with the protein crystal.     

Biophysics of underwater hearing in the clawed frog,Xenopus laevis

Anesthetized clawed frogs (Xenopus laevis) were stimulated with underwater sound and the tympanic disk vibrations were studied using laser vibrometry. The tympanic disk velocities ranged from 0.01 to 0.5 mm/s (at a sound pressure of 2 Pa) in the frequency range of 0.4–4 kHz and were 20–40 dB higher than those of the surrounding tissue. The frequency response of the disk had two peaks, in the range of 0.6–1.1 kHz and 1.6–2.2 kHz, respectively. The first peak corresponded to the peak vibrations of the body wall overlying the lung. The second peak matched model predictions of the pulsations of the air bubble in the middle ear cavity. Filling the middle ear cavity with water lowered the disk vibrations by 10–30 dB in the frequency range of 0.5–3 kHz.Inflating the lungs shifted the low-frequency peak downwards, but did not change the high-frequency peak. Thus, the disk vibrations in the frequency range of the mating call (main energy at 1.7–1.9 kHz) were mainly caused by pulsations of the air in the middle ear cavity; sound transmission via the lungs was more important at low frequencies (below 1 kHz). Furthermore, the low-frequency peak could be reversibly reduced in amplitude by loading the larynx with metal or tissue glue. This shows that the sound-induced vibrations of the lungs are probably coupled to the middle ear cavities via the larynx. Also, anatomical observations show that the two middle ear cavities and the larynx are connected in an air-filled recess in submerged animals.This arrangement is unique to pipid frogs and may be a structural adaptation to connect all the air spaces of the frog and improve low-frequency underwater hearing. Another function of the recess may be to allow cross-talk between the two middle ear cavities. Thus, the ear might be directional. Our pilot experiments show up to 10 dB difference between ipsi- and contralateral stimulus directions in a narrow frequency range around 2 kHz.     

A Physical Picture of Protein Dynamics and Conformational Changes

A physical model is reviewed which explains different aspects of protein dynamics consistently. At low temperatures, the molecules are frozen in conformational substates. Their average energy is 3/2RT. Solid-state vibrations occur on a time scale of femtoseconds to nanoseconds. Above a characteristic temperature, often called the dynamical transition temperature, slow modes of motions can be observed occurring on a time scale between about 140 and 1 ns. These motions are overdamped, quasidiffusive, and involve collective motions of segments of the size of an α-helix. Molecules performing these types of motion are in the “flexible state”. This state is reached by thermal activation. It is shown that these motions are essential for conformational relaxation. Based on this picture, a new approach is proposed to understand conformational changes. It connects structural fluctuations and conformational transitions.     

Green Fluorescent Protein-based Expression Screening of Membrane Proteins in Escherichia coli

The production of recombinant membrane proteins for structural and functional studies remains technically challenging due to low levels of expression and the inherent instability of many membrane proteins once solubilized in detergents. A protocol is described that combines ligation independent cloning of membrane proteins as GFP fusions with expression in Escherichia coli detected by GFP fluorescence. This enables the construction and expression screening of multiple membrane protein/variants to identify candidates suitable for further investment of time and effort. The GFP reporter is used in a primary screen of expression by visualizing GFP fluorescence following SDS polyacrylamide gel electrophoresis (SDS-PAGE). Membrane proteins that show both a high expression level with minimum degradation as indicated by the absence of free GFP, are selected for a secondary screen. These constructs are scaled and a total membrane fraction prepared and solubilized in four different detergents. Following ultracentrifugation to remove detergent-insoluble material, lysates are analyzed by fluorescence detection size exclusion chromatography (FSEC). Monitoring the size exclusion profile by GFP fluorescence provides information about the mono-dispersity and integrity of the membrane proteins in different detergents. Protein: detergent combinations that elute with a symmetrical peak with little or no free GFP and minimum aggregation are candidates for subsequent purification. Using the above methodology, the heterologous expression in E. coli of SED (shape, elongation, division, and sporulation) proteins from 47 different species of bacteria was analyzed. These proteins typically have ten transmembrane domains and are essential for cell division. The results show that the production of the SEDs orthologues in E. coli was highly variable with respect to the expression levels and integrity of the GFP fusion proteins. The experiment identified a subset for further investigation.     

Lattice vibrations in crystalline L-alanine

Zwitterionic L -alanine forms crystals containing strong hydrogen-bonding and methyl-methyl interactions. Well-defined low-frequency lattice vibrations exist in the crystals involving correlated intermolecular motions on the picosecond timescale. A characterization of these vibrations is expected to provide useful information on the nature of nonbonded interactions in peptides and proteins. We examine some of the vibrations using coherent inelastic neutron scattering and computer simulation techniques. The neutron scattering measurements are used to determine phonon dispersion relations for the acoustic and some low-frequency optic modes in the crystal. There is evidence for interaction between the two lowest frequency optical phonons and the longitudinal acoustic mode. The velocity of sound is anisotropic and can be correlated with the hydrogen-bonding arrangement in the crystal. Corresponding phonon dispersion relations are derived from normal mode analyses of the crystal using the program CHARMM. Although some calculated vibrational frequencies are somewhat too high, the form of the calculated dispersion relations are in good agreement with experiment. © 1993 John Wiley & Sons, Inc.     

Local dynamics of DNA probed with optical absorption spectroscopy of bound ethidium bromide.

We have studied the local dynamics of calf thymus double-helical DNA by means of an "optical labeling" technique. The study has been performed by measuring the visible absorption band of the cationic dye ethidium bromide, both free in solution and bound to DNA, in the temperature interval 360-30 K and in two different solvent conditions. The temperature dependence of the absorption line shape has been analyzed within the framework of the vibronic coupling theory, to extract information on the dynamic properties of the system; comparison of the thermal behavior of the absorption band of free and DNA-bound ethidium bromide gave information on the local dynamics of the double helix in the proximity of the chromophore. For the dye free in solution, large spectral heterogeneity and coupling to a "bath" of low-frequency (soft) modes is observed; moreover, anharmonic motions become evident at suitably high temperatures. The average frequency of the soft modes and the amplitude of anharmonic motions depend upon solvent composition. For the DNA-bound dye, at low temperatures, heterogeneity is decreased, the average frequency of the soft modes is increased, and anharmonic motions are hindered. However, a new dynamic regime characterized by a large increase in anharmonic motions is observed at temperatures higher than approximately 280 K. The DNA double helix therefore appears to provide, at low temperatures, a rather rigid environment for the bound chromophore, in which conformational heterogeneity is reduced and low-frequency motions (both harmonic vibrations and anharmonic contributions) are hindered. The system becomes anharmonic at approximately 180 K; however, above approximately 280 K, anharmonicity starts to increase much more rapidly than for the dye free in solution; this can be attributed to the onset of wobbling of the dye in its intercalation site, which is likely connected with the onset of (functionally relevant) DNA motions, involving local opening/unwinding of the double helix. As shown by parallel measurements of the melting curves, these motions precede the melting of the double helix and depend upon solvent composition much more than does the melting itself.     

Filtering molecular dynamics trajectories to reveal low-frequency collective motions: phospholipase A2

A novel method for analysing molecular dynamics trajectories has been developed, which filters out high frequencies using digital signal processing techniques and facilitates focusing on the low-frequency collective motions of proteins. These motions involve low energy slow motions, which lead to important biological phenomena such as domain closure and allosteric effects in enzymes. The filtering method treats each of the atomic trajectories obtained from the molecular dynamics simulation as a "signal". The trajectories of each of the atoms in the system (or any subset of interest) are Fourier transformed to the frequency domain, a filtering function is applied and then an inverse transformation back to the time domain yields the filtered trajectory. The filtering method has been used to study the dynamics of the enzyme phospholipase A2. In the filtered trajectory, all the high frequency bond and valence angle vibrations were eliminated, leaving only low-frequency motion, mainly fluctuations in torsions and conformational transitions. Analysis of this trajectory revealed interesting motions of the protein, including concerted movements of helices, and changes in shape of the active site cavity. Unlike normal mode analysis, which has been used to study the motion of proteins, this method does not require converged minimizations or diagonalization of a matrix of second derivatives. In addition, anharmonicity, multiple minima and conformational transitions are treated explicitly. Thus, the filtering method avoids most of the approximations implicit in other investigations of the dynamic behaviour of large systems.     

Sound Waves Effectively Assist Tobramycin in Elimination of Pseudomonas aeruginosa Biofilms In vitro

Microbial biofilms are highly refractory to antimicrobials. The aim of this study was to investigate the use of low-frequency vibration therapy (20–20 kHz) on antibiotic-mediated Pseudomonas aeruginosa biofilm eradication. In screening studies, low-frequency vibrations were applied on model biofilm compositions to identify conditions in which surface standing waves were observed. Alginate surface tension and viscosity were also measured. The effect of vibration on P. aeruginosa biofilms was studied using a standard biofilm assay. Subminimal inhibitory concentrations (sub-MIC) of tobramycin (5 μg/ml) were added to biofilms 3 h prior, during, and immediately after vibration and quantitatively assessed by (2,3-bis-(2-methoxy-4-nitro-5-sulfophenyl)-2H-tetrazolium-5-carboxanilide) reduction assay (XTT) and, qualitatively, by confocal laser scanning microscopy (CLSM). The standing waves occurred at frequencies <1,000 Hz. Biofilms vibrated without sub-MIC tobramycin showed a significantly reduced metabolism compared to untreated controls (p < 0.05). Biofilms treated with tobramycin and vibrated simultaneously (450, 530, 610, and 650 Hz), or vibrated (450 and 650 Hz) then treated with tobramycin subsequently, or vibrated (610Hz, 650Hz) after 3 h of tobramycin treatment showed significantly lower metabolism compared to P. aeruginosa biofilm treated with tobramycin alone (p < 0.05). CLSM imaging further confirmed these findings. Low frequency vibrations assisted tobramycin in killing P. aeruginosa biofilms at sub-MIC. Thus, sound waves together with antibiotics are a promising approach in eliminating pathogenic biofilms.KEY WORDS: alginate, biofilm, Pseudomonas, tobramycin, vibration     

Global collective motions in the mammalian and bacterial respiratory complex I

The respiratory complex I is an enzyme responsible for the conversion of chemical energy into an electrochemical proton motive force across the membrane. Despite extensive studies, the mechanism by which the activity of this enormous, ca. 1?MDa, redox-coupled proton pump is regulated still remains unclear. Recent structural studies (Zhu et al., Nature 2016; Fiedorczuk et al., Nature 2016) resolved complex I in different conformations connected to the active-to-deactive (A/D) transition that regulate complex I activity in several species. Based on anisotropic network models (ANM) and principal component analysis (PCA), we identify here transitions between experimentally resolved structures of the mammalian complex I as low-frequency collective motions of the enzyme, highlighting similarities and differences between the bacterial and mammalian enzymes. Despite the reduced complexity of the smaller bacterial enzyme, our results suggest that the global dynamics of complex I is overall conserved. We further probe how the supernumerary subunits could be involved in the modulation of the A/D-transition, and show that in particular the 42?kDa and B13 subunits affect the global motions of the mammalian enzyme.     

Internal DNA Modes Below 25 cm−1: A Resonance Raman Spectroscopy Observation

Abstract

The first resonance Raman scattering observation of the low-frequency (LF) region (below 40 up to 12 cm^?1) of DNA motions is presented. Since the concentration of the studied DNA solution was very low (1 mg/ml), the spectra features reflect internal vibrations of the macro-molecule. The decomposition of the spectra into Lorentzians clearly indicate three intrahelical DNA modes: the corresponding peaks are located at the frequencies 16,19, and 23 (±1) cm^?1. This result is in agreement with our quasi-continuity model of the LF B-form DNA dynamics (V. Lisy, P. Miskovsky and P. Schreiber, J. Biomol. Struct. Dyn. 13, 707 (1996)). The fit of the experimental frequencies to the theory, using the Genetic Algorithms approach, allowed us to make some conclusions about the model force constants which could be found by independent conformational energy calculations. Possible positions of five lowest-frequency DNA peaks, predicted by the model, are discussed.     

Influence of SigB inactivation on Corynebacterium glutamicum protein secretion

The non-essential Corynebacterium glutamicum sigma factor, sigB, modulates global gene expression during the transition from exponential growth to the stationary phase. Utilizing a signal peptide derived from C. glutamicum R CgR_0949, a sigB disruption mutant able to secrete 3- to 5-fold more green fluorescence protein (GFP) and α-amylase than the wild type strain was isolated. The signal peptide selectively enabled the mutant to produce greater amounts of both proteins, which were in turn secreted in culture medium in greater quantities than previously acknowledged. A peak GFP productivity of 2.8 g/l was attained, representing the highest GFP productivity reported in C. glutamicum to date. CgR_0949 signal sequence length (30 residues), type (Tat) or the target protein identity (GFP or α-amylase) had no measurable effect on the magnitude of the protein accumulation and consequent secretion. It therefore follows that actual experimentation remains the fastest way to identify suitable signal sequences in C. glutamicum. More secretion studies may reveal even greater secretion productivity by C. glutamicum and consequently present an attractive avenue to further enhance the utility of C. glutamicum as an industrial workhorse.     

Low-frequency Raman scattering study of six nucleosides

Raman spectroscopy was used to study the low-frequency (?200?cm^?1) vibrations in crystalline samples of six naturally occurring nucleosides: deoxythymidine (dT), deoxycytidine (dC), deoxyadenosine (dA), uridine (rU), cytidine (rC), and adenosine (rA). Such low-frequency vibrations are important for biological processes in which the conformation of a nucleic acid molecule changes. These experiments also provide a test for the low-frequency vibrational modes of dT, dC, and dA predicted by Shishkin et al.     

Dual color microscopic imagery of cells expressing the green fluorescent protein and a red-shifted variant

The green fluorescent protein (GFP) from the jellyfish, Aequorea victoria, has become a versatile reporter for monitoring gene expression and protein localization in a variety of cells and organisms. GFP emits bright green light (λ_max = 510 nm) when excited with ultraviolet (UV) or blue light (λ_max = 395 nm, minor peak at 470 nm). The chromophore in GFP is intrinsic to the primary structure of the protein, and fluorescence from GFP does not require additional gene products, substrates or other factors. GFP fluorescence is stable, species-independent and can be monitored noninvasively using the techniques of fluorescence microscopy and flow cytometry [Chalfie et al., Science 263 (1994) 802–805; Stearns, Curr. Biol. 5 (1995) 262–264]. The protein appears to undergo an autocatalytic reaction to create the fluorophore [Heim et al., Proc. Natl. Acad. Sci. USA 91 (1994) 12501–12504] in a process involving cyclization of a Tyr⁶⁶ aa residue. Recently [Delagrave et al., Bio/Technology 13 (1995) 151–154], a combinatorial mutagenic strategy was targeted at aa 64 through 69, which spans the chromophore of A. victoria GFP, yielding a number of different mutants with redshifted fluorescence excitation spectra. One of these, RSGFP4, retains the characteristic green emission spectra (λ_max = 505 nm), but has a single excitation peak (λ_max = 490 nm). The fluorescence properties of RSGFP4 are similar to those of another naturally occurring GFP from the sea pansy, Renilla reniformis [Ward and Cormier, Photobiochem. Photobiol. 27 (1978) 389–396]. In the present study, we demonstrate by fluorescence microscopy that selective excitation of A. victoria GFP and RSGFP4 allows for spectral separation of each fluorescent signal, and provides the means to image these signals independently in a mixed population of bacteria or mammalian cells.     

Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the “boson peak” observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure. Received: 12 October 1998 / Revised version: 19 February 1999 / Accepted: 1 March 1999     

Water-coupled low-frequency modes of myoglobin and lysozyme observed by inelastic neutron scattering.

Conformational changes of proteins often involve the relative motion of rigid structural domains. Normal mode analysis and molecular dynamics simulations of small globular proteins predict delocalized vibrations with frequencies below 20 cm(-1), which may be overdamped in solution due to solvent friction. In search of these modes, we have studied deuterium-exchanged myoglobin and lysozyme using inelastic neutron scattering in the low-frequency range at full and low hydration to modify the degree of damping. At room temperature, the hydrated samples exhibit a more pronounced quasielastic spectrum due to diffusive motions than the dehydrated samples. The analysis of the corresponding lineshapes suggests that water modifies mainly the amplitude, but not the characteristic time of fast protein motions. At low temperatures, in contrast, the dehydrated samples exhibit larger motional amplitudes than the hydrated ones. The excess scattering, culminating at 16 cm(-1), is suggested to reflect water-coupled librations of polar side chains that are depressed in the hydrated system by strong intermolecular hydrogen bonding. Both myoglobin and lysozyme exhibit ultra-low-frequency modes below 10 cm(-1) in the dry state, possibly related to the breathing modes predicted by harmonic analysis.