Simulations of impulsive laser scattering of biological protein assemblies: application to M13 bacteriophage

We develop a theoretical framework, based on a bond-polarizability model, for simulating the impulsive force experienced on a protein or an assembly of proteins from a pulsed light source by coupling the laser electric field to an atomic distortion. The mechanism is impulsive stimulated Raman scattering (ISRS) where mechanical distortions produce variation in the electronic polarization through atomic displacements similar to vibrational Raman scattering. The magnitude of the impulsive force is determined from the empirical two-body bond-polarizability model and the intensity of the incident light. We apply the method to the M13 bacteriophage protein capsid system by performing several classical molecular-dynamics simulations that include the additional impulsive laser scattering force at various light intensities and pulse widths. The results of the molecular-dynamics simulations are then qualitatively interpreted with a simple harmonic oscillator model driven by ISRS. The intensity of light required to produce damage to the capsid in the simulations was found to be far higher than what was found in recent pulsed laser scattering experiments of M13 phage, suggesting that the observed inactivation of viruses with ultrashort laser pulses involves processes and/or mechanisms not taken into account in the present simulations.