Skip to content

Commit 1bd7f3a

Browse files
committed
Added documentation for restraints
1 parent 789d298 commit 1bd7f3a

File tree

1 file changed

+16
-4
lines changed

1 file changed

+16
-4
lines changed

docs/README.md

+16-4
Original file line numberDiff line numberDiff line change
@@ -104,22 +104,34 @@ After the execution of `lightdock_setup` script, several files and directories w
104104

105105
### 2.2. Distance restraints (Experimental)
106106

107-
In version `0.5.1`, distance restraints on the receptor surface have been implemented. This new feature works in the following way:
107+
In version `0.5.2`, distance restraints on both receptor and ligand have been implemented. This new feature works in the following way:
108108

109109
From a `restraints_file` file containing the following information:
110110

111111
```
112112
R A.SER.150
113113
R A.TYR.151
114+
L B.LYS.2
115+
L B.ARG.3
114116
```
115117

116-
Where `R` means for `Receptor` (at the moment only residue restraints on the receptor are accepted) and then the rest of the line is a residue identifier in the format `Chain.Residue_Name.Residue_Number` in the original PDB file numeration.
118+
Where `R` means for `Receptor` and `L` means for `Ligand`, and then the rest of the line is a residue identifier in the format `Chain.Residue_Name.Residue_Number` in the original PDB file numeration.
117119

118-
For each of these residues, only the closest swarms are considered during `lightdock_setup`. At the moment, only the **10 closest swarms** for each residue are considered. If some of the swarms overlap, then only one copy of that swarm is used.
120+
For each of residue restraints specified for the receptor, only the closest swarms are considered during `lightdock_setup`. At the moment, only the **10 closest swarms** for each residue are considered. If some of the swarms overlap, then only one copy of that swarm is used.
121+
122+
Then, the simulation will try to optimize both energy and restraints satisfied.
123+
124+
The following scoring functions are prepared to be used with restraints:
125+
126+
* [DFIRE](../lightdock/scoring/dfire)
127+
* [FastDFIRE](../lightdock/scoring/fastdfire)
128+
* [CPyDock](../lightdock/scoring/cpydock)
129+
130+
Once the simulation has ended, the script `lgd_filter_restraints.py` should be used in order to remove predictions which don't satisfy the given restraints.
119131

120132
### 2.3. Tips and tricks
121133

122-
- As a general rule of thumb, the receptor structure is the bigger molecule and the ligand the smaller. With bigger and smaller, the metric used is the longest diameter of the minimum ellipsoid containing the molecule. A script called `lgd_calculate_diameter.py` can be found in `$LIGHTDOCK_HOME/bin/support` path in order to calculate an approximation of that metric.
134+
- As a general rule of thumb, the receptor structure is the bigger molecule and the ligand the smaller. For this concept of size, the metric used is the longest diameter of the minimum ellipsoid containing the molecule. A script called `lgd_calculate_diameter.py` can be found in `$LIGHTDOCK_HOME/bin/support` path in order to calculate an approximation of that metric.
123135

124136
- If the `init` directory exists, `lightdock_setup` makes use of the information contained in that directory.
125137

0 commit comments

Comments
 (0)