Title: Relaxation of protein structures using FIRSTFRODA
1Relaxation of protein structures using FIRST/FRODA
- Dan Farrell
- Arizona State University
2Overview
- Steric overlap in protein X-ray crystal
structures - New relaxation approach FIRST/FRODA
- Results
- Future Enhancements
3(No Transcript)
4(No Transcript)
5(No Transcript)
6Why are relaxed structures important to
FIRST/FRODA?
- FIRST
- Rigid cluster decomposition based on input
geometry - FRODA
- Locked-in steric overlap causes instability
- Cannot fully reach a target that has steric
overlap
7Why are there steric overlaps in crystal
structures?
- Hydrogen placement
- Side chain and backbone freedom
- Resolution limitations
8RelaxationTraditional Method
- Molecular Dynamics--Minimize potential energy
- Bond stretching
- Bond bending
- Dihedrals
- Sterics
- Coulomb
9Some (partisan) comments on MD versus FRODA
- MD will do nanosecs. Biological functionality
takes place millisecs. - Thus a factor 106 to go for MD. Moores law says
that computer power doubles every 2 years
therefore must wait 25 years for MD to give
answers. - Phenomenological potentials in MD dubious,
especially away from the native state. Not the
gold standard - Water makes and breaks h-bonds on nanosecs, so
huge problem for MD. For much longer times, water
flows to fill space, somewhat like a superfluid. - FRODA gives large RMSDs that correspond to
millisecs, and uses only the most important, and
non-controversial, parts of the potential. Water
can be ignored as it is diffusing so quickly.
10(No Transcript)
11(No Transcript)
12FIRST and FRODA
13FRODA Initial State
14FRODA Throw atoms
15FRODA Cycle 1 fit ghosts to atoms
16FRODA Cycle 1 move each atom along its average
mismatch vector
17FRODA Cycle 2 fit ghosts to atoms
18FRODA Cycle 2 move each atom along its average
mismatch vector
19FRODA Initial State
20FRODA no throw
21FRODA Cycle 1 fit ghosts to atoms
22FRODA Cycle 1 move each atom along its average
mismatch vector
23FRODA Cycle 2 fit ghosts to atoms
24FRODA Cycle 2 move each atom along its average
mismatch vector
25FRODA van der Waals radii
- Radii decrease
- when oppositely charged atoms interact
- donor/acceptor in potentially H-bonding geometry
- polar hydrogens reduced to 0.25 Å (except when
interacting with another positively charged atom) - Repulsive mismatch
- proportional to overlap
- drops to zero when overlap is small
26Relaxation with FRODA
- Dont count H-bonds and hydrophobic tethers as
rigid constraints. - Dont throw atoms
- Allow steric repulsion
- Allow attraction (under development)
- Keep C-alphas fixed
27(No Transcript)
28(No Transcript)
29Relaxing the PDB
- PDB Select List 1500 crystal structures, all
better than 3.0 Å resolution and 30 R-factor. - Richardson clashscore Number of atomic overlaps
4.0 Å per 1000 atoms.
30Future Enhancements
- New stickiness routines to replace current
attraction routines - Ideal ghost templates
31Applications
- Large structures 106 atoms
- NMR
- CryoEM
32Acknowledgements
- Stephen Wells
- Mike Thorpe
- Brandon Hespenheide
- Scott Menor