Observation of the low frequency vibrational modes of bacteriophage M13 in water by Raman spectroscopy
© Tsen et al; licensee BioMed Central Ltd. 2006
Received: 23 August 2006
Accepted: 22 September 2006
Published: 22 September 2006
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 [7–11]. 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 . 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  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 . 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
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.
We model the long fiber as a solid cylindrical tube of elastic continuous material, as shown in Figure 6c. Following Balandin et al.  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 e ik z Z . 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 (r, θ) with no z dependence (e ik z Z = 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 , σ rθ , 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 , 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 
where , describe the polarization of the incident and scattered light. The Raman tensor Δ 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  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 () and perpendicular () 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
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 and , 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(nθ)or sin(nθ), 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,
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,
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  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.
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.
- 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.123287View ArticleGoogle Scholar
- 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.124806View ArticleGoogle Scholar
- Kamins TI, Williams RS: Lithographic positioning of self-assembled Ge islands on Si (001). Appl Phys Lett 1997, 71: 1201-1203. 10.1063/1.119625View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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/010View ArticleGoogle Scholar
- 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.125681View ArticleGoogle Scholar
- 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-7View ArticleGoogle Scholar
- 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.031View ArticleGoogle Scholar
- 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.1092740View ArticlePubMedGoogle Scholar
- 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/nl0342545View ArticleGoogle Scholar
- 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.200304376View ArticleGoogle Scholar
- 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
- 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
- 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.2557View ArticlePubMedGoogle Scholar
- 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.048View ArticlePubMedGoogle Scholar
- Balandin A, Fonoberov VA: Vibrational modes ofNano-Template Viruses. Journal of Biomedical Nanotechnology 2005, 1: 90-95. 10.1166/jbn.2005.005View ArticleGoogle Scholar
- Graff KF: Wave motion in elastic solids. Ohio State University Press, New York; 1991.Google Scholar
- Landau LD, Lifshitz EM: Theory of Elasticity. third edition. Pergamon Press, London; 1986.Google Scholar
- 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-7View ArticleGoogle Scholar
- Yu PY, Cardona M: Fundamentals of Semiconductors – Physics and Materials Properties. 2nd edition. Springer-Verlag, Berlin; 1999:362-371.View ArticleGoogle Scholar
- 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.580View ArticleGoogle Scholar
- Dykeman EC, Tsen KT, Sankey OFto be published
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.