Tight junctions (TJs) are key players in determining tissue-specific paracellular permeability across epithelial and endothelial membranes. Claudin proteins, the primary determinants of TJs structure and functionality, assemble in paracellular spaces to form channels and pores that are charge and size selective. Here, using molecular dynamics (MD) simulations, we elucidate the molecular assembly of claudin-3 and claudin-5 proteins of blood-brain barrier TJs. Despite having a high degree of sequence and structural similarity, these two claudins form different types of cis-interactions. Molecular docking of the observed cis-interfaces into trans-forms revealed two putative pore models that were also observed in the self-assembly simulations. The observed pore structures (pore I and II) have pore-lining residues that have been previously reported in the literature. The pore I model is consistent with a previously reported claudin-15 model. The pore II model, also consistent with biochemical results, has not been reported previously. Further analysis using in silico site-directed mutations provide convincing support for the validity of the pore II model. Using steered MD and umbrella sampling, we computed the transport properties of water and α-d-glucose through pore II. The study offers new insight into the selectivity of blood-brain barrier TJs.

The blood-brain barrier (BBB) constituted by claudin-5 tight junctions is critical in maintaining the homeostasis of the central nervous system, but this highly selective molecular interface is an impediment for therapeutic interventions in neurodegenerative and neurological diseases. Therapeutic strategies that can exploit the paracellular transport remain elusive due to lack of molecular insights of the tight junction assembly. This study focuses on analyzing the membrane driven cis interactions of claudin-5 proteins in the formation of the BBB tight junctions. We have adopted a synergistic approach employing in silico multiscale dynamics and in vitro cross-linking experiments to study the claudin-5 interactions. Long time scale simulations of claudin-5 monomers, in seven different lipid compositions, show formation of cis dimers that subsequently aggregate into strands. In vitro formaldehyde cross-linking studies also conclusively show that cis-interacting claudin-5 dimers cross-link with short methylene spacers. Using this synergistic approach, we have identified five unique dimer interfaces in our simulations that correlate with the cross-linking experiments, four of which are mediated by transmembrane (TM) helices and the other mediated by extracellular loops (ECL). Potential of mean force calculations of these five dimers revealed that the TM mediated interfaces, which can have distinctive leucine zipper interactions in some cases, are more stable than the ECL mediated interface. Additionally, simulations show that claudin-5 dimerization is significantly influenced by the lipid microenvironment. This study captures the fundamental interactions responsible for the BBB tight junction assembly and offers a framework for extending this work to other tight junctions found in the body.

No abstract

The cell envelope of Gram-negative bacteria contains a lipopolysaccharide (LPS) rich outer membrane that acts as the first line of defense for bacterial cells in adverse physical and chemical environments. The LPS macromolecule has a negatively charged oligosaccharide domain that acts as an ionic brush, limiting the permeability of charged chemical agents through the membrane. Besides the LPS, the outer membrane has radially extending O-antigen polysaccharide chains and β-barrel membrane proteins that make the bacterial membrane physiologically unique compared to phospholipid cell membranes. Elucidating the interplay of these contributing macromolecular components and their role in the integrity of the bacterial outer membrane remains a challenge. To bridge the gap in our current understanding of the Gram-negative bacterial membrane, we have developed a coarse grained force field for outer membrane that is computationally affordable for simulating dynamical process over physiologically relevant time scales. The force field was benchmarked against available experimental and atomistic simulations data for properties such as membrane thickness, density profiles of the residues, area per lipid, gel to liquid-crystalline phase transition temperatures, and order parameters. More than 17 membrane compositions were studied with a combined simulation time of over 100 μs. A comparison of simulated structural and dynamical properties with corresponding experimental data shows that the developed force field reproduces the overall physiology of LPS rich membranes. The affordability of the developed model for long time scale simulations can be instrumental in determining the mechanistic aspects of the antimicrobial action of chemical agents as well as assist in designing antimicrobial peptides with enhanced outer membrane permeation properties.

Tight junction (TJ) protein assembly controls permeability across epithelial and endothelial cells; thus, biochemical interactions that control the TJ assembly have physiological and biomedical significance. In this work, we employed multiscale simulations to probe the TJ self-assembly of five classic claudins (-1, -2, -4, -15, and -19). Claudin proteins assembled into dimeric and occasionally trimeric interfaces that subsequently formed larger polymeric strands. Using orientation-angle analysis to decompose polymeric strands, we found that individual claudins prefer certain dimer interfaces to others. Despite variations in the exact dimer populations observed in individual claudins, there appears to be an overall conformational uniformity in the type of dimeric interactions formed by the claudin family of proteins. A detailed structural characterization of the trimeric assemblies revealed that they could be putative receptors for trimeric Clostridium perfringens enterotoxin. Full characterization of the claudin-2 dimer interface revealed a cysteine cross-linkable interaction, which could be assembled into a symmetric pore of 7.4 Å average diameter. We extended the analysis of pore structure to other classic claudins and found that the distribution of polar residues lining the pore volume varied considerably between the barrier- and pore-forming claudins, potentially delineating the functionality in classic claudins.