Skip to content

Commit 54dcfa6

Browse files
committed
New tutorial chapter
1 parent e42737d commit 54dcfa6

1 file changed

Lines changed: 36 additions & 0 deletions

File tree

structure/crystal-contacts.md

Lines changed: 36 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,36 @@
1+
# How to calculate all crystal contacts in a PDB structure
2+
3+
## Why crystal contacts?
4+
5+
A protein structure is determined by X-ray diffraction by producing a crystal - an infinite lattice of molecules - of the protein. Thus the end result of the diffraction experiment is a crystal lattice and not just a single molecule. However the PDB file only contains the coordinates of the Asymmetric Unit, defined as the minimum unit needed to reconstruct the full crystal using symmetry operators.
6+
7+
[here](http://www.wwpdb.org/news/news_2013.html#22-May-2013)
8+
9+
10+
11+
## Getting the set of unique contacts in the crystal lattice
12+
13+
This code snippet will produce a list of all non-redundant interfaces present in the crystal lattice of PDB entry [1SMT](http://www.rcsb.org/pdb/explore.do?structureId=1SMT):
14+
15+
```java
16+
AtomCache cache = new AtomCache();
17+
18+
StructureIO.setAtomCache(cache);
19+
20+
Structure structure = StructureIO.getStructure("1SMT");
21+
22+
CrystalBuilder cb = new CrystalBuilder(structure);
23+
24+
StructureInterfaceList interfaces = cb.getUniqueInterfaces(6);
25+
26+
interfaces.calcAsas(3000, 1, -1);
27+
28+
// now interfaces are sorted by areas, we can get the largest interface in the crystal and look at its area
29+
interfaces.get(1).getTotalArea();
30+
31+
```
32+
33+
An interface is defined here as any 2 chains with at least a pair of atoms within the given distance cutoff (6 A in the example above)
34+
35+
36+

0 commit comments

Comments
 (0)