CSULA Seminar: February 3, 2015
National Institutes of Health (NIH)
National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK)
Laboratory of Chemical Physics (LCP), Theoretical Biophysical Chemistry (TBC)
How do we predict phase separations of protein solutions?
How do we make predictions about intrinsically disordered proteins
given their large conformational landscape?
*Support provided by the Intramural Research Division of the NIDDK, NIH.
GSIGAASMEF CFDVFKELKV HHANENIFYC PIAIMSALAM VYLGAKDSTR TQINKVVRFD KLPGFGDEIE AQCGTSVNVH
SSLRDILNQI TKPNDVYSFS LASRLYAEER YPILPEYLQC VKELYRGGLE PINFQTAADQ ARELINSWVE SQTNGIIRNV
LQPSSVDSQT AMVLVNAIVF KGLWEKAFKD EDTQAMPFRV TEQESKPVQM MYQIGLFRVA SMASEKMKIL ELPFASGTMS
MLVLLPDEVS GLEQLESIIN FEKLTEWTSS NVMEERKIKV YLPRMKMEEK YNLTSVLMAM GITDVFSSSA NLSGISSAES
LKISQAVHAA HAEINEAGRE VVGGAEAGVD AASVSEEFRA DHPFLFCIKH IATNAVLFFG RCVSP
Predict structure from sequence
Native structure, folding pathways, ...
How do we predict phase separations of protein solutions?
How do we model many protein-protein interactions?
Can we predict aggregates from experimental structure?
Need a way of validating model.
An equation of state expanded in powers of density
is the pairwise interaction of two moleculesFor rotationally invariant molecules*
Goal: Develop a realistic pair potential for virial calculation.
Start with the crystallized PDB Structure
e.g. Human Serum Albumin PDB:1A06
Solve for with the Adaptive Poisson-Boltzmann Solver (APBS),
Typically (in the absence of ions), and .
Decompose the field, determine a region of excluded volume +
spherical harmonic decomposition for large distances.
Calculate the non-ideality of a protein molecule after including
both the excluded volume and electrostatics.
Predict the second-virial coefficient as a function of pH values, protein concentrations, binary mixtures, and salt concentrations.
Ongoing research: Use the model in higher-order simulations
to predict phase behavior via Gibb's ensembles.
How do we make predictions about intrinsically disordered proteins
given their large conformational landscape?
Goal: Develop a model for IDP interactions.
MJ contact energy reproduces MF energy
MF Energy distributions: Physically reasonable
Pairwise decomposition of protein complex energy; Binding affinity
Contact matrix is not symmetric
Specificity score: Define "decoys" as weakly bound
structures in protein network.
What's next? Add structure to mean field calculations.
Lattices may be optimal for IDP's, they can reproduce native-energies but quickly sample extended conformational space.