The carboxyl terminus of Neph family members binds to the PDZ domain protein zonula occludens-1

The PSD95/Dlg/ZO-1 (PDZ) domain-containing protein zonula occludens-1 (ZO-1) selectively localizes to the cytoplasmic basis of the slit diaphragm, a specialized cell-cell contact in between glomerular podocytes necessary to prevent the loss of protein in the urine. However, the function of ZO-1 at the slit diaphragm has remained elusive. Deletion of Neph1, a slit diaphragm protein of the immunoglobulin superfamily with a cytoplasmic PDZ binding site, causes proteinuria in mice. We demonstrate now that Neph1 binds ZO-1. This interaction was mediated by the first PDZ domain of ZO-1 and involved the conserved PDZ domain binding motif present in the carboxyl terminus of the three known Neph family members. Furthermore, Neph1 co-immunoprecipitates with ZO-1 from lysates of mouse kidneys, demonstrating that this interaction occurs in vivo. Both deletion of the PDZ binding motif of Neph1 as well as threonine-to-glutamate mutation of the threonine within the binding motif abrogated binding of ZO-1, suggesting that phosphorylation may regulate this interaction. ZO-1 binding was associated with a strong increase in tyrosine phosphorylation of the cytoplasmic tail of Neph1 and dramatically accelerated the ability of Neph1 to induce signal transduction. Thus, our data suggest that ZO-1 may organize Neph proteins and recruit signal transduction components to the slit diaphragm of podocytes.     

Comparative structural analysis of the Erbin PDZ domain and the first PDZ domain of ZO-1. Insights into determinants of PDZ domain specificity

We report a structural comparison of the first PDZ domain of ZO-1 (ZO1-PDZ1) and the PDZ domain of Erbin (Erbin-PDZ). Although the binding profile of Erbin-PDZ is extremely specific ([D/E][T/S]WV(COOH)), that of ZO1-PDZ1 is similar ([R/K/S/T][T/S][W/Y][V/I/L](COOH)) but broadened by increased promiscuity for three of the last four ligand residues. Consequently, the biological function of ZO-1 is also broadened, as it interacts with both tight and adherens junction proteins, whereas Erbin is restricted to adherens junctions. Structural analyses reveal that the differences in specificity can be accounted for by two key differences in primary sequence. A reduction in the size of the hydrophobic residue at the base of the site(0) pocket enables ZO1-PDZ1 to accommodate larger C-terminal residues. A single additional difference alters the specificity of both site(-1) and site(-3). In ZO1-PDZ1, an Asp residue makes favorable interactions with both Tyr(-1) and Lys/Arg(-3). In contrast, Erbin-PDZ contains an Arg at the equivalent position, and this side chain cannot accommodate either Tyr(-1) or Lys/Arg(-3) but, instead, interacts favorably with Glu/Asp(-3). We propose a model for ligand recognition that accounts for interactions extending across the entire binding site but that highlights several key specificity switches within the PDZ domain fold.     

Ischemic injury to kidney induces glomerular podocyte effacement and dissociation of slit diaphragm proteins Neph1 and ZO-1

Glomerular injury is often characterized by the effacement of podocytes, loss of slit diaphragms, and proteinuria. Renal ischemia or the loss of blood flow to the kidneys has been widely associated with tubular and endothelial injury but rarely has been shown to induce podocyte damage and disruption of the slit diaphragm. In this study, we have used an in vivo rat ischemic model to demonstrate that renal ischemia induces podocyte effacement with loss of slit diaphragm and proteinuria. Biochemical analysis of the ischemic glomerulus shows that ischemia induces rapid loss of interaction between slit diaphragm junctional proteins Neph1 and ZO-1. To further understand the effect of ischemia on molecular interactions between slit diaphragm proteins, a cell culture model was employed to study the binding between Neph1 and ZO-1. Under physiologic conditions, Neph1 co-localized with ZO-1 at cell-cell contacts in cultured human podocytes. Induction of injury by ATP depletion resulted in rapid loss of Neph1 and ZO-1 binding and redistribution of Neph1 and ZO-1 proteins from cell membrane to the cytoplasm. Recovery resulted in increased Neph1 tyrosine phosphorylation, restoring Neph1 and ZO-1 binding and their localization at the cell membrane. We further demonstrate that tyrosine phosphorylation of Neph1 mediated by Fyn results in significantly increased Neph1 and ZO-1 binding, suggesting a critical role for Neph1 tyrosine phosphorylation in reorganizing the Neph1-ZO-1 complex. This study documents that renal ischemia induces dynamic changes in the molecular interactions between slit diaphragm proteins, leading to podocyte damage and proteinuria.     

Small and Wide Angle X-Ray Scattering Studies of Biological Macromolecules in Solution

In this paper, Small and Wide Angle X-ray Scattering (SWAXS) analysis of macromolecules is demonstrated through experimentation. SWAXS is a technique where X-rays are elastically scattered by an inhomogeneous sample in the nm-range at small angles (typically 0.1 - 5°) and wide angles (typically > 5°). This technique provides information about the shape, size, and distribution of macromolecules, characteristic distances of partially ordered materials, pore sizes, and surface-to-volume ratio. Small Angle X-ray Scattering (SAXS) is capable of delivering structural information of macromolecules between 1 and 200 nm, whereas Wide Angle X-ray Scattering (WAXS) can resolve even smaller Bragg spacing of samples between 0.33 nm and 0.49 nm based on the specific system setup and detector. The spacing is determined from Bragg''s law and is dependent on the wavelength and incident angle.In a SWAXS experiment, the materials can be solid or liquid and may contain solid, liquid or gaseous domains (so-called particles) of the same or another material in any combination. SWAXS applications are very broad and include colloids of all types: metals, composites, cement, oil, polymers, plastics, proteins, foods, and pharmaceuticals. For solid samples, the thickness is limited to approximately 5 mm.Usage of a lab-based SWAXS instrument is detailed in this paper. With the available software (e.g., GNOM-ATSAS 2.3 package by D. Svergun EMBL-Hamburg and EasySWAXS software) for the SWAXS system, an experiment can be conducted to determine certain parameters of interest for the given sample. One example of a biological macromolecule experiment is the analysis of 2 wt% lysozyme in a water-based aqueous buffer which can be chosen and prepared through numerous methods. The preparation of the sample follows the guidelines below in the Preparation of the Sample section. Through SWAXS experimentation, important structural parameters of lysozyme, e.g. the radius of gyration, can be analyzed.     

The Src homology 3 domain is required for junctional adhesion molecule binding to the third PDZ domain of the scaffolding protein ZO-1

Tight junctions are cell-cell contacts that regulate the paracellular flux of solutes and prevent pathogen entry across cell layers. The assembly and permeability of this barrier are dependent on the zonula occludens (ZO) membrane-associated guanylate kinase (MAGUK) proteins ZO-1, -2, and -3. MAGUK proteins are characterized by a core motif of protein-binding domains that include a PDZ domain, a Src homology 3 (SH3) domain, and a region of homology to guanylate kinase (GUK); the structure of this core motif has never been determined for any MAGUK. To better understand how ZO proteins organize the assembly of protein complexes we have crystallized the entire PDZ3-SH3-GUK core motif of ZO-1. We have also crystallized this core motif in complex with the cytoplasmic tail of the ZO-1 PDZ3 ligand, junctional adhesion molecule A (JAM-A) to determine how the activity of different domains is coordinated. Our study shows a new feature for PDZ class II ligand binding that implicates the two highly conserved Phe(-2) and Ser(-3) residues of JAM. Our x-ray structures and NMR experiments also show for the first time a role for adjacent domains in the binding of ligands to PDZ domains in the MAGUK proteins family.     

Self-Assembly of MinE on the Membrane Underlies Formation of the MinE Ring to Sustain Function of the Escherichia coli Min System

The pole-to-pole oscillation of the Min proteins in Escherichia coli results in the inhibition of aberrant polar division, thus facilitating placement of the division septum at the midcell. MinE of the Min system forms a ring-like structure that plays a critical role in triggering the oscillation cycle. However, the mechanism underlying the formation of the MinE ring remains unclear. This study demonstrates that MinE self-assembles into fibrillar structures on the supported lipid bilayer. The MinD-interacting domain of MinE shows amyloidogenic properties, providing a possible mechanism for self-assembly of MinE. Supporting the idea, mutations in residues Ile-24 and Ile-25 of the MinD-interacting domain affect fibril formation, membrane binding ability of MinE and MinD, and subcellular localization of three Min proteins. Additional mutations in residues Ile-72 and Ile-74 suggest a role of the C-terminal domain of MinE in regulating the folding propensity of the MinD-interacting domain for different molecular interactions. The study suggests a self-assembly mechanism that may underlie the ring-like structure formed by MinE-GFP observed in vivo.     

Domain-swapped dimerization of the second PDZ domain of ZO2 may provide a structural basis for the polymerization of claudins

Zonula occludens proteins (ZOs), including ZO1/2/3, are tight junction-associated proteins. Each of them contains three PDZ domains. It has been demonstrated that ZO1 can form either homodimers or heterodimers with ZO2 or ZO3 through the second PDZ domain. However, the underlying structural basis is not well understood. In this study, the solution structure of the second PDZ domain of ZO2 (ZO2-PDZ2) was determined using NMR spectroscopy. The results revealed a novel dimerization mode for PDZ domains via three-dimensional domain swapping, which can be generalized to homodimers of ZO1-PDZ2 or ZO3-PDZ2 and heterodimers of ZO1-PDZ2/ZO2-PDZ2 or ZO1-PDZ2/ZO3-PDZ2 due to high conservation between PDZ2 domains in ZO proteins. Furthermore, GST pulldown experiments and immunoprecipitation studies demonstrated that interactions between ZO1-PDZ2 and ZO2-PDZ2 and their self-associations indeed exist both in vitro and in vivo. Chemical cross-linking and dynamic laser light scattering experiments revealed that both ZO1-PDZ2 and ZO2-PDZ2 can form oligomers in solution. This PDZ domain-mediated oligomerization of ZOs may provide a structural basis for the polymerization of claudins, namely the formation of tight junctions.     

Neph1, a component of the kidney slit diaphragm, is tyrosine-phosphorylated by the Src family tyrosine kinase and modulates intracellular signaling by binding to Grb2

There are several lines of evidence that the podocyte slit diaphragm (SD), which serves as a structural framework for the filtration barrier in kidney glomerulus, also plays an essential role as a signaling platform. Several SD components including nephrin and TRPC6 are known to be phosphorylated by a Src family tyrosine kinase (SFK), Fyn. Here we have characterized Neph1, another SD component, as a novel substrate of SFK. Fyn interacts with and phosphorylates the cytoplasmic domain of Neph1 in vitro and in intact cells. Peptide mass fingerprinting and site-directed mutagenesis identified several tyrosine phosphorylation sites. In pull-down assays using rat glomerular lysates, Neph1 but not nephrin specifically binds to adaptor protein Grb2 and tyrosine kinase Csk in a phosphorylation-dependent manner. Both tyrosine 637 and 638 of Neph1 are crucial for Neph1-Grb2 binding. Phosphorylation of tyrosine 637 is significantly up-regulated in in vivo models of podocyte injury. Furthermore, Neph1 attenuates ERK activation elicited by Fyn, and this inhibitory effect requires the intact binding motif for the Grb2 SH2 domain. Our results shown here demonstrate that Neph1 is a novel in vivo substrate of SFK and suggest that Neph1 modulates ERK signaling through phosphorylation-dependent interaction with Grb2. Thus, SFK orchestrates a wide spectrum of protein-protein interactions and intracellular signaling networks at SD through tyrosine phosphorylation.     

The Structure of the Ternary Complex of Krev Interaction Trapped 1 (KRIT1) Bound to Both the Rap1 GTPase and the Heart of Glass (HEG1) Cytoplasmic Tail

Loss of function mutation in Krev interaction trapped 1 (KRIT1) causes autosomal dominant familial cerebral cavernous malformations and disrupts cardiovascular development. The biological function of KRIT1 requires that its FERM (band 4.1, ezrin, radixin, moesin) domain physically interact with both the small GTPase Rap1 and the cytoplasmic tail of the Heart of glass (HEG1) membrane anchor. In this study, we show that the KRIT1 FERM domain can bind both Rap1 and HEG1 simultaneously, and we solved the crystal structure of the KRIT1-Rap1-HEG1 ternary complex. Rap1 binds on the surface of the F1 and F2 subdomains, in an interaction that leaves its Switch II region accessible to other potential effectors. HEG1 binds in a hydrophobic pocket at the KRIT1 F1 and F3 interface, and there is no overlap with the Rap1-binding site. Indeed, the affinity of KRIT1 or the KRIT1-Rap1 complex for HEG1 is comparable (K_d = 1.2 and 0.96 μm, respectively) showing that there is no competition between the two sites. Furthermore, analysis of this structure revealed a specific ionic interaction between the F2 lobe of KRIT1 and Rap1 that could explain the remarkable Rap1 specificity of KRIT1. This structural insight enabled design of KRIT1(K570I), a mutant that binds Rap1 with 8-fold lower affinity and exhibits increased binding to HRas. These data show that HEG1 can recruit the Rap1-KRIT complex to the plasma membrane where Rap1''s Switch II region remains accessible and reveals an important determinant of KRIT1''s specificity for Rap1.     

The BAR domain protein PICK1 regulates cell recognition and morphogenesis by interacting with Neph proteins

Neph proteins are evolutionarily conserved membrane proteins of the immunoglobulin superfamily that control the formation of specific intercellular contacts. Cell recognition through these proteins is essential in diverse cellular contexts such as patterning of the compound eye in Drosophila melanogaster, neuronal connectivity in Caenorhabditis elegans, and the formation of the kidney filtration barrier in mammals. Here we identify the PDZ and BAR domain protein PICK1 (protein interacting with C-kinase 1) as a Neph-interacting protein. Binding required dimerization of PICK1, was dependent on PDZ domain protein interactions, and mediated stabilization of Neph1 at the plasma membrane. Moreover, protein kinase C (PKCα) activity facilitated the interaction through releasing Neph proteins from their binding to the multidomain scaffolding protein zonula occludens 1 (ZO-1), another PDZ domain protein. In Drosophila, the Neph homologue Roughest is essential for sorting of interommatidial precursor cells and patterning of the compound eye. RNA interference-mediated knockdown of PICK1 in the Drosophila eye imaginal disc caused a Roughest destabilization at the plasma membrane and a phenotype that resembled rst mutation. These data indicate that Neph proteins and PICK1 synergistically regulate cell recognition and contact formation.     

Crystal structure of N-glycosylated human glypican-1 core protein: structure of two loops evolutionarily conserved in vertebrate glypican-1

Glypicans are a family of cell-surface proteoglycans that regulate Wnt, hedgehog, bone morphogenetic protein, and fibroblast growth factor signaling. Loss-of-function mutations in glypican core proteins and in glycosaminoglycan-synthesizing enzymes have revealed that glypican core proteins and their glycosaminoglycan chains are important in shaping animal development. Glypican core proteins consist of a stable α-helical domain containing 14 conserved Cys residues followed by a glycosaminoglycan attachment domain that becomes exclusively substituted with heparan sulfate (HS) and presumably adopts a random coil conformation. Removal of the α-helical domain results in almost exclusive addition of the glycosaminoglycan chondroitin sulfate, suggesting that factors in the α-helical domain promote assembly of HS. Glypican-1 is involved in brain development and is one of six members of the vertebrate family of glypicans. We expressed and crystallized N-glycosylated human glypican-1 lacking HS and N-glycosylated glypican-1 lacking the HS attachment domain. The crystal structure of glypican-1 was solved using crystals of selenomethionine-labeled glypican-1 core protein lacking the HS domain. No additional electron density was observed for crystals of glypican-1 containing the HS attachment domain, and CD spectra of the two protein species were highly similar. The crystal structure of N-glycosylated human glypican-1 core protein at 2.5 Å, the first crystal structure of a vertebrate glypican, reveals the complete disulfide bond arrangement of the conserved Cys residues, and it also extends the structural knowledge of glypicans for one α-helix and two long loops. Importantly, the loops are evolutionarily conserved in vertebrate glypican-1, and one of them is involved in glycosaminoglycan class determination.     

Structural changes in the carboxyl terminus of the gap junction protein connexin43 indicates signaling between binding domains for c-Src and zonula occludens-1

Regulation of cell-cell communication by the gap junction protein connexin43 can be modulated by a variety of connexin-associating proteins. In particular, c-Src can disrupt the connexin43 (Cx43)-zonula occludens-1 (ZO-1) interaction, leading to down-regulation of gap junction intercellular communication. The binding sites for ZO-1 and c-Src correspond to widely separated Cx43 domains (approximately 100 residues apart); however, little is known about the structural modifications that may allow information to be transferred over this distance. Here, we have characterized the structure of the connexin43 carboxyl-terminal domain (Cx43CT) to assess its ability to interact with domains from ZO-1 and c-Src. NMR data indicate that the Cx43CT exists primarily as an elongated random coil, with two regions of alpha-helical structure. NMR titration experiments determined that the ZO-1 PDZ-2 domain affected the last 19 Cx43CT residues, a region larger than that reported to be required for Cx43CT-ZO-1 binding. The c-Src SH3 domain affected Cx43CT residues Lys-264-Lys-287, Ser-306-Glu-316, His-331-Phe-337, Leu-356-Val-359, and Ala-367-Ser-372. Only region Lys-264-Lys-287 contains the residues previously reported to act as an SH3 binding domain. The specificity of these interactions was verified by peptide competition experiments. Finally, we demonstrated that the SH3 domain could partially displace the Cx43CT-PDZ-2 complex. These studies represent the first structural characterization of a connexin domain when integrated in a multimolecular complex. Furthermore, we demonstrate that the structural characteristics of a disordered Cx43CT are advantageous for signaling between different binding partners that may be important in describing the mechanism of channel closure or internalization in response to pathophysiological stimuli.     

ZO-3, a Novel Member of the MAGUK Protein Family Found at the Tight Junction, Interacts with ZO-1 and Occludin

A 130-kD protein that coimmunoprecipitates with the tight junction protein ZO-1 was bulk purified from Madin-Darby canine kidney (MDCK) cells and subjected to partial endopeptidase digestion and amino acid sequencing. A resulting 19–amino acid sequence provided the basis for screening canine cDNA libraries. Five overlapping clones contained a single open reading frame of 2,694 bp coding for a protein of 898 amino acids with a predicted molecular mass of 98,414 daltons. Sequence analysis showed that this protein contains three PSD-95/SAP90, discs-large, ZO-1 (PDZ) domains, a src homology (SH3) domain, and a region similar to guanylate kinase, making it homologous to ZO-1, ZO-2, the discs large tumor suppressor gene product of Drosophila, and other members of the MAGUK family of proteins. Like ZO-1 and ZO-2, the novel protein contains a COOH-terminal acidic domain and a basic region between the first and second PDZ domains. Unlike ZO-1 and ZO-2, this protein displays a proline-rich region between PDZ2 and PDZ3 and apparently contains no alternatively spliced domain. MDCK cells stably transfected with an epitope-tagged construct expressed the exogenous polypeptide at an apparent molecular mass of ~130 kD. Moreover, this protein colocalized with ZO-1 at tight junctions by immunofluorescence and immunoelectron microscopy. In vitro affinity analyses demonstrated that recombinant 130-kD protein directly interacts with ZO-1 and the cytoplasmic domain of occludin, but not with ZO-2. We propose that this protein be named ZO-3.     

Structural Basis of the Binding of Merlin FERM Domain to the E3 Ubiquitin Ligase Substrate Adaptor DCAF1

The tumor suppressor gene Nf2 product, Merlin, plays vital roles in controlling proper development of organ sizes by specifically binding to a large number of target proteins localized both in cytoplasm and nuclei. The FERM domain of Merlin is chiefly responsible for its binding to target proteins, although the molecular basis governing these interactions are poorly understood due to lack of structural information. Here, we report the crystal structure of the Merlin FERM domain in complex with its binding domain derived from the E3 ubiquitin ligase substrate adaptor DCAF1 (also known as VPRBP). Unlike target binding modes found in ERM proteins, the Merlin-FERM binding domain of DCAF1 folds as a β-hairpin and binds to the α1/β5-groove of the F3 lobe of Merlin-FERM via extensive hydrophobic interactions. In addition to providing the first structural glimpse of a Merlin-FERM·target complex, the structure of the Merlin·DCAF1 complex is likely to be valuable for understanding the interactions of Merlin with its binding partners other than DCAF1.     

Nephrin and Neph1 co-localize at the podocyte foot process intercellular junction and form cis hetero-oligomers

Glomerular visceral epithelial cells (podocytes) appear to play a central role in maintaining the selective filtration barrier of the renal glomerulus. While the immunoglobulin superfamily member Nephrin was proposed to act as a cell adhesion molecule at the podocyte intercellular junction necessary for maintaining glomerular perm selectivity, the Nephrin ligand has not been identified. The existence of a new subfamily of Nephrin-like molecules including Neph1 was recently described. Genetic deletion of Nephrin or Neph1 resulted in similar phenotypes of podocyte foot process effacement and proteinuria. The subcellular localization of Neph1 and the possibility that Nephrin and Neph1 interact was investigated. Polyclonal antiserum for Neph1 was raised and characterized. Neph1 migrated as a 90-kDa protein on SDS-PAGE under reducing conditions. Neph1 was identified in a glomerular and podocyte-specific distribution in adult rat kidney. Like Nephrin and Podocin, Neph1 was enriched in Triton X-100 detergent-resistant membrane fractions. Consistent with this observation, immunogold electron microscopy demonstrated that Neph1 localized exclusively to lateral margins of podocyte foot processes at the insertion of the slit diaphragm. Neph1 and Nephrin participate in a direct cis-interaction involving their cytoplasmic domains. In addition, interactions between the extracellular domain of Nephrin and itself and between the extracellular domain of Nephrin and that of Neph1 were detected. Neph1 did not interact via a homophilic interaction. These observations suggest that Nephrin and Neph1 form a hetero-oligomeric receptor complex in the plane of the membrane that might interact across the foot process intercellular junction through interactions between Nephrin with itself and Neph1.     

In vitro protein complex formation with cytoskeleton-anchoring domain of occludin identified by limited proteolysis

Occludin is an essential membrane protein component of cellular tight junctions, participating in both cell-cell adhesion in the paracellular space and anchoring of the junctional complex to the cytoskeleton. The latter function is accomplished through binding of the C-terminal cytoplasmic region to scaffolding proteins that mediate binding to cytoskeletal actin. We isolated a structural domain from both the bacterial-expressed C-terminal cytoplasmic region of human occludin and native cellular occludin, extracted from epithelial (Madin-Darby canine kidney) or endothelial (human brain) cells, by limited proteolysis with trypsin. This human occludin domain contains the last 119 amino acids as identified by N-terminal sequencing and peptide mass fingerprinting using matrix-assisted laser desorption ionization-time of flight mass spectrometry. Based on the sequence and secondary structure prediction, this domain contains 4 of 5 alpha-helices in the C-terminal region and is linked to the fourth membrane-spanning region by a loosely structured tethering polypeptide. Comparison of circular dichroism spectra of recombinant proteins corresponding to the entire C-terminal region versus only the binding domain region also supports the interpretation that the helical structural elements are concentrated in that domain. Co-immunoprecipitation of this domain with ZO-2 demonstrated preservation of the specificity of the scaffolding protein-binding function, and binding studies with immobilized ZO-2 suggest the presence of multiple ZO-2 binding sites in this domain. These results provide a basis for development of a structural model of the ZO-binding site that can be used to investigate regulation of tight junction anchoring by intracellular signaling events.     

Molecular dynamics simulation of human immunodeficiency virus protein U (Vpu) in lipid/water Langmuir monolayer

Virus protein U (Vpu) is an accessory membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). Various NMR and CD studies have shown that the transmembrane domain of Vpu has a helical conformation and that the cytoplasmic domain adopts the helix-loop-helix-turn motif. This 3.5-ns molecular dynamics (MD) simulation of Vpu in a lipid/membrane environment has fully reproduced these structural characteristics. Membrane propensities of two amphipathic helices in the cytoplasmic domain are further compared here to understand better their complicated orientational behavior known from experiment. This study first reveals that the highly conserved loop region in the cytoplasmic domain can be closely associated with the membrane surface. It is known from the simulation that Vpu is associated with 34 lipids in this Langmuir monolayer. The lipids that are located between the Vpu transmembrane helix and the first helix in the cytoplasmic domain are pushed up by Vpu. These elevated lipids have increased P-N tilt angles for the head groups but unchanged acyl-chain tilt angles compared with lipids that do not interact with Vpu. This study verifies the significance of applying MD simulation in refining protein structure and revealing detailed protein-lipid interaction in membrane/water environment. Figure XZ view of a snapshot of Vpu/DLGPC/water system after 3.5 ns NP(N)gamma T MD simulation. Coloring scheme: Vpu, red; C, green; H, pink; N, blue; O, orange; P, magenta; water, light blue