In the correct context, molecular dynamics calculations can provide valuable information for project development, saving time and money. Molecular dynamics can be useful for assessing the energetic stability of proteins and identifying available conformations and, therefore, analyzing available binding pockets. However, the calculations are computationally expensive and require careful analysis to avoid misleading results. In order to get the most insight from molecular dynamics it’s important to have background understanding of the system of interest and of protein movements in general.
Proteins are highly dynamic systems. It‘s easy to forget that a single crystal structure is a static snapshot of a single conformation from an ensemble of available conformational states. In addition, there’s no guarantee that the crystallographically observed conformation is representative of the bioactive conformation rather than an artefact from the crystal growing process or due to crystal packing forces. In the same way the ligand binding mode, if present, may not be truly representative of the biologically relevant conformation.
Due to the complexity of protein and ligand movements, it would not be possible to generate all the available biologically relevant conformational states from the crystallographically observed conformational state(s). Even when there are a series of related crystal structures in differing conformational states, this still does not give the complete picture of all the available protein conformations or indeed any information on the protein energy landscape, so the energy barriers and population states of the observed conformations remain unknown. This is when molecular dynamics can help fill in the missing pieces of the conformational puzzle and protein energy landscape.
Figure 1: Crystal structures 1sw2 (top right) and 1sw5 (top left) modeled in Flare™ showing the protein movement observed upon ligand binding for the osmoprotection protein ProX. The movement is in the relative position of the two domains rather than a restructuring of either domain.
Molecular dynamics calculations help to map the energy landscape by solvating the protein with explicit water molecules then adding energy into the system to allow it to move, generating an ensemble of protein structures. During this process snapshots of the protein are taken at predefined time steps, calculating a series of protein structures which are likely to be representative of the biologically relevant protein conformations. Analysis of these conformations from a molecular dynamics calculation can enable the assessment of protein movements and the identification of different lower energy conformational states for the system. Because there are explicit water molecules in the system, it’s possible to identify features like bridging water molecules and those water molecules which are held rigid by the protein and, therefore, should be treated as structural features.
Like all molecular mechanics methods, the quality of the simulation depends on many factors, but two key determining factors are the quality of the force field used and the quality of the starting system, therefore, correct preparation the protein system is critical.
Currently a significant use for molecular dynamics is in the field of Free Energy Perturbation calculations (FEP). The current fashion for using pure molecular dynamics calculations is to run multiple simulations to generate a collection of energy landscapes for exploring the available conformations. This type of approach requires a large CPU or GPU farm to distribute hundreds or thousands of discrete calculations. Generating a large number of structures from a collection of discrete molecular dynamics calculations allows for more confidence in the obtained results and less susceptibility to the numerical noise that can occur in a single calculation. This approach can provide a greater understanding of the available conformations but requires extensive time, computational resource and monetary commitment. The availability of cloud computing has opened the ability to run this type of large-scale molecular dynamics calculation to a larger market, but the time and resource required to fully analyze the calculation is still extensive.
Often in a series of crystal structures for a protein with a set of closely related ligands bound, the only movement observed is in those residues directly interacting or adjacent to the bound ligand. These series of crystal structures help us to assess the side-chain mobility in the active site and can assist with improving design of the current ligand series but are of less use when trying to exploit unexplored regions of the active site to develop novel series.
Systems where the active site is rigid, with only side chain movements and minimal backbone movements, can be relatively simple to investigate using a constrained molecular dynamics calculation. In these calculations only the side chains of the residues in the active site and slightly more distance co-operative residues are allowed to move. The backbone atoms and remote side chains can be either restrained or completely fixed. This allows a rapid investigation of the side chains which are known to adapt upon ligand binding. This type of calculation can also once again provide useful information on water molecules observed in or around the active site in the crystal structure, allowing for the identification of those water molecules which are stable and should be considered as structural, as opposed to those which are transient.
Figure 2: 1fax, 1mtw, 1x7a, 3lc5, 4yzu, 4z0k and 5tnt crystal structures for Factor IXa, modeled in Flare , show the overall structural conservation but the more subtle protein movement around the active site to accommodate the different ligand classes observed.
A more complex example is when there is movement in the secondary or tertiary structure of the protein and the ligand occupies a different conformational state of the protein from that available in the crystal structures. The typical example of this is in kinase systems with the DFG loop in or out conformation, but even with these reasonable well represented systems in the PDB database the crystal structures are not always available in the correct conformation states.
Figure 3: DFG In and Out conformations for ABL kinase, modeled in Flare , using crystal structures 2gqg and 3qrk.
For relatively well understood and highly represented crystallographic systems like kinases it’s fairly easy to find a suitable template from which to construct a homology model of the correct conformational state. However, for less well represented systems it’s not always possible to have this wealth of structural information available.
The investigation of this type of system firstly requires an understanding of the system to appreciate the potential protein movements and conformational states available. Comparison to related family members can help to elucidate conserved protein movements and assist with this type of analysis. These identified conformational states can then be further investigated by the employment of molecular dynamics calculations to assess protein stability and the energy profile of these different conformational states and transitions between them.
Molecular dynamics is also an invaluable tool when performing homology modeling as homology models tend to be in a high energy state when first created. The regions around loop insertions and deletions and side chains generally tend to be in sub-optimum conformational states and local energy minimisation only moves these groups into the nearest energy minimum. Molecular dynamics calculations, when used properly, can be used to smooth out the energy surface and let the whole system relax from a higher energy state to a more stable lower energy state, which can be more representative of the biologically relevant system. The stability of the homology model during the molecular dynamics calculation also gives an indication of the overall quality of the homology model. Those homology models which are of lower quality tend not to survive molecular dynamics processing with their secondary or tertiary structure intact.
We help you reach your next milestone faster and more cost effectivelyContact us for a free confidential discussion