Obtaining Protein Association Energy Landscape for Integral Membrane Proteins

Biophysics of claudin proteins in tight junction architecture: Three decades of progress

Tight junctions are cell-cell adhesion complexes that act as gatekeepers of the paracellular space. Formed by several transmembrane proteins, the claudin family performs the primary gate-keeping function. The claudin proteins form charge and size-selective diffusion barriers to maintain homeostasis across endothelial and epithelial tissue. Of the 27 known claudins in mammals, some are known to seal the paracellular space, while others provide selective permeability. The differences in permeability arise due to the varying expression levels of claudins in each tissue. The tight junctions are observed as strands in freeze-fracture electron monographs; however, at the molecular level, tight junction strands form when multiple claudin proteins assemble laterally (cis assembly) within a cell and head-on (trans assembly) with claudins of the adjacent cell in a zipper-like architecture, closing the gap between the neighboring cells. The disruption of tight junctions caused by changing claudin expression levels or mutations can lead to diseases. Therefore, knowledge of the molecular architecture of the tight junctions and how that is tied to tissue-specific function is critical for fighting diseases. Here, we review the current understanding of the tight junctions accrued over the last three decades from experimental and computational biophysics perspectives.

Combined Computational-Biochemical Approach Offers an Accelerated Path to Membrane Protein Solubilization

Membrane proteins are difficult to isolate and purify due to their dependence on the surrounding lipid membrane for structural stability. Detergents are often used to solubilize these proteins, with this approach requiring a careful balance between protein solubilization and denaturation. Determining which detergent is most appropriate for a given protein has largely been done empirically through screening, which requires large amounts of membrane protein and associated resources. Here, we describe an alternative to conventional detergent screening using a computational modeling approach to identify the most likely candidate detergents for solubilizing a protein of interest. We demonstrate our approach using ghrelin O-acyltransferase (GOAT), a member of the membrane-bound O-acyltransferase family of integral membrane enzymes that has not been solubilized or purified in active form. A computationally derived GOAT structural model provides the only structural information required for this approach. Using computational analysis of detergent ability to penetrate phospholipid bilayers and stabilize the GOAT structure, a panel of common detergents were rank-ordered for their proposed ability to solubilize GOAT. The simulations were performed at all-atom resolution for a combined simulation time of 24 μs. Independently, we biologically screened these detergents for their solubilization of fluorescently tagged GOAT constructs. We found computational prediction of protein structural stabilization was the better predictor of detergent solubilization ability, but neither approach was effective for predicting detergents that would support GOAT enzymatic function. The current rapid expansion of membrane protein computational models lacking experimental structural information and our computational detergent screening approach can greatly improve the efficiency of membrane protein detergent solubilization, supporting downstream functional and structural studies.

Claudin-23 reshapes epithelial tight junction architecture to regulate barrier function

Claudin family tight junction proteins form charge- and size-selective paracellular channels that regulate epithelial barrier function. In the gastrointestinal tract, barrier heterogeneity is attributed to differential claudin expression. Here, we show that claudin-23 (CLDN23) is enriched in luminal intestinal epithelial cells where it strengthens the epithelial barrier. Complementary approaches reveal that CLDN23 regulates paracellular ion and macromolecule permeability by associating with CLDN3 and CLDN4 and regulating their distribution in tight junctions. Computational modeling suggests that CLDN23 forms heteromeric and heterotypic complexes with CLDN3 and CLDN4 that have unique pore architecture and overall net charge. These computational simulation analyses further suggest that pore properties are interaction-dependent, since differently organized complexes with the same claudin stoichiometry form pores with unique architecture. Our findings provide insight into tight junction organization and propose a model whereby different claudins combine to form multiple distinct complexes that modify epithelial barrier function by altering tight junction structure.

The impact of pathogenic and artificial mutations on Claudin-5 selectivity from molecular dynamics simulations

Tight-junctions (TJs) are multi-protein complexes between adjacent endothelial or epithelial cells. In the blood-brain-barrier (BBB), they seal the paracellular space and the Claudin-5 (Cldn5) protein forms their backbone. Despite the fundamental role in brain homeostasis, little is known on Cldn5-based TJ assemblies. Different structural models were suggested, with Cldn5 protomers generating paracellular pores that restrict the passage of ions and small molecules. Recently, the first Cldn5 pathogenic mutation, G60R, was identified and shown to induce Cl^--selective channels and Na⁺ barriers in BBB TJs, providing an excellent opportunity to validate the structural models. Here, we used molecular dynamics to study the permeation of ions and water through two distinct G60R-Cldn5 paracellular architectures. Only the so-called Pore I reproduces the functional modification observed in experiments, displaying a free energy (FE) minimum for Cl^- and a barrier for Na⁺ consistent with anionic selectivity. We also studied the artificial Q57D and Q63D mutations in the constriction region, Q57 being conserved in Cldns except for cation permeable homologs. In both cases, we obtain FE profiles consistent with facilitated passage of cations. Our calculations provide the first in-silico description of a Cldn5 pathogenic mutation, further assessing the TJ Pore I model and yielding new insight on BBB's paracellular selectivity.

Computational Assessment of Different Structural Models for Claudin-5 Complexes in Blood-Brain Barrier Tight Junctions

The blood-brain barrier (BBB) strictly regulates the exchange of ions and molecules between the blood and the central nervous system. Tight junctions (TJs) are multimeric structures that control the transport through the paracellular spaces between the adjacent brain endothelial cells of the BBB. Claudin-5 (Cldn5) proteins are essential for TJ formation and assemble into multiprotein complexes via cis-interactions within the same cell membrane and trans-interactions across two contiguous cells. Despite the relevant biological function of Cldn5 proteins and their role as targets of brain drug delivery strategies, the molecular details of their assembly within TJs are still unclear. Two different structural models have been recently introduced, in which Cldn5 dimers belonging to opposite cells join to generate paracellular pores. However, a comparison of these models in terms of ionic transport features is still lacking. In this work, we used molecular dynamics simulations and free energy (FE) calculations to assess the two Cldn5 pore models and investigate the thermodynamic properties of water and physiological ions permeating through them. Despite different FE profiles, both structures present single/multiple FE barriers to ionic permeation, while being permissive to water flux. These results reveal that both models are compatible with the physiological role of Cldn5 TJ strands. By identifying the protein-protein surface at the core of TJ Cldn5 assemblies, our computational investigation provides a basis for the rational design of synthetic peptides and other molecules capable of opening paracellular pores in the BBB.

The Ism between Endothelial Cilia and Endothelial Nanotubules Is an Evolving Concept in the Genesis of the BBB

The blood–brain barrier (BBB) is fundamental in maintaining central nervous system (CNS) homeostasis by regulating the chemical environment of the underlying brain parenchyma. Brain endothelial cells (BECs) constitute the anatomical and functional basis of the BBB. Communication between adjacent BECs is critical for establishing BBB integrity, and knowledge of its nanoscopic landscape will contribute to our understanding of how juxtaposed zones of tight-junction protein interactions between BECs are aligned. The review discusses and critiques types of nanostructures contributing to the process of BBB genesis. We further critically evaluate earlier findings in light of novel high-resolution electron microscopy descriptions of nanoscopic tubules. One such phenotypic structure is BEC cytoplasmic projections, which, early in the literature, is postulated as brain capillary endothelial cilia, and is evaluated and compared to the recently discovered nanotubules (NTs) formed in the paracellular spaces between BECs during barrier-genesis. The review attempts to elucidate a myriad of unique topographical ultrastructures that have been reported to be associated with the development of the BBB, viz., structures ranging from cilia to BEC tunneling nanotubules (TUNTs) and BEC tethering nanotubules (TENTs).

Engineering of 2D nanomaterials to trap and kill SARS-CoV-2: a new insight from multi-microsecond atomistic simulations

In late 2019, coronavirus disease 2019 (COVID-19) was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Spike protein is one of the surface proteins of SARS-CoV-2 that is essential for its infectious function. Therefore, it received lots of attention for the preparation of antiviral drugs, vaccines, and diagnostic tools. In the current study, we use computational methods of chemistry and biology to study the interaction between spike protein and its receptor in the body, angiotensin-I-converting enzyme-2 (ACE2). Additionally, the possible interaction of two-dimensional (2D) nanomaterials, including graphene, bismuthene, phosphorene, p-doped graphene, and functionalized p-doped graphene, with spike protein is investigated. The functionalized p-doped graphene nanomaterials were found to interfere with spike protein better than the other tested nanomaterials. In addition, the interaction of the proposed nanomaterials with the main protease (M^pro) of SARS-CoV-2 was studied. Functionalized p-doped graphene nanomaterials showed more capacity to prevent the activity of M^pro. These 2D nanomaterials efficiently reduce the transmissibility and infectivity of SARS-CoV-2 by both the deformation of the spike protein and inhibiting the M^pro. The results suggest the potential use of 2D nanomaterials in a variety of prophylactic approaches, such as masks or surface coatings, and would deserve further studies in the coming years.

All-Atom Molecular Dynamics Elucidating Molecular Mechanisms of Single-Transmembrane Model Peptide Dimerization in a Lipid Bilayer

Protein-protein interactions between transmembrane helices are essential elements for membrane protein structures and functions. To understand the effects of peptide sequences and lipid compositions on these interactions, single-molecule experiments using model systems comprising artificial peptides and membranes have been extensively performed. However, their dynamic behavior at the atomic level remains largely unclear. In this study, we applied the all-atom molecular dynamics (MD) method to simulate the interactions of single-transmembrane helical peptide dimers in membrane environments, which has previously been analyzed by single-molecule experiments. The simulations were performed with two peptides (Ala- and Leu-based artificially designed peptides, termed "host peptide", and the host peptide added with the GXXXG motif, termed "GXXXG peptide"), two membranes (pure-POPC and POPC mixed with 30% cholesterols), and two dimer directions (parallel and antiparallel), consistent with those in the previous experiment. As a result, the MD simulations with parallel dimers reproduced the experimentally observed tendency that introducing cholesterols weakened the interactions in the GXXXG dimer and facilitated those in the host dimer. Our simulation suggested that the host dimer formed hydrogen bonds but the GXXXG dimer did not. However, some discrepancies were also observed between the experiments and simulations. Limitations in the space and time scales of simulations restrict the large-scale undulation and peristaltic motions of the membranes, resulting in differences in lateral pressure profiles. This effect could cause a discrepancy in the rotation angles of helices against the membrane normal.

Model of Environmental Membrane Field for Transmembrane Proteins

The water environment determines the activity of biological processes. The role of such an environment interpreted in the form of an external field expressed by the 3D Gaussian distribution in the fuzzy oil drop model directs the folding process towards the generation of a centrally located hydrophobic core with the simultaneous exposure of polar residues on the surface. In addition to proteins soluble in the water environment, there is a significant group of membrane proteins that act as receptors or channels, including ion channels in particular. The change of the polar (water) environment into a highly hydrophobic (membrane) environment is quite radical, resulting in a different hydrophobicity distribution within the membrane protein. Modification of the notation of the force field expressing the presence of the hydrophobic environment has been proposed in this work. A modified fuzzy oil drop model with its adaptation to membrane proteins was used to interpret the structure of membrane proteins–mechanosensitive channel. The modified model was also used to describe the so-called negative cases—i.e., for water-soluble proteins with a clear distribution consistent with the fuzzy oil drop model.

Molecular architecture and assembly of the tight junction backbone

The functional and structural concept of tight junctions has developed after discovery of claudin and TAMP proteins. Many of these proteins contribute to epi- and endothelial barrier but some, in contrast, form paracellular channels. Claudins form the backbone of tight junction (TJ) strands whereas other proteins regulate TJ dynamics. The current joined double-row model of TJ strands and channels is crucially based on the linear alignment of claudin-15 in the crystal. Molecular dynamics simulations, protein docking, mutagenesis, cellular TJ reconstitution, and electron microscopy studies largely support stability and functionality of the model. Here, we summarize in silico and in vitro data about TJ strand assembly including comparison of claudin crystal structures and alternative models. Sequence comparisons, experimental and structural data substantiate differentiation of classic and non-classic claudins differing in motifs related to strand assembly. Classic claudins seem to share a similar mechanism of strand formation. Interface variations likely contribute to TJ strand flexibility. Combined in vitro/in silico studies are expected to elucidate mechanistic keys determining TJ regulation.

Computational Nanoscopy of Tight Junctions at the Blood-Brain Barrier Interface

The selectivity of the blood-brain barrier (BBB) is primarily maintained by tight junctions (TJs), which act as gatekeepers of the paracellular space by blocking blood-borne toxins, drugs, and pathogens from entering the brain. The BBB presents a significant challenge in designing neurotherapeutics, so a comprehensive understanding of the TJ architecture can aid in the design of novel therapeutics. Unraveling the intricacies of TJs with conventional experimental techniques alone is challenging, but recently developed computational tools can provide a valuable molecular-level understanding of TJ architecture. We employed the computational methods toolkit to investigate claudin-5, a highly expressed TJ protein at the BBB interface. Our approach started with the prediction of claudin-5 structure, evaluation of stable dimer conformations and nanoscale assemblies, followed by the impact of lipid environments, and posttranslational modifications on these claudin-5 assemblies. These led to the study of TJ pores and barriers and finally understanding of ion and small molecule transport through the TJs. Some of these in silico, molecular-level findings, will need to be corroborated by future experiments. The resulting understanding can be advantageous towards the eventual goal of drug delivery across the BBB. This review provides key insights gleaned from a series of state-of-the-art nanoscale simulations (or computational nanoscopy studies) performed on the TJ architecture.