Thank you all for your comments! Let me give you more details about my case:<br><br>I want to study the interaction modes of known inhibitors with the monomeric state of an IDP (~100 aa). IDPs only assume secondary structure when in complex with their partners, hence there are no crystal structures of the protein as a monomer. The only experimental data that exist are NMR and CD data of a close homologue (~78% sequence identity) which show that the monomer is partially folded in the following simplified pattern:<br>
<br>CHHHHHHHHCCCCCCCChhHHHhhhC<br><br>where "h" denotes transient helix, "H" relative rigid helix, and "C" random coil. The whole IDP is ~100 aa, but we have indications that the inhibitors bind to the central disordered 40 aa. Therefore I was thinking of running two simulations, one for the whole ~100 aa using secondary restraints wherever applicable, and a longer one for the central 40 amino acids. The problem with the whole-IDP simulation is that the unfolded ~100 aa will occupy substantially bigger space than the central 40 aa, therefore the dimensions of the box will be larger and will include much more solvent. This will make the simulation cumbersome and I won't be able to do enough sampling for the whole IDP.<br>
<br>In both simulations I am thinking of keeping the desired amino acids from the crystal structure of the IDP in complex with its partner, carry out REMD to unfold it, and then add an inhibitor to see where it binds. <br>
<br><br>You comments about the proposed protocol, force field and water model to use, will be highly appreciated!<br><br>thanks,<br>Thomas<br><br><br><br><div class="gmail_quote">On 8 June 2011 22:37, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im"><br>
<br>
Da-Wei Li wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
I really don't think you can get adequate sampling for IDPs that have 40 residues, using full atomic MD.<br>
<br>
</blockquote>
<br></div>
I disagree. Perhaps brute force MD would not accomplish the task (unless you have considerable resources and don't want your answers very quickly, but even then...), but there are certainly a number of techniques, including implicit solvent, REMD, generalized Hamiltonian replica exchange, and fancy umbrella sampling and/or free energy settings that can vastly improve sampling for such systems.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
<br>
On Wed, Jun 8, 2011 at 3:25 PM, Michael Daily <<a href="mailto:mdaily.work@gmail.com" target="_blank">mdaily.work@gmail.com</a> <mailto:<a href="mailto:mdaily.work@gmail.com" target="_blank">mdaily.work@gmail.com</a>>> wrote:<br>
<br>
Do you have some experimental data to compare to your IDP<br>
simulations, like X-ray scattering or some such? I'd imagine that<br>
IDP simulations with either forcefield would only be qualitatively<br>
accurate given that the forcefields are calibrated, as you say, on<br>
rigid proteins and small molecules.<br>
<br>
On 6/8/11 8:00 AM, Thomas Evangelidis wrote:<br>
</div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
Dear Prof van der Spoel and GROMACS users,<br>
<br>
I have read the article "Scrutinizing Molecular Mechanics Force<br>
Fields..." where it is demonstrated that AMBER99sb yields protein<br>
conformations that are in better agreement with residual dipolar<br>
coupling, J-coupling and NOE data, compared with other force<br>
fields. However, both proteins used for this benchmark study<br>
(ubiquitin and protein G) are rather rigid, so I was wondering if<br>
there is a similar analysis for flexible proteins/peptides. I want<br>
to simulate a few intrinsically disordered proteins/peptides of<br>
length between 40 and 100 aa and would like to know what would be<br>
the best choice of force field to use. Any experience or knowledge<br>
on this matter would be greatly appreciated!<br>
<br>
thanks in advance,<br>
Thomas<br>
<br>
<br>
<br>
On 27 May 2011 19:01, David van der Spoel <<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a><br></div><div><div></div><div class="h5">
<mailto:<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a>>> wrote:<br>
<br>
On 2011-05-27 17.50, simon sham wrote:<br>
<br>
Hi,<br>
I have recently done two simulations on a protein at high<br>
temperature<br>
near its melting temperature. At the beginning I used<br>
CHARMM forcefield,<br>
and then OPLSAA to double check the results. There are<br>
some differences<br>
in the structures between the forcefield used. I know<br>
different<br>
forcefields can give different results. All parameters in the<br>
simulations were the same except the forcefield. Is there<br>
anyway I can<br>
tell which forcefield gives more reliable results?<br>
<br>
Thanks for the insights,<br>
<br>
Simon<br>
<br>
You might want to check<br>
<br>
Oliver Lange, David van der Spoel and Bert de Groot:<br>
Scrutinizing Molecular Mechanics Force Fields on the<br>
Microsecond Timescale With NMR Data Biophys. J. 99 pp. 647-655<br>
(2010)<br>
<br>
where we compare a number of FFs to NMR data.<br>
<br>
</div></div></blockquote></blockquote></blockquote></div>