Home >
Community >
Why is it necessary to add hydrogen and delete water before protein-ligand docking?
Upvote
24
Downvote
+ Docking
+ Bioinformatics
+ Biochemistry
Posted by
Luis Becerra
Why is it necessary to add hydrogen and delete water before protein-ligand docking?
Hydrogens
Most crystallographic structures are X-ray and the hydrogens are absent. In the case of neutron diffraction models, there may be many missing high b-factor ones (example in PDB:5MOP). As a consequence models are reprotonated. Also missing carbons and missing loop residues are added for some in silico applications too (not docking).
However, it is always good to manually tweak the protonation. In the case of enzymes with protanatable side chains in the active site, some reading is required to get the charges right: by default aspartates will be unprotonated, but they may be protonated (glutamic acid) in the active site and it is important to get it right (cf.this post about it). Some times its the ligand parameterisation may give the physiologically incorrect protonation (cf.this post).
Waters
Crystallographic waters can be important catalytically and structurally, but can be pushed out and actually it is a good thing: you have a hit that take its place. Polar bonds are stronger than hydrophobic ones and most types of π bonding.
Waters in most docking software is implicit and is modelled as a parameter in the force-fields —think of it as a ghostly watery haze. So a left behind water, not only will get in the way of the docked small molecule, it may for some programs cause issues as they will try to dock it. A notable exception is WaterDock mod for AutoDock, but the jury is out on whether it is better. Water is polarisable in reality in a very complex way (protonated (OH3+) or deprotonated (OH-) with very different charges) and most water models do not cut it and instead do not account for this.
The way implicit waters are computed differs in some applications, which is also something to look into.
Most crystallographic structures are X-ray and the hydrogens are absent. In the case of neutron diffraction models, there may be many missing high b-factor ones (example in PDB:5MOP). As a consequence models are reprotonated. Also missing carbons and missing loop residues are added for some in silico applications too (not docking).
However, it is always good to manually tweak the protonation. In the case of enzymes with protanatable side chains in the active site, some reading is required to get the charges right: by default aspartates will be unprotonated, but they may be protonated (glutamic acid) in the active site and it is important to get it right (cf.this post about it). Some times its the ligand parameterisation may give the physiologically incorrect protonation (cf.this post).
Waters
Crystallographic waters can be important catalytically and structurally, but can be pushed out and actually it is a good thing: you have a hit that take its place. Polar bonds are stronger than hydrophobic ones and most types of π bonding.
Waters in most docking software is implicit and is modelled as a parameter in the force-fields —think of it as a ghostly watery haze. So a left behind water, not only will get in the way of the docked small molecule, it may for some programs cause issues as they will try to dock it. A notable exception is WaterDock mod for AutoDock, but the jury is out on whether it is better. Water is polarisable in reality in a very complex way (protonated (OH3+) or deprotonated (OH-) with very different charges) and most water models do not cut it and instead do not account for this.The way implicit waters are computed differs in some applications, which is also something to look into.
My answer is mainly based on molecular dynamics,but in this context it shouldn't make much of a difference.
The main reason has to do with the force-field-water model used during the simulation. Examples of water models are SPC, TIP3P or TIP4P; same molecule (water) but different topology.
A pdb usually comes from different sources and it can contain molecules that are not standard and therefore not recognised by the model of your simulation. Usually, once the structure is loaded, your system is re-built so that it conforms to the model topology. There are exception though, for example if your water is part of a functional active site. In that case you can keep it, but its structures (topology of the water molecule) has to be recognised by the model.
Hydrogens are usually added to "saturate" and stabilise the bonds.
My answer is mainly based on molecular dynamics,but in this context it shouldn't make much of a difference.
The main reason has to do with the force-field-water model used during the simulation. Examples of water models are SPC, TIP3P or TIP4P; same molecule (water) but different topology.
A pdb usually comes from different sources and it can contain molecules that are not standard and therefore not recognised by the model of your simulation. Usually, once the structure is loaded, your system is re-built so that it conforms to the model topology. There are exception though, for example if your water is part of a functional active site. In that case you can keep it, but its structures (topology of the water molecule) has to be recognised by the model.
Hydrogens are usually added to "saturate" and stabilise the bonds.
Totally correct for MD. However, a wee note: most ligand docking program use implicit water (except special case such as WaterDock), which means that TIP3P waters or similar are not added, but a fudge factor is used to account for the dynamic but stabilising polar properties of water.More
Upvote
VOTE
Downvote
The molecular docking programs usually tread the target structure as rigid. It does not know which atom belongs to the receptor and which is bulk water. So, if you keep the water in the structure the program cannot properly place the ligand there.
The molecular docking programs usually tread the target structure as rigid. It does not know which atom belongs to the receptor and which is bulk water. So, if you keep the water in the structure the program cannot properly place the ligand there.
Hydrogens
Most crystallographic structures are X-ray and the hydrogens are absent. In the case of neutron diffraction models, there may be many missing high b-factor ones (example in PDB:5MOP). As a consequence models are reprotonated. Also missing carbons and missing loop residues are added for some in silico applications too (not docking).
However, it is always good to manually tweak the protonation. In the case of enzymes with protanatable side chains in the active site, some reading is required to get the charges right: by default aspartates will be unprotonated, but they may be protonated (glutamic acid) in the active site and it is important to get it right (cf. this post about it). Some times its the ligand parameterisation may give the physiologically incorrect protonation (cf. this post).
Waters
Crystallographic waters can be important catalytically and structurally, but can be pushed out and actually it is a good thing: you have a hit that take its place. Polar bonds are stronger than hydrophobic ones and most types of π bonding.
Waters in most docking software is
implicitand is modelled as a parameter in the force-fields —think of it as a ghostly watery haze. So a left behind water, not only will get in the way of the docked small molecule, it may for some programs cause issues as they will try to dock it. A notable exception is WaterDock mod for AutoDock, but the jury is out on whether it is better. Water is polarisable in reality in a very complex way (protonated (OH3+) or deprotonated (OH-) with very different charges) and most water models do not cut it and instead do not account for this. The way implicit waters are computed differs in some applications, which is also something to look into.Hydrogens
Most crystallographic structures are X-ray and the hydrogens are absent. In the case of neutron diffraction models, there may be many missing high b-factor ones (example in PDB:5MOP). As a consequence models are reprotonated. Also missing carbons and missing loop residues are added for some in silico applications too (not docking).
However, it is always good to manually tweak the protonation. In the case of enzymes with protanatable side chains in the active site, some reading is required to get the charges right: by default aspartates will be unprotonated, but they may be protonated (glutamic acid) in the active site and it is important to get it right (cf. this post about it). Some times its the ligand parameterisation may give the physiologically incorrect protonation (cf. this post).
Waters
Crystallographic waters can be important catalytically and structurally, but can be pushed out and actually it is a good thing: you have a hit that take its place. Polar bonds are stronger than hydrophobic ones and most types of π bonding.
Waters in most docking software is
implicitand is modelled as a parameter in the force-fields —think of it as a ghostly watery haze. So a left behind water, not only will get in the way of the docked small molecule, it may for some programs cause issues as they will try to dock it. A notable exception is WaterDock mod for AutoDock, but the jury is out on whether it is better. Water is polarisable in reality in a very complex way (protonated (OH3+) or deprotonated (OH-) with very different charges) and most water models do not cut it and instead do not account for this.The way implicit waters are computed differs in some applications, which is also something to look into.More
VOTE
My answer is mainly based on molecular dynamics,but in this context it shouldn't make much of a difference.
The main reason has to do with the
force-field-watermodel used during the simulation. Examples of water models areSPC,TIP3PorTIP4P; same molecule (water) but different topology.A
pdbusually comes from different sources and it can contain molecules that are not standard and therefore not recognised by the model of your simulation. Usually, once the structure is loaded, your system isre-builtso that it conforms to the model topology. There are exception though, for example if your water is part of a functional active site. In that case you can keep it, but its structures (topology of the water molecule) has to be recognised by the model.Hydrogens are usually added to "saturate" and stabilise the bonds.
My answer is mainly based on molecular dynamics,but in this context it shouldn't make much of a difference.
The main reason has to do with the
force-field-watermodel used during the simulation. Examples of water models areSPC,TIP3PorTIP4P; same molecule (water) but different topology.A
pdbusually comes from different sources and it can contain molecules that are not standard and therefore not recognised by the model of your simulation. Usually, once the structure is loaded, your system isre-builtso that it conforms to the model topology. There are exception though, for example if your water is part of a functional active site. In that case you can keep it, but its structures (topology of the water molecule) has to be recognised by the model.Hydrogens are usually added to "saturate" and stabilise the bonds.
More
VOTE
VOTE
The molecular docking programs usually tread the target structure as rigid. It does not know which atom belongs to the receptor and which is bulk water. So, if you keep the water in the structure the program cannot properly place the ligand there.
The molecular docking programs usually tread the target structure as rigid. It does not know which atom belongs to the receptor and which is bulk water. So, if you keep the water in the structure the program cannot properly place the ligand there.
More
VOTE