Skip to main content

Observation of the low frequency vibrational modes of bacteriophage M13 in water by Raman spectroscopy



Recently, a technique which departs radically from conventional approaches has been proposed. This novel technique utilizes biological objects such as viruses as nano-templates for the fabrication of nanostructure elements. For example, rod-shaped viruses such as the M13 phage and tobacco mosaic virus have been successfully used as biological templates for the synthesis of semiconductor and metallic nanowires.

Results and discussion

Low wave number (≤ 20 cm-1) acoustic vibrations of the M13 phage have been studied using Raman spectroscopy. The experimental results are compared with theoretical calculations based on an elastic continuum model and appropriate Raman selection rules derived from a bond polarizability model. The observed Raman mode has been shown to belong to one of the Raman-active axial torsion modes of the M13 phage protein coat.


It is expected that the detection and characterization of this low frequency vibrational mode can be used for applications in nanotechnology such as for monitoring the process of virus functionalization and self-assembly. For example, the differences in Raman spectra can be used to monitor the coating of virus with some other materials and nano-assembly process, such as attaching a carbon nanotube or quantum dots.


Although viruses are generally regarded as etiologic agents of disease, they can also be made useful, and the concept of using them as a template to build uniform semiconductor nanostructures seems to be increasingly plausible. Among the crucial issues associated with molecular-beam epitaxy (MBE) self-assembly for applications in nanotechnology is the variation in component size [1, 2]. The production of identical or nearly identical structures at nanoscale is desirable; however, this is extremely difficult to accomplish with conventional techniques. To date, MBE-grown nanowires, nanorods, and quantum dots, which are proposed to be elements in future nanoelectronic circuits, have been characterized by size dispersion despite numerous efforts towards the control of size and shape [3, 4]. Although the method of self-assembly by employment of chemically generated templates [5, 6] decreases to some extent the size dispersion of the grown structures, the results are still not satisfactory.

Recently, a technique which departs radically from conventional approaches has been proposed [711]. This novel technique utilizes biological objects such as viruses as nano-templates for the fabrication of nanostructure elements. For example, rod-shaped viruses such as the M13 phage and tobacco mosaic virus (TMV) have been successfully used as biological templates for the synthesis of semiconductor and metallic nanowires [7, 10, 11]. Furthermore, genetically modified TMV and M13 phage have been shown to be useful for the self-assembly of nanomaterials into liquid crystals, films and fibers [8, 9]. It is therefore very likely that genetically programmed viruses will play an important role in developing the next generation of optoelectronic devices and nanoelectronic circuits.

For monitoring the abovementioned self-assembly processes, an in-situ, non-destructive technique is desirable. Raman spectroscopy has been shown to be a non-invasive technique in material research. To the authors' knowledge, previous studies of viruses using Raman spectroscopy have focused only on the high frequency (large wave number) regions (≥ 600 cm-1) where the internal virus composition, i.e. localized vibrations of multiply bonded or electron-rich groups in proteins, was studied [12]. In this paper, we report the first observation of low wave number (≤ 20 cm-1) acoustic vibrations of the M13 phage using Raman spectroscopy. The observed vibrations are compared with theoretical calculations based on an elastic continuum model and appropriate Raman intensities and selection rules derived from a continuum limit of the bond polarizability model. The observed Raman mode has been shown to belong to one of the Raman-active axial torsion modes of the M13 phage protein coat. Because of the sensitivity of these frequencies upon environments, it is expected that the detection of this low frequency vibrational mode can be used to monitor and help to control the process of virus functionalization, By virus functionalization we mean virus being used to make useful things such as nanostructure and nanoelectronic devices; for instance, when coating viruses with different materials, attaching viruses to quantum dots and carbon nanotubes, and forming multiple superstructures.

Samples and experimental technique

The M13 phage samples in water solution used in this work were prepared as follows: To propagate the M13 phage, an overnight culture of the host cell, JM101, was diluted 20-fold into 125-ml flasks containing 20 ml of LB medium. When the culture reached 0.5 of optical density at 550 nm (OD550), the phage was added at a multiplicity of infection of 20 and further grown until stationary phase (ca. 12 h postinfection). Crude phage suspensions were prepared by centrifugation (10,000 × g, 15 min) of the culture to remove the cells and passing the supernatants through a membrane filter (0.45-μm pore size). To concentrate the phage titer, the filtrated supernatant was precipitated by 0.25 M NaCl and 2.5% polyethylene glycol 6000 for 4 hr on ice. The pellet collected by centrifugation (13500 × g for 15 min at 4°C) was dissolved in 100 μl of distilled water. To determine the phage titer, a double-layer bioassay [13] was performed on an LB agar plate.

To prevent heating of the samples by laser irradiation during the Raman scattering experiments, the second harmonic output of a cw mode-locked YAlG laser was used as an excitation source [14]. The laser, which has photon energy of 2.34 eV, was operated at a repetition rate of 76 MHz and has a pulse width of about 70 ps. 90°-scattering geometry was employed. The Raman scattered signal was collected and analyzed by a standard computer-controlled Raman system which included a double spectrometer and a photomultiplier tube with associated photon counting electronics. The spectrometer had a spectral resolution of about 1.0 cm-1. All the data reported here were taken at T = 300K.

Experimental results, theoretical model calculations and discussions

A typical Raman scattering spectrum taken for M13 phages at 1020 pfu/ml and in the spectral range between 2 and 20 cm-1 is shown in Fig. 1 (the solid circles). The experiments have been repeated at least 30 times and the data have been found to be reproducible. They were an average of at least 30 experiments. The statistic error (mostly from that of photon counting) associated with the data is minimal. The distinctive feature of the spectrum is a broad shoulder around 8.5 cm-1 sitting on top of a background. To rule out both instrumental artifacts and the possibility of contributions from the DNA within the phages, we repeated the experiments with M13 phages without protein coats, i.e. with only the single-stranded M13 phage DNA at the same concentration in water. We notice that from the dimension of M13 phage and the weight of a single stranded DNA inside the phage, one can determine the equivalent weight of single stranded DNA for M13 phages with 1020 and 1021 pfu/ml. It is this equivalent weight of single stranded DNA that was used in the experiments. The results are shown as a solid line in Fig. 1. Comparison of the two spectra shows that there is indeed a broad structure at 8.5 cm-1 associated with scattering of light from the M13 phage protein coats. The remaining background is due to imperfections in the rejection of elastic light by the spectrometer. The actual low frequency acoustic vibrational mode signal from the phages is obtained by the subtraction of these two spectra. Figure 2 shows the resulting structure after the subtraction. The broad peak has been found to center around 8.5 cm-1 and is asymmetric. It has a full-width-at-half-maximum (FWHM) of about 5.0 cm-1. Since the spectral resolution of our Raman system is about 1.0 cm-1, we conclude that the relatively broad Raman peak observed here likely resulted from inhomogeneous broadening.

Figure 1
figure 1

Raman scattering spectrum of M13 phages in water for a concentration of 1020 pfu/ml with (solid circles) and without (solid curve) protein coating. The solid curve represents background signal resulting from the imperfection of rejection of elastic scattering of light by spectrometer.

Figure 2
figure 2

The low frequency vibrational mode of M13 with a concentration of 1020 pfu/ml obtained after the subtraction of the background signal.

In order to further test our interpretations, we carried out similar Raman experiments with M13 phages at the higher concentration of 1021 pfu/ml in water solution. The result is shown in Fig. 3. Again, the solid squares and solid line refer to the results for M13 phages with and without protein coats, respectively. Fig. 4 shows the resulting structure after the subtraction. Similar to Fig. 2, the observed broad Raman peak is centered around 8.5 cm-1 and is asymmetric; however, the signal-to-noise ratio is significantly improved. This further confirms that what we have observed is the low frequency acoustic vibrational mode associated with the M13 phage protein coat.

Figure 3
figure 3

Raman scattering spectrum of M13 phages in water for a concentration of 1021 pfu/ml with (solid circles) and without (solid curve) protein coating. The solid curve corresponds to background signal due to the imperfection of rejection of elastic scattering of light by spectrometer.

Figure 4
figure 4

The low frequency vibrational mode of M13 with a concentration of 1021 pfu/ml obtained after the subtraction of the background signal.

Fig. 5 shows a comparison of the observed broad structures for the two concentrations studied, 1020 pfu/ml and 1021 pfu/ml. The data are properly normalized. We find that the integrated area under the peaks scales very well with M13 phage concentration, as expected.

Figure 5
figure 5

Raman spectra of low frequency vibrational mode of M13 phages for two concentrations, as indicated. The Raman intensity has been found to scale with theconcentration of the phages, as expected.

To obtain better insight into the character of the observed low frequency vibrational Raman mode of M13 phages, we performed theoretical calculations in which an elastic continuum model is assumed and appropriate Raman selection rules were derived based on a continuum extension of the bond polarizability model.

Experimental Raman spectroscopy is able to observe a low frequency mode of a cylindrical type viral particle; however, it is unable to provide a qualitative picture of the type of atomic displacements involved in the vibration. To understand what types of modes are being excited by the Raman spectroscopy experiment, we turn to a theoretical analysis of the problem.

There are several routes that can be followed in developing a theory capable of predicting Raman spectral lines, including atomistic models which give full atomic detail, to a continuum model which averages over the local geometry on the atomic scale and expands its vision to a larger length scale. Low frequency modes being probed here through Raman scattering have long wavelengths and thus are not sensitive to atomic details. Rather they are the result of the average global properties of the phage. Our goal here is to provide qualitative insights and a semi-quantitative understanding concerning the types of vibrational modes accessible through Raman scattering of a cylindrical phage. Thus a continuum approach, when used cautiously, is a method of choice due to its relative simplicity and ease of interpretation.

In building a continuum model, it is important to first understand the structure of the M13 bacteriophage. M13 is composed of individual coat proteins bundled together so that in cross section the bundle has approximately the shape of a long circular cylinder (>860 nm), or tube [15]. Each individual coat protein in the bundle is comprised of 50 amino acids and has the conformation of an α-helix roughly aligned with the tube axis. Figure 6 shows a small segment of the filament. In cross section, there are 10 individual proteins packed together. The coat proteins are packed in a partially helical fashion so that they form the shape of a cone. The narrow neck at the top of one cone fits inside the large base of the cone above it. One end of each cone contains acidic residues while the other end is basic. This produces strong cone-cone interactions which increase the stability of the assembly. The central region of the cone is apolar and enhances the interactions between α-helices of a single cone. The stacking of thousands of such cones gives the phage its long worm-like structure. Inside the core of the tube resides the ssDNA and water.

Figure 6
figure 6

Diagrams of the structure of M13 phage. (a) A molecular model containing 10 α-helix coat proteins forming a cone of the structure. (b) A cross section of the phage composed of α-helix coat proteins. (c) The elastic continuum model of M13; a circular cylinder of outer radius "b" and inner radius "a."

We model the long fiber as a solid cylindrical tube of elastic continuous material, as shown in Figure 6c. Following Balandin et al. [16] the viral protein coat has an inner radius "a" and outer radius "b". The inner and outer radii of the M13 shell are set to a = 1.9 nm and b = 3.4 nm, respectively, which were determined from the atomistic structure of the capsid by averaging over distances between atoms in the structure. Elastic wave theory [17, 18] is used to arrive at an analytic expression that describes the frequencies and displacement patterns of each vibrational mode of the cylindrical shell. Cylindrical coordinates (r, θ, z) are used where z is the direction along the axis of the cylinder. In general, the vibrational modes depend on the wave vector kZ varying as eikzZ. We assume that the viral particle is longer than the wavelength of the light; Raman scattering is then only produced by waves with k z ≈ 0. The displacement vectors of the vibrational modes are then u MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWG1bqDgaWcaaaa@2E31@ (r, θ) with no z dependence (eikzZ= 1), and with components ur (radial direction), uθ (angular direction) and uZ (axial direction). The radial dependence is a linear combination of Bessel functions of the 1st and 2nd kinds, and the angular dependence is sinusoidal being proportional to sin(nθ) and cos(nθ) where n is an integer, n = [0,∞]. Since the net radial force is zero on the surfaces of the shell, the radial stress components σ rr , σ , and σ rz must vanish at the inner and outer radii "a" and "b." This gives six conditions relating the amplitudes that produce a 6 × 6 determinant17 which must vanish at the allowed frequencies ω. There are just two parameters in the theory, and they are the longitudinal and transverse sound speeds, c l and c t . We use c l = 1817 m/s and c t = 915 m/s, the measured sound speeds in lysozymes [19], which is similar to the measured speeds of sound in different protein crystals such as ribonuclease and hemoglobin.

Elastic wave theory provides a complete description of the possible frequencies and their displacement patterns of vibrational modes which are then used to determine their Raman scattering intensity. The Raman intensity is proportional to the Raman tensor Δ χ as [20]

I | d ^ s · ( Δ χ ) · d ^ i | 2 Eq . 1 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqscqGHDisTdaabdiqaaiqbdsgaKzaajaWaaSbaaSqaaiabdohaZbqabaGccqWIpM+zdaqadiqaaiabfs5aeHGaciqb=D8aJzaanaaacaGLOaGaayzkaaGaeS4JPFMafmizaqMbaKaadaWgaaWcbaGaemyAaKgabeaaaOGaay5bSlaawIa7amaaCaaaleqabaGaeGOmaidaaOGaaCzcaiaaxMaacqqGfbqrcqqGXbqCcqqGUaGlcqqGGaaicqqGXaqmcqGGSaalaaa@4A60@

where d ^ i MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWGKbazgaqcamaaBaaaleaacqWGPbqAaeqaaaaa@2F94@ , d ^ s MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWGKbazgaqcamaaBaaaleaacqWGZbWCaeqaaaaa@2FA8@ describe the polarization of the incident and scattered light. The Raman tensor Δ χ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWFhpWygaqdaaaa@2E7E@ depends on how the polarizability (due to the electromagnetic field of the laser) of the material changes with displacements produced by the vibrational mode.

For the computation of the Raman tensor we generalize the bond polarizability model [21] to a continuum-like model. We assume that the material is composed of bonds between atoms in which each bond has a derivative of the polarizability tensor with components parallel to the bond direction ( α ' MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFXoqydaqhaaWcbaGaeSyjIafabaGaei4jaCcaaaaa@307D@ ) and perpendicular ( α ' MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFXoqydaqhaaWcbaGaeyyPI4fabaGaei4jaCcaaaaa@3106@ ) due to changes in bond length. Each bond is assumed to occupy a volume V b and all bonds are assumed to have the same polarizability derivates. The bonds are then averaged over all random orientations. The spirit of the model is to produce a long length-scale theory that averages over many atoms but is derived from a tested atomistic model. The derived change of the polarizability tensor due to the strain field is

Δ α = α s V b ( U 1 3 I T r ( U ) ) + α c V b T r ( U ) I . Eq . 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqqHuoariiGacuWFXoqygaqdaiabg2da9maalaaabaGae8xSde2aaSbaaSqaaiabdohaZbqabaaakeaacqWGwbGvdaWgaaWcbaGaemOyaigabeaaaaGcdaqadiqaaiqbdwfavzaanaGaeyOeI0YaaSaaaeaacqaIXaqmaeaacqaIZaWmaaGafmysaKKba0aacqGHflY1cqWGubavcqWGYbGCdaqadiqaaiqbdwfavzaanaaacaGLOaGaayzkaaaacaGLOaGaayzkaaGaey4kaSYaaSaaaeaacqWFXoqydaWgaaWcbaGaem4yamgabeaaaOqaaiabdAfawnaaBaaaleaacqWGIbGyaeqaaaaakiabdsfaujabdkhaYnaabmGabaGafmyvauLba0aaaiaawIcacaGLPaaacuWGjbqsgaqdaiabc6caUiaaxMaacaWLjaGaeeyrauKaeeyCaeNaeeOla4IaeeiiaaIaeeOmaiJaeiOla4caaa@5AF8@

The parameters α s and α c are the polarizability changes of the material due to shear strain and compressional strain, respectively. The parameters α s and α c are related to α ' MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFXoqydaqhaaWcbaGaeSyjIafabaGaei4jaCcaaaaa@307D@ and α ' MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFXoqydaqhaaWcbaGaeyyPI4fabaGaei4jaCcaaaaa@3106@ , but we treat them as parameters. Our calculation only predicts relative (not absolute) intensities, so only the ratio α c /α s is important. We set this ratio to 1/2, which models carbon-carbon bonds well. Only the n = 0 modes depend on this ratio (see Eqs. 3a-c). Details of the model are described in Ref. 22. Integrating Eq. 2 over the volume of the cylindrical viral shell gives the total susceptibly tensor of the viral coat. Since the polarizabilty per unit volume, Eq. 2, is dependant on either cos()or sin(), the integral over the angular dependence implies that the total susceptibility tensor in Cartesian coordinates is non zero for only n = 0,1 and 2 modes of the cylindrical shell.

The results are,

I ( 2 λ + μ ) [ 1 3 + 2 α c α s ] 2 Γ 0 2 ,  for n = 0 ; ( Eq . 3A ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqscqGHDisTdaqadiqaaiabikdaYGGaciab=T7aSjabgUcaRiab=X7aTbGaayjkaiaawMcaamaadmGabaWaaSaaaeaacqaIXaqmaeaacqaIZaWmaaGaey4kaSIaeGOmaiZaaSaaaeaacqWFXoqydaWgaaWcbaGaem4yamgabeaaaOqaaiab=f7aHnaaBaaaleaacqWGZbWCaeqaaaaaaOGaay5waiaaw2faamaaCaaaleqabaGaeGOmaidaaOGaeu4KdC0aa0baaSqaaiabgcdaWaqaaiabgkdaYaaakiabcYcaSiabbccaGiabbAgaMjabb+gaVjabbkhaYjabbccaGiabb6gaUjabg2da9iabicdaWiabcUda7iaaxMaacaWLjaWaaeWaceaacqqGfbqrcqqGXbqCcqqGUaGlcqqGGaaicqqGZaWmcqqGbbqqaiaawIcacaGLPaaaaaa@5AAF@

I μ Γ 1 2 ,  for n = 1 ; ( Eq 3B ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqscqGHDisTiiGacqWF8oqBcqGHflY1cqqHtoWrdaqhaaWcbaGaeGymaedabaGaeGOmaidaaOGaeiilaWIaeeiiaaIaeeOzayMaee4Ba8MaeeOCaiNaeeiiaaIaeeOBa4Maeyypa0JaeGymaeJaei4oaSJaaCzcaiaaxMaadaqadiqaaiabbweafjabbghaXjabbccaGiabbodaZiabbkeacbGaayjkaiaawMcaaaaa@49BF@

I ( 2 λ + μ ) Γ 2 2 4 ,  for n = 2. ( Eq . 3C ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqscqGHDisTdaqadiqaaiabikdaYGGaciab=T7aSjabgUcaRiab=X7aTbGaayjkaiaawMcaamaalaaabaGaeu4KdC0aa0baaSqaaiabgkdaYaqaaiabgkdaYaaaaOqaaiabisda0aaacqGGSaalcqqGGaaicqqGMbGzcqqGVbWBcqqGYbGCcqqGGaaicqqGUbGBcqGH9aqpcqaIYaGmcqGGUaGlcaWLjaGaaCzcamaabmGabaGaeeyrauKaeeyCaeNaeeOla4IaeeiiaaIaee4mamJaee4qameacaGLOaGaayzkaaaaaa@4E54@

Here λ and μ are functions of the scattering angle θ, λ = 4 - 6 sin2 (θ) μ = 14 - sin2 (θ), and Γn are numerically evaluated integrals over the displacement pattern of the modes. The experiments were performed at a scattering angle of 90°. The Γ integrals are over the displacement pattern of the vibration. For example, the n = 0 integral is,

Γ 0 = a b ω 2 c l 2 r d r [ A J n ( ω r / c l ) + B Y n ( ω r / c l ) ] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqqHtoWrdaWgaaWcbaGaeyimaadabeaakiabg2da9maapehabaWaaSaaaeaaiiGacqWFjpWDdaahaaWcbeqaaiabikdaYaaaaOqaaiabdogaJnaaDaaaleaacqWGSbaBaeaacqaIYaGmaaaaaaqaaiabdggaHbqaaiabdkgaIbqdcqGHRiI8aOGaemOCaiNaemizaqMaemOCaiNaei4waSLaemyqaeKaemOsaO0aaSbaaSqaaiabd6gaUbqabaGcdaqadiqaaiab=L8a3jabdkhaYjabc+caViabdogaJjabdYgaSbGaayjkaiaawMcaaiabgUcaRiabdkeacjabdMfaznaaBaaaleaacqWGUbGBaeqaaOWaaeWaceaaiiaacqGFjpWDcqWGYbGCcqGGVaWlcqWGJbWycqWGSbaBaiaawIcacaGLPaaacqGGDbqxaaa@5C00@

where Jn and Yn are Bessel functions.

Applying this theoretical foundation to the M13 bacteriophage, we find several peaks in the range of frequency range of 0–25 cm-1. Figure 7 shows the modes and their Raman intensities (Gaussians broadened by 1 cm-1). Figure 7a shows the radial models, Figure 7b shows the axial models, and all modes are shown together in Figure 7c. It is important to separate the modes by their character because damping of the two types of modes is expected to be quite different. The radial modes will be highly damped since the atomic displacements of the virus move in or out radially either pushing against the water on the exterior or against the water/DNA solution in the interior cavity of the virus. Axial modes produce displacements of the atoms that are parallel to the axis of the virus and hence should not interact strongly with water. The motion relative to water of the axial modes will be shear motion (rather than compressional as for radial modes) and damping will be marginal. The modes marked with a star (*) in Fig. 7c are axial. We hypothesize that only the axial modes (those with a star) are observable in the Raman experiment because the other (radial) modes are strongly damped.

Figure 7
figure 7

Theoretical Raman Intensity of (a) radial modes; (b) axial modes; (c) all modes. Radial modes are expected to be highly damped by the solution. In (c), the undamped axial modes are indicated by a star (*).

There are two main Raman intensities predicted in Figure 7c – one at 10.5 cm-1 that is axial, and one at 20.4 cm-1 that is radial. The mode patterns of these two modes are shown in Figure 8a and 8b. The first main peak at 10.5 cm-1 is the largest and corresponds to an n = 1 axial mode of vibration for the cylindrical shell. In this mode each α-helix that comprises the protein coat undergoes axial shearing. Circling about the virus by an angle θ produces a modulation of the displacement by cos(θ) (or sin(θ), as it is doubly degenerate). The mode pattern for one plane of atoms in a plane of the virus is similar to that of a drum-head. The second major peak at 20.4 cm-1 corresponds to radial compression/expansion of the virus. The mode has n = 0, so that there is no dependence on the angle θ.

Figure 8
figure 8

Cross sectional displacement patterns for three low frequency modes with significant Raman intensity. (a) An axial mode at 10.5 cm-1. This mode we identify as the 8.5 cm-1 mode observed in experiment. (b) A radial mode at 20.4 cm-1, expected to be damped by solution. (c) An axial mode at 30.6 cm-1.

Our hypothesis for the observed Raman scattering profile is shown in Figure 9. In this figure we have only included the axial modes since they should resonate with little damping. We have included a Gaussian broadening of 5 cm-1, consistent with inhomogeneous broadening observed in the experiment. The low frequency vibrational mode observed in our Raman measurements (~8.5 cm-1) agrees remarkably well with the theoretical prediction of a single peak (~10.5 cm-1) produced by an axial modes.

Figure 9
figure 9

Final theoretical Raman intensity prediction after considering the damping effects of the solvent and the validity of the elastic continuum model.

There are many other modes of higher frequency that are predicted by the continuum model. We do not investigate frequencies higher than 25 cm-1 here because either (i) these modes correspond to n>2 and hence are not active in Raman spectroscopy, or (ii) the modes have internal wavelengths of a few Ångstroms which means that they will not be predicted reliably in a continuum model. An example of such a mode is shown in Figure 8c. This mode is just over 30 cm-1 and has an internal wavelength of about 10Å. A continuum model averages over many atoms (distances of several bondlengths) and this mode is not expected to be computed accurately. Future work will develop the model further to include atomistic details and coupling of the modes with the solvent. Such extensions will further test the limits of the continuum theory, and probe deeper into our hypothesis concerning the relative importance of damping of the various modes.

In summary, the theory has identified the low frequency modes of vibration of the M13 phage in solution and qualitatively and semi-quantitatively reproduced the observed Raman spectrum. This indicates that Raman spectroscopy is a useful, non-destructive technique for probing the low frequency modes of the virus. These modes are global, being affected by the overall conformational properties of the virus, and will be useful in providing information on how viruses are functionalized for applications in nanotechnology.

We note that Balandin's article [16] calculated the lowest vibrational modes of M13 and TMV viruses in water and air. The authors also pointed out the possible observation of these modes by micro-Raman spectroscopy. They neither (1) performed the micro-Raman experiments; nor (2) calculated/predicted the Raman intensity in a quantitative way. On the other hand, our present work is the first experimental results reported in the literature so far for the low vibrational modes of M13 phages. In addition, we have used bond-polarizability model to calculate/predict the relative strength of various Raman active modes in M13 phages. By using this theoretical model, we were able to not only explain the frequency of the mode observed in our Raman experiments but also identify it to be the axial torsion mode of the protein coat of M13 phages. Therefore, our current contributions significantly advance the work of Balandin's.


Low frequency vibrational modes of the M13 phage have been studied by Raman spectroscopy. The observed vibrational mode at 8.5 cm-1 agrees well with theoretical model calculations that are based on an elastic continuum model and a bond polarizability model with Raman selection rules. With the help of theoretical data, we have determined that the observed Raman mode corresponds to an axial (drum-head-like) vibrational mode of the M13 phage protein coat which suffers minimal damping in the exterior solvent. The differences in Raman spectra can be used to monitor the coating of virus with some other materials and nano-assembly process, such as attaching a carbon nanotube or quantum dots. Our results suggest that Raman spectroscopy is a feasible, non-destructive technique for probing the process of virus functionalization. Bacteriophage M13 is one of the simplest paradigms for viral capsids, and the characterization of its low frequency vibrational modes is a significant step towards the use of these modes for functionalization studies.


  1. Liu JL, Wu WG, Balandin A, Jin G, Wang KL: Intersubband Absorption in boron-doped multiple Ge quantum dot. Appl Phys Lett 1999, 74: 185-187. 10.1063/1.123287

    Article  CAS  Google Scholar 

  2. Liu JL, Wu WG, Balandin A, Lin G, Luo YH, Thomas SG, Lu Y, Wang KL: Observation of inter-sub-level transitions in modulation-doped Ge quantum dots. Appl Phys Lett 1999, 75: 1745-1747. 10.1063/1.124806

    Article  CAS  Google Scholar 

  3. Kamins TI, Williams RS: Lithographic positioning of self-assembled Ge islands on Si (001). Appl Phys Lett 1997, 71: 1201-1203. 10.1063/1.119625

    Article  CAS  Google Scholar 

  4. Balandin A, Jin G, Wang KL: Issues of practical realization of a quantum dot register for quantum computing. J Electronics Materials 2000, 29: 549-553.

    Article  CAS  Google Scholar 

  5. Bandyopadhyay S, Miller AE, Chang HC, Banerjee G, Yue DE, Ricker RE, Jones S, Eastman JA, Chandrasekhar M: Electrochemically assembled quasi-periodic quantum dot arrays. Nanotechnology 1996, 7: 360-371. 10.1088/0957-4484/7/4/010

    Article  CAS  Google Scholar 

  6. Balandin A, Wang KL, Kouklin N, Bandyopadhyay S: Raman spectroscopy of electrochemically self-assembled CdS quantum dots. Appl Phys Lett 2000, 76: 137-139. 10.1063/1.125681

    Article  CAS  Google Scholar 

  7. Shenton W, Douglas , Young TM, Stubbs G, Mann S: Inorganic-organic nanotube composites from template mineralzation of tobacco mosaic virus. Adv Mater 1999, 11: 253-256. 10.1002/(SICI)1521-4095(199903)11:3<253::AID-ADMA253>3.0.CO;2-7

    Article  CAS  Google Scholar 

  8. Flynn CE, Lee SW, Peelle BR, Belcher AM: Viruses asvehicles for growth, organization and assembly of materials. Acta Materialia 2003, 51: 5867-5880. 10.1016/j.actamat.2003.08.031

    Article  CAS  Google Scholar 

  9. Mao C, Solis DJ, Reiss BD, Kottmann SD, Sweeney RY, Hayhurst A, Georgiou G, Iverson B, Belcher AM: Virus-based toolkit for the directed synthesis of magnetic and semiconductoring nanowires. Science 2004, 303: 213-217. 10.1126/science.1092740

    Article  CAS  PubMed  Google Scholar 

  10. Knez M, Bittner AM, Boes F, Wege C, Jeske H, Maiss E, Kern K: Biotemplate of 3-nm nickel and cobalt nanowires. Nano letters 2003, 3: 1079-1082. 10.1021/nl0342545

    Article  CAS  Google Scholar 

  11. Knez M, Sumser M, Bittner AM, Wege C, Jeske H, Martin TP, Kern K: Spatially selective nucleation of metal clusters on the tobacco mosaic viruses. Adv Funct Mater 2004, 14: 116-124. 10.1002/adfm.200304376

    Article  CAS  Google Scholar 

  12. Tuma R, Thomas GJ: Raman Spectroscopy of Viruses, in Handbook of Vibrational Spectroscopoy. Edited by: Chalmers JM, Griffiths PR. John Wileys & Sons, Chichester; 2002.

    Google Scholar 

  13. Eisenstark A: Bacteriophage techniques. In Methods in virology. Volume 1. Edited by: Maramorosch K, Koprowski H. Academic Press, New York, N.Y; 1967:449-524.

    Google Scholar 

  14. Tsen KT, Wald KR, Ruf Tobias, Yu PY, Morkoc H: Electron-optical phonon interactions in ultrathin GaAs-AlAs multiple quantum well structures. Phys Rev Lett 1991, 67: 2557-2560. 10.1103/PhysRevLett.67.2557

    Article  CAS  PubMed  Google Scholar 

  15. Marvin DA, Welsh LC, Symmons MF, Scott WRP, Strauss SK: Molecular structure and fd (f1, M13) filamentous bacteriophage refined with respect to X-ray fibre diffraction and solid-state NMR data supports specific models of phage assembly at the bacterial membrane. J Mol Biol 2006, 355: 294-309. 10.1016/j.jmb.2005.10.048

    Article  CAS  PubMed  Google Scholar 

  16. Balandin A, Fonoberov VA: Vibrational modes ofNano-Template Viruses. Journal of Biomedical Nanotechnology 2005, 1: 90-95. 10.1166/jbn.2005.005

    Article  CAS  Google Scholar 

  17. Graff KF: Wave motion in elastic solids. Ohio State University Press, New York; 1991.

    Google Scholar 

  18. Landau LD, Lifshitz EM: Theory of Elasticity. third edition. Pergamon Press, London; 1986.

    Google Scholar 

  19. Tachibana M, Kojima K, Ikuyama R, Kobayashi Y, Ataka M: Sound velocity and dynamic elastic constants of lysozyme single crystals. Chem Phys Lett 2000, 332: 259-264. 10.1016/S0009-2614(00)01267-7

    Article  CAS  Google Scholar 

  20. Yu PY, Cardona M: Fundamentals of Semiconductors – Physics and Materials Properties. 2nd edition. Springer-Verlag, Berlin; 1999:362-371.

    Book  Google Scholar 

  21. Go S, Bilz H, Cardona M: Bond Charge, BondPolarizability, and Phonon Spectra in Semiconductors. Phys Rev Lett 1975, 34: 580-583. 10.1103/PhysRevLett.34.580

    Article  CAS  Google Scholar 

  22. Dykeman EC, Tsen KT, Sankey OFto be published

Download references


This work is supported in part by the National Science Foundation under Grant No. DMR-0305147 and by DOD RAD II STO C. The opinions or assertions contained herein are the private views of the authors and are not to be construed as official or reflecting the views of the US Department of the Army, the Uniformed Services University of the Health Sciences, or the US Department of Defense.

Author information

Authors and Affiliations


Corresponding author

Correspondence to KT Tsen.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Tsen, K., Dykeman, E.C., Sankey, O.F. et al. Observation of the low frequency vibrational modes of bacteriophage M13 in water by Raman spectroscopy. Virol J 3, 79 (2006).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: