D-I-TASSER Q96FL8-D1

D-I-TASSER results for Q96FL8-D1

[Click result.zip to download all results on this page]

Input Sequence in FASTA format

>Q96FL8 (397 residues)
MEAPEEPAPVRGGPEATLEVRGSRCLRLSAFREELRALLVLAGPAFLVQLMVFLISFISS
VFCGHLGKLELDAVTLAIAVINVTGVSVGFGLSSACDTLISQTYGSQNLKHVGVILQRSA
LVLLLCCFPCWALFLNTQHILLLFRQDPDVSRLTQTYVTIFIPALPATFLYMLQVKYLLN
QGIVLPQIVTGVAANLVNALANYLFLHQLHLGVIGSALANLISQYTLALLLFLYILGKKL
HQATWGGWSLECLQDWASFLRLAIPSMLMLCMEWWAYEVGSFLSGILNWKKACQQAQVHA
NLKVNNVPRSGNSALPQDPLHPGCPENLEGILTNDVGKTGEPQSDQQMRQEEPLPEHPQD
GAKLSRKQLVLRRGLLLLGVFLILLVGILVRFYVRIQ

Predicted Secondary Structure

Sequence	20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 380 \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| MEAPEEPAPVRGGPEATLEVRGSRCLRLSAFREELRALLVLAGPAFLVQLMVFLISFISSVFCGHLGKLELDAVTLAIAVINVTGVSVGFGLSSACDTLISQTYGSQNLKHVGVILQRSALVLLLCCFPCWALFLNTQHILLLFRQDPDVSRLTQTYVTIFIPALPATFLYMLQVKYLLNQGIVLPQIVTGVAANLVNALANYLFLHQLHLGVIGSALANLISQYTLALLLFLYILGKKLHQATWGGWSLECLQDWASFLRLAIPSMLMLCMEWWAYEVGSFLSGILNWKKACQQAQVHANLKVNNVPRSGNSALPQDPLHPGCPENLEGILTNDVGKTGEPQSDQQMRQEEPLPEHPQDGAKLSRKQLVLRRGLLLLGVFLILLVGILVRFYVRIQ
Prediction	CCCCCCCCCCCCCCCCCHHHHCCCCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHCCCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHCCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHSSSSSSCC
Confidence	9986667666789864322211223444799999999999999999999999999999999972669367899999999999999999999999799999998658986999999999999999999999999999999999969699999999999999999899999999999999981872799999999999999999999860165511379999999999999999999865641543278548999999999999999999999999999999999981994999999999999999999999999999999999852362565567646999389999988765444545346888999999999899999999999826304677529
	H:Helix; S:Strand; C:Coil

Predicted Solvent Accessibility

Sequence	20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360 380 \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| \| MEAPEEPAPVRGGPEATLEVRGSRCLRLSAFREELRALLVLAGPAFLVQLMVFLISFISSVFCGHLGKLELDAVTLAIAVINVTGVSVGFGLSSACDTLISQTYGSQNLKHVGVILQRSALVLLLCCFPCWALFLNTQHILLLFRQDPDVSRLTQTYVTIFIPALPATFLYMLQVKYLLNQGIVLPQIVTGVAANLVNALANYLFLHQLHLGVIGSALANLISQYTLALLLFLYILGKKLHQATWGGWSLECLQDWASFLRLAIPSMLMLCMEWWAYEVGSFLSGILNWKKACQQAQVHANLKVNNVPRSGNSALPQDPLHPGCPENLEGILTNDVGKTGEPQSDQQMRQEEPLPEHPQDGAKLSRKQLVLRRGLLLLGVFLILLVGILVRFYVRIQ
Prediction	8656556541455664435444553141430140012003001111101322333233100000101531000000031323333313321311000000000101532630020011001002221222121213022003004034300510120000001212010102002100201430300000012112201310020122131101000100000122112001200021431450132011400530230000010002001212321122001012044320001001100112122030211000010011011320100001002366362034003202200113333130102101022001001011112021201010218
	Values range from 0 (buried residue) to 8 (highly exposed residue)

Predicted Contact, Hydrogen and Distance Map Used in D-I-TASSER simulation

Contact Map

Distance Map

Hydrogen bond networks

D-I-TASSER simulation is guided by the consensus contact map (left figure), distance map (middle figure) and Hydrogen bond network (right figure) derived based on confidence scores of AttentionPotential. In the contact, distance map and hydrogen bond networks, the axes mark the residue index along the sequence. For the contact map, each dot represents a residue pair with predicted contact, while for the distance map and hydrogen bond network, a color scale represents a distance of 1-20+ angstroms or a angle of 0-180 degree.

Top 10 threading templates used by D-I-TASSER

Rank

PDB
hit

ID1

ID2

Cov

Norm.
Zscore

Download
alignment

                  20                  40                  60                  80                 100                 120                 140                 160                 180                 200                 220                 240                 260                 280                 300                 320                 340                 360                 380

                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |                   |

SS
Seq

CCCCCCCCCCCCCCCCCHHHHCCCCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHCCCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCCCHHHHHHHHHCCCCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHSSSSSSCC
MEAPEEPAPVRGGPEATLEVRGSRCLRLSAFREELRALLVLAGPAFLVQLMVFLISFISSVFCGHLGKLELDAVTLAIAVINVTGVSVGFGLSSACDTLISQTYGSQNLKHVGVILQRSALVLLLCCFPCWALFLNTQHILLLFRQDPDVSRLTQTYVTIFIPALPATFLYMLQVKYLLNQGIVLPQIVTGVAANLVNALANYLFLHQLHLGVIGSALANLISQYTLALLLFLYILGKKLHQATWGGWSLECLQDWASFLRLAIPSMLMLCMEWWAYEVGSFLSGILNWKKACQQAQVHANLKVNNVPRSGNSALPQDPLHPGCPENLEGILTNDVGKTGEPQSDQQMRQEEPLPEHPQDGAKLSRKQLVLRRGLLLLGVFLILLVGILVRFYVRIQ

5xjjA

0.29

0.19

5.91

1.00

DEthreader

--------------------------SG-QLTAELKRVTRLAAPMATVTIAQYLLPVISVMVAGHNGELQLSGVALATSFTNVTGFSIMYGLVGALETLCGQAYGAKQYEKIGTYTYSAIASNIPICFIISIIWFYIENILISLGQDPDISRIAGSYAFWLIPVLFAQAIVIPLTRFLLTQGLVLPLLYTAVTTLLFHVFVCWVFVLVFVLGSNGPAMATSVSFWFYAVILSCYVRFSSCEKTRGFVSE-DFVSCVKQFFQYGVPSAAMICLEWWLFELLILCSGLPKLETS------------------------------------N------------------------------------AF--N---------------------------

5xjjA1

0.31

0.20

6.09

2.03

SPARKS-K

---------------------------SGQLTAELKRVTRLAAPMATVTIAQYLLPVISVMVAGHNGELQLSGVALATSFTNVTGFSIMYGLVGALETLCGQAYGAKQYEKIGTYTYSAIASNIPICFIISIIWFYIENILISLGQDPDISRIAGSYAFWLIPVLFAQAIVIPLTRFLLTQGLVLPLLYTAVTTLLFHVFVCWVFVLVFVLGSNGPAMATSVSFWFYAVILSCYVRFSSSCEKTRGFVSEDFVSCVKQFFQYGVPSAAMICLEWWLFELLILCSGL---------------------------------------------------------------------------------------------------------------

5y50A

0.23

0.21

6.71

1.08

MapAlign

-----------------------------GFLRETKKLSYIAGAMIAVNSSMYVLQVISIMMVGHLGELFLSSTAIAVSFCSVTGFSVVFGLASALETLCGQANGAKQYEKLGVHTYTGIVSLFLVCIPLSLLWTYIGDILSLIGQDAMVAQEAGKFATWLIPALFGYATLQPLVRFFQAQSLILPLVMSSVSSLCIHIVLCWSLVFKFGLGSLGAAIAIGVSYWLNVTVLGLYMTFSSSCSKSRATISMSLFEGMGEFFRFGIPSASMICLEWWSFEFLVLLSGKLEASVLSVCLSTQSSLYQIPESLGAAASTRVANELGAGNPKQARMAVYTAMVITGVESIMVGAIVFGARNVFGYLFSSETEVVDYVAPLLSLSVIFDALHAALSGVARGSG

5y50A

0.22

0.20

6.44

0.67

CEthreader

----------------------------DGFLRETKKLSYIAGAMIAVNSSMYVLQVISIMMVGHLGELFLSSTAIAVSFCSVTGFSVVFGLASALETLCGQANGAKQYEKLGVHTYTGIVSLFLVCIPLSLLWTYIGDILSLIGQDAMVAQEAGKFATWLIPALFGYATLQPLVRFFQAQSLILPLVMSSVSSLCIHIVLCWSLVFKFGLGSLGAAIAIGVSYWLNVTVLGLYMTFSSSCSKSRATISMSLFEGMGEFFRFGIPSASMICLEWWSFEFLVLLSGILPNPKLEASVLSVCLSTQSSLYQIPESLGAAASTRVANELGAGNPKQARMAVYTAMVITGVESIMVGAIVFGARNVFGYLFSSETEVVDYVKSMAPLLSLSVIFDALHAAL

5xjjA

0.25

0.23

7.03

1.45

MUSTER

---------------------------SGQLTAELKRVTRLAAPMATVTIAQYLLPVISVMVAGHNGELQLSGVALATSFTNVTGFSIMYGLVGALETLCGQAYGAKQYEKIGTYTYSAIASNIPICFIISIIWFYIENILISLGQDPDISRIAGSYAFWLIPVLFAQAIVIPLTRFLLTQGLVLPLLYTAVTTLLFHVFVCWVFVLVFVLGSNGPAMATSVSFWFYAVILSCYVRFSSSCEKTRGFVSEDFVSCVKQFFQYGVPSAAMICLEWWLFELLILCSGLLSN-PKLETSVLSICLTTETLHYVISSGVAA---------AVSTRVSNNLGAGNPQVARVSVLAGLCLWLRNIIGYAFSNSKEVVDYVADLSPLLCLSFILDGFTAVLNGV

5xjjA1

0.31

0.20

6.09

2.45

HHsearch

---------------------------SGQLTAELKRVTRLAAPMATVTIAQYLLPVISVMVAGHNGELQLSGVALATSFTNVTGFSIMYGLVGALETLCGQAYGAKQYEKIGTYTYSAIASNIPICFIISIIWFYIENILISLGQDPDISRIAGSYAFWLIPVLFAQAIVIPLTRFLLTQGLVLPLLYTAVTTLLFHVFVCWVFVLVFVLGSNGPAMATSVSFWFYAVILSCYVRFSSSCEKTRGFVSEDFVSCVKQFFQYGVPSAAMICLEWWLFELLILCSGL---------------------------------------------------------------------------------------------------------------

5y50A

0.24

0.21

6.42

2.65

FFAS-3D

-----------------------------GFLRETKKLSYIAGAMIAVNSSMYVLQVISIMMVGHLGELFLSSTAIAVSFCSVTGFSVVFGLASALETLCGQANGAKQYEKLGVHTYTGIVSLFLVCIPLSLLWTYIGDILSLIGQDAMVAQEAGKFATWLIPALFGYATLQPLVRFFQAQSLILPLVMSSVSSLCIHIVLCWSLVFKFGLGSLGAAIAIGVSYWLNVTVLGLYMTFSSSCSKSRATISMSLFEGMGEFFRFGIPSASMICLEWWSFEFLVLLSGILPKLEASVLSVCLSTQSS------------LYQIPESLGAAASTRVANELGAGNPKQARMAV-------------------YTAMVITGVESIMVGAIVFNVFGYLFSS--

5xjjA

0.24

0.22

6.85

1.75

EigenThreader

---------------------------SGQLTAELKRVTRLAAPMATVTIAQYLLPVISVMVAGHNGELQLSGVALATSFTNVTGFSIMYGLVGALETLCGQAYGAKQYEKIGTYTYSAIASNIPICFIISIIWFYIENILISLGQDPDISRIAGSYAFWLIPVLFAQAIVIPLTRFLLTQGLVLPLLYTAVTTLLFHVFVCWVFVLVFVLGSNGPAMATSVSFWFYAVILSCYVRFSSSCEKTRGFVSEDFVSCVKQFFQYGVPSAAMICLEWWLFELLILCSGLLSWLVESAFFSILLFTFRNIIGYAFSNSKEVVDYVADLSPLLCLSFILDGFTAVLNGVARGSGWQHIGAWNNIFSYYLVGAPVGVYLAFRGVVIGSTVQATVLAIVTASMN

5y50A

0.27

0.21

6.49

1.85

CNFpred

------------------------------FLRETKKLSYIAGAMIAVNSSMYVLQVISIMMVGHLGELFLSSTAIAVSFCSVTGFSVVFGLASALETLCGQANGAKQYEKLGVHTYTGIVSLFLVCIPLSLLWTYIGDILSLIGQDAMVAQEAGKFATWLIPALFGYATLQPLVRFFQAQSLILPLVMSSVSSLCIHIVLCWSLVFKFGLGSLGAAIAIGVSYWLNVTVLGLYMTFSSSCSKSRATISMSLFEGMGEFFRFGIPSASMICLEWWSFEFLVLLSGILP-NPKLEASVLSVCLSTQSS----------------------------------------------------------LYQIPESLGAAASTRVANELARMAVYTAMVIT

5y50A

0.27

0.19

5.71

1.00

DEthreader

---------------------------D-GFLRETKKLSYIAGAMIAVNSSMYVLQVISIMMVGHLGELFLSSTAIAVSFCSVTGFSVVFGLASALETLCGQANGAKQYEKLGVHTYTGIVSLFLVCIPLSLLWTYIGDILSLIGQDAMVAQEAGKFATWLIPALFGYATLQPLVRFFQAQSLILPLVMSSVSSLCIHIVLCWSLVFKFGLGSLGAAIAIGVSYWLNVTVLGLYMTFSSCSSRATISMS--LFEGMGEFFRFGIPSASMICLEWWSFEFLVLSGPKLEASV-----------------------------------SS-------------------------------------AFKMRG--L-----------------------

(a)	ID1 is the number of template residues identical to query divided by number of aligned residues.
(b)	ID2 is the number of template residues identical to query divided by query sequence length.
(c)	Cov is equal the number of aligned template residues divided by query sequence length.
(d)	Norm. Zscore is the normalized Z-score of the threading alignments. A Normalized Z-score >1 means a good alignment and is highlighted in bold.
(e)	Download alignment lists the threading program used to identify the template, and provide the 3D structure of aligned regions of threading templates (threading[1-10].pdb.gz).
(f)	Template residues identical to query sequence are highlighted in color.

Top 1 final models from D-I-TASSER

Click to view	Rank^a	Download	Estimated TM-score^b
	1	model1.pdb.gz	0.6953

(a)

D-I-TASSER simulations generate a large ensemble of structural conformations, i.e. decoys. These decoys are clustered by SPICKER based on pairwise structure similarity to report up to five final models from the five largest clusters. Models are ranked in descending order of cluster size. If the simulations converge well, it is possible to have less than 5 models generated, which is usually an indication of good model quality.

(b)

The model confidence is quantitatified by estimated TM-score (eTM-score), calculated based on significance of threading template alignments, contact map satisfaction rate, mean absolute error between distance of model and distance of AttentionPotential, and convergence of D-I-TASSER simulations. eTM-score is typically in the range of [0, 1], with higher eTM-score signifies higher model confidence.

Proteins with similar structure

Top 10 structural analogs in PDB (as identified by TM-align)

Rank	PDB Hit	TM-score	RMSD^a	IDEN^a	Cov.	Download Alignment
1	5y50A	0.80	2.94	0.235	0.892	model1_5y50A.pdb.gz
2	5xjjA	0.79	3.16	0.239	0.897	model1_5xjjA.pdb.gz
3	4lz6A	0.76	3.59	0.128	0.892	model1_4lz6A.pdb.gz
4	3mktA	0.75	3.49	0.146	0.884	model1_3mktA.pdb.gz
5	3vvnA	0.74	3.97	0.127	0.892	model1_3vvnA.pdb.gz
6	4hukA	0.73	3.72	0.136	0.869	model1_4hukA.pdb.gz
7	6idpA	0.73	4.04	0.137	0.884	model1_6idpA.pdb.gz
8	6fv6A	0.70	3.64	0.151	0.824	model1_6fv6A.pdb.gz
9	6cc4A	0.70	4.72	0.111	0.909	model1_6cc4A.pdb.gz
10	4z3nA	0.70	4.02	0.127	0.846	model1_4z3nA.pdb.gz

(a)	Query structure is shown in cartoon, while the structural analog is displayed using backbone trace.
(b)	Ranking of proteins is based on TM-score of the structural alignment between the query structure and known structures in the PDB library.
(c)	RMSD^a is the RMSD between residues that are structurally aligned by TM-align.
(d)	IDEN^a is the percentage sequence identity in the structurally aligned region.
(e)	Cov. represents the coverage of the alignment by TM-align and is equal to the number of structurally aligned residues divided by length of the query protein.

Predicted Gene Ontology (GO) Terms

Molecular Function (MF)

GO term	Cscore^GO	Name
GO:0015238	1.00	drug transmembrane transporter activity

Download full result of the above consensus prediction.

Click the graph to show a high resolution version.

(a)

Cscore^GO is the confidence score of predicted GO terms. Cscore^GO values range in between [0-1]; where a higher value indicates a better confidence in predicting the function using the template.

(b)

The graph shows the predicted terms within the Gene Ontology hierachy for Molecular Function. Confidently predicted terms are color coded by Cscore^GO:

[0.4,0.5)

[0.5,0.6)

[0.6,0.7)

[0.7,0.8)

[0.8,0.9)

[0.9,1.0]

Biological Process (BP)

GO term	Cscore^GO	Name
GO:0006855	1.00	drug transmembrane transport

Download full result of the above consensus prediction.

Click the graph to show a high resolution version.

(a)

Cscore^GO is the confidence score of predicted GO terms. Cscore^GO values range in between [0-1]; where a higher value indicates a better confidence in predicting the function using the template.

(b)

The graph shows the predicted terms within the Gene Ontology hierachy for Biological Process. Confidently predicted terms are color coded by Cscore^GO:

[0.4,0.5)

[0.5,0.6)

[0.6,0.7)

[0.7,0.8)

[0.8,0.9)

[0.9,1.0]

Cellular Component (CC)

GO term	Cscore^GO	Name
GO:0044464	0.75	cell part
GO:0016020	0.75	membrane

Download full result of the above consensus prediction.

Click the graph to show a high resolution version.

(a)

Cscore^GO is the confidence score of predicted GO terms. Cscore^GO values range in between [0-1]; where a higher value indicates a better confidence in predicting the function using the template.

(b)

The graph shows the predicted terms within the Gene Ontology hierachy for Cellular Component. Confidently predicted terms are color coded by Cscore^GO:

[0.4,0.5)

[0.5,0.6)

[0.6,0.7)

[0.7,0.8)

[0.8,0.9)

[0.9,1.0]

Predicted Enzyme Commission (EC) Numbers

Top 5 enzyme homologs in PDB

Rank	Cscore^EC	PDB Hit	TM-score	RMSD^a	IDEN^a	Cov.	EC Number	Predicted Active Site Residues
1	0.060	2nyaF	0.360	6.99	0.032	0.602	1.7.99.4	NA
2	0.060	1t3tA	0.359	6.57	0.051	0.574	6.3.5.3	NA
3	0.060	1vncA	0.376	6.66	0.030	0.597	1.11.1.10	76
4	0.060	3hxxA	0.356	6.82	0.047	0.577	6.1.1.7	NA
5	0.060	1ti2A	0.344	7.01	0.030	0.582	1.97.1.2	NA

	Click on the radio buttons to visualize predicted active site residues.
(a)	Cscore^EC is the confidence score for the Enzyme Commission (EC) number prediction. Cscore^EC values range in between [0-1]; where a higher score indicates a more reliable EC number prediction.
(b)	TM-score is a measure of global structural similarity between query and template protein.
(c)	RMSD^a is the RMSD between residues that are structurally aligned by TM-align.
(d)	IDEN^a is the percentage sequence identity in the structurally aligned region.
(e)	Cov. represents the coverage of global structural alignment and is equal to the number of structurally aligned residues divided by length of the query protein.

Predicted Ligand Binding Sites

Template proteins with similar binding site:

Rank	Cscore^LB	PDB Hit	TM-score	RMSD^a	IDEN^a	Cov.	BS-score	Lig. Name	Download Complex	Predicted binding site residues
1	0.01	1ce8C	0.316	6.55	0.035	0.501	0.20	ADP	complex1.pdb.gz	65,66,67,74,75,76,78,214,215
2	0.01	1bxrA	0.366	6.60	0.050	0.589	0.24	ORN	complex2.pdb.gz	68,284,285
3	0.01	1a9xC	0.310	6.63	0.042	0.491	0.21	ORN	complex3.pdb.gz	70,211,215,216
4	0.01	1ce8A	0.297	6.96	0.036	0.494	0.25	IMP	complex4.pdb.gz	64,68,69,71,291
5	0.01	1a9x2	0.310	6.63	0.042	0.491	0.21	III	complex5.pdb.gz	63,64,67
6	0.01	1t36A	0.330	7.46	0.038	0.572	0.21	U	complex6.pdb.gz	67,74,75,150,151,154
7	0.01	3bb0A	0.401	6.36	0.045	0.625	0.23	PO3	complex7.pdb.gz	60,74,75,76,152
8	0.01	1mo8A	0.177	5.68	0.021	0.259	0.18	ATP	complex8.pdb.gz	59,68,74,75
9	0.01	3n23C	0.374	6.72	0.035	0.607	0.15	OBN	complex9.pdb.gz	69,72,73

	Click on the radio buttons to visualize predicted binding site and residues.
(a)	Cscore^LB is the confidence score of predicted binding site. Cscore^LB values range in between [0-1]; where a higher score indicates a more reliable ligand-binding site prediction.
(b)	BS-score is a measure of local similarity (sequence & structure) between template binding site and predicted binding site in the query structure. Based on large scale benchmarking analysis, we have observed that a BS-score >1 reflects a significant local match between the predicted and template binding site.
(c)	TM-score is a measure of global structural similarity between query and template protein.
(d)	RMSD^a the RMSD between residues that are structurally aligned by TM-align.
(e)	IDEN^a is the percentage sequence identity in the structurally aligned region.
(f)	Cov. represents the coverage of global structural alignment and is equal to the number of structurally aligned residues divided by length of the query protein.

[Click result.zip to download all results on this page]

References:
1.	Wei Zheng, Chengxin Zhang, Yang Li, Robin Pearce, Eric W. Bell, Yang Zhang. Folding non-homology proteins by coupling deep-learning contact maps with I-TASSER assembly simulations. In preparation, 2020.