OmegaFold icon indicating copy to clipboard operation
OmegaFold copied to clipboard

reproduce result of omegafold and find a different result

Open pengzhangzhi opened this issue 2 years ago • 8 comments

Hi. I am trying to reproduce the result shown in your paper, specifically, the result in Fig.2, as shown in the figure below. image

However, the CDRH3 RMSD is 3.6Å in our test, which is significantly larger than your 0.38 Å. The visualized structure between the native and predicted shows the deviation in CDRH3 region, as circled in the red line. 49ca1d074119b4e6c9b193101634cd6

I bet there are something wrong with our program. To ensure fair reproduce, I hope you can help me find out what is went wrong.

We use the omegafold from this repo and infer the structure of 7kpj_2_C, this is the heavy chain of the 7KPJ. We append the predicted and native pdb to this issue. pdb.zip

pengzhangzhi avatar Aug 10 '22 02:08 pengzhangzhi

let me know if you want know more detail of our reproduction.

pengzhangzhi avatar Aug 10 '22 02:08 pengzhangzhi

Hi,

Thanks a lot for your reproduction and reply!

The light chain part will influence CDR H3, so in paper we test our model with another ensemble model with complex modeling ability.

For that model we train on multi-chain data and inference with light and heavy chain together. This helps loop modeling.

We will release more ensemble code and model soon. Stay tuned!

mooninrain avatar Aug 10 '22 09:08 mooninrain

Thanks! That explains a lot. Does the ensemble model consist of multiple model or just one with more sophisticated training?

pengzhangzhi avatar Aug 10 '22 16:08 pengzhangzhi

Thanks! We will release multiple models.

There are several ways to encourage divergent predictions, including changing the model architecture, training curriculum and data sampling algorithm. So these ensembles can be helpful sometimes.

mooninrain avatar Aug 10 '22 16:08 mooninrain

Hi. Another question, Do you evaluate the RMSD on full-atoms or just backbone? Some methods in the antibody structure prediction only do backbone evaluation.

pengzhangzhi avatar Aug 16 '22 10:08 pengzhangzhi

You can predcit the antibody complexes by adding a long linker like 'GGGGGGGGGGGGGGGGGGGGGGGGGGGGGG' between different chains now.

ijayden-lung avatar Aug 25 '22 02:08 ijayden-lung

GGGGGGGGGGGGGGGGGGGGGGGGGGGGGG

So that you would know the partition of each chain by use the linker as a sign right?

pengzhangzhi avatar Aug 25 '22 02:08 pengzhangzhi

Sure.

ijayden-lung avatar Aug 26 '22 02:08 ijayden-lung