Craniofacial anomalies in a murine model of heterozygous fibroblast growth factor 10 gene mutation.

OBJECTIVE: Dysregulation of Fibroblast Growth Factor 10 (FGF10), a member of the family of Fibroblast Growth Factor (FGF) proteins, has been implicated in craniofacial and dental anomalies, including craniosynostosis, cleft palate, and Lacrimo-Auriculo-Dento-Digital Syndrome. The aim of this murine study was to assess the craniofacial and dental phenotypes associated with a heterozygous FGF10 gene (FGF10+/- ) mutation at skeletal maturity. METHODS: Skulls of 40 skeletally mature mice, comprising two genotypes (heterozygous FGF10+/- mutation, n = 22; wildtype, n = 18) and two sexes (male, n = 23; female, n = 17), were subjected to micro-computed tomography. Landmark-based linear dimensions were measured for the cranial vault, maxilla, mandible, and first molar teeth. Multivariate analysis of variance was performed to assess whether there were significant differences in the craniofacial and dental structures between genotypes and sexes. RESULTS: The craniomaxillary skeleton and the first molar teeth were smaller in the FGF10+/- mice (P < .05), but the mandible was unaffected. Sex did not have a significant effect on these structures (P > .05). Cranial sutural defects were noted in 5/22 (22.7%) mutant versus 2/18 (11.1%) wildtype mice, and cleft palate in only one (4.5%) mutant mouse. None of the mice displayed craniosynostosis, expansive bony lesions, bifid condyles, or impacted teeth. CONCLUSION: The FGF10+/- mutation was associated with craniomaxillary skeletal hypoplasia that probably arose from deficient (delayed) intramembranous ossification of the sutured bones. Overall, the skeletal and dental data suggest that the FGF10 gene plays an important role in the aetiology of craniofacial dysmorphology and malocclusion.

The Role of Fibroblast Growth Factor 10 Signaling in Duodenal Atresia.

INTRODUCTION: Duodenal atresia (DA) is a congenital bowel obstruction requiring major surgery in the first week of life. Three morphological phenotypes are described, reflecting increasing degrees of obstruction and discontinuity of the duodenum. The cause of DA is not known. Tandler's original "solid cord" hypothesis conflicts with recent biological evidence, and is unable to account for differing DA types. In humans, a genetic etiology is supported by the association between Trisomy 21 and DA, and reports of familial inheritance patterns. Interruption of FGF10/FGFR2b signaling is the best demonstrated genetic link to DA in mice, with 35-75% of homozygous knockout embryos developing DA. PURPOSE: This review examines the current evidence surrounding the etiology of DA. We focus on research regarding FGF10/FGFR2b signaling and its role in duodenal and other intestinal atresia. Further, we outline planned future research in this area, that we consider necessary to validate and better understand this murine model in order to successfully translate this research into clinical practice. CONCLUSION: Determining the etiology of DA in humans is a clinical and scientific imperative. Fgf10/Fgfr2b murine models represent current science's best key to unlocking this mystery. However, further research is required to understand the complex role of FGF10/FGFR2b signaling in DA development. Such complexity is expected, given the lethality of their associated defects makes ubiquitous interruption of either Fgf10 or Fgfr2b genes an unlikely cause of DA in humans. Rather, local or tissue-specific mutation in Fgf10, Fgfr2b, or their downstream targets, is the hypothesized basis of DA etiology.

Clinical data and Pediatric Quality of Life Inventory (PedsQL™) scores for children with duodenal atresia.

This article presents raw data obtained from a prospectively collected database of children with duodenal atresia at tertiary pediatric surgery hospital. For all potential participants, pertinent demographic, clinical and operative data was obtained from the database. Potential participants were then contacted and invited to complete a Pediatric Quality of Life Inventory (PedsQL™) 4.0 core score and gastrointestinal module questionnaires. Participant's response to each item in the questionnaires is provided, as well as their calculated health related quality of life scores. Data has the potential to be reused in future studies examining quality of life in duodenal atresia, paediatric gastrointestinal conditions, surgical neonatal conditions and children with trisomy 21. Further analysis and discussion is contained in related research article titled "Quality of life outcomes in children born with duodenal atresia" [1].

Quality of life outcomes in children born with duodenal atresia.

PURPOSE: The aim of this study was to determine long term quality of life (QoL) outcome for children who underwent surgery for duodenal atresia (DA). METHODS: Patients were identified from a prospective database of neonatal DA cases managed at a tertiary pediatric surgical centre. The QoL was measured using the validated PedsQL™ 4.0 core score and PedsQL™ gastrointestinal module; higher score equates to better QoL. Participants' scores were compared to published control cohorts, age-matching the core score. Trisomy 21 was identified a priori as a possible confounder, informing subgroup analyses for children with and without trisomy 21. RESULTS: Fifty-five families were invited to participate, with 38 surveys returned (39% male; median age 6.7y, range 2.7-17.3y). Seven participants had trisomy 21. There were no differences in QoL measures between all DA participants and controls. The PedsQL™ core score was significantly lower for DA participants with trisomy 21, but there was no accompanying difference in PedsQL™ gastrointestinal score. CONCLUSIONS: Children undergoing DA surgery in the neonatal period typically grow up to have a QoL comparable to a healthy population. Children with DA and trisomy 21 were more likely to have reduced overall QoL, albeit without an associated difference in gastrointestinal QoL score. LEVEL OF EVIDENCE: Prognosis study - level II (prospective cohort study).

FGF10 and the Mystery of Duodenal Atresia in Humans.

Background: Duodenal atresia (DA) is a congenital obstruction of the duodenum, which affects 1 in 7000 pregnancies and requires major surgery in the 1st days of life. Three morphological DA types are described. In humans, the association between DA and Down syndrome suggests an underlying, albeit elusive, genetic etiology. In mice, interruption of fibroblast growth factor 10 (Fgf10) gene signaling results in DA in 30-50% of embryos, supporting a genetic etiology. This study aims to validate the spectrum of DA in two novel strains of Fgf10 knock-out mice, in preparation for future and translational research. Methods: Two novel CRISPR Fgf10 knock-out mouse strains were derived and embryos generated by heterozygous plug-mating. E15.5-E19.5 embryos were genotyped with respect to Fgf10 and micro-dissected to determine the presence and type of DA. Results: One twenty seven embryos (32 wild-type, 34 heterozygous, 61 null) were analyzed. No wild-type or heterozygous embryos had DA. However, 74% of Fgf10 null embryos had DA (49% type 1, 18% type 2, and 33% type 3). Conclusion: Our CRISPR-derived strains showed higher penetrance of DA due to single-gene deletion of Fgf10 in mice than previously reported. Further, the DA type distribution in these mice more closely reiterated that observed in humans. Future experiments will document RNA and protein expression of FGF10 and its key downstream signaling targets in normal and atretic duodenum. This includes exploitation of modern, high-fidelity developmental tools, e.g., Fgf10 flox/+-tomatoflox/flox mice.

Activation of the Notch Signaling Pathway In Vivo Elicits Changes in CSL Nuclear Dynamics.

A key feature of Notch signaling is that it directs immediate changes in transcription via the DNA-binding factor CSL, switching it from repression to activation. How Notch generates both a sensitive and accurate response-in the absence of any amplification step-remains to be elucidated. To address this question, we developed real-time analysis of CSL dynamics including single-molecule tracking in vivo. In Notch-OFF nuclei, a small proportion of CSL molecules transiently binds DNA, while in Notch-ON conditions CSL recruitment increases dramatically at target loci, where complexes have longer dwell times conferred by the Notch co-activator Mastermind. Surprisingly, recruitment of CSL-related corepressors also increases in Notch-ON conditions, revealing that Notch induces cooperative or "assisted" loading by promoting local increase in chromatin accessibility. Thus, in vivo Notch activity triggers changes in CSL dwell times and chromatin accessibility, which we propose confer sensitivity to small input changes and facilitate timely shut-down.

Optimization and Validation of Efficient Models for Predicting Polythiophene Self-Assembly.

We develop an optimized force-field for poly(3-hexylthiophene) (P3HT) and demonstrate its utility for predicting thermodynamic self-assembly. In particular, we consider short oligomer chains, model electrostatics and solvent implicitly, and coarsely model solvent evaporation. We quantify the performance of our model to determine what the optimal system sizes are for exploring self-assembly at combinations of state variables. We perform molecular dynamics simulations to predict the self-assembly of P3HT at ∼350 combinations of temperature and solvent quality. Our structural calculations predict that the highest degrees of order are obtained with good solvents just below the melting temperature. We find our model produces the most accurate structural predictions to date, as measured by agreement with grazing incident X-ray scattering experiments.

Tying Together Multiscale Calculations for Charge Transport in P3HT: Structural Descriptors, Morphology, and Tie-Chains.

Evaluating new, promising organic molecules to make next-generation organic optoelectronic devices necessitates the evaluation of charge carrier transport performance through the semi-conducting medium. In this work, we utilize quantum chemical calculations (QCC) and kinetic Monte Carlo (KMC) simulations to predict the zero-field hole mobilities of ∼100 morphologies of the benchmark polymer poly(3-hexylthiophene), with varying simulation volume, structural order, and chain-length polydispersity. Morphologies with monodisperse chains were generated previously using an optimized molecular dynamics force-field and represent a spectrum of nanostructured order. We discover that a combined consideration of backbone clustering and system-wide disorder arising from side-chain conformations are correlated with hole mobility. Furthermore, we show that strongly interconnected thiophene backbones are required for efficient charge transport. This definitively shows the role "tie-chains" play in enabling mobile charges in P3HT. By marrying QCC and KMC over multiple length- and time-scales, we demonstrate that it is now possible to routinely probe the relationship between molecular nanostructure and device performance.

Conformational behavior and stacking interactions of contorted polycyclic aromatics.

We present a systematic computational analysis of the conformations and stacking interactions of a set of 18 saddle-shaped, contorted polycyclic aromatic compounds at the B97-D3M(BJ)/TZV(2d,2p)//B97-D/TZV(2d,2p) level of theory. These doubly-concave systems offer a means of tuning the strength of stacking interactions through variations in molecular curvature, and understanding the intermolecular non-covalent interactions exhibited by these systems will aid the design of contorted polycyclic systems with precisely defined packing in the solid state. Computations reveal wide variations in both the nature of the low-lying conformations and the stacking affinities of these systems. In particular, the introduction of both thiophene rings around the periphery of these systems and the incorporation of B and N atoms into the coronene core can greatly enhance their tendency to form strongly stacked dimers. Overall, these data provide a reminder that curvature does not always lead to stronger stacking interactions.

Enhanced Computational Sampling of Perylene and Perylothiophene Packing with Rigid-Body Models.

Molecular simulations have the potential to advance the understanding of how the structure of organic materials can be engineered through the choice of chemical components but are limited by computational costs. The computational costs can be significantly lowered through the use of modeling approximations that capture the relevant features of a system, while lowering algorithmic complexity or by decreasing the degrees of freedom that must be integrated. Such methods include coarse-graining techniques, approximating long-range electrostatics with short-range potentials, and the use of rigid bodies to replace flexible bonded constraints between atoms. To understand whether and to what degree these techniques can be leveraged to enhance the understanding of planar organic molecules, we investigate the morphologies predicted by molecular dynamic simulations using simplified molecular models of perylene and perylothiophene. Approximately, 10 000 wall-clock hours of graphics processing unit-accelerated simulations are performed using both rigid and flexible models to test their efficiency and predictive capability with the two chemistries. We characterize the 1191 resulting morphologies using simulated X-ray diffraction and cluster analysis to distinguish structural transitions, summarized by four phase diagrams. We find that the morphologies generated by the rigid model of perylene and perylothiophene match with those generated by the flexible model. We find that ordered, hexagonally packed columnar phases are thermodynamically favored over a wide range of densities and temperatures for both molecules, in qualitative agreement with experiments. Furthermore, we find the rigid model to be more computationally efficient for both molecules, providing more samples per second and shorter times to equilibrium. Owing to the structural accuracy and improved computational efficiency of modeling polyaromatic groups as rigid bodies, we recommend this modeling choice for enhancing the sampling in polyaromatic molecular simulations.

Study of Plasmodium falciparum DHHC palmitoyl transferases identifies a role for PfDHHC9 in gametocytogenesis.

Palmitoylation is the post-translational reversible addition of the acyl moiety, palmitate, to cysteine residues of proteins and is involved in regulating protein trafficking, localization, stability and function. The Aspartate-Histidine-Histidine-Cysteine (DHHC) protein family, named for their highly conserved DHHC signature motif, is thought to be responsible for catalysing protein palmitoylation. Palmitoylation is widespread in all eukaryotes, including the malaria parasite, Plasmodium falciparum, where over 400 palmitoylated proteins are present in the asexual intraerythrocytic schizont stage parasites, including proteins involved in key aspects of parasite maturation and development. The P. falciparum genome includes 12 proteins containing the conserved DHHC motif. In this study, we adapted a palmitoyl-transferase activity assay for use with P. falciparum proteins and demonstrated for the first time that P. falciparum DHHC proteins are responsible for the palmitoylation of P. falciparum substrates. This assay also reveals that multiple DHHCs are capable of palmitoylating the same substrate, indicating functional redundancy at least in vitro. To test whether functional redundancy also exists in vivo, we investigated the endogenous localization and essentiality of a subset of schizont-expressed PfDHHC proteins. Individual PfDHHC proteins localized to distinct organelles, including parasite-specific organelles such as the rhoptries and inner membrane complex. Knock-out studies identified individual DHHCs that may be essential for blood-stage growth and others that were functionally redundant in the blood stages but may have functions in other stages of parasite development. Supporting this hypothesis, disruption of PfDHHC9 had no effect on blood-stage growth but reduced the formation of gametocytes, suggesting that this protein could be exploited as a transmission-blocking target. The localization and stage-specific expression of the DHHC proteins may be important for regulating their substrate specificity and thus may provide a path for inhibitor development.

Protease-Activated Receptor 4 Variant p.Tyr157Cys Reduces Platelet Functional Responses and Alters Receptor Trafficking.

OBJECTIVE: Protease-activated receptor 4 (PAR4) is a key regulator of platelet reactivity and is encoded by F2RL3, which has abundant rare missense variants. We aimed to provide proof of principle that rare F2LR3 variants potentially affect platelet reactivity and responsiveness to PAR1 antagonist drugs and to explore underlying molecular mechanisms. APPROACH AND RESULTS: We identified 6 rare F2RL3 missense variants in 236 cardiac patients, of which the variant causing a tyrosine 157 to cysteine substitution (Y157C) was predicted computationally to have the greatest effect on PAR4 structure. Y157C platelets from 3 cases showed reduced responses to PAR4-activating peptide and to α-thrombin compared with controls, but no reduction in responses to PAR1-activating peptide. Pretreatment with the PAR1 antagonist vorapaxar caused lower residual α-thrombin responses in Y157C platelets than in controls, indicating greater platelet inhibition. HEK293 cells transfected with a PAR4 Y157C expression construct had reduced PAR4 functional responses, unchanged total PAR4 expression but reduced surface expression. PAR4 Y157C was partially retained in the endoplasmic reticulum and displayed an expression pattern consistent with defective N-glycosylation. Mutagenesis of Y322, which is the putative hydrogen bond partner of Y157, also reduced PAR4 surface expression in HEK293 cells. CONCLUSIONS: Reduced PAR4 responses associated with Y157C result from aberrant anterograde surface receptor trafficking, in part, because of disrupted intramolecular hydrogen bonding. Characterization of PAR4 Y157C establishes that rare F2RL3 variants have the potential to markedly alter platelet PAR4 reactivity particularly after exposure to therapeutic PAR1 antagonists.

A versatile strategy for rapid conditional genome engineering using loxP sites in a small synthetic intron in Plasmodium falciparum.

Conditional genome engineering in the human malaria pathogen Plasmodium falciparum remains highly challenging. Here we describe a strategy for facile and rapid functional analysis of genes using an approach based on the Cre/lox system and tailored for organisms with short and few introns. Our method allows the conditional, site-specific removal of genomic sequences of essential and non-essential genes by placing loxP sites into a short synthetic intron to produce a module (loxPint) can be placed anywhere in open reading frames without compromising protein expression. When duplicated, the loxPint module serves as an intragenic recombineering point that can be used for the fusion of gene elements to reporters or the conditional introduction of point mutations. We demonstrate the robustness and versatility of the system by targeting the P. falciparum merozoite surface protein 1 gene (msp1), which has previously proven refractory to genetic interrogation, and the parasite exported kinase FIKK10.1.

Processing of Plasmodium falciparum Merozoite Surface Protein MSP1 Activates a Spectrin-Binding Function Enabling Parasite Egress from RBCs.

The malaria parasite Plasmodium falciparum replicates within erythrocytes, producing progeny merozoites that are released from infected cells via a poorly understood process called egress. The most abundant merozoite surface protein, MSP1, is synthesized as a large precursor that undergoes proteolytic maturation by the parasite protease SUB1 just prior to egress. The function of MSP1 and its processing are unknown. Here we show that SUB1-mediated processing of MSP1 is important for parasite viability. Processing modifies the secondary structure of MSP1 and activates its capacity to bind spectrin, a molecular scaffold protein that is the major component of the host erythrocyte cytoskeleton. Parasites expressing an inefficiently processed MSP1 mutant show delayed egress, and merozoites lacking surface-bound MSP1 display a severe egress defect. Our results indicate that interactions between SUB1-processed merozoite surface MSP1 and the spectrin network of the erythrocyte cytoskeleton facilitate host erythrocyte rupture to enable parasite egress.

Diversity and impact of rare variants in genes encoding the platelet G protein-coupled receptors.

Platelet responses to activating agonists are influenced by common population variants within or near G protein-coupled receptor (GPCR) genes that affect receptor activity. However, the impact of rare GPCR gene variants is unknown. We describe the rare single nucleotide variants (SNVs) in the coding and splice regions of 18 GPCR genes in 7,595 exomes from the 1,000-genomes and Exome Sequencing Project databases and in 31 cases with inherited platelet function disorders (IPFDs). In the population databases, the GPCR gene target regions contained 740 SNVs (318 synonymous, 410 missense, 7 stop gain and 6 splice region) of which 70 % had global minor allele frequency (MAF) < 0.05 %. Functional annotation using six computational algorithms, experimental evidence and structural data identified 156/740 (21 %) SNVs as potentially damaging to GPCR function, most commonly in regions encoding the transmembrane and C-terminal intracellular receptor domains. In 31 index cases with IPFDs (Gi-pathway defect n=15; secretion defect n=11; thromboxane pathway defect n=3 and complex defect n=2) there were 256 SNVs in the target regions of 15 stimulatory platelet GPCRs (34 unique; 12 with MAF< 1 % and 22 with MAF≥ 1 %). These included rare variants predicting R122H, P258T and V207A substitutions in the P2Y12 receptor that were annotated as potentially damaging, but only partially explained the platelet function defects in each case. Our data highlight that potentially damaging variants in platelet GPCR genes have low individual frequencies, but are collectively abundant in the population. Potentially damaging variants are also present in pedigrees with IPFDs and may contribute to complex laboratory phenotypes.

Are hot charge transfer states the primary cause of efficient free-charge generation in polymer:fullerene organic photovoltaic devices? A kinetic Monte Carlo study.

Kinetic Monte Carlo simulations are used to examine the effect of high-energy, 'hot' delocalised charge transfer (HCT) states for donor:acceptor and mixed:aggregate blends, the latter relating to polymer:fullerene photovoltaic devices. Increased fullerene aggregation is shown to enhance charge generation and short-circuit device current - largely due to the increased production of HCT states at the aggregate interface. However, the instances where HCT states are predicted to give internal quantum efficiencies in the region of 50% do not correspond to HCT delocalisation or electron mobility measured in experiments. These data therefore suggest that HCT states are not the primary cause of high quantum efficiencies in some polymer:fullerene OPVs. Instead it is argued that HCT states are responsible for the fast charge generation seen in spectroscopy, but that regional variation in energy levels are the cause of long-term, efficient free-charge generation.

A novel thromboxane A2 receptor N42S variant results in reduced surface expression and platelet dysfunction.

A small number of thromboxane receptor variants have been described in patients with a bleeding history that result in platelet dysfunction. We have identified a patient with a history of significant bleeding, who expresses a novel heterozygous thromboxane receptor variant that predicts an asparagine to serine substitution (N42S). This asparagine is conserved across all class A GPCRs, suggesting a vital role for receptor structure and function.We investigated the functional consequences of the TP receptor heterozygous N42S substitution by performing platelet function studies on platelet-rich plasma taken from the patient and healthy controls. We investigated the N42S mutation by expressing the wild-type (WT) and mutant receptor in human embryonic kidney (HEK) cells. Aggregation studies showed an ablation of arachidonic acid responses in the patient, whilst there was right-ward shift of the U46619 concentration response curve (CRC). Thromboxane generation was unaffected. Calcium mobilisation studies in cells lines showed a rightward shift of the U46619 CRC in N42S-expressing cells compared to WT. Radioligand binding studies revealed a reduction in BMax in platelets taken from the patient and in N42S-expressing cells, whilst cell studies confirmed poor surface expression. We have identified a novel thromboxane receptor variant, N42S, which results in platelet dysfunction due to reduced surface expression. It is associated with a significant bleeding history in the patient in whom it was identified. This is the first description of a naturally occurring variant that results in the substitution of this highly conserved residue and confirms the importance of this residue for correct GPCR function.

Characterization of multiple platelet activation pathways in patients with bleeding as a high-throughput screening option: use of 96-well Optimul assay.

Up to 1% of the population have mild bleeding disorders, but these remain poorly characterized, particularly with regard to the roles of platelets. We have compared the usefulness of Optimul, a 96-well plate-based assay of 7 distinct pathways of platelet activation to characterize inherited platelet defects in comparison with light transmission aggregometry (LTA). Using Optimul and LTA, concentration-response curves were generated for arachidonic acid, ADP, collagen, epinephrine, Thrombin receptor activating-peptide, U46619, and ristocetin in samples from (1) healthy volunteers (n = 50), (2) healthy volunteers treated with antiplatelet agents in vitro (n = 10), and (3) patients with bleeding of unknown origin (n = 65). The assays gave concordant results in 82% of cases (κ = 0.62, P < .0001). Normal platelet function results were particularly predictive (sensitivity, 94%; negative predictive value, 91%), whereas a positive result was not always substantiated by LTA (specificity, 67%; positive predictive value, 77%). The Optimul assay was significantly more sensitive at characterizing defects in the thromboxane pathway, which presented with normal responses with LTA. The Optimul assay is sensitive to mild platelet defects, could be used as a rapid screening assay in patients presenting with bleeding symptoms, and detects changes in platelet function more readily than LTA. This trial was registered at www.isrctn.org as #ISRCTN 77951167.

Modulation of global low-frequency motions underlies allosteric regulation: demonstration in CRP/FNR family transcription factors.

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distinct site. There is growing evidence that allosteric cooperativity can be communicated by modulation of protein dynamics without conformational change. The mechanisms, however, for communicating dynamic fluctuations between sites are debated. We provide a foundational theory for how allostery can occur as a function of low-frequency dynamics without a change in structure. We have generated coarse-grained models that describe the protein backbone motions of the CRP/FNR family transcription factors, CAP of Escherichia coli and GlxR of Corynebacterium glutamicum. The latter we demonstrate as a new exemplar for allostery without conformation change. We observe that binding the first molecule of cAMP ligand is correlated with modulation of the global normal modes and negative cooperativity for binding the second cAMP ligand without a change in mean structure. The theory makes key experimental predictions that are tested through an analysis of variant proteins by structural biology and isothermal calorimetry. Quantifying allostery as a free energy landscape revealed a protein "design space" that identified the inter- and intramolecular regulatory parameters that frame CRP/FNR family allostery. Furthermore, through analyzing CAP variants from diverse species, we demonstrate an evolutionary selection pressure to conserve residues crucial for allosteric control. This finding provides a link between the position of CRP/FNR transcription factors within the allosteric free energy landscapes and evolutionary selection pressures. Our study therefore reveals significant features of the mechanistic basis for allostery. Changes in low-frequency dynamics correlate with allosteric effects on ligand binding without the requirement for a defined spatial pathway. In addition to evolving suitable three-dimensional structures, CRP/FNR family transcription factors have been selected to occupy a dynamic space that fine-tunes biological activity and thus establishes the means to engineer allosteric mechanisms driven by low-frequency dynamics.

A member of the Plasmodium falciparum PHIST family binds to the erythrocyte cytoskeleton component band 4.1.

BACKGROUND: Plasmodium falciparum parasites export more than 400 proteins into the cytosol of their host erythrocytes. These exported proteins catalyse the formation of knobs on the erythrocyte plasma membrane and an overall increase in erythrocyte rigidity, presumably by modulating the endogenous erythrocyte cytoskeleton. In uninfected erythrocytes, Band 4.1 (4.1R) plays a key role in regulating erythrocyte shape by interacting with multiple proteins through the three lobes of its cloverleaf-shaped N-terminal domain. In P. falciparum-infected erythrocytes, the C-lobe of 4.1R interacts with the P. falciparum protein mature parasite-infected erythrocyte surface antigen (MESA), but it is not currently known whether other P. falciparum proteins bind to other lobes of the 4.1R N-terminal domain. METHODS: In order to identify novel 4.1R interacting proteins, a yeast two-hybrid screen was performed with a fragment of 4.1R containing both the N- and α-lobes. Positive interactions were confirmed and investigated using site-directed mutagenesis, and antibodies were raised against the interacting partner to characterise it's expression and distribution in P. falciparum infected erythrocytes. RESULTS: Yeast two-hybrid screening identified a positive interaction between the 4.1R N- and α-lobes and PF3D7_0402000. PF3D7_0402000 is a member of a large family of exported proteins that share a domain of unknown function, the PHIST domain. Domain mapping and site-directed mutagenesis established that it is the PHIST domain of PF3D7_0402000 that interacts with 4.1R. Native PF3D7_0402000 is localized at the parasitophorous vacuole membrane (PVM), and colocalizes with a subpopulation of 4.1R. DISCUSSION: The function of the majority of P. falciparum exported proteins, including most members of the PHIST family, is unknown, and in only a handful of cases has a direct interaction between P. falciparum-exported proteins and components of the erythrocyte cytoskeleton been established. The interaction between 4.1R and PF3D7_0402000, and localization of PF3D7_0402000 with a sub-population of 4.1R at the PVM could indicate a role in modulating PVM structure. Further investigation into the mechanisms for 4.1R recruitment is needed. CONCLUSION: PF3D7_0402000 was identified as a new binding partner for the major erythrocyte cytoskeletal protein, 4.1R. This interaction is consistent with a growing body of literature that suggests the PHIST family members function by interacting directly with erythrocyte proteins.