IUCrJ

Human cytochrome P450 2C9 (CYP2C9) is a hepatic microsomal enzyme involved in the oxidative metabolism of clinically important drugs, but the structural organization of its oligomeric assemblies outside crystallographic packing environments remains poorly understood. Here, we report the cryo-EM structure of human CYP2C9 determined under aqueous, membrane-free conditions at 3.31 Angstrom resolution. The structure reveals a C2-symmetric hexameric assembly organized as a dimer of trimers. Individual protomers retain the conserved P450 fold and heme-binding architecture observed in previously reported crystal structures, indicating that assembly formation does not substantially perturb the catalytic core. The hexamer is stabilized by defined intra-trimer interfaces involving the N-terminal region and residues around Trp212 and Phe482, together with inter-trimer interfaces involving Leu71 and the 220-227 loop. These interfaces are distinct from the crystal packing contacts observed in CYP2C9 crystal structures, demonstrating that the assembly is not a simple recapitulation of crystallographic packing. Notably, the inter-trimer interface is located near the FG-loop-containing surface previously implicated in membrane association. This suggests that the observed hexamer may represent a membrane-free association of two trimers through membrane-related surfaces, whereas the trimeric arrangement itself may be compatible with membrane-associated organization. The structure therefore provides a framework for investigating how trimer formation, membrane interaction and local conformational changes in the FG-loop region may influence CYP2C9 function.

Cryo-electron tomography (cryo-ET) enables near-native visualization of subcellular architectures, yet applying it to moderately thick, multilayered tissues such as the retina is hampered by inadequate vitrification and inaccurate depth-targeting. Here, we developed PLCT, an integrated approach combining modified high-pressure freezing, cryo-ultramicrotome trimming, and plasma-based cryo-FIB milling to overcome these barriers. PLCT reliably vitrified <100 m retinal strips with minimal ice artifacts, navigates precisely to the outer plexiform layer using morphological landmarks, and produces high-quality lamellae suitable for high-resolution cryo-ET. Subtomogram averaging (STA) analysis identified microtubules at 16.33 [A] within retinal horizontal cell processes. Importantly, STA also resolved a 10-nm-diameter filamentous structure at 24.81 [A] in the same processes, featuring six peripheral strands surrounding an elongated central density with continuous intervening cavities, an architecture consistent with intermediate filaments. Together with its native localization and immunoreactivity, these features collectively identify the filaments as neurofilaments. Separately, 3D reconstruction of synaptic ribbons uncovered a previously unrecognized "mahjong tile"-like fine ultrastructure. These results demonstrate that PLCT-produced lamellae are of sufficient quality to support structural analysis in native tissue. Although demonstrated on retinal photoreceptor synapses as a proof-of-principle, PLCT is inherently generalizable, with its depth-navigation and vitrification strategies directly applicable to any multilayered tissues. This work establishes PLCT as a robust, reproducible platform for depth-resolved in situ cryo-ET of multilayered tissues.

Ammonia oxidation initiates nitrification and is closely linked to microbial N2O production. Ammonia monooxygenase (AMO) catalyzes the first and rate-limiting step of nitrification and is widespread across evolutionarily distinct ammonia-oxidizing archaea (AOA) and bacteria (AOB). The ocean is the largest biome for AOA and AOB, which have distinct ecological niches and markedly different sensitivities to nitrification inhibitors. However, the lack of archaeal AMO structures and inhibitor-bound AMO complexes has hindered mechanistic understanding of the architectural, catalytic, and inhibitory divergence between these two enzyme systems. Here, we report high-resolution cryo-electron microscopy (cryo-EM) structures of marine archaeal AMO captured in active and inactivated states within its native membrane environment, together with inhibitor-bound structures of estuarine bacterial AMO. Archaeal AMO forms an unexpected cup-shaped homotrimer composed of eight subunits per protomer and exhibits substantial architectural divergence from bacterial AMO. Integrated structural, biochemical, kinetic, and computational analyses reveal distinct periplasmic architectures, copper-center organization, and hydrophobic channels between archaeal and bacterial AMOs for ammonium acquisition, catalysis and inhibitor response. These findings provide a structural and mechanistic framework for understanding how archaeal and bacterial AMOs have diverged to distinct ammonia-oxidizing strategies and inhibitor susceptibilities across environmentally important ammonia oxidizers.

This paper presents and critically reviews the results of original and some literature based experimental studies conducted by the authors last years on the structural organization of DNA in dormant (starvation stress), anabiotic dormant (4 HR treatment) E. coli cells, as well as the K12 {Delta}dps strain, which lacks the Dps protein (Dps null E. coli). The experimental data includes small-angle synchrotron radiation diffraction (SAXS) and transmission electron microscopy (TEM) data. Synchrotron radiation diffraction experiments on K12{Delta}dps cells allowed us to conclude that peaks at 44.3, 22.1, and 14.8 angstrom resolutions are associated exclusively with ordered DNA organization. Peaks at 44.3, 22.1, and 14.8 angstrom resolutions are also observed for samples of dormant (starvation stress) cells and anabiotically dormant cells. Therefore, this ordered DNA organization also applies to samples of dormant and anabiotically dormant cells. A model is proposed that considers the ordered DNA organization in the cell as a cholesteric liquid crystal. The powder diffraction pattern calculated based on this model is compared with experimental small angle X ray scattering (SAXS) data obtained on Dps-null cell samples. The model completely reproduces the key features of the experimental diffraction pattern from Dps-null cell samples. Accordingly, the cholesteric liquid crystal model corresponds to DNA packaging in dormant and anabiotically dormant cells. Cholesteric liquid crystal ordering should be further considered in all models of cellular DNA packaging. To address the question of which structural organization of DNA predominates in the cell: the cholesteric liquid crystal or nanocrystalline or whether they coexist and fully manifest themselves under different external conditions, it is necessary to utilize the latest methodological advances in structural analysis.

Deep learning-based protein structure prediction methods that leverage evolutionary information from multiple sequence alignments (MSAs), exemplified by AlphaFold2, have achieved remarkable accuracy. However, existing methods still struggle to predict challenging proteins, particularly those with novel folds or limited evolutionary information, and to recover alternative conformational states. Here we show that structure prediction models trained under different MSA-depth distributions corresponding to different levels of evolutionary information exhibit complementary generalization behaviors, and that a model trained on a mixture of these distributions can combine their complementary generalization strengths. Building on this insight, we developed ProtMonomer, a deep learning framework trained on MSA-depth distributions representing a broad range of evolutionary information levels to improve structure prediction. Across benchmarks comprising CASP15 targets, non-redundant experimentally determined structures, orphan proteins, and short peptides, ProtMonomer performed comparably to or better than leading methods, including AlphaFold2 and AlphaFold3, with particularly strong performance on challenging targets. For fold-switching proteins, ProtMonomer also recovered alternative conformational states more accurately than AlphaFold2 and AlphaFold3 across diverse homologous sequence sampling strategies. In addition to improving predictive accuracy, ProtMonomer substantially reduced inference cost through an efficient architecture, enabling high-throughput applications. Together, these findings provide insights into the generalization of evolution-informed structure prediction models and support ProtMonomer as an accurate and efficient framework for protein structure prediction.

Spatial-omics workflows enable molecular analysis within tissue spatial context. Despite the prognostic value of tissue stiffness, these approaches have not incorporated direct, mechanical measurements. This omission reflects several challenges, including sample requirements, low throughput, specialized equipment, and complex data registration. Here, we introduce mechanoMaST (mechanics mapped to spatial transcriptomics), the first workflow to combine absolute mechanical measurements with spatial-omics. It pairs atomic force microscopy-based nanoindentation stiffness maps with spatial transcriptomics maps from adjacent tissue cryosections. The two modalities are then computationally co-registered to enable direct spatial correlation at 100 um resolution, with mapping accuracy quantified through error propagation, providing ground-truth mechanical data directly linked to spatial gene expression. We demonstrate mechanoMaST in human colorectal cancer liver metastasis, generating a spatial resource from 10 patients and revealing a four-gene stiffness signature. mechanoMaST is readily adaptable to other tissues across development and disease, and extendable to additional spatial-omics modalities in adjacent sections.

Despite the well-known health risks of neutron exposures, key gaps remain in understanding neutron-induced molecular responses and identifying reliable biodosimetric markers that distinguish neutrons from photon exposure. We provide the first genome-wide analysis of the human blood transcriptional response to an accelerator-derived fission-like spectrum of neutrons versus photons, evaluating transcriptomic relative biological effectiveness (RBE) and radiation quality-discriminating gene signatures. Whole blood from healthy donors was irradiated ex vivo with X-rays (140 kV, 0-4 Gy, n = 3) or neutrons (0.1-8 MeV, 0-1 Gy, n = 2), incubated for 6 h or 24 h, and processed for RNA sequencing from peripheral blood mononuclear cells (PBMCs). Neutrons were markedly more potent than X-rays at inducing differentially expressed genes (DEGs) at equal doses, showing a peak response 6 h post-irradiation followed by a decline. In contrast, X-rays caused a continuous increase in DEGs up to 24 h (neutrons vs. X-rays at 1 Gy: 1,449 vs. 121 DEGs at 6 h; 996 vs. 621 DEGs at 24 h). A universal p53-centered 34-gene signature, including FDXR, EDA2R, GADD45A, and ZMAT3, showed highly monotonic dose responses (Spearman correlation coefficient {approx} 1) across donors, radiation qualities, and timepoints. Additionally, difference-in-differences analysis identified radiation quality-discriminating genes only at 6 h, with transcriptional convergence observed by 24 h, suggesting a very narrow time window for biodosimetric differentiation. We identified a neutron-specific gene signature driven by cGAS-STING-NF-{kappa}B signaling (RELB, NFKB1, C3, MALAT1) and suppression of B-cell and myeloid identity genes (IGHD, TCL1A, CLEC7A, TLR2), defining a biologically coherent neutron quality index with distinct immunomodulatory effects. For the first time, we assessed neutron RBEs at the gene, pathway, and global transcriptomic levels in a human blood model, reporting a global transcriptomic neutron RBE of 1.30 (95% CI: 1.14-1.49) at 6 h and 1.21 (95% CI: 1.14-1.28) at 24 h, providing a valuable basis for biodosimetry in mixed-field exposure scenarios. Our findings advance the mechanistic understanding of neutron radiation responses and support the development of biodosimetric approaches for mixed-field exposure scenarios.

Disulfide-rich peptides possess exceptional potency and selectivity but are often limited by the instability and synthetic challenges associated with native disulfide bonds. Here, we report the design, synthesis, pharmacological evaluation, and structural characterisation of triazole-based peptidomimetics of the -GI conotoxin, a selective antagonist of the muscle-type nicotinic acetylcholine receptor (nAChR). A series of 1,4- and 1,5-disubstituted triazole analogues were prepared entirely on resin using CuAAC and RuAAC chemistry to replace the native Cys3/13 disulfide bridge. Functional evaluation against human muscle nAChRs revealed that 1,5-triazole analogues retained low-nanomolar potency, with the lead mimetic exhibiting activity comparable to native -GI. Cryo-electron microscopy of the lead compound bound to the muscle-type nAChR provided the first structure of a disulfide-isostere peptidomimetic in complex with a membrane receptor. The structure demonstrates that the 1,5-triazole reproduces the native peptide fold with high fidelity while contributing receptor-facing interactions not available to the native disulfide bridge. Molecular dynamics simulations further revealed conserved hydration networks and similar conformational sampling between the native peptide and lead mimetic. Together, these findings establish triazoles as effective disulfide surrogates and provide a structural framework for the rational design of stabilised conotoxin therapeutics.

Many biological processes rely on mechanical forces, with protein molecules acting as key mediators. Understanding how proteins respond to mechanical stress is essential for conditions including cardiomyopathy and muscular dystrophy. Natural proteins such as dystrophin and utrophin are composed of heterogeneous folding domains with distinct mechanical properties; deciphering domain-level behavior provides insights into disease mechanisms and informs therapeutic strategies. Single-molecule force spectroscopy (SMFS) enables probing the mechanical properties of entire proteins, yet current approaches struggle to identify heterogeneous folding domains, particularly without prior knowledge. Here, we present the first automated framework to identify heterogeneous folding domains in SMFS data, applying both existing clustering methods and a novel physics-aware deep clustering architecture, LatentUnfold. LatentUnfold learns complementary latent representations from force magnitude and the force-extension physical relationship through dual autoencoders, jointly optimized for clustering assignments. We apply our framework to experimental SMFS data collected from a synthetic two-domain protein (ddFLN4-Titin I27) as well as natural protein constructs of dystrophin and utrophin, with Monte Carlo simulated datasets serving as controlled validation. For the synthetic protein, we recover mechanical properties consistent with previously reported values for each domain. For the natural proteins, we uncover two mechanically distinct domain populations - corresponding to the N-terminal domain and spectrin-like repeats - with differences in both unfolding force and contour length increase, and reveal different unfolding order between them for the first time. This work enables domain-level biological inference, overcoming prior limitations that relied on averaging and overlooked heterogeneity, thus advancing the understanding of mechanical behavior in protein unfolding.

UCH-L1 is a monomeric deubiquitinating enzyme whose native structure embeds a shallow $5_2$ knot located near the N-terminus, placing the knotted topology in direct proximity to both the substrate-binding pocket and the catalytic site. While our previous work established that N-terminal integrity is critical for catalytic activity, the energetic cost of unknotting and its structural consequences remained unquantified. Here, we combine steered molecular dynamics with an umbrella sampling scheme to generate topologically modified variants of UCH-L1 and, for the first time, reconstruct the free-energy profile of UCH-L1 unknotting. The potential of mean force reveals a steep energetic barrier to knot disruption, consistent with knotting being a late, rate-limiting folding step that is effectively locked in once the native structure is established. Long unbiased MD simulations of fully unknotted variants in both apo and holo states show that knot removal increases local flexibility at the N-terminus without inducing significant global structural destabilization. Binding energy calculations indicate that the unknotted variant binds to ubiquitin less tightly than the wild-type ($\sim$-62~vs~$\sim$-76~kcal/mol), suggesting that topological integrity contributes to substrate affinity. Together, these results show that the $5_2$ knot in UCH-L1 is not a passive structural feature but a functional element that fine-tunes folding kinetics and contributes to substrate binding efficiency.

The navigational prowess of migratory birds is thought to arise from light-dependent radical-pair chemistry in cryptochrome 4 (CRY4), yet the slow structural transitions that couple photochemistry to signaling remain elusive. Here, we combine temperature-controlled steady-state UV-visible spectroscopy and hydrogen-deuterium exchange mass spectrometry (HDX-MS) to elucidate the photochemical and conformational dynamics of pigeon CRY4 (ClCRY4). Steady-state measurements at 5-25 {degrees}C reveal that lower temperatures slow FAD photoreduction and prolong the FAD neutral semiquinone signaling state. This occurs without a solvent kinetic isotope effect, implicating a conformational change rather than proton transfer as the rate determining step in FAD neutral semiquinone formation. Simultaneous HDX-MS under blue-light exposure identifies protection near the FAD-binding site and C-terminal region. To enhance sensitivity, we developed a pump-probe HDX-MS approach at 10 {degrees}C. This reveals eight peptides (within the phosphate-binding loop, protrusion motif, electron-transfer-chain loops and C-terminal tail) that exhibit rapid ([&le;]10 s) and sustained light-induced protection, delineating early conformational rearrangements as a prerequisite for FAD neutral semiquinone accumulation. The findings of slower onset HDX protection as well as a bimodal pattern of deuterium uptake in the phosphate-binding loop further identify a local redistribution of conformational substates on the time scale of the accumulation of the signaling species. Site specific mutagenesis within the CTT supports the findings, which lead to a model in which blue light triggers rapid clamping down of protein near the two regions of spin pair separation, followed by a rate limiting closure of a surface loop. The resolution of time-dependent structural transitions that follow photoactivation of CRY4 resolves the interface between quantum radical-pair formation and classical conformational changes, while providing an enhanced structural framework for the molecular events that underlie avian magnetoreception.

Life began with the emergence of a molecule that could replicate its own genetic material, a task plausibly mediated by an RNA-dependent RNA polymerase ribozyme. Here, we present the structure of such a polymerase ribozyme, bound to RNA substrates comprising the template, primer, and nucleoside triphosphate (NTP) analog. The structure reveals how directed evolution shaped flanking elements around a highly conserved catalytic core derived from the ancestral class I ligase ribozyme. Each element serves as a functional module, positioning the primer-template duplex and incoming NTP within the active site of the enzyme. This emergent domain organization is remarkably similar to the "right hand" configuration of polymerase proteins, suggesting a common functional form for copying nucleic acids, regardless of biopolymer catalyst.

Adhesion of cells is the key factor determining functioning of multicellular organisms. Viscoelastic properties of cells can be studied by multiple methods. However, attractiveness of cells or extracellular matrix without the elastic component (dispersive adhesion) is not accessible. We present an extension of force spectrometry technology: the Multivalent Adhesive Probe Atomic Force Microscopy (MAPA) that delivers dispersive adhesion maps of live cells and biosurfaces, and identifies differences unresolved by viscoelastic probing.

High-quality experimental datasets that link protein-ligand structures with binding affinity data are essential for developing and evaluating structure-based machine learning methods. To help address this need, we established OpenBind as an open-science initiative to generate large-scale experimental datasets for structure-based AI and molecular discovery. Here, we describe the first public OpenBind release, which, to the best of our knowledge, is the largest public single-target experimental structure-affinity dataset. The dataset focuses on enteroviral 2A protease, comprising 925 crystallographic binding events from 699 compounds and associated affinity measurements for 601 compounds. It combines structures from an initial fragment screen and follow-on molecules, together with affinity data, linking experimentally determined protein-ligand binding modes to biophysical measurements within a coherent antiviral discovery campaign. We used this dataset to evaluate protein-ligand structure prediction, binding-affinity prediction, and virtual screening using representative structure-based methods, including docking and cofolding. This exposed several challenges that are central to practical structure-based modelling: docking performance depends strongly on binding-pocket conformation, poses are difficult to rank, and structure-based affinity prediction remains challenging. Fine-tuning OpenFold3-p2 on the fragment-screen structures substantially improved pose prediction and virtual screening for related follow-on compounds, demonstrating how early-stage experimental structures can support target-specific model adaptation.

Neutron scattering is a powerful technique to investigate atomic structures and molecular dynamics of proteins at the nano-scale. When it comes to dynamics, incoherent and coherent scattering respectively provide information on the single and collective dynamics of nuclei. In proteins, hydrogen has the highest incoherent cross-section, and it is common practice to overlook the contribution of coherent terms stemming from all nuclei. However, the fast collective dynamics of heavier nuclei could also be studied if coherent scattering and incoherent scattering were experimentally separated. The recent advent of polarized neutron spectroscopy with sufficient flux and energy resolution has made it possible, and opens new perspectives to investigate the relative importance of coherent scattering and the information it provides on biological samples. The present study reports on the use of polarized quasi-elastic neutron scattering (QENS) and the application of a minimalistic model adapted to both individual and collective dynamics. Using a perdeuterated green fluorescent protein as a model globular protein, the study provides an interpretation of the dynamical parameters obtained with QENS, and a comparative study of the Elastic Coherent and Incoherent Scattering Factor. Based on both experiments and calculations, we discuss the relative importance of distinct and self components of coherent scattering, which is often wrongly assumed to be representative of collective dynamics only. The results highlight the current impediments rendering complicated a straightforward analysis of fast collective dynamics in hydrated protein samples.

Linker chemistry and conformation are central determinants of PROTAC activity, shaping ternary-complex geometry, cooperativity, target-lysine presentation and cellular permeability. Existing linker generators often lack explicit control over linker flexibility, require predefined attachment sites and linker lengths, or produce structures that demand substantial geometric correction, limiting their utility in practical PROTAC design. Here we introduce FlexiTAC, a Bayesian flow network that jointly generates linker atom types and coordinates from the warhead and E3-ligase-ligand contexts. We also assemble PROTAC-3D, a quality-controlled collection of 63,554 component-resolved PROTAC structures for model training, and PROTAC-Bench, which covers molecular quality, fragment preservation, geometric fidelity, conformational stability, fragment awareness, rediscovery and sampling efficiency. Compared to the best 3D baseline models, FlexiTAC improves validity by 12.0-12.7% and achieves the highest PoseBusters pass rate of 79.5%-80.0%. A differentiable guidance module shifted generated linkers along a conformational ensemble-derived rigidity axis without retraining the generator. In silico case studies further show that the model can accept crystal-derived, redocked or predicted structural inputs. Together, FlexiTAC, PROTAC-3D and PROTAC-Bench establish an integrated and reproducible framework for data-driven PROTAC linker design, combining controllable structure-conditioned generation with standardized training data and evaluation protocols. This framework expands the linker chemical and conformational space accessible to computational exploration, provides a foundation for future method development and enables the systematic generation of structure-conditioned linker designs with tunable conformational flexibility.

Pseudomonas aeruginosa is an opportunistic human pathogen capable of infecting a wide range of tissues and organs. Its persistence during chronic infection is strongly associated with biofilm formation, which depends on extracellular polysaccharides such as alginate. The HtrA-like periplasmic serine protease MucD is a key regulator of bacterial virulence, stress response, and alginate production, yet its molecular mechanism has remained largely unclear. Here, we discovered the alginate acetylation and export proteins AlgX and AlgK as MucD substrates, and characterized their degradation by mass spectrometry and bioinformatic analysis. We further determined the cryo-EM structure of MucD bound to an AlgK-derived substrate peptide, offering atomic insights into MucD oligomerization assembly, substrate recognition, and specificity. Together with structure-guided mutagenesis and biochemical assays, our results revealed that MucD proteolytic activity is governed by an equilibrium between a resting 12-mer and an active trimer. Crucially, we demonstrate that MucD represses alginate biosynthesis post-translationally, in addition to its previously implicated role in transcriptional regulation. These findings define a distinct activation mechanism and regulatory function for MucD and provide new insight into bacterial HtrA-like serine proteases.

Infectious bronchitis virus (IBV) is a member of the Gammacoronavirus genus responsible for respiratory illness and weakened eggshells in infected chickens, adversely impacting the poultry industry. Escaping innate immune detection during infection is crucial for coronavirus proliferation in the host. The production of double-stranded RNA during coronavirus replication triggers innate immune sensors to create an antiviral state within infected cells. To counter this response, coronaviruses employ nonstructural protein 15 (nsp15) endoribonuclease to degrade double-stranded RNA. Here, we use cryo-electron microscopy and biochemistry to characterize IBV nsp15 interactions with RNA. While the overall structure and active site of IBV nsp15 strongly resemble previous studies of nsp15 from other coronaviral genera, we note that double-stranded RNA contacts several non-conserved residues peripheral to the enzyme active site. Our data show that these residue positions can have strong impacts on RNA cleavage suggesting unique solutions for RNA engagement across coronavirus species. We also demonstrate a preference for IBV nsp15 to cleave double-stranded RNA over single-stranded RNA and observe nsp15 hexamers with two double-stranded RNAs bound simultaneously. This study reinforces the need to study diverse coronavirus species to identify distinct viral enzyme characteristics.

We present the solution structures of HIV-1 proteins NC(p7)1-55 corresponding to the full-length NC(p7) and mature p6. The studies were carried in water and, to mimic the membrane, in micellar DPC (Dodecylphosphocholine) conditions. Our results unravel for the first time the structure adopted by the N-terminal amino acids of the free NC(p7)1-55, with the formation of a small helix spanning residues F6 to R10. Our NMR and Fluorescence Anisotropy data disclose an interaction between NC(p7)1-55 and p6 both in water and DPC, with respective Kd of 2.5mM and 370 mM at 23{degrees}C. The interaction is thus strengthened in lipidic conditions. Protein p6 stabilizes the N-terminus of NC(p7)1-55 while increasing at the same time the dynamic of the first zinc finger. Although the entire p6 sequence is involved in the interaction, we show that its C-terminal region is particularly sensitive to the presence of NC(p7)1-55, with a propensity of forming a a helix ranging from amino acids S111 to F116. This study brings experimental evidence of a direct protein-protein interaction between p6 and the N-terminal region of NC(p7)1-55. We further show that such interaction is readily accommodated within the NC(p15) framework and hypothesize that it may facilitate the selective assembly of assembly of the viral genomic RNA (gRNA) in the cell.

The cluster of differentiation 36 (CD36) is a membrane protein with broad physiological roles in health and disease, and its function is regulated in part by phosphorylation. Experimental evidence shows that phosphorylation of Thr92 reduces CD36 affinity for thrombospondin-1 (TSP-1), binding of which initiates antiangiogenic signaling, whereas phosphorylation of Ser237 decreases CD36-mediated fatty acid uptake, with implications for energy metabolism. However, the only available crystal structure of CD36 lacks phosphorylation, and the molecular mechanisms by which phosphorylation regulates CD36 function remain largely unknown. This study provides an atomically detailed computational characterization of CD36 in unphosphorylated and dual phosphorylated states, using molecular dynamics simulations with a total sampling time of 30 microseconds in combination with Markov state models. We present, to our knowledge, the first evidence of a cryptic pocket on CD36 surface that is formed by phosphorylation. This cryptic surface pocket and a loop spanning residues 121-131 form a high affinity binding site for TSP-1 derived ligands, shifting their binding away from the canonical site. We propose that this altered binding provides a molecular basis for the disruption of antiangiogenic signaling upon CD36 phosphorylation. Additionally, our data indicate that, phosphorylation increases helicity and compaction within the helix-loop region spanning residues 296-331, narrowing one of the entrances to the internal cavity and reducing its overall volume. These conformational changes provide a potential mechanistic explanation for the decrease in fatty acid uptake upon CD36 phosphorylation. Our findings provide structural insights that may inform the future design of CD36 modulators and emphasize the importance of targeting phosphorylation induced CD36 conformations in angiogenic and metabolic diseases.