"Development of an atomistic model of the entire PCV2 virus using a molecular dynamics simulation has allowed us to model the structure and dynamics of the PCV2 virus test system over a microsecond timescale."
Want HECtime on ARCHER?
Interested in how to apply for time on ARCHER like these researchers did?
You can find out more here
State-of-the-art computers can model very large biomolecular systems (e.g. organelles, cellular membranes, entire viruses) at atomistic resolution using highly efficient supercomputers and complementary numerical methods (‘Molecular Dynamics’ (MD)) that can simulate molecular systems of millions of atoms in size; and knowledge of the average 3D structural coordinates for almost all atoms in the system derived from high resolution experimental techniques, e.g. X-ray crystallography and cryo-electron microscopy. All-atom MD computer models can provide information about, and unique insights into, complex biomolecular systems which is impossible to obtain by any other means; e.g., the reconstruction of the parts of the molecular system that cannot be revealed experimentally (typically flexible or variable parts of the structure), and an understanding of dynamics at physiological temperature. Simulation of viruses is most realistic because, in contrast to other large biomolecular systems, viruses are self-contained biological units, which exist in isolation from the rest of the organism (although they cannot reproduce in isolation).
Computer models of cellular organelles, for example, consider only part of the system that interacts with the rest of the cell via complicated, poorly understood mechanisms making comparison with the experimental results difficult. For viruses, interaction with aqueous solution is the only external force that defines the structure and dynamics and modelling of water is well developed in MD. The main difficulty in computer modelling of entire viruses is the number of atoms in the model, which is very large (at the edge of feasibility of such simulations). This is because the surrounding water should be included in the system explicitly as water molecules, (not as a continuum), and the number of water molecules can reach 95% of the number of atoms in the whole system. It has been recognised that the water surrounding biomolecules plays the major role in controlling the dynamics, emphasising certain motions whilst restricting others.
Antimicrobial resistance is a growing problem for many bacterial infections, including those produced by E. coli. For E. coli, a well-studied phage is the (+)ssRNA phage MS2. Its genome of 3569 base pairs and a small icosahedral capsid (~27nm in diameter) makes it a suitable candidate for all-atom MD simulations. MS2 has been used as a carrier for drug molecules and as such could be used in non-conventional phage therapy to tackle antimicrobial resistance. MS2 is not the best candidate for phage therapy, but it has all relevant properties, is well studied and relatively small. We will use it as a proof-of-concept system to develop methodologies which can be applied to other phages.
Enabled formulation of recommendations for applying the approach to larger viruses for effective phage therapy method developments.
Fostered collaboration with Odessa University (Ukraine), Centro Nacional de Biotecnologia (Spain).
Allowed training of PhD students.
The goal of this project was to develop an atomistic model of the entire MS2 virus. The developed model can be used to understand how alterations to the genome enhance their potential antimicrobial resistance properties as phage therapy agents. Packaging the genome inside the capsid problem will also need to be considered. The genome structure is more difficult to measure experimentally, in contrast to the capsid. Even though significant work has been done, unconventional approaches will be needed.
ARCHER resources have been used for classical MD simulations of a very large biological system where all atoms of the system were included. We have completed the simulations of a different virus, PCV2, which is smaller than the target virus MS2. Our current work concentrates on MS2.
RIKEN Quantitative Biology Center - Japan
Prof. Makoto Taiji group
City College of New York - US
Prof. Reza Khayat group
Aston University - UK
Dr Dmitry Nerukh
Dr Michael Stich
This study made use of a HECBioSim ARCHER project allocation EP/L000253/1