Next Article in Journal
Neutrophil Maturation and Survival Is Controlled by IFN-Dependent Regulation of NAMPT Signaling
Next Article in Special Issue
Tilivalline- and Tilimycin-Independent Effects of Klebsiella oxytoca on Tight Junction-Mediated Intestinal Barrier Impairment
Previous Article in Journal
Interspecies Outer Membrane Vesicles (OMVs) Modulate the Sensitivity of Pathogenic Bacteria and Pathogenic Yeasts to Cationic Peptides and Serum Complement
Previous Article in Special Issue
Mouse Models of Human Claudin-Associated Disorders: Benefits and Limitations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

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

Department of Biomedical and Chemical Engineering, Syracuse University, Syracuse, NY 13244, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(22), 5583; https://doi.org/10.3390/ijms20225583
Submission received: 15 October 2019 / Revised: 5 November 2019 / Accepted: 6 November 2019 / Published: 8 November 2019
(This article belongs to the Special Issue The Tight Junction and Its Proteins: More Than Just a Barrier)

Abstract

:
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.

Graphical Abstract

1. The Blood–Brain Barrier

The blood–brain barrier (BBB) is a selective permeability barrier that impedes the influx of potentially harmful blood-borne chemicals into the brain [1,2,3,4]. Although discovered by Ehrlich in 1885, the term blood–brain barrier was coined almost four decades later in the 1920s, after a series of reports by Lewandowsky, Goldmann, Stern, and Gautier [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. Their collective work in vertebrate brain showed that tracers and toxins were unable to permeate the blood–brain interface; for example, trypan blue dye injected into the brain of a mouse remained confined to the brain, whereas the same dye injected outside the brain stained the entire body except the brain (Figure 1a). This permeability barrier is considered an evolutionary adaptation in vertebrates, and is vital for maintaining homeostasis in the CNS [1,5,21,22,23]. Unlike other organs, the brain lacks energy reserves [24,25,26]; therefore, for normal function, it requires a constant supply of nutrients from the bloodstream [27,28]. The brain’s microvasculature regulates the transfer of essential nutrients while protecting the brain against small chemical imbalances, which can lead to serious disease pathologies [29,30,31,32,33].
A critical component of the BBB is a specialized physical barrier comprised of tight junctions (TJs) present along the apical side of endothelial cells [1,23,34]. These cell-to-cell adhesion structures act as gatekeepers of the paracellular pathway (Figure 1b) by regulating the passive diffusion of molecules and ions into the brain [23,35,36]. The establishment of TJs at the BBB leads to a high electric resistance (>1000–3000 Ω cm2) across the interface [37,38]. The presence of TJs also leads to cell polarization whereby the apical and basolateral sides of the endothelia are isolated [2,35]. The polarized endothelia control transcellular permeability by specific expression of transporters and receptor proteins in the apical side of the membrane (Figure 1b). In the event of disease pathologies, TJs are compromised, leading to the loss of barrier selectivity in both the transcellular and paracellular pathways [33,35].
Although the BBB is vital for neurological functions, its selective permeability is a deterrent to drug delivery into the brain [31,39,40,41]. The BBB obstructs effective treatment of neurodegenerative diseases such as Alzheimer’s and Parkinson’s, because therapeutics are unable to cross the barrier [39]. To develop new therapeutics, in vivo animal models provide an experimental environment that mimics the intricacies of human physiology. However, high cost and lack of complete translation of disease pathologies from animal models to humans make in vivo studies less effective for drug development. In vitro models are therefore being used early in a drug’s development to investigate its permeability across the BBB [42,43,44]. Drug permeability is typically studied using cell cultures grown on permeable membranes in a transwell, where the transendothelial/epithelial electrical resistance (TEER) is used to measure the change in the electrical resistance across the cellular monolayer (Figure 1c) [37,45]. The TEER measurements evaluate the permeability and integrity of the endothelial cell layer [45,46], but the static cell culture cannot mimic the fluid characteristics of the BBB. More recently, microfluidic organ-on-chip in vitro models have been employed for the assessment of TJ barrier properties. Efforts are underway to develop BBB-on-a-chip that can recreate the human BBB using cells derived from the brain endothelial cells [47,48]. Additionally, in silico models of the BBB TJs have been developed to reliably compute the permeability of ions and small molecules across the paracellular barrier [49,50,51,52]. With advances in computer algorithms and hardware, strategies are being developed to design novel therapeutics that can cross the BBB. A concerted in vivo, in vitro, and in silico effort is underway to unravel the complexities of the BBB.

2. Tight Junctions: Role of Claudin-5 in the BBB

Claudin family of membrane proteins establishes TJs in endothelial and epithelial cells [53,54,55]. Claudins are expressed in vascular endothelial cells and all known epithelial cells throughout the body [56,57,58]. To date, 27 members of the claudin family are known to be functionally expressed in mammals [59,60]. They are classified into classic claudins (1–10, 14, 15, 17, and 19) and non-classic claudins (11–13, 16, 18, and 20–27) based on their sequence similarity and structural homology [59,61,62]. Most tissues express multiple members of the claudin family at the TJs [56,63,64,65,66,67]. Claudin proteins (~20–30 kDa) fold into a four-transmembrane helix bundle (TM1–4) with two extracellular loops (ECL1–2) and an intracellular loop (ICL; Figure 2) [60,62,68]. The lateral interactions mediated by TM and ECL domains of claudins within the membrane of the same cell are termed cis interactions. Cis interactions between multiple claudins occur in membranes of expressing cells. Subsequently, the claudins undergo trans assembly via head-on interactions of their ECL loops and form the macromolecular TJ assembly with adjacent claudin expressing cells. The TJ strands (often observable under freeze-fracture microscopy) are multi-component structures with claudins as the primary component alongside other TJ-associated proteins such as occludin, tricellulin, junctional adhesion molecules (JAMs), and zonula occludens (ZO-1,2) [53,60,61].
Claudin-5 is the key TJ protein at the BBB with an expression level that is much higher than claudin-1, -3, and -12 [69,70]. Knockout experiments have shown that claudin-5 is responsible for controlling the paracellular permeability of molecules up to ~800 Da [37,38,71]. In vitro biochemical analyses have suggested that claudin-5′s ECL2 loop strongly contributes to the trans interactions [72,73]. Mutations to conserved residues in ECL2 lead to a marked increase in paracellular permeability because of compromised size selectivity [59,74]. Studies also show that claudin-5 knockout mice die within 10 h of birth, although their TJs remain impermeable to molecules >800 Da in size [71]. In contrast a double knockdown study of both claudin-5 and occludin showed enhanced permeability of 3–10 kDa tracers indicating that occludin serves a structural role in the BBB TJs [75,76].

3. Need for Computational Modeling

In silico models complement in vitro and in vivo approaches used to investigate biological systems. Computer models augment our ability to analyze aspects of biological systems that are too complex to be investigated by conventional experiments alone [77,78,79,80]. Computational approaches provide a molecular-level basis for experimentally observed biomacromolecular assembly, selectivity, dynamics, and macroscopic properties.
For example, X-ray crystallography of a membrane protein provides essential structural information about its conformation at low temperatures [81,82]. However, proteins participate in dynamic processes and have thermal motions at physiological temperatures that cannot be resolved by current experimental techniques [83,84]. In silico methods can readily build upon the crystal structure data to provide information about molecular displacements, conformational dynamics, and reorientation frequencies of proteins [80]. The resulting models can also provide insights into the intramolecular and intermolecular interactions of a protein in its native environment, including other proteins, lipids, ions, and solvents. Another area where computational methods have been transformative is in predicting the three-dimensional structure of membrane proteins. The hydrophobic nature of membrane proteins makes them experimentally unwieldy for crystallization and X-ray structure determination [83,84,85]. In recent years, advances in protein sequencing have enabled the utilization of sequence co-evolution information to infer and predict the native fold(s) of a protein to a high-level of accuracy that can be validated by experimental biochemical assays [86,87,88,89,90].
Advances in computational hardware also have facilitated the study of complex biological systems. In the last decade, molecular simulations of biological systems with millions of interacting particles over microsecond to millisecond timescales have been reported [64,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116]. These simulations provide unprecedented insights into the system under investigation and bridge nanoscale data to macroscale experimental results.

4. Computational Toolkit

Computational methods, in the form of in-house computer programs, freely-distributed and commercial software suites, and online servers, offer an extensive toolkit that can be used to investigate different aspects of a biological system. Based on the desired qualitative or quantitative output, an array of computational models can be employed. Figure 3 summarizes a hierarchical computational approach relevant to the study of membrane proteins, using claudin-5 as an example.

4.1. Protein Structure Prediction

Crystal structures of murine and human claudins mClaudin-15 (4P79) [68], mClaudin-19 (3X29) [117], hClaudin-4 (5B2G) [118], and mClaudin-3 (6AKE) [119] are available. Most of the currently available claudin structures have missing loop and c-terminal domains, due to crystallographic artifacts. To study a claudin that has not yet been structurally resolved, the first step is to construct its homology model [120,121,122]. Reliable homology models can be derived if the target and template structures have at least 30% residue identity spanning the entire length of the domain of interest [120,123]. Since most claudins share ~30% sequence identity, this technique can reliably model structures of the other members of the claudin protein family [59,62,124]. Multiple tools are available online that can build homology models; [125,126] among them the I-TASSER [127,128,129] and RaptorX [130,131] servers are well established in the structural biology field (Figure 3). Both servers take the primary sequence of the target protein as input and give as output the lowest energy 3D-modeled structure of the protein. In building homology models care should be taken in modeling the nonconserved domains. For example, the flexible c-terminal domains of claudin family of proteins are challenging to model because they are intrinsically disordered and adopt nonspecific secondary structures [132]. Homology modeling methods are designed to provide a correct global structure but do not optimize side-chain orientations or flexible-loop conformations [121,132,133,134]. Such details can be obtained by loop modeling techniques such as ModLoop and FoldX [135,136]. Once the 3D-structure has been computationally determined, the next step is to embed the claudin in the desired lipid membrane and optimize the structure using molecular dynamics (MD) [137] simulation at atomistic resolution.

4.2. Protein Dynamics

Simulations of biological systems are performed at multiple resolutions—atomistic, coarse-grained, and mesoscale [137,138,139]. In atomistic resolution, the motion of each atom is evaluated as a function of time, whereas, in the coarse-grain (CG) resolution, a group of atoms in an amino acid is mapped to a bead to evaluate its dynamics. Mesoscale models take coarse graining up another notch; clusters of amino acids are mapped to a superbead, and simple harmonic oscillators represent the interaction between those superbeads [140,141,142]. Since each resolution has merits and pitfalls, multiscale modeling approaches are increasingly being used to take advantage of the benefits available at each resolution [143]. Multiscale simulations are vital to interpreting experimental results and guiding new lines of experimental investigation [80,144,145,146,147,148].
The classical atomistic MD approach is to simulate each atom of the system using pairwise interatomic potentials to estimate the interaction strength of each atom with its remaining atoms in the system. The interaction energy of each atom in a system is modeled using a molecular force field, which comprises of both bonded (bond stretch, angle bend, dihedral rotation, and torsional) and non-bonded (electrostatic and van der Waal’s) interaction terms (Figure 4). There are many force fields available in the literature for proteins, such as CHARMM [137], AMBER [149], GROMOS [150], and OPLS [151] among others. Each force field has a unique set of parameters that are optimized to reproduce experimental reference data. Therefore, selection of a force field is an important decision before starting an MD simulation for any system.
The all-atom simulations are restricted to small time-steps (1–4 fs) to capture the intra- and intermolecular motion of atoms over time. Atomistic simulations result in detailed structural, thermodynamic, and kinetic data that can be compared to experimental data. A typical biomolecular system comprises of millions of atoms, embedded in ionic media to mimic the physiological conditions, simulating biological processes over micro- or millisecond timescales. Although atomic- level molecular modeling provides specific interatomic interactions, its practical application is limited by computational resources.
Setting up membrane protein systems for molecular simulations is challenging when it comes to embedding a protein in multicomponent lipid membranes. To streamline the process, the recent development of publicly hosted interactive web-based servers (such as charmm-gui) [152] is enabling the generation of input files required for membrane protein simulations. Among the existing simulation packages are popular MD engines such as GROMACS [153,154], NAMD [155], AMBER [156,157], CHARMM [152,158,159], and visualization tools such as VMD [160], Pymol [161], and YASARA (Figure 3) [162,163,164].
The use of CG simulations has increased in the last 10 years with the development of several force fields for a variety of biomolecules. Popular CG force fields include MARTINI [165,166,167,168,169], SIRAH [170,171], PRIME [172,173], and PACE [174,175]. Among these, MARTINI is the most popular for membrane protein simulations [165,166,176]. The MARTINI force field maps ~4 nonhydrogenic atoms into one bead; for example, a glycine amino acid residue is mapped as one bead [165,166,168]. By mapping a small number of atoms into a bead, CG resolution reduces the number of particles in a simulation and allows the system to be sampled less frequently with a larger time-step (20–40 fs). This approach enables simulation timescales that are 2 to 3 orders of magnitude larger than atomistic simulations. Thus, the advantages of CG simulations have spurred interest in simulations of larger sizes (hundreds of nanometers) and longer timescales (tens of microseconds) [177,178,179,180]. Others and we have rigorously demonstrated the ability of MD simulations to capture residue-level precision in membrane protein interactions [91,94,97,108,110,115].

4.3. Protein–Protein Interactions

4.3.1. Self-Assembly Simulations

Molecular self-assembly simulation is a powerful technique that mimics interactions among molecules over time [94,98,145,181]. These simulations capture the spontaneous organization of molecules dictated by intra- and intermolecular nonbonded interactions (van der Waals forces, electrostatics, hydrogen bonds, hydrophobic, and hydrophilic interactions) under near-equilibrium conditions. Self-assembly simulations are routinely used to study complex membrane protein–protein interactions. These simulations, however, have several limitations including high computational cost. Membrane protein assembly is particularly computationally expensive in part because protein diffusion requires much longer equilibration time than the relaxation of membrane lipids. Moreover, the simulations need to be performed over a sufficiently long time to obtain reliable data, although ensuring convergence of such simulations is a challenge. There is no guarantee that a 10–100 µs timescale is long enough for proteins to assemble [182,183]. Employing multiscale simulations mitigates the challenges of long timescale simulation.

4.3.2. Molecular Docking

Molecular docking, another popular method of studying protein associations, is a highly effective method to examine and rank interactions based on rigorous algorithms that evaluate empirical energetics data. Docking approaches are simple geometry-based methods where tens of thousands of protein–protein interfaces are generated via geometric rigid body motions followed by a scoring scheme where interfaces are ranked based on the energetics [184,185]. The energy function that is used in the scoring algorithm is often empirically modeled or derived from a knowledge-based method. Docking approaches tend to be fast because generating geometries of proteins is a simple task. Depending upon the type of protein–protein interaction, docking methods can give realistic results. Popular docking methods include Cluspro [186,187], Rosetta [188], HADDOCK [189,190], and ZDOCK [191] among others [192,193,194,195]. Cluspro has robust reproducibility, and it is the current state of the art in protein–protein docking. Rosetta is a more advanced software suite of protein modeling and analysis tools. The Rosetta docking protocol has been optimized for membrane proteins, and one can define membrane domain boundaries to guide the docking. However, despite the ease and speed of docking methods, static docking is not optimal for macromolecular assemblies that involve multiple proteins [196]. Moreover, for membrane proteins, the presence of the lipid environment significantly influences proteins’ association patterns, but none of the currently available docking methods offers a realistic lipid environment for protein docking [197,198,199,200,201,202].

4.3.3. Protein Association Energy Landscape

An alternate strategy to study membrane protein association is to employ ensemble-based sampling methods. These methods use hundreds to a few thousand starting structures that each run independently to yield an unbiased sampling of the protein–protein association landscape. One such ensemble-based method, docking assay for transmembrane (DAFT) components [203] uses the population density of the resulting protein association conformations as a proxy for stability. Other ensemble-based methods use amino acid contact-based stability as a metric [204,205]; however, none of the available methods provides the association energies needed to determine the relative stability of protein–protein conformations obtained in the ensemble.
In our work, we used a combined stochastic sampling and MD approach to sample pairwise protein–protein interactions [206]. Taking advantage of in-plane diffusion of membrane proteins relative to each other in the lipid bilayer, we defined a two-dimensional rotational space of the interacting proteins. We adopted a stochastic sampling technique to generate thousands of seed geometries that uniformly sample the rotational space around a pair of interacting proteins. These seed geometries were then propagated independently to equilibrium using standard MD. The interaction energy associated with each conformation visited in the dynamics was computed to generate the protein association energy landscape (PANEL) [206]. By placing the proteins within the interacting distance of the seed geometries, we eliminated the need for protein diffusion, thereby reducing the simulation times for each system and boosting the PANEL performance by two orders of magnitude over classic self-assembly simulations. The PANEL method also provides an absolute measure of interaction energy for any pair of interacting proteins in dimeric or even in trimeric association [206]. We used the PANEL method to report the landscape of claudin-5 dimers (Figure 3), corroborating our previous self-assembly simulation results [115]. Furthermore, the robustness, computational affordability, and user-friendly implementation of the PANEL method makes it attractive for application to any membrane protein.

4.4. Pore Characterization and Paracellular Transport

The molecular architecture of the TJ pores and barriers are the determinants of the charge and size-selective transport of ions and molecules. Analysis of the pore’s structural features is important to understand TJ selectivity. The structural characteristics of TJ pores (such as length, diameter, volume, and cross-sectional area) provide valuable insight into the size selectivity of the pore. Additionally, identification of pore-lining residues in the channel walls and the bottleneck regions of the pore offer significant clues about their charge selectivity. Several software tools, such as CAVER [207], MOLE [208], HOLE [209], MOLAXIS [210,211], and POREWALKER [212], (each based on slightly different algorithms) provide structural analysis of protein pores (Figure 3).
Paracellular transport involves calculation of the free energy profiles of the translocation of a substrate through the pore. This calculated energy directly correlates to the permeability of the substrate across the pore [213,214]. The free energy profile can be obtained from the advanced sampling techniques of umbrella sampling [215] or metadynamics [216,217,218]. In the umbrella sampling technique, where multiple frames or conformations of a solute are generated along the length of the pore, an external force-constant is applied on the solute to restrain its position. Based on the amount of force applied, the free energy profile along the transport path can be evaluated [215,219]. An alternative method to determine the free energy profile is the recently introduced metadynamics approach where the system is biased gradually such that the ion translocation is made feasible within the affordable timescale of MD simulation. Metadynamics is available in GROMACS through the PLUMED plugin [220].
A simplistic equilibrium MD simulation can also be set up with different ion concentrations to observe the number of ion translocation events that occur during the simulation. This method may provide insights about preferred ion transport under unbiased, equilibrium conditions; however, the limitation of the method is the timesale in which an ion transport event will occur may not be within the simulaion time frame. To accelerate the ion transport event, computational electrophysiology has emerged as a method that involves biasing the molecular system by creating localized charge imbalances in system, which induces ion transport [221,222]. The small charge imbalances provide biologically realistic electrochemical gradients across the pore and facilitate ion transport.

5. Computational Nanoscopy of Claudin-5 in BBB Tight Junction Architecture

Claudin proteins embedded in the lipid bilayer undergo cis assembly via their transmembrane helical domains, or by the extracellular loops, or by both [94,115]. In general, protein–protein cis interactions are mediated by a multicomponent mixture of saturated and unsaturated lipids with variety of charged head groups and lipid tail lengths [91,106,201,202,223]. The helical bundle of membrane proteins has a rigid tertiary structure that provides stability to the protein to reside in the membrane environment. This rigidity in turn affects the local lipid environment and manifests as a hydrophobic mismatch between the protein and the lipids [198,224,225]. Hydrophobic mismatch is energetically unfavorable and serves as a driving force for the embedded protein to find interactions that would minimize its hydrophobic mismatch with the lipids and undergo cis assembly [103,226,227,228,229,230].
There are multiple aspects of claudin cis and trans assembly that influence the TJ architecture. The cis-assembled claudins can provide the template for the trans interactions [97,231,232]; therefore, the study of cis interfaces is vital to understand the TJ architecture. Besides, claudins are reported to be post-translationally modified (PTM) with phosphorylation and palmitoylation. Of these PTMs, studies suggest that palmitoylation alters recruitment of claudin proteins in the TJs and TJ permeability [91,233,234]. Since most aspartate–histidine–histidine–cysteine (DHHC) enzymes that assist in palmitoylation [235] are thought to be localized to the endoplasmic reticulum (ER) [198,236], it is likely that claudin cis interactions that are influenced by palmitoylation are localized in the ER. However, DHHC enzymes have been reported in other organelles and in the plasma membrane [237], which makes both the location of claudin-5 palmitoylation in the secretory pathway and the cis dimerization ambiguous. It is not surprising that the complexity of the BBB TJ architecture has been a deterrent in unraveling the mystery of the size and charge selectivity of the TJ. A list of basic questions that have confounded experts in the TJ research field include:
(1)
How do claudin proteins form contiguous strands? Do the claudin proteins have preferred cis conformations?
(2)
Where along the secretory pathway do claudin-5 proteins undergo cis assembly?
(3)
How do posttranslational modifications of claudin-5 influence the cis interactions?
(4)
How are TJ pores and barriers formed?
(5)
How do BBB TJs get their size and change selectivity?
To provide a new perspective to these long-standing challenges, we developed tools for in silico investigation of the TJ claudin-5 proteins. In a series of publications over the past few years [91,94,97,115,206], we have reported on the various aspects of claudin-5 BBB TJ architecture. Here, we reviewed the outcomes of these five lines of inquiry (Figure 5).

5.1. Claudin-5 Forms Symmetric and Asymmetric Dimer Conformations that Form Contiguous Strands

The cis interactions among claudins, determined by computational methods, corroborate results of biochemical assays [115]. The claudin-5 cis interactions were studied using both the self-assembly and the PANEL approaches [115,206].
The self-assembly simulations, performed in CG resolution, suggest that claudin-5 monomers embedded in a simple mixture of saturated and unsaturated lipids readily form dimers (within 400–500 ns) [115]. These dimers subsequently (in 5–10 µs) form contiguous strands and high-order assemblies (Figure 6a). For example, in a system composed of 64 claudin-5 monomers, a contiguous strand of 36 monomers was observed with a ring-like pentamer (Figure 6b), which match the size description of ~10 nm particles observed under freeze fracture [238,239,240]. The claudin strands present at the end of the 10 µs simulation were extracted and reverse mapped to atomistic resolution for conformational and structural analysis (Figure 6c). A pair-wise conformational analysis of the self-assembled claudins revealed four, most prevalent conformations (labeled dimer A−D) that concurrently exist to form the strands. Several repeat simulations have been performed that all show formation of dimers A−D although in varying populations. Since it is practically impossible to achieve convergence in self-assembly simulations due to finite timescales, we were unable to provide converged relative populations of each of the four dimers. However, these simulations indicate that a dimer is the smallest stable unit of the self-assembled strand and that the cis assembly can occur in the absence of trans interacting partners from the adjacent cell. In an experimental study, claudin-2 cis dimerization was also observed independent of cell-surface expression and intercellular cross-linking [241].
To delve deeper into the stability of the claudin-5 dimers, we developed the PANEL approach [206]. The PANEL plot for claudin-5 provided a comprehensive sampling of dimer conformations and a clear visualization of the entire energy landscape as a function of the orientation angles (Figure 7). The claudin-5 dimer interaction energies range from 0 (least stable) to approximately −1500 kJ mol−1 (most stable) [206]. The PANEL contour plot shows multiple low energy conformations. Conformations along the diagonal (Figure 7b) represent symmetric dimers, while those off the diagonal are asymmetric dimers. It was reassuring that the key symmetric Dimers B, C, and D, also observed in self-assembly and in experiments, are shown in the minimum energy PANEL plot. Dimer B is a TM3:TM3-mediated dimer, while Dimer C is mediated by TM4:TM4, and Dimer D, the symmetric ECL1-mediated interaction, also proposed by Suzuki et al. [242]. There are also several asymmetric dimers (Dimer A), which are marked as Dimers A1, A2, and A3 in the PANEL plot. Dimer A1 corresponds to the linear arrangement observed in the claudin-15 crystal structure [242]. The relative orientation of the claudins in each of the Dimers A to D (Figure 7b) is described in our earlier work [115,206]. Although the presence of claudin dimers and higher order complexes has been reported in several experimental studies [73,240,241,243,244], isolating intact dimers or higher order complexes still remains a challenge. Use of computational approach provides a direct look into the structural details of these claudin interfaces.

5.2. Primary Cis Interactions Occur Early in the Secretory Pathway

Like other TJ proteins, claudins are secreted in the endoplasmic reticulum (ER) and reach the plasma membrane via the Golgi complex [231] Along this secretory pathway, the membrane composition changes drastically [245,246,247] To mimic the complexity of cell organelles encountered by a typical membrane protein during its secretory pathway, we performed claudin-5 self-assembly in three model membrane compositions: ER, cholesterol-enriched endoplasmic reticulum (ERc), and plasma membrane (PM) [91]. The ER membrane model constituted a binary mixture of saturated and unsaturated lipids in 1:1 ratio; the ERc was a 2:2:1 ternary mixture of saturated lipid, unsaturated lipid, and cholesterol; and the PM model comprised ten lipid components including gangliosides, sphingomyelin, and cholesterol [91].
Our findings indicate that claudin-5 monomers are more likely to form dimers in simple membranes like the ER where saturated and unsaturated phospholipids are the principal components. As claudins move to the Golgi complex, trans-Golgi vesicles, and the plasma membrane, the higher cholesterol concentrations in these membranes limit claudin-5 assembly. We observed that an increase in membrane complexity was negatively correlated with the length of claudin-5 strands (Figure 8a) [91].
A comparison of PANEL plots for claudin-5 dimers in a simple unsaturated di-oleoyl-phosphatidylcholine (DOPC) lipid membrane and ERc membrane shows the impact of the lipid composition on dimer stability (Figure 8a). The effect is quantified by taking the difference in interaction energies (ΔPANEL) between corresponding geometries (θ, θ’) on the two plots. A vast majority of the dimer interactions were destabilized in ERc membrane, although a few conformations did show more favorable interactions. In particular, Dimer B conformation (mediated via TM3:TM3 interaction) is stabilized. Overall, our results suggest that dimerization is more favorable in conditions existing in the early stages of the secretory pathway.

5.3. Posttranslational Modifications Alter Relative Stability of cis Dimers

Several studies have indicated the significance of palmitoylation in the recruitment of specific claudins to TJ and the consequent effect on the TJ permeability [233,234]. Since palmitoylation and depalmitoylation are highly labile reactions, the evaluation of these modifications on claudin assembly is experimentally challenging [237]. However, in silico methods provide the necessary spatial and temporal resolutions required to study these posttranslational modifications.
The presence of palmitoyl chains at the protein–lipid interface introduces variations in interfacial properties, effecting both the extent of dimerization as well as the distribution of specific dimer interfaces. Our results show that palmitoylation increases cholesterol and saturated lipid accumulation around the claudins and enhances cis interactions to the Dimer B interface while disrupting Dimer C largely due to steric factors (Figure 9a,b) [91]. The impact of palmitoylation on claudin-5 cis dimers is evident from the ΔPANEL landscape plotted as the difference between palmitoylated claudin (claudin-5P) and non-palmitoylated claudin-5 PANEL plots (Figure 9c). The ΔPANEL plot shows the stabilization of Dimer B and the destabilization of Dimer C regions. The change in dimer stability upon palmitoylation indicates that posttranslational modification can be a potential strategy to control paracellular permeability [91,97,233,234,237,248,249].

5.4. Cis Dimers Act as Precursors to Form Trans Interfaces

Trans interactions are formed when claudins on adjacent cells interact head-on to form the macromolecular TJ assembly. We adopted both CG- and atomistic-level simulations to decipher the nature of trans interfaces. Numerous challenges were posed both on conceptual and computational resource fronts while designing the trans assembly simulations.
To mimic paracellular space, our initial studies of trans interactions used two membrane patches with multiple (64) equidistant claudin-5 proteins [97]. They were placed head-on, such that the ECL loops from either membrane could interact within the gap between the two membranes (Figure 10). Sampling 20 µs of trans self-assembly took about 12 days on a supercomputer using 20 high performance computational nodes. Despite the intensive computation, the simulations represented only a small membrane patch that contained 128 claudin monomers across the two membranes. In contrast to a single membrane simulation, trans assembly simulations were unable to capture a few cis dimeric interactions that were prevalent in the single membrane self-assembly simulations. The simulations were intentionally designed so that a diffusing claudin-5 monomer would encounter a trans interacting partner before encountering a cis interacting partner; this ensured generation of more trans interaction variants. The absence of cis dimers rendered this a monomer-driven, trans interaction assembly [97]. Multiple pore-like assemblies were observed in these simulations showing ECL contacts stabilized by hydrogen-bonds and salt bridges. These trans pairs also strongly resembled pore models suggested in previous studies [72,73,74,250]. However, this model represented a very low fraction of the actual trans interactions, which span the entire perimeter of a cell and involve tens of millions of protein monomers.
From the insights obtained from the cis and trans self-assembly simulations a dimer-guided trans interaction study was carried out using molecular docking approaches. This method yielded highly promising pore-forming trans interfaces. Moreover, the formation of cis dimers before trans interactions has been suggested in several findings [231,241,251,252]. Based on our work, we suggest that the cis dimers act as precursors for forming trans interfaces [91,97,115]. The fact that cis assembly simulation also self-assembles into a dimer in shorter timescales, lends support to that cis interactions may act as the template to guide the trans assembly.
The findings from these simulations were able to fill significant gaps in our understanding of the structure–function relationships between claudin monomers. Although several questions still need to be answered, and further investigations need to be performed to achieve a comprehensive understanding of the TJ architecture, our current understanding has been made possible by the an array of techniques in the computational toolkit.

5.5. Sub-Nanometer Wide Claudin-5 TJ Pores, Lined by Charged Residue Determine Pore Selectivity

For claudin-5, two pore models, mediated by different cis dimers, have been reported based on trans interaction simulations and mutation analyses of cis dimers [97]. The pore I model is a face-to-face stacking of the dimer D interface. This pore resembles the proposed claudin-15 pore [242]. The pore II model was based on the Dimer B interactions observed in claudin-5. These models were generated using docking methods in combination with symmetric refinement to arrive at the pore structures based on Dimers D and B (Figure 11). Fundamental pore features, such as diameter, length of the pore, pore bottleneck diameter, and the pore lining residues, were characterized.
Pore characteristics provide a general idea of the nature of the pore. However, evaluating the energetics of molecular transport along these pores is essential to further enhance our understanding of paracellular selectivity of the BBB’s TJs. Using advanced sampling methods, translocation of water and glucose was reported across claudin-5 pore model formed by Dimer B. These simulations showed that the claudin-5 Dimer B pore allows water to translocate without energy barriers but had a higher energy barrier for glucose translocation compared with the GLUT1 glucose transporter [97,253,254]. Further extension of this work revealed that claudin-2 is selective to sodium ions with multiple observations of ion translocation under equilibrium conditions. Similar approaches are gaining attention within the field (e.g., claudin-15 pores probed using MD simulations) [255,256,257]. The more comprehensive undertaking of looking at ion transport in pores formed by different claudins is still in its infancy, but several computational groups are developing the essential methods needed to pursue this line of inquiry. In the meantime, the proposed pore structures present opportunities for further refinement to obtain better insights about their structural as well as functional aspects.

6. Additional Factors that Influence Tight Junctions

Computational methods can efficiently study cis and trans interactions of claudins at the molecular level. Claudins and other TJ proteins exist in a chemical environment that varies in type and amount of proteins and lipids, hydrophobicity, protein lipidation states, and other aspects [63,200,225,231,246,247]. These physiological factors contribute to the final TJ architecture. Moreover, TJs are known to be dynamic entities that undergo constant structural change [258,259]. Thus, to understand paracellular transport in its entirety, other factors need to be considered. A few key aspects important to TJ formation are discussed below.

6.1. Higher Order Claudin Aggregates

Several studies have reported the presence of claudin aggregates (in addition to strands) [61,118], suggesting that higher-order aggregates such as trimers, tetramers, or hexamers may be formed. Clostridium perfringens enterotoxin, which causes food poisoning, has been reported to form higher-order aggregates upon interacting with claudins. We proposed that the trimeric form of the enterotoxin could potentially bind to the trimeric form of claudins in the intestinal epithelium [94,117]. In another recent report of claudin-11 TJs, the carbon replica method revealed an interesting view of consistent double-stranded structures [260]. These reports emphasize that in order to fully understand the architecture of the TJ strands, it is important to consider the higher-order interfaces in addition to dimers.

6.2. Involvement of Other TJ Proteins: Key to TJ Structural Stability

In addition to claudins, proteins from approximately 80 other families are known to be involved in the TJ protein complex assembly including occludins and JAMs [53,54,60]. Occludins and JAMs contribute to the structural stability of TJs. Although these proteins are not known to affect directly TJs’ selective permeability or barrier function, their role in maintaining the structural integrity of the TJ strongly suggests they form interfaces with claudins in order to support the strand. This may lead to masking of some claudins’ interfacing surfaces while binding to others, thereby redistributing the relative populations of dimer interfaces. The study of occludins and JAMs is crucial to provide a complete perspective of TJs due to their significant presence within the TJ assembly.

6.3. Cytoskeletal Involvement

Claudins and other TJ membrane proteins are dynamically regulated by scaffolding proteins (zonula occludens-1/2/3) and actin cytoskeleton [53,231]. Although studying claudins individually provides insights about their interface forming behavior, thereby helping us to predict the architecture and nature of the corresponding TJ strands, continuous regulation by the cytoskeleton may induce several factors such as localized redistribution of claudin populations, or regulation of biophysical properties such as protein tilting, membrane anchoring, and so on. Since the cytoskeletal regulation tends to be highly dynamic and specific to the needs of the cell, deciphering their roles remains one of the greatest challenges.

7. Future Work

Future work on the BBB research will be directed towards the development of strategies to enable and enhance drug delivery to the brain. A fundamental understanding of the molecular construct that guards the blood–brain interface has been obtained through computational tools, which have enabled the design of in silico models to explore and the understanding of the architecture of this interface. One future research goal is to identify small molecules that can act as TJ modulators for use in screening claudin-5 trans interactions, reducing specificity, increasing pore size, and enhancing pore permeability in a controlled and reversible manner. Having the molecular-level pore structure is foundational to this endeavor. We envisioned that in vitro studies would be designed in the future to validate the predictions of the computational models. Other future work includes pursuing an integrated computational and in vitro/in vivo study of multiple drug pathways for neurotherapeutics, evaluating molecular drug interactions, obtaining the dynamics of transport of potential drug molecules with the vision of quantifying drug permeability through the BBB.

Author Contributions

N.R., F.J.I., and S.N.; conceptualization, S.N., N.R., and F.J.I.; methodology, S.N., N.R., and F.J.I.; software, N.R. and F.J.I.; validation, N.R. and F.J.I.; formal analysis, N.R. and F.J.I.; investigation, N.R. and F.J.I.; data curation, N.R. and F.J.I.; writing—original draft preparation, N.R. and F.J.I.; writing—review and editing, S.N., N.R., and F.J.I.; visualization, S.N., N.R., and F.J.I.; supervision, S.N.; project administration, S.N.; funding acquisition, S.N.

Funding

This research was funded by the CAREER CBET-1453312 grant from the National Science Foundation. Computational resources for the work were provided by the following groups—Information and Technology Services at Syracuse University; Anton 2, provided by the Pittsburgh Supercomputing Center (PSC) through Grant R01GM116961 from the National Institutes of Health (NIH); Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1053575.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

BBBBlood–brain barrier
CGCoarse grained
Claudin-5PPalmitoylated claudin
DAFTDocking assay for transmembrane
DHHCAspartate-Histidine-Histidine-Cysteine
DOPC1,2-dioleoyl-sn-glycero-3-phosphocholine
ECLExtracellular loop
EREndoplasmic Reticulum
ERcCholesterol enriched Endoplasmic Reticulum
ICLIntracellular loop
JAMSJunctional adhesion molecules
MDMolecular dynamics
PANELProtein Association Energy Landscape
TEERTransendothelial/epithelial electrical resistance
TJTight junction
TMTransmembrane loop
ZOZonula Occuldens

References

  1. Abbott, N.J. Comparative Physiology of the Blood-Brain Barrier. In Physiology and Pharmacology of the Blood-Brain Barrier; Bradbury, M.W.B., Ed.; Springer: Berlin/Heidelberg, Germany, 1992; pp. 371–396. [Google Scholar]
  2. Bauer, H.C.; Krizbai, I.A.; Bauer, H.; Traweger, A. “You Shall Not Pass”-tight junctions of the blood brain barrier. Front. Neurosci. 2014, 8, 392. [Google Scholar] [CrossRef] [PubMed]
  3. Begley, D.J.; Brightman, M.W. Structural and functional aspects of the blood-brain barrier. In Peptide Transport and Delivery into the Central Nervous System; Prokai, L., Prokai-Tatrai, K., Eds.; Birkhäuser Basel: Basel, Switzerland, 2003; pp. 39–78. [Google Scholar]
  4. Daneman, R.; Prat, A. The blood-brain barrier. Cold Spring Harb. Perspect. Biol. 2015, 7, a020412. [Google Scholar] [CrossRef] [PubMed]
  5. Saunders, N.R.; Dreifuss, J.-J.; Dziegielewska, K.M.; Johansson, P.A.; Habgood, M.D.; Møllgård, K.; Bauer, H.-C. The rights and wrongs of blood-brain barrier permeability studies: A walk through 100 years of history. Front. Neurosci. 2014, 8, 404. [Google Scholar] [CrossRef] [PubMed]
  6. Lewandowsky, M. Zur Lehre von der Cerebrospinalflüssgkeit. Z. Clin. Med. 1900, 40, 480–494. [Google Scholar]
  7. Goldmann, E.E. Die äussere und innere Sekretion des gesunden und kranken Organismus im Lichte der ‘vitalen Färbung. Beiträg Klin. Chir. 1909, 64, 192–265. [Google Scholar]
  8. Goldmann, E.E. Vitalfärbung am Zentralnervensyatem. Beitrag zur Physio-Pathologie des plexus chorioideus und der Hirnhäute. Abh. Preuss Akad. Wiss. Phys. Math. Kl. 1913, 1, 1–60. [Google Scholar]
  9. Ehrlich, P. The Relations Existing between Chemical Constitution, Distribution and Pharmacological Action; Collected Studies on Immunity New York: John Wiley & Sons, translated by C Bolduana from Ch XXXIV of Gesammelte Arbeiten zur Immunitätsforschung; Ehrlich, P., Ed.; Hirschwald: Berlin, Germany, 1906. [Google Scholar]
  10. Ehrlich, P. Das Sauerstoffbedürfnis des Organismus. In Eine Farbenanalytische Studie; Hirschwald: Berlin, Germany, 1885. [Google Scholar]
  11. Stern, L.; Gautier, R. II.—Les Rapports Entre Le Liquide Céphalo-Rachidien Et Les éléments Nerveux De L’axe Cerebrospinal. Arch. Int. Physiol. 1922, 17, 391–448. [Google Scholar] [CrossRef]
  12. Stern, L.; Gautier, R. III.—Rapports Entre Le Liquide Céphalo-Rachidien Des Espaces Ventriculaires Et Celui Des Espaces Sous-Arachnoïdiens. Arch. Int. Physiol. 1923, 20, 403–436. [Google Scholar] [CrossRef]
  13. Stern, L.; Gautier, R. Passage simultané des substances dans le liquide céphalo-rachidien et dans les centres nerveux. Rcr D. Ia Soc. De Phys. Et D’hist. Natur. De Genève 1918, 35, 58–60. [Google Scholar]
  14. Stern, L.; Gautier, R. Le passage dans le liquide céphalo-rachidien de substances introduites dans la circulation et leur action sur le système nerveux central chez les différentes espèces animales. Rcr D. Ia Soc. De Phys. Et D’hist. Natur. De Genève 1918, 35, 91–94. [Google Scholar]
  15. Stern, L.; Gautier, R. Recherches sur le liquide céphalo-rachidien. 1. Les rapports entre le liquide céphalo-rachidien et la circulation sanguine. Arch. Int. Physiol. 1921, 17, 138–192. [Google Scholar] [CrossRef]
  16. Roux, E.; Borrel, A. Tétanos cérébral et immunité contre le tétanus. Ann. Inst. Pasteur 1898, 12, 225–239. [Google Scholar]
  17. Bouffard, G. Injection des couleurs de benzidine aux animaux normaux. Ann. D. L’lnst. Pasteur. Paris 1906, 20, 539–548. [Google Scholar]
  18. Weed, L.H. An anatomical consideration of the cerebro-spinal fluid. Anat. Res. 1917, 461–496. [Google Scholar] [CrossRef]
  19. Weed, L.H. The development of the cerebrospinal fluid spaces in pig and in man. Contrib. Embryol. Carnegie Inst. 1917, 5, 3–116. [Google Scholar]
  20. Wislocki, G.B. Experimental studies on fetal absorption. I. The vitality stained fetus. Embryol. Carnegie Inst. 1920, 11, 45–60. [Google Scholar]
  21. Cohen, H.; Davies, S. The Morphology and Permeability of the Roof of the Fourth Ventricle in some Mammalian Embryos. J. Anat. 1938, 72, 430–455. [Google Scholar]
  22. Friedemann, U. Blood-brain barrier. Physiol. Rev. 1942, 22, 125–145. [Google Scholar] [CrossRef]
  23. Wolburg, H.; Lippoldt, A. Tight junctions of the blood–brain barrier: Development, composition and regulation. Vasc. Pharmacol. 2002, 38, 323–337. [Google Scholar] [CrossRef]
  24. Duelli, R.; Kuschinsky, W. Brain Glucose Transporters: Relationship to Local Energy Demand. Physiology 2001, 16, 71–76. [Google Scholar] [CrossRef]
  25. Howarth, C.; Gleeson, P.; Attwell, D. Updated energy budgets for neural computation in the neocortex and cerebellum. J. Cereb. Blood Flow Metab. 2012, 32, 1222–1232. [Google Scholar] [CrossRef] [PubMed]
  26. Harris, J.; Jolivet, R.; Attwell, D. Synaptic Energy Use and Supply. Neuron 2012, 75, 762–777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Mergenthaler, P.; Lindauer, U.; Dienel, G.A.; Meisel, A. Sugar for the brain: The role of glucose in physiological and pathological brain function. Trends Neurosci. 2013, 36, 587–597. [Google Scholar] [CrossRef] [PubMed]
  28. Patching, S.G. Glucose Transporters at the Blood-Brain Barrier: Function, Regulation and Gateways for Drug Delivery. Mol. Neurobiol. 2017, 54, 1046–1077. [Google Scholar] [CrossRef]
  29. Abbott, N.J. Inflammatory Mediators and Modulation of Blood–Brain Barrier Permeability. Cell. Mol. Neurobiol. 2000, 20, 131–147. [Google Scholar] [CrossRef]
  30. Huber, J.D.; Egleton, R.D.; Davis, T.P. Molecular physiology and pathophysiology of tight junctions in the blood–brain barrier. Trends Neurosci. 2001, 24, 719–725. [Google Scholar] [CrossRef]
  31. Pardridge, W.M. Drug and Gene Delivery to the Brain: The Vascular Route. Neuron 2002, 36, 555–558. [Google Scholar] [CrossRef]
  32. Ballabh, P.; Braun, A.; Nedergaard, M. The blood–brain barrier: An overview: Structure, regulation, and clinical implications. Neurobiol. Dis. 2004, 16, 1–13. [Google Scholar] [CrossRef]
  33. Hawkins, B.T.; Davis, T.P. The Blood-Brain Barrier/Neurovascular Unit in Health and Disease. Pharmacol. Rev. 2005, 57, 173. [Google Scholar] [CrossRef]
  34. Kniesel, U.; Wolburg, H. Tight Junctions of the Blood–Brain Barrier. Cell. Mol. Neurobiol. 2000, 20, 57–76. [Google Scholar] [CrossRef]
  35. Luissint, A.-C.; Artus, C.; Glacial, F.; Ganeshamoorthy, K.; Couraud, P.-O. Tight junctions at the blood brain barrier: Physiological architecture and disease-associated dysregulation. Fluids Barriers CNS 2012, 9, 23. [Google Scholar] [CrossRef] [PubMed]
  36. Greene, C.; Campbell, M. Tight junction modulation of the blood brain barrier: CNS delivery of small molecules. Tissue Barriers 2016, 4, e1138017. [Google Scholar] [CrossRef] [PubMed]
  37. Butt, A.M.; Jones, H.C.; Abbott, N.J. Electrical resistance across the blood-brain barrier in anaesthetized rats: A developmental study. J. Physiol. 1990, 429, 47–62. [Google Scholar] [CrossRef] [PubMed]
  38. Abbott, N.J.; Patabendige, A.A.K.; Dolman, D.E.M.; Yusof, S.R.; Begley, D.J. Structure and function of the blood–brain barrier. Neurobiol. Dis. 2010, 37, 13–25. [Google Scholar] [CrossRef] [PubMed]
  39. Pardridge, W.M. The blood-brain barrier: Bottleneck in brain drug development. NeuroRX 2005, 2, 3–14. [Google Scholar] [CrossRef]
  40. Pardridge, W.M. Blood–brain barrier delivery. Drug Discov. Today 2007, 12, 54–61. [Google Scholar] [CrossRef]
  41. Pardridge, W.M. Drug Transport across the Blood–Brain Barrier. J. Cereb. Blood Flow Metab. 2012, 32, 1959–1972. [Google Scholar] [CrossRef]
  42. Cucullo, L.; Bennani-Baiti, B.; Rapp, E.; Janigro, D. Drug delivery and in vitro models of the blood–brain barrier. Curr. Opin. Drug Discov. Dev. 2005, 8, 89–99. [Google Scholar]
  43. Wilhelm, I.; Krizbai, I. In Vitro Models of the Blood-Brain Barrier for the Study of Drug Delivery to the Brain. Mol. Pharm. 2014, 11, 1949–1963. [Google Scholar] [CrossRef]
  44. Heymans, M.; Sevin, E.; Gosselet, F.; Lundquist, S.; Culot, M. Mimicking brain tissue binding in an in vitro model of the blood-brain barrier illustrates differences between in vitro and in vivo methods for assessing the rate of brain penetration. Eur. J. Pharm. Biopharm. 2018, 127, 453–461. [Google Scholar] [CrossRef]
  45. Srinivasan, B.; Kolli, A.R.; Esch, M.B.; Abaci, H.E.; Shuler, M.L.; Hickman, J.J. TEER Measurement Techniques for In Vitro Barrier Model Systems. J. Lab. Autom. 2015, 20, 107–126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Santaguida, S.; Janigro, D.; Hossain, M.; Oby, E.; Rapp, E.; Cucullo, L. Side by side comparison between dynamic versus static models of blood–brain barrier in vitro: A permeability study. Brain Res. 2006, 1109, 1–13. [Google Scholar] [CrossRef] [PubMed]
  47. Achyuta, A.K.H.; Conway, A.J.; Crouse, R.B.; Bannister, E.C.; Lee, R.N.; Katnik, C.P.; Behensky, A.A.; Cuevas, J.; Sundaram, S.S. A modular approach to create a neurovascular unit-on-a-chip. Lab. A Chip 2013, 13, 542–553. [Google Scholar] [CrossRef] [PubMed]
  48. Wang, Y.I.; Abaci, H.E.; Shuler, M.L. Microfluidic blood-brain barrier model provides in vivo-like barrier properties for drug permeability screening. Biotechnol. Bioeng. 2017, 114, 184–194. [Google Scholar] [CrossRef] [PubMed]
  49. Bennion, B.J.; Be, N.A.; McNerney, M.W.; Lao, V.; Carlson, E.M.; Valdez, C.A.; Malfatti, M.A.; Enright, H.A.; Nguyen, T.H.; Lightstone, F.C.; et al. Predicting a Drug’s Membrane Permeability: A Computational Model Validated With in Vitro Permeability Assay Data. J. Phys. Chem. B 2017, 121, 5228–5237. [Google Scholar] [CrossRef] [PubMed]
  50. Dickson, C.J.; Hornak, V.; Pearlstein, R.A.; Duca, J.S. Structure–Kinetic Relationships of Passive Membrane Permeation from Multiscale Modeling. J. Am. Chem. Soc. 2017, 139, 442–452. [Google Scholar] [CrossRef]
  51. Carpenter, T.S.; Kirshner, D.A.; Lau, E.Y.; Wong, S.E.; Nilmeier, J.P.; Lightstone, F.C. A method to predict blood-brain barrier permeability of drug-like compounds using molecular dynamics simulations. Biophys. J. 2014, 107, 630–641. [Google Scholar] [CrossRef]
  52. Geldenhuys, W.J.; Mohammad, A.S.; Adkins, C.E.; Lockman, P.R. Molecular determinants of blood-brain barrier permeation. Ther. Deliv. 2015, 6, 961–971. [Google Scholar] [CrossRef]
  53. González-Mariscal, L.; Betanzos, A.; Nava, P.; Jaramillo, B.E. Tight junction proteins. Prog. Biophys. Mol. Biol. 2003, 81, 1–44. [Google Scholar] [CrossRef]
  54. Tsukita, S.; Furuse, M.; Itoh, M. Multifunctional strands in tight junctions. Nat. Rev. Mol. Cell Biol. 2001, 2, 285–293. [Google Scholar] [CrossRef]
  55. Morita, K.; Furuse, M.; Fujimoto, K.; Tsukita, S. Claudin multigene family encoding four-transmembrane domain protein components of tight junction strands. Proc. Natl. Acad. Sci. USA 1999, 96, 511. [Google Scholar] [CrossRef] [PubMed]
  56. Schulzke, J.D.; Ploeger, S.; Amasheh, M.; Fromm, A.; Zeissig, S.; Troeger, H.; Richter, J.; Bojarski, C.; Schumann, M.; Fromm, M. Epithelial Tight Junctions in Intestinal Inflammation. Ann. N. Y. Acad. Sci. 2009, 1165, 294–300. [Google Scholar] [CrossRef] [PubMed]
  57. Morita, K.; Sasaki, H.; Furuse, M.; Tsukita, S. Endothelial Claudin. J. Cell Biol. 1999, 147, 185. [Google Scholar] [CrossRef] [PubMed]
  58. Furuse, M.; Hata, M.; Furuse, K.; Yoshida, Y.; Haratake, A.; Sugitani, Y.; Noda, T.; Kubo, A.; Tsukita, S. Claudin-based tight junctions are crucial for the mammalian epidermal barrier. J. Cell Biol. 2002, 156, 1099. [Google Scholar] [CrossRef]
  59. Günzel, D.; Yu, A.S.L. Claudins and the Modulation of Tight Junction Permeability. Physiol. Rev. 2013, 93, 525–569. [Google Scholar] [CrossRef] [Green Version]
  60. Günzel, D.; Fromm, M. Claudins and Other Tight Junction Proteins. Compr. Physiol. 2012, 2, 1819–1852. [Google Scholar]
  61. Van Itallie, C.M.; Anderson, J.M. Claudin interactions in and out of the tight junction. Tissue Barriers 2013, 1, e25247. [Google Scholar] [CrossRef] [Green Version]
  62. Krause, G.; Winkler, L.; Mueller, S.L.; Haseloff, R.F.; Piontek, J.; Blasig, I.E. Structure and function of claudins. Biochim. Biophys. Acta (BBA) Biomembr. 2008, 1778, 631–645. [Google Scholar] [CrossRef] [Green Version]
  63. Buckley, A.; Turner, J.R. Cell Biology of Tight Junction Barrier Regulation and Mucosal Disease. Cold Spring Harb. Perspect. Biol. 2018, 10, a029314. [Google Scholar] [CrossRef]
  64. Garcia-Hernandez, V.; Quiros, M.; Nusrat, A. Intestinal epithelial claudins: Expression and regulation in homeostasis and inflammation. Ann. N. Y. Acad. Sci. 2017, 1397, 66–79. [Google Scholar] [CrossRef]
  65. Hou, J.; Rajagopal, M.; Yu, A.S.L. Claudins and the kidney. Annu. Rev. Physiol. 2013, 75, 479–501. [Google Scholar] [CrossRef] [PubMed]
  66. Hewitt, K.J.; Agarwal, R.; Morin, P.J. The claudin gene family: Expression in normal and neoplastic tissues. BMC Cancer 2006, 6, 186. [Google Scholar] [CrossRef] [PubMed]
  67. Wang, H.; Yang, X. The expression patterns of tight junction protein claudin-1, -3, and -4 in human gastric neoplasms and adjacent non-neoplastic tissues. Int. J. Clin. Exp. Pathol. 2015, 8, 881–887. [Google Scholar] [PubMed]
  68. Suzuki, H.; Nishizawa, T.; Tani, K.; Yamazaki, Y.; Tamura, A.; Ishitani, R.; Dohmae, N.; Tsukita, S.; Nureki, O.; Fujiyoshi, Y. Crystal Structure of a Claudin Provides Insight into the Architecture of Tight Junctions. Science 2014, 344, 304. [Google Scholar] [CrossRef] [PubMed]
  69. Ohtsuki, S.; Yamaguchi, H.; Katsukura, Y.; Asashima, T.; Terasaki, T. mRNA expression levels of tight junction protein genes in mouse brain capillary endothelial cells highly purified by magnetic cell sorting. J. Neurochem. 2008, 104, 147–154. [Google Scholar] [CrossRef]
  70. Daneman, R.; Zhou, L.; Agalliu, D.; Cahoy, J.D.; Kaushal, A.; Barres, B.A. The Mouse Blood-Brain Barrier Transcriptome: A New Resource for Understanding the Development and Function of Brain Endothelial Cells. PLoS ONE 2010, 5, e13741. [Google Scholar] [CrossRef]
  71. Nitta, T.; Hata, M.; Gotoh, S.; Seo, Y.; Sasaki, H.; Hashimoto, N.; Furuse, M.; Tsukita, S. Size-selective loosening of the blood-brain barrier in claudin-5-deficient mice. J. Cell Biol. 2003, 161, 653–660. [Google Scholar] [CrossRef]
  72. Krause, G.; Winkler, L.; Piehl, C.; Blasig, I.; Piontek, J.; Müller, S.L. Structure and Function of Extracellular Claudin Domains. Mol. Biol. Cell 2009, 1165, 4333–4346. [Google Scholar] [CrossRef]
  73. Rossa, J.; Ploeger, C.; Vorreiter, F.; Saleh, T.; Protze, J.; Günzel, D.; Wolburg, H.; Krause, G.; Piontek, J. Claudin-3 and Claudin-5 Protein Folding and Assembly into the Tight Junction Are Controlled by Non-conserved Residues in the Transmembrane 3 (TM3) and Extracellular Loop 2 (ECL2) Segments. J. Biol. Chem. 2014, 289, 7641–7653. [Google Scholar] [CrossRef] [Green Version]
  74. Piehl, C.; Piontek, J.; Cording, J.; Wolburg, H.; Blasig, I.E. Participation of the second extracellular loop of claudin-5 in paracellular tightening against ions, small and large molecules. Cell. Mol. Life Sci. 2010, 67, 2131–2140. [Google Scholar] [CrossRef]
  75. Keaney, J.; Walsh, D.M.; O’Malley, T.; Hudson, N.; Crosbie, D.E.; Loftus, T.; Sheehan, F.; McDaid, J.; Humphries, M.M.; Callanan, J.J.; et al. Autoregulated paracellular clearance of amyloid-β across the blood-brain barrier. Sci. Adv. 2015, 1, e1500472. [Google Scholar] [CrossRef] [PubMed]
  76. Greene, C.; Hanley, N.; Campbell, M. Claudin-5: Gatekeeper of neurological function. Fluids Barriers CNS 2019, 16, 3. [Google Scholar] [CrossRef] [PubMed]
  77. Di Ventura, B.; Lemerle, C.; Michalodimitrakis, K.; Serrano, L. From in vivo to in silico biology and back. Nature 2006, 443, 527–533. [Google Scholar] [CrossRef] [PubMed]
  78. Ekins, S.; Mestres, J.; Testa, B. In silico pharmacology for drug discovery: Methods for virtual ligand screening and profiling. Br. J. Pharmacol. 2007, 152, 9–20. [Google Scholar] [CrossRef]
  79. McInnes, C. Virtual screening strategies in drug discovery. Curr. Opin. Chem. Biol. 2007, 11, 494–502. [Google Scholar] [CrossRef]
  80. Huggins, D.J.; Biggin, P.C.; Dämgen, M.A.; Essex, J.W.; Harris, S.A.; Henchman, R.H.; Khalid, S.; Kuzmanic, A.; Laughton, C.A.; Michel, J.; et al. Biomolecular simulations: From dynamics and mechanisms to computational assays of biological activity. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2019, 9, e1393. [Google Scholar] [CrossRef]
  81. Hope, H. Crystallography of Biological Macromolecules at Ultra-Low Temperature. Annu. Rev. Biophys. Biophys. Chem. 1990, 19, 107–126. [Google Scholar] [CrossRef]
  82. Smyth, M.S.; Martin, J.H. x ray crystallography. Mol. Pathol. 2000, 53, 8–14. [Google Scholar] [CrossRef]
  83. Carpenter, E.P.; Beis, K.; Cameron, A.D.; Iwata, S. Overcoming the challenges of membrane protein crystallography. Curr. Opin. Struct. Biol. 2008, 18, 581–586. [Google Scholar] [CrossRef] [Green Version]
  84. Caffrey, M. Crystallizing Membrane Proteins for Structure Determination: Use of Lipidic Mesophases. Annu. Rev. Biophys. 2009, 38, 29–51. [Google Scholar] [CrossRef]
  85. Seddon, A.M.; Curnow, P.; Booth, P.J. Membrane proteins, lipids and detergents: Not just a soap opera. Biochim. Biophys. Acta (BBA) Biomembr. 2004, 1666, 105–117. [Google Scholar] [CrossRef] [PubMed]
  86. Dimmic, M.W.; Hubisz, M.J.; Bustamante, C.D.; Nielsen, R. Detecting coevolving amino acid sites using Bayesian mutational mapping. Bioinformatics 2005, 21, i126–i135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Lui, S.; Tiana, G. The network of stabilizing contacts in proteins studied by coevolutionary data. J. Chem. Phys. 2013, 139, 155103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  88. Jana, B.; Morcos, F.; Onuchic, J.N. From structure to function: The convergence of structure based models and co-evolutionary information. Phys. Chem. Chem. Phys. 2014, 16, 6496–6507. [Google Scholar] [CrossRef] [PubMed]
  89. Wang, Y.; Barth, P. Evolutionary-guided de novo structure prediction of self-associated transmembrane helical proteins with near-atomic accuracy. Nat. Commun. 2015, 6, 7196. [Google Scholar] [CrossRef] [PubMed]
  90. Teixeira, P.L.; Mendenhall, J.L.; Heinze, S.; Weiner, B.; Skwark, M.J.; Meiler, J. Membrane protein contact and structure prediction using co-evolution in conjunction with machine learning. PLoS ONE 2017, 12, e0177866. [Google Scholar] [CrossRef]
  91. Rajagopal, N.; Irudayanathan, F.J.; Nangia, S. Palmitoylation of Claudin-5 Proteins Influences Their Lipid Domain Affinity and Tight Junction Assembly at the Blood–Brain Barrier Interface. J. Phys. Chem. B 2019, 123, 983–993. [Google Scholar] [CrossRef]
  92. Murthy, A.C.; Dignon, G.L.; Kan, Y.; Zerze, G.H.; Parekh, S.H.; Mittal, J.; Fawzi, N.L. Molecular interactions underlying liquid−liquid phase separation of the FUS low-complexity domain. Nat. Struct. Mol. Biol. 2019, 26, 637–648. [Google Scholar] [CrossRef]
  93. Zheng, M.; Zhao, J.; Cui, C.; Fu, Z.; Li, X.; Liu, X.; Ding, X.; Tan, X.; Li, F.; Luo, X.; et al. Computational chemical biology and drug design: Facilitating protein structure, function, and modulation studies. Med. Res. Rev. 2018, 38, 914–950. [Google Scholar] [CrossRef]
  94. Irudayanathan, F.J.; Wang, X.; Wang, N.; Willsey, S.R.; Seddon, I.A.; Nangia, S. Self-Assembly Simulations of Classic Claudins—Insights into the Pore Structure, Selectivity, and Higher Order Complexes. J. Phys. Chem. B 2018, 122, 7463–7474. [Google Scholar] [CrossRef]
  95. Ibsen, K.N.; Ma, H.; Banerjee, A.; Tanner, E.E.L.; Nangia, S.; Mitragotri, S. Mechanism of Antibacterial Activity of Choline-Based Ionic Liquids (CAGE). ACS Biomater. Sci. Eng. 2018, 4, 2370–2379. [Google Scholar] [CrossRef]
  96. Sengupta, D.; Prasanna, X.; Mohole, M.; Chattopadhyay, A. Exploring GPCR–Lipid Interactions by Molecular Dynamics Simulations: Excitements, Challenges, and the Way Forward. J. Phys. Chem. B 2018, 122, 5727–5737. [Google Scholar] [CrossRef] [PubMed]
  97. Irudayanathan, F.J.; Wang, N.; Wang, X.; Nangia, S. Architecture of the paracellular channels formed by claudins of the blood–brain barrier tight junctions. Ann. N. Y. Acad. Sci. 2017, 1405, 131–146. [Google Scholar] [CrossRef] [PubMed]
  98. Grime, J.M.A.; Dama, J.F.; Ganser-Pornillos, B.K.; Woodward, C.L.; Jensen, G.J.; Yeager, M.; Voth, G.A. Coarse-grained simulation reveals key features of HIV-1 capsid self-assembly. Nat. Commun. 2016, 7, 11568. [Google Scholar] [CrossRef] [PubMed]
  99. Ingólfsson, H.I.; Arnarez, C.; Periole, X.; Marrink, S.J. Computational ‘microscopy’ of cellular membranes. J. Cell Sci. 2016, 129, 257. [Google Scholar] [CrossRef] [PubMed]
  100. Reddy, T.; Sansom, M.S. The Role of the Membrane in the Structure and Biophysical Robustness of the Dengue Virion Envelope. Structure 2016, 24, 375–382. [Google Scholar] [CrossRef] [Green Version]
  101. Arnarez, C.; Marrink, S.J.; Periole, X. Molecular mechanism of cardiolipin-mediated assembly of respiratory chain supercomplexes. Chem. Sci. 2016, 7, 4435–4443. [Google Scholar] [CrossRef] [Green Version]
  102. Holdbrook, D.A.; Huber, R.G.; Piggot, T.J.; Bond, P.J.; Khalid, S. Dynamics of Crowded Vesicles: Local and Global Responses to Membrane Composition. PLoS ONE 2016, 11, e0156963. [Google Scholar] [CrossRef]
  103. Castillo, N.; Monticelli, L.; Barnoud, J.; Tieleman, D.P. Free energy of WALP23 dimer association in DMPC, DPPC, and DOPC bilayers. Chem. Phys. Lipids 2013, 169, 95–105. [Google Scholar] [CrossRef] [Green Version]
  104. Bennett, W.F.D.; MacCallum, J.L.; Hinner, M.J.; Marrink, S.J.; Tieleman, D.P. Molecular View of Cholesterol Flip-Flop and Chemical Potential in Different Membrane Environments. J. Am. Chem. Soc. 2009, 131, 12714–12720. [Google Scholar] [CrossRef] [Green Version]
  105. Periole, X.; Knepp, A.M.; Sakmar, T.P.; Marrink, S.J.; Huber, T. Structural Determinants of the Supramolecular Organization of G Protein-Coupled Receptors in Bilayers. J. Am. Chem. Soc. 2012, 134, 10959–10965. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Schäfer, L.V.; de Jong, D.H.; Holt, A.; Rzepiela, A.J.; de Vries, A.H.; Poolman, B.; Killian, J.A.; Marrink, S.J. Lipid packing drives the segregation of transmembrane helices into disordered lipid domains in model membranes. Proc. Natl. Acad. Sci. USA 2011, 108, 1343. [Google Scholar] [CrossRef] [PubMed]
  107. Duncan, A.L.; Reddy, T.; Koldsø, H.; Hélie, J.; Fowler, P.W.; Chavent, M.; Sansom, M.S.P. Protein crowding and lipid complexity influence the nanoscale dynamic organization of ion channels in cell membranes. Sci. Rep. 2017, 7, 16647. [Google Scholar] [CrossRef] [PubMed]
  108. Koldsø, H.; Sansom, M.S.P. Organization and Dynamics of Receptor Proteins in a Plasma Membrane. J. Am. Chem. Soc. 2015, 137, 14694–14704. [Google Scholar] [CrossRef] [Green Version]
  109. Shorthouse, D.; Hedger, G.; Koldsø, H.; Sansom, M.S.P. Molecular simulations of glycolipids: Towards mammalian cell membrane models. Biochimie 2016, 120, 105–109. [Google Scholar] [CrossRef]
  110. Koldsø, H.; Reddy, T.; Fowler, P.W.; Duncan, A.L.; Sansom, M.S.P. Membrane Compartmentalization Reducing the Mobility of Lipids and Proteins within a Model Plasma Membrane. J. Phys. Chem. B 2016, 120, 8873–8881. [Google Scholar] [CrossRef] [Green Version]
  111. Duncan, A.L.; Song, W.; Sansom, M.S.P. Lipid-Dependent Regulation of Ion Channels and G Protein–Coupled Receptors: Insights from Structures and Simulations. Annu. Rev. Pharmacol. Toxicol. 2019. [Google Scholar] [CrossRef]
  112. Koldsø, H.; Shorthouse, D.; Hélie, J.; Sansom, M.S.P. Lipid Clustering Correlates with Membrane Curvature as Revealed by Molecular Simulations of Complex Lipid Bilayers. PLoS Comput. Biol. 2014, 10, e1003911. [Google Scholar] [CrossRef]
  113. Duncan, A.L.; Bandurka, M.A.R.; Chavent, M.G.; Rassam, P.; Song, W.; Birkholz, O.; Helie, J.; Reddy, T.; Beliaev, D.; Hambly, B.; et al. How Nanoscale Protein Interactions Determine the Mesoscale Dynamic Organisation of Membrane Proteins. Biophys. J. 2019, 116, 365a. [Google Scholar] [CrossRef] [Green Version]
  114. Lin, X.; Gorfe, A.A.; Levental, I. Protein Partitioning into Ordered Membrane Domains: Insights from Simulations. Biophys. J. 2018, 114, 1936–1944. [Google Scholar] [CrossRef] [Green Version]
  115. Irudayanathan, F.J.; Trasatti, J.P.; Karande, P.; Nangia, S. Molecular Architecture of the Blood Brain Barrier Tight Junction Proteins–A Synergistic Computational and In Vitro Approach. J. Phys. Chem. B 2016, 120, 77–88. [Google Scholar] [CrossRef] [PubMed]
  116. Ma, H.; Irudayanathan, F.J.; Jiang, W.; Nangia, S. Simulating Gram-Negative Bacterial Outer Membrane: A Coarse Grain Model. J. Phys. Chem. B 2015, 119, 14668–14682. [Google Scholar] [CrossRef] [PubMed]
  117. Saitoh, Y.; Suzuki, H.; Tani, K.; Nishikawa, K.; Irie, K.; Ogura, Y.; Tamura, A.; Tsukita, S.; Fujiyoshi, Y. Structural insight into tight junction disassembly by Clostridium perfringens enterotoxin. Science 2015, 347, 775. [Google Scholar] [CrossRef] [PubMed]
  118. Shinoda, T.; Shinya, N.; Ito, K.; Ohsawa, N.; Terada, T.; Hirata, K.; Kawano, Y.; Yamamoto, M.; Kimura-Someya, T.; Yokoyama, S.; et al. Structural basis for disruption of claudin assembly in tight junctions by an enterotoxin. Sci. Rep. 2016, 6, 33632. [Google Scholar] [CrossRef] [PubMed]
  119. Nakamura, S.; Irie, K.; Tanaka, H.; Nishikawa, K.; Suzuki, H.; Saitoh, Y.; Tamura, A.; Tsukita, S.; Fujiyoshi, Y. Morphologic determinant of tight junctions revealed by claudin-3 structures. Nat. Commun. 2019, 10, 816. [Google Scholar] [CrossRef] [PubMed]
  120. Xiang, Z. Advances in homology protein structure modeling. Curr. Protein Pept. Sci. 2006, 7, 217–227. [Google Scholar] [CrossRef]
  121. Shen, M.-y.; Sali, A. Statistical potential for assessment and prediction of protein structures. Protein Sci. 2006, 15, 2507–2524. [Google Scholar] [CrossRef] [Green Version]
  122. Kopp, J.; Schwede, T. Automated protein structure homology modeling: A progress report. Pharmacogenomics 2004, 5, 405–416. [Google Scholar] [CrossRef]
  123. Fiser, A. Template-based protein structure modeling. Methods Mol. Biol. 2010, 673, 73–94. [Google Scholar]
  124. Mineta, K.; Yamamoto, Y.; Yamazaki, Y.; Tanaka, H.; Tada, Y.; Saito, K.; Tamura, A.; Igarashi, M.; Endo, T.; Takeuchi, K.; et al. Predicted expansion of the claudin multigene family. FEBS Lett. 2011, 585, 606–612. [Google Scholar] [CrossRef] [Green Version]
  125. Kelley, L.A.; Mezulis, S.; Yates, C.M.; Wass, M.N.; Sternberg, M.J.E. The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 2015, 10, 845. [Google Scholar] [CrossRef] [PubMed]
  126. Schwede, T.; Kopp, J.; Guex, N.; Peitsch, M.C. SWISS-MODEL: An automated protein homology-modeling server. Nucleic Acids Res. 2003, 31, 3381–3385. [Google Scholar] [CrossRef] [PubMed]
  127. Roy, A.; Kucukural, A.; Zhang, Y. I-TASSER: A unified platform for automated protein structure and function prediction. Nat. Protoc. 2010, 5, 725–738. [Google Scholar] [CrossRef] [PubMed]
  128. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein structure and function prediction. Nat. Methods 2014, 12, 7. [Google Scholar] [CrossRef] [PubMed]
  129. Yang, J.; Zhang, Y. I-TASSER server: New development for protein structure and function predictions. Nucleic Acids Res. 2015, 43, W174–W181. [Google Scholar] [CrossRef]
  130. Källberg, M.; Wang, H.; Wang, S.; Peng, J.; Wang, Z.; Lu, H.; Xu, J. Template-based protein structure modeling using the RaptorX web server. Nat. Protoc. 2012, 7, 1511–1522. [Google Scholar] [CrossRef] [Green Version]
  131. Peng, J.; Xu, J. RaptorX: Exploiting structure information for protein alignment by statistical inference. Proteins 2011, 79 (Suppl. 10), 161–171. [Google Scholar] [CrossRef]
  132. Fiser, A.; Do, R.K.G.; Šali, A. Modeling of loops in protein structures. Protein Sci. 2000, 9, 1753–1773. [Google Scholar] [CrossRef] [Green Version]
  133. Baker, D.; Sali, A. Protein Structure Prediction and Structural Genomics. Science 2001, 294, 93. [Google Scholar] [CrossRef]
  134. Wiltgen, M. Algorithms for Structure Comparison and Analysis: Homology Modelling of Proteins. In Encyclopedia of Bioinformatics and Computational Biology; Ranganathan, S., Gribskov, M., Nakai, K., Schönbach, C., Eds.; Academic Press: Oxford, UK, 2019; pp. 38–61. [Google Scholar] [CrossRef]
  135. Schymkowitz, J.; Borg, J.; Stricher, F.; Nys, R.; Rousseau, F.; Serrano, L. The FoldX web server: An online force field. Nucleic Acids Res. 2005, 33, W382–W388. [Google Scholar] [CrossRef]
  136. Fiser, A.; Sali, A. ModLoop: Automated modeling of loops in protein structures. Bioinformatics 2003, 19, 2500–2501. [Google Scholar] [CrossRef] [PubMed]
  137. MacKerell, A.D.; Bashford, D.; Bellott, M.; Dunbrack, R.L.; Evanseck, J.D.; Field, M.J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S.; et al. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586–3616. [Google Scholar] [CrossRef] [PubMed]
  138. Wang, H.; Junghans, C.; Kremer, K. Comparative atomistic and coarse-grained study of water: What do we lose by coarse-graining? Eur. Phys. J. E 2009, 28, 221–229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  139. Kmiecik, S.; Gront, D.; Kolinski, M.; Wieteska, L.; Dawid, A.E.; Kolinski, A. Coarse-Grained Protein Models and Their Applications. Chem. Rev. 2016, 116, 7898–7936. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  140. Arkhipov, A.; Freddolino, P.L.; Imada, K.; Namba, K.; Schulten, K. Coarse-Grained Molecular Dynamics Simulations of a Rotating Bacterial Flagellum. Biophys. J. 2006, 91, 4589–4597. [Google Scholar] [CrossRef] [Green Version]
  141. Shih, A.Y.; Arkhipov, A.; Freddolino, P.L.; Schulten, K. Coarse Grained Protein−Lipid Model with Application to Lipoprotein Particles. J. Phys. Chem. B 2006, 110, 3674–3684. [Google Scholar] [CrossRef]
  142. Akhmatskaya, E.; Reich, S. Meso-GSHMC: A stochastic algorithm for meso-scale constant temperature simulations. Procedia Comput. Sci. 2011, 4, 1353–1362. [Google Scholar] [CrossRef] [Green Version]
  143. de Jong, D.H.; Periole, X.; Marrink, S.J. Dimerization of Amino Acid Side Chains: Lessons from the Comparison of Different Force Fields. J. Chem. Theory Comput. 2012, 8, 1003–1014. [Google Scholar] [CrossRef] [Green Version]
  144. Ayton, G.S.; Voth, G.A. Systematic multiscale simulation of membrane protein systems. Curr. Opin. Struct. Biol. 2009, 19, 138–144. [Google Scholar] [CrossRef] [Green Version]
  145. Yuan, C.; Li, S.; Zou, Q.; Ren, Y.; Yan, X. Multiscale simulations for understanding the evolution and mechanism of hierarchical peptide self-assembly. Phys. Chem. Chem. Phys. 2017, 19, 23614–23631. [Google Scholar] [CrossRef]
  146. Singh, N.; Li, W. Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications. Int. J. Mol. Sci. 2019, 20, 3774. [Google Scholar] [CrossRef] [PubMed]
  147. Jiang, W.; Luo, J.; Nangia, S. Multiscale Approach to Investigate Self-Assembly of Telodendrimer Based Nanocarriers for Anticancer Drug Delivery. Langmuir 2015, 31, 4270–4280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  148. Jiang, W.; Wang, X.; Guo, D.; Luo, J.; Nangia, S. Drug-Specific Design of Telodendrimer Architecture for Effective Doxorubicin Encapsulation. J. Phys. Chem. B 2016, 120, 9766–9777. [Google Scholar] [CrossRef] [PubMed]
  149. Cornell, W.D.; Cieplak, P.; Bayly, C.I.; Gould, I.R.; Merz, K.M.; Ferguson, D.M.; Spellmeyer, D.C.; Fox, T.; Caldwell, J.W.; Kollman, P.A. A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules. J. Am. Chem. Soc. 1995, 117, 5179–5197. [Google Scholar] [CrossRef] [Green Version]
  150. Oostenbrink, C.; Villa, A.; Mark, A.E.; Van Gunsteren, W.F. A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force-field parameter sets 53A5 and 53A6. J. Comput. Chem. 2004, 25, 1656–1676. [Google Scholar] [CrossRef] [PubMed]
  151. Jorgensen, W.L.; Maxwell, D.S.; Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. [Google Scholar] [CrossRef]
  152. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  153. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef]
  154. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef]
  155. Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kalé, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef] [Green Version]
  156. Thompson, J.M.T.; Grindon, C.; Harris, S.; Evans, T.; Novik, K.; Coveney, P.; Laughton, C. Large-scale molecular dynamics simulation of DNA: Implementation and validation of the AMBER98 force field in LAMMPS. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2004, 362, 1373–1386. [Google Scholar]
  157. Case, D.A.; Cheatham Iii, T.E.; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M., Jr.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem. 2005, 26, 1668–1688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  158. Wu, E.L.; Cheng, X.; Jo, S.; Rui, H.; Song, K.C.; Dávila-Contreras, E.M.; Qi, Y.; Lee, J.; Monje-Galvan, V.; Venable, R.M.; et al. CHARMM-GUI Membrane Builder toward realistic biological membrane simulations. J. Comput. Chem. 2014, 35, 1997–2004. [Google Scholar] [CrossRef] [PubMed]
  159. Brooks, B.R.; Brooks Iii, C.L.; Mackerell Jr, A.D.; Nilsson, L.; Petrella, R.J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; et al. CHARMM: The biomolecular simulation program. J. Comput. Chem. 2009, 30, 1545–1614. [Google Scholar] [CrossRef] [PubMed]
  160. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  161. The PyMOL Molecular Graphics System; Version 2.0; Schrödinger, LLC: New York, NY, USA, 2015.
  162. Land, H.; Humble, M.S. YASARA: A Tool to Obtain Structural Guidance in Biocatalytic Investigations. In Protein Engineering: Methods and Protocols; Bornscheuer, U.T., Höhne, M., Eds.; Springer New York: New York, NY, USA, 2018; pp. 43–67. [Google Scholar] [CrossRef]
  163. Krieger, E.; Vriend, G. YASARA View—molecular graphics for all devices—From smartphones to workstations. Bioinformatics 2014, 30, 2981–2982. [Google Scholar] [CrossRef] [PubMed]
  164. Krieger, E.; Koraimann, G.; Vriend, G. Increasing the precision of comparative models with YASARA NOVA—A self-parameterizing force field. Proteins Struct. Funct. Bioinform. 2002, 47, 393–402. [Google Scholar] [CrossRef]
  165. Marrink, S.J.; Risselada, H.J.; Yefimov, S.; Tieleman, D.P.; de Vries, A.H. The MARTINI Force Field: Coarse Grained Model for Biomolecular Simulations. J. Phys. Chem. B 2007, 111, 7812–7824. [Google Scholar] [CrossRef] [Green Version]
  166. Monticelli, L.; Kandasamy, S.K.; Periole, X.; Larson, R.G.; Tieleman, D.P.; Marrink, S.-J. The MARTINI Coarse-Grained Force Field: Extension to Proteins. J. Chem. Theory Comput. 2008, 4, 819–834. [Google Scholar] [CrossRef]
  167. Marrink, S.J.; Tieleman, D.P. Perspective on the Martini model. Chem. Soc. Rev. 2013, 42, 6801–6822. [Google Scholar] [CrossRef] [Green Version]
  168. Periole, X.; Marrink, S.-J. The Martini Coarse-Grained Force Field. In Biomolecular Simulations: Methods and Protocols; Monticelli, L., Salonen, E., Eds.; Humana Press: Totowa, NJ, USA, 2013; pp. 533–565. [Google Scholar]
  169. de Jong, D.H.; Singh, G.; Bennett, W.F.D.; Arnarez, C.; Wassenaar, T.A.; Schäfer, L.V.; Periole, X.; Tieleman, D.P.; Marrink, S.J. Improved Parameters for the Martini Coarse-Grained Protein Force Field. J. Chem. Theory Comput. 2013, 9, 687–697. [Google Scholar] [CrossRef] [PubMed]
  170. Darré, L.; Machado, M.R.; Brandner, A.F.; González, H.C.; Ferreira, S.; Pantano, S. SIRAH: A Structurally Unbiased Coarse-Grained Force Field for Proteins with Aqueous Solvation and Long-Range Electrostatics. J. Chem. Theory Comput. 2015, 11, 723–739. [Google Scholar] [CrossRef] [PubMed]
  171. Machado, M.R.; Pantano, S. SIRAH tools: Mapping, backmapping and visualization of coarse-grained models. Bioinformatics 2016, 32, 1568–1570. [Google Scholar] [CrossRef] [PubMed]
  172. Kar, P.; Gopal, S.M.; Cheng, Y.-M.; Predeus, A.; Feig, M. PRIMO: A Transferable Coarse-grained Force Field for Proteins. J. Chem. Theory Comput. 2013, 9, 3769–3788. [Google Scholar] [CrossRef] [PubMed]
  173. Cheon, M.; Chang, I.; Hall, C.K. Extending the PRIME model for protein aggregation to all 20 amino acids. Proteins 2010, 78, 2950–2960. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  174. Han, W.; Wan, C.-K.; Wu, Y.-D. PACE Force Field for Protein Simulations. 2. Folding Simulations of Peptides. J. Chem. Theory Comput. 2010, 6, 3390–3402. [Google Scholar] [CrossRef]
  175. Han, W.; Wan, C.-K.; Jiang, F.; Wu, Y.-D. PACE Force Field for Protein Simulations. 1. Full Parameterization of Version 1 and Verification. J. Chem. Theory Comput. 2010, 6, 3373–3389. [Google Scholar] [CrossRef]
  176. Uusitalo, J.J.; Ingólfsson, H.I.; Marrink, S.J.; Faustino, I. Martini Coarse-Grained Force Field: Extension to RNA. Biophys. J. 2017, 113, 246–256. [Google Scholar] [CrossRef] [Green Version]
  177. Bernardi, R.C.; Melo, M.C.R.; Schulten, K. Enhanced sampling techniques in molecular dynamics simulations of biological systems. Biochim. Et Biophys. Acta (BBA) Gen. Subj. 2015, 1850, 872–877. [Google Scholar] [CrossRef]
  178. Perilla, J.R.; Goh, B.C.; Cassidy, C.K.; Liu, B.; Bernardi, R.C.; Rudack, T.; Yu, H.; Wu, Z.; Schulten, K. Molecular dynamics simulations of large macromolecular complexes. Curr. Opin. Struct. Biol. 2015, 31, 64–74. [Google Scholar] [CrossRef] [Green Version]
  179. Perilla, J.R.; Schulten, K. Physical properties of the HIV-1 capsid from all-atom molecular dynamics simulations. Nat. Commun. 2017, 8, 15959. [Google Scholar] [CrossRef] [PubMed]
  180. Xiao, G.; Ren, M.; Hong, H. 50 million atoms scale molecular dynamics modelling on a single consumer graphics card. Adv. Eng. Softw. 2018, 124, 66–72. [Google Scholar] [CrossRef]
  181. Orsi, M. 15-Molecular simulation of self-assembly. In Self-Assembling Biomaterials; Azevedo, H.S., da Silva, R.M.P., Eds.; Woodhead Publishing: Cambridge, MA, USA, 2018; pp. 305–318. [Google Scholar] [CrossRef]
  182. Ben-Nissan, G.; Sharon, M. Capturing protein structural kinetics by mass spectrometry. Chem. Soc. Rev. 2011, 40, 3627–3637. [Google Scholar] [CrossRef] [PubMed]
  183. Sekhar, A.; Vallurupalli, P.; Kay, L.E. Defining a length scale for millisecond-timescale protein conformational exchange. Proc. Natl. Acad. Sci. USA 2013, 110, 11391. [Google Scholar] [CrossRef]
  184. Agrawal, P.; Singh, H.; Srivastava, H.K.; Singh, S.; Kishore, G.; Raghava, G.P.S. Benchmarking of different molecular docking methods for protein-peptide docking. BMC Bioinform. 2019, 19, 426. [Google Scholar] [CrossRef]
  185. Pagadala, N.S.; Syed, K.; Tuszynski, J. Software for molecular docking: A review. Biophys. Rev. 2017, 9, 91–102. [Google Scholar] [CrossRef]
  186. Kozakov, D.; Hall, D.R.; Xia, B.; Porter, K.A.; Padhorny, D.; Yueh, C.; Beglov, D.; Vajda, S. The ClusPro web server for protein–protein docking. Nat. Protoc. 2017, 12, 255. [Google Scholar] [CrossRef]
  187. Comeau, S.R.; Gatchell, D.W.; Vajda, S.; Camacho, C.J. ClusPro: A fully automated algorithm for protein–protein docking. Nucleic Acids Res. 2004, 32, W96–W99. [Google Scholar] [CrossRef]
  188. Lyskov, S.; Gray, J.J. The RosettaDock server for local protein-protein docking. Nucleic Acids Res. 2008, 36, W233–W238. [Google Scholar] [CrossRef] [Green Version]
  189. de Vries, S.J.; van Dijk, M.; Bonvin, A.M.J.J. The HADDOCK web server for data-driven biomolecular docking. Nat. Protoc. 2010, 5, 883. [Google Scholar] [CrossRef]
  190. Dominguez, C.; Boelens, R.; Bonvin, A.M.J.J. HADDOCK: A Protein−Protein Docking Approach Based on Biochemical or Biophysical Information. J. Am. Chem. Soc. 2003, 125, 1731–1737. [Google Scholar] [CrossRef] [PubMed]
  191. Chen, R.; Li, L.; Weng, Z. ZDOCK: An initial-stage protein-docking algorithm. Proteins Struct. Funct. Bioinform. 2003, 52, 80–87. [Google Scholar] [CrossRef] [PubMed]
  192. Liu, Y.; Zhao, L.; Li, W.; Zhao, D.; Song, M.; Yang, Y. FIPSDock: A new molecular docking technique driven by fully informed swarm optimization algorithm. J. Comput. Chem. 2013, 34, 67–75. [Google Scholar] [CrossRef] [PubMed]
  193. Gagnon, J.K.; Law, S.M.; Brooks Iii, C.L. Flexible CDOCKER: Development and application of a pseudo-explicit structure-based docking method within CHARMM. J. Comput. Chem. 2016, 37, 753–762. [Google Scholar] [CrossRef] [PubMed]
  194. Pons, C.; Garzon, J.I.; Lopéz-Blanco, J.R.; Fernandez-Recio, J.; Kovacs, J.; Chacon, P.; Abagyan, R. FRODOCK: A new approach for fast rotational protein–protein docking. Bioinformatics 2009, 25, 2544–2551. [Google Scholar]
  195. Liu, M.; Wang, S. MCDOCK: A Monte Carlo simulation approach to the molecular docking problem. J. Comput. Aided Mol. Des. 1999, 13, 435–451. [Google Scholar] [CrossRef]
  196. Yellapu, N.K. Molecular Modelling, Dynamics, and Docking of Membrane Proteins: Still a Challenge. In Applied Case Studies and Solutions in Molecular Docking-Based Drug Design; IGI Global: Hershey, PA, USA, 2016; pp. 186–208. [Google Scholar]
  197. Burke, K.A.; Yates, E.A.; Legleiter, J. Biophysical insights into how surfaces, including lipid membranes, modulate protein aggregation related to neurodegeneration. Front. Neurol 2013, 4, 17. [Google Scholar] [CrossRef]
  198. Bogdanov, M.; Dowhan, W.; Vitrac, H. Lipids and topological rules governing membrane protein assembly. Biochim. Biophys. Acta 2014, 1843, 1475–1488. [Google Scholar] [CrossRef]
  199. Bogdanov, M.; Mileykovskaya, E.; Dowhan, W. Lipids in the assembly of membrane proteins and organization of protein supercomplexes: Implications for lipid-linked disorders. Subcell. Biochem. 2008, 49, 197–239. [Google Scholar]
  200. Ho, C.-C. Chapter 7—Membranes for Bioseparations. In Bioprocessing for Value-Added Products from Renewable Resources; Yang, S.-T., Ed.; Elsevier: Amsterdam, The Netherlands, 2007; pp. 163–183. [Google Scholar] [CrossRef]
  201. Page, R.C.; Li, C.; Hu, J.; Gao, F.P.; Cross, T.A. Lipid bilayers: An essential environment for the understanding of membrane proteins. Magn. Reson. Chem. 2007, 45, S2–S11. [Google Scholar] [CrossRef]
  202. Sengupta, D.; Marrink, S.J. Lipid-mediated interactions tune the association of glycophorin A helix and its disruptive mutants in membranes. Phys. Chem. Chem. Phys. 2010, 12, 12987–12996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  203. Wassenaar, T.A.; Pluhackova, K.; Moussatova, A.; Sengupta, D.; Marrink, S.J.; Tieleman, D.P.; Böckmann, R.A. High-Throughput Simulations of Dimer and Trimer Assembly of Membrane Proteins. The DAFT Approach. J. Chem. Theory Comput. 2015, 11, 2278–2291. [Google Scholar] [CrossRef] [PubMed]
  204. Altwaijry, N.A.; Baron, M.; Wright, D.W.; Coveney, P.V.; Townsend-Nicholson, A. An Ensemble-Based Protocol for the Computational Prediction of Helix–Helix Interactions in G Protein-Coupled Receptors using Coarse-Grained Molecular Dynamics. J. Chem. Theory Comput. 2017, 13, 2254–2270. [Google Scholar] [CrossRef] [PubMed]
  205. Hall, B.A.; Halim, K.B.A.; Buyan, A.; Emmanouil, B.; Sansom, M.S.P. Sidekick for Membrane Simulations: Automated Ensemble Molecular Dynamics Simulations of Transmembrane Helices. J. Chem. Theory Comput. 2014, 10, 2165–2175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  206. Rajagopal, N.; Nangia, S. Obtaining Protein Association Energy Landscape (PANEL) for Integral Membrane Proteins. J. Chem. Theory Comput. 2019. [Google Scholar] [CrossRef] [PubMed]
  207. Petrek, M.; Otyepka, M.; Banás, P.; Kosinová, P.; Koca, J.; Damborský, J. CAVER: A new tool to explore routes from protein clefts, pockets and cavities. BMC Bioinform. 2006, 7, 316. [Google Scholar] [CrossRef]
  208. Petřek, M.; Košinová, P.; Koča, J.; Otyepka, M. MOLE: A Voronoi Diagram-Based Explorer of Molecular Channels, Pores, and Tunnels. Structure 2007, 15, 1357–1363. [Google Scholar] [CrossRef] [Green Version]
  209. Smart, O.S.; Neduvelil, J.G.; Wang, X.; Wallace, B.A.; Sansom, M.S. HOLE: A program for the analysis of the pore dimensions of ion channel structural models. J. Mol Graph. 1996, 14, 354–360. [Google Scholar] [CrossRef]
  210. Yaffe, E.; Fishelovitch, D.; Wolfson, H.J.; Halperin, D.; Nussinov, R. MolAxis: A server for identification of channels in macromolecules. Nucleic Acids Res. 2008, 36, W210–W215. [Google Scholar] [CrossRef]
  211. Yaffe, E.; Fishelovitch, D.; Wolfson, H.J.; Halperin, D.; Nussinov, R. MolAxis: Efficient and accurate identification of channels in macromolecules. Proteins 2008, 73, 72–86. [Google Scholar] [CrossRef] [Green Version]
  212. Pellegrini-Calace, M.; Maiwald, T.; Thornton, J.M. PoreWalker: A Novel Tool for the Identification and Characterization of Channels in Transmembrane Proteins from Their Three-Dimensional Structure. PLoS Comput. Biol. 2009, 5, e1000440. [Google Scholar] [CrossRef] [PubMed]
  213. Peter, C.; Hummer, G. Ion Transport through Membrane-Spanning Nanopores Studied by Molecular Dynamics Simulations and Continuum Electrostatics Calculations. Biophys. J. 2005, 89, 2222–2234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  214. Hummer, G.; Pratt, L.R.; García, A.E. Molecular Theories and Simulation of Ions and Polar Molecules in Water. J. Phys. Chem. A 1998, 102, 7885–7895. [Google Scholar] [CrossRef] [Green Version]
  215. Torrie, G.M.; Valleau, J.P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187–199. [Google Scholar] [CrossRef]
  216. Barducci, A.; Bonomi, M.; Parrinello, M. Metadynamics. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2011, 1, 826–843. [Google Scholar] [CrossRef]
  217. Laio, A.; Gervasio, F.L. Metadynamics: A method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science. Rep. Prog. Phys. 2008, 71, 126601. [Google Scholar] [CrossRef]
  218. Barducci, A.; Bussi, G.; Parrinello, M. Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method. Phys. Rev. Lett. 2008, 100, 020603. [Google Scholar] [CrossRef] [Green Version]
  219. Kumar, S.; Rosenberg, J.M.; Bouzida, D.; Swendsen, R.H.; Kollman, P.A. THE weighted histogram analysis method for free-energy calculations on biomolecules. I. The method. J. Comput. Chem. 1992, 13, 1011–1021. [Google Scholar] [CrossRef]
  220. Bonomi, M.; Branduardi, D.; Bussi, G.; Camilloni, C.; Provasi, D.; Raiteri, P.; Donadio, D.; Marinelli, F.; Pietrucci, F.; Broglia, R.A.; et al. PLUMED: A portable plugin for free-energy calculations with molecular dynamics. Comput. Phys. Commun. 2009, 180, 1961–1972. [Google Scholar] [CrossRef] [Green Version]
  221. Kutzner, C.; Grubmüller, H.; de Groot, B.L.; Zachariae, U. Computational electrophysiology: The molecular dynamics of ion channel permeation and selectivity in atomistic detail. Biophys. J. 2011, 101, 809–817. [Google Scholar] [CrossRef]
  222. Kutzner, C.; Köpfer, D.A.; Machtens, J.-P.; de Groot, B.L.; Song, C.; Zachariae, U. Insights into the function of ion channels by computational electrophysiology simulations. Biochim. Biophys. Acta (BBA) Biomembr. 2016, 1858, 1741–1752. [Google Scholar] [CrossRef] [Green Version]
  223. Epand, R.M. Lipid polymorphism and protein–lipid interactions. Biochim. Biophys. Acta (BBA) Rev. Biomembr. 1998, 1376, 353–368. [Google Scholar] [CrossRef]
  224. Killian, J.A. Hydrophobic mismatch between proteins and lipids in membranes. Biochim. Biophys. Acta (BBA) Rev. Biomembr. 1998, 1376, 401–416. [Google Scholar] [CrossRef]
  225. de Jesus, A.J.; Allen, T.W. The determinants of hydrophobic mismatch response for transmembrane helices. Biochim. Biophys. Acta (BBA) Biomembr. 2013, 1828, 851–863. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  226. Parton, D.L.; Klingelhoefer, J.W.; Sansom, M.S.P. Aggregation of model membrane proteins, modulated by hydrophobic mismatch, membrane curvature, and protein class. Biophys. J. 2011, 101, 691–699. [Google Scholar] [CrossRef] [PubMed]
  227. Bavi, O.; Vossoughi, M.; Naghdabadi, R.; Jamali, Y. The Combined Effect of Hydrophobic Mismatch and Bilayer Local Bending on the Regulation of Mechanosensitive Ion Channels. PLoS ONE 2016, 11, e0150578. [Google Scholar] [CrossRef] [PubMed]
  228. Webb, R.J.; East, J.M.; Sharma, R.P.; Lee, A.G. Hydrophobic Mismatch and the Incorporation of Peptides into Lipid Bilayers: A Possible Mechanism for Retention in the Golgi. Biochemistry 1998, 37, 673–679. [Google Scholar] [CrossRef]
  229. Milovanovic, D.; Honigmann, A.; Koike, S.; Göttfert, F.; Pähler, G.; Junius, M.; Müllar, S.; Diederichsen, U.; Janshoff, A.; Grubmüller, H.; et al. Hydrophobic mismatch sorts SNARE proteins into distinct membrane domains. Nat. Commun. 2015, 6, 5984. [Google Scholar] [CrossRef] [Green Version]
  230. Fowler, P.W.; Williamson, J.J.; Sansom, M.S.P.; Olmsted, P.D. Roles of Interleaflet Coupling and Hydrophobic Mismatch in Lipid Membrane Phase-Separation Kinetics. J. Am. Chem. Soc. 2016, 138, 11633–11642. [Google Scholar] [CrossRef]
  231. Koval, M. Differential pathways of claudin oligomerization and integration into tight junctions. Tissue Barriers 2013, 1, e24518. [Google Scholar] [CrossRef] [Green Version]
  232. Koval, M. Claudins—Key Pieces in the Tight Junction Puzzle. Cell Commun. Adhes. 2006, 13, 127–138. [Google Scholar] [CrossRef] [PubMed]
  233. Heiler, S.; Mu, W.; Zöller, M.; Thuma, F. The importance of claudin-7 palmitoylation on membrane subdomain localization and metastasis-promoting activities. Cell Commun. Signal. 2015, 13. [Google Scholar] [CrossRef] [PubMed]
  234. Van Itallie, C.M.; Gambling, T.M.; Carson, J.L.; Anderson, J.M. Palmitoylation of claudins is required for efficient tight-junction localization. J. Cell Sci. 2005, 118, 1427–1436. [Google Scholar] [CrossRef] [Green Version]
  235. Fukata, Y.; Bredt, D.S.; Fukata, M. Protein Palmitoylation by DHHC Protein Family, the Dynamic Synapse: Molecular Methods in Ionotropic Receptor Biology; CRC Press/Taylor & Francis: Boca Raton, FL, USA, 2006. [Google Scholar]
  236. Gorleku, O.A.; Barns, A.-M.; Prescott, G.R.; Greaves, J.; Chamberlain, L.H. Endoplasmic Reticulum Localization of DHHC Palmitoyltransferases Mediated by Lysine-based Sorting Signals. J. Biol. Chem. 2011, 286, 39573–39584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  237. Chamberlain, L.H.; Shipston, M.J. The Physiology of Protein S-acylation. Physiol. Rev. 2015, 95, 341–376. [Google Scholar] [CrossRef]
  238. Hou, J. Chapter 2—Paracellular Channel Formation. In The Paracellular Channel; Hou, J., Ed.; Academic Press: Cambridge, MA, USA, 2019; pp. 9–27. [Google Scholar] [CrossRef]
  239. Hou, J.; Renigunta, A.; Konrad, M.; Gomes, A.S.; Schneeberger, E.E.; Paul, D.L.; Waldegger, S.; Goodenough, D.A. Claudin-16 and claudin-19 interact and form a cation-selective tight junction complex. J. Clin. Investig. 2008, 118, 619–628. [Google Scholar] [CrossRef] [Green Version]
  240. Mitic, L.L.; Unger, V.M.; Anderson, J.M. Expression, solubilization, and biochemical characterization of the tight junction transmembrane protein claudin-4. Protein Sci. 2003, 12, 218–227. [Google Scholar] [CrossRef] [Green Version]
  241. Van Itallie, C.M.; Mitic, L.L.; Anderson, J.M. Claudin-2 Forms Homodimers and Is a Component of a High Molecular Weight Protein Complex. J. Biol. Chem. 2011, 286, 3442–3450. [Google Scholar] [CrossRef] [Green Version]
  242. Suzuki, H.; Tani, K.; Tamura, A.; Tsukita, S.; Fujiyoshi, Y. Model for the architecture of claudin-based paracellular ion channels through tight junctions. J. Mol. Biol. 2015, 427, 291–297. [Google Scholar] [CrossRef]
  243. Coyne, C.B.; Gambling, T.M.; Boucher, R.C.; Carson, J.L.; Johnson, L.G. Role of claudin interactions in airway tight junctional permeability. Am. J. Physiol. Lung Cell. Mol. Physiol. 2003, 285, L1166–L1178. [Google Scholar] [CrossRef]
  244. Angelow, S.; Yu, A.S.L. Structure-Function Studies of Claudin Extracellular Domains by Cysteine-scanning Mutagenesis. J. Biol. Chem. 2009, 284, 29205–29217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  245. Anbazhagan, V.; Schneider, D. The membrane environment modulates self-association of the human GpA TM domain—Implications for membrane protein folding and transmembrane signaling. Biochim. Biophys. Acta (BBA) Biomembr. 2010, 1798, 1899–1907. [Google Scholar] [CrossRef] [PubMed]
  246. van Meer, G.; Voelker, D.R.; Feigenson, G.W. Membrane lipids: Where they are and how they behave. Nat. Rev. Mol. Cell Biol. 2008, 9, 112. [Google Scholar] [CrossRef] [PubMed]
  247. Von Heijne, G. The membrane protein universe: What’s out there and why bother? J. Intern. Med. 2007, 261, 543–557. [Google Scholar] [CrossRef]
  248. Blaskovic, S.; Blanc, M.; Goot, F.G. What does S-palmitoylation do to membrane proteins? FEBS J. 2013, 280, 2766–2774. [Google Scholar] [CrossRef] [Green Version]
  249. Linder, M.E.; Deschenes, R.J. Palmitoylation: Policing protein stability and traffic. Nat. Rev. Mol. Cell Biol. 2007, 8, 74–84. [Google Scholar] [CrossRef]
  250. Rossa, J.; Lorenz, D.; Ringling, M.; Veshnyakova, A.; Piontek, J. Overexpression of claudin-5 but not claudin-3 induces formation of trans-interaction–dependent multilamellar bodies. Ann. N. Y. Acad. Sci. 2012, 1257, 59–66. [Google Scholar] [CrossRef]
  251. Piontek, J.; Winkler, L.; Wolburg, H.; Müller, S.L.; Zuleger, N.; Piehl, C.; Wiesner, B.; Krause, G.; Blasig, I.E. Formation of tight junction: Determinants of homophilic interaction between classic claudins. FASEB J. 2007, 22, 146–158. [Google Scholar] [CrossRef]
  252. Kausalya, P.J.; Amasheh, S.; Günzel, D.; Wurps, H.; Müller, D.; Fromm, M.; Hunziker, W. Disease-associated mutations affect intracellular traffic and paracellular Mg2+ transport function of Claudin-16. J. Clin Investig. 2006, 116, 878–891. [Google Scholar] [CrossRef]
  253. Galochkina, T.; Ng Fuk Chong, M.; Challali, L.; Abbar, S.; Etchebest, C. New insights into GluT1 mechanics during glucose transfer. Sci. Rep. 2019, 9, 998. [Google Scholar] [CrossRef]
  254. Pan, Y.; Zhang, Y.; Gongpan, P.; Zhang, Q.; Huang, S.; Wang, B.; Xu, B.; Shan, Y.; Xiong, W.; Li, G.; et al. Single glucose molecule transport process revealed by force tracing and molecular dynamics simulations. Nanoscale Horiz. 2018, 3, 517–524. [Google Scholar] [CrossRef]
  255. Samanta, P.; Wang, Y.; Fuladi, S.; Zou, J.; Li, Y.; Shen, L.; Weber, C.; Khalili-Araghi, F. Molecular determination of claudin-15 organization and channel selectivity. J. Gen. Physiol. 2018, 150, 949. [Google Scholar] [CrossRef] [PubMed]
  256. Alberini, G.; Benfenati, F.; Maragliano, L. Molecular Dynamics Simulations of Ion Selectivity in a Claudin-15 Paracellular Channel. J. Phys. Chem. B 2018, 122, 10783–10792. [Google Scholar] [CrossRef] [PubMed]
  257. Alberini, G.; Benfenati, F.; Maragliano, L. A refined model of claudin-15 tight junction paracellular architecture by molecular dynamics simulations. PLoS ONE 2017, 12, e0184190. [Google Scholar] [CrossRef] [PubMed]
  258. Weber, C.R.; Turner, J.R. Dynamic modeling of the tight junction pore pathway. Ann. N. Y. Acad. Sci. 2017, 1397, 209–218. [Google Scholar] [CrossRef]
  259. Shen, L.; Weber, C.R.; Turner, J.R. The tight junction protein complex undergoes rapid and continuous molecular remodeling at steady state. J. Cell Biol. 2008, 181, 683. [Google Scholar] [CrossRef]
  260. Krystofiak, E.S.; Heymann, J.B.; Kachar, B. Carbon replicas reveal double stranded structure of tight junctions in phase-contrast electron microscopy. Commun. Biol. 2019, 2, 98. [Google Scholar] [CrossRef]
Figure 1. The blood–brain barrier. (a) Discovery of the blood–brain barrier (BBB), (b) pathways of molecular transport across the neurovascular endothelial barrier: transcellular (arrow through the cell) and paracellular (arrow through space between the two cells) and (c) in vitro characterization of the paracellular barrier strength of movement of therapeutic molecules (downward arrows) using transepithelial electric resistance (TEER) model.
Figure 1. The blood–brain barrier. (a) Discovery of the blood–brain barrier (BBB), (b) pathways of molecular transport across the neurovascular endothelial barrier: transcellular (arrow through the cell) and paracellular (arrow through space between the two cells) and (c) in vitro characterization of the paracellular barrier strength of movement of therapeutic molecules (downward arrows) using transepithelial electric resistance (TEER) model.
Ijms 20 05583 g001
Figure 2. Homology modeled claudin-5 structure shows four transmembrane helices (TM1−4), extracellular loops (ECL1 and ECL2), and intracellular loop (ICL) in cartoon representation colored from N-terminal (blue) to truncated C-terminal (red). The membrane spanning region lies within the horizontal (red and blue) planes.
Figure 2. Homology modeled claudin-5 structure shows four transmembrane helices (TM1−4), extracellular loops (ECL1 and ECL2), and intracellular loop (ICL) in cartoon representation colored from N-terminal (blue) to truncated C-terminal (red). The membrane spanning region lies within the horizontal (red and blue) planes.
Ijms 20 05583 g002
Figure 3. An overview of computational approaches and tools used to study the tight junctions at the blood–brain barrier. The stepwise research strategy included: obtaining protein sequence from online libraries; determining protein structure; equilibrating the protein in membrane lipids; investigating cis assembly (dimer and higher-order); investigating trans assembly for predicting tight junction barrier and pores; evaluating transport properties on ions and molecules; and testing tight junction modulators for enhanced drug delivery. The computational approaches (homology modeling, atomistic and coarse-grained molecular dynamics self-assembly simulations, PANEL, etc.) and tools (online libraries, webservers, force-fields, and software suites) associated with each step of our research strategy are also provided.
Figure 3. An overview of computational approaches and tools used to study the tight junctions at the blood–brain barrier. The stepwise research strategy included: obtaining protein sequence from online libraries; determining protein structure; equilibrating the protein in membrane lipids; investigating cis assembly (dimer and higher-order); investigating trans assembly for predicting tight junction barrier and pores; evaluating transport properties on ions and molecules; and testing tight junction modulators for enhanced drug delivery. The computational approaches (homology modeling, atomistic and coarse-grained molecular dynamics self-assembly simulations, PANEL, etc.) and tools (online libraries, webservers, force-fields, and software suites) associated with each step of our research strategy are also provided.
Ijms 20 05583 g003
Figure 4. Typical form of a molecular force-field constituted by bonded (bond, angle, dihedral, and torsion) interactions and nonbonded (electrostatic and van der Waal’s) interaction energy terms. The force filed parameter set include: equilibrium bond length ( l 0 ) , bond angle ( θ 0 ) , dihedral angle ( ω 0 ) , torsion angle ( φ 0 ) , and their respective force constants, k l , k θ , k ω , and k φ , as well as charge ( q i )   on each atom i , dielectric constant ( ϵ 0 ) , strength of dispersion interactions ( ϵ i j ) , and contact distance ( σ i j ) between the atoms i and j. The interatomic distance is represented by r i j .
Figure 4. Typical form of a molecular force-field constituted by bonded (bond, angle, dihedral, and torsion) interactions and nonbonded (electrostatic and van der Waal’s) interaction energy terms. The force filed parameter set include: equilibrium bond length ( l 0 ) , bond angle ( θ 0 ) , dihedral angle ( ω 0 ) , torsion angle ( φ 0 ) , and their respective force constants, k l , k θ , k ω , and k φ , as well as charge ( q i )   on each atom i , dielectric constant ( ϵ 0 ) , strength of dispersion interactions ( ϵ i j ) , and contact distance ( σ i j ) between the atoms i and j. The interatomic distance is represented by r i j .
Ijms 20 05583 g004
Figure 5. A schematic view of claudin-5′s cis and trans assembly along the secretory pathway. The cis dimerization of claudin-5 proteins in endoplasmic membrane (ER) and Golgi complex membrane can result in multiple conformations that are categorized as dimers A, B, C, and D [94,97,115]; posttranslational modification can influence the relative stability of these dimers [91]; the trans interaction of C dimer pair results in a paracellular barrier, whereas trans interaction of B–B or D–D dimer pairs result in paracellular pore.
Figure 5. A schematic view of claudin-5′s cis and trans assembly along the secretory pathway. The cis dimerization of claudin-5 proteins in endoplasmic membrane (ER) and Golgi complex membrane can result in multiple conformations that are categorized as dimers A, B, C, and D [94,97,115]; posttranslational modification can influence the relative stability of these dimers [91]; the trans interaction of C dimer pair results in a paracellular barrier, whereas trans interaction of B–B or D–D dimer pairs result in paracellular pore.
Ijms 20 05583 g005
Figure 6. Claudin-5 strand formation in coarse-grain molecular dynamics (CG MD) self-assembly simulations. (a) Snapshots (top view) of claudin-5 (purple beads) assembly in lipid (cyan beads) membrane at t = 0, 5, and 10 µs, (b) close up view of a self-assembled claudin-5 pentamer, and (c) reversed-mapped strand showing dimer A, C, and D conformations.
Figure 6. Claudin-5 strand formation in coarse-grain molecular dynamics (CG MD) self-assembly simulations. (a) Snapshots (top view) of claudin-5 (purple beads) assembly in lipid (cyan beads) membrane at t = 0, 5, and 10 µs, (b) close up view of a self-assembled claudin-5 pentamer, and (c) reversed-mapped strand showing dimer A, C, and D conformations.
Ijms 20 05583 g006
Figure 7. Claudin-5 cis dimerization. (a) Key cis dimers formed in claudin-5, reverse mapped to atomistic resolutions, and (b) PANEL plot showing location of low energy dimers A1, A2, and A3 (blue circles); and dimers B, C, and D (black circles).
Figure 7. Claudin-5 cis dimerization. (a) Key cis dimers formed in claudin-5, reverse mapped to atomistic resolutions, and (b) PANEL plot showing location of low energy dimers A1, A2, and A3 (blue circles); and dimers B, C, and D (black circles).
Ijms 20 05583 g007
Figure 8. Effects of lipid complexity on claudin-5 dimerization. (a) Snapshots (top view) of self-assembled claudins (green) embedded in ER, cholesterol-enriched endoplasmic reticulum (ERc), and plasma membrane (PM) membranes with saturated lipids (red), unsaturated (blue), and cholesterol (white). (b) PANEL plots generated for claudin-5 dimers embedded in symmetric DOPC bilayer and ERc membrane. The ΔPANEL plot shows stabilization (blue) and destabilization (red) of the dimer geometries due to change in lipid environment DOPC to ERc. The stabilized dimer B location is marked (circle, black) on the ΔPANEL plot.
Figure 8. Effects of lipid complexity on claudin-5 dimerization. (a) Snapshots (top view) of self-assembled claudins (green) embedded in ER, cholesterol-enriched endoplasmic reticulum (ERc), and plasma membrane (PM) membranes with saturated lipids (red), unsaturated (blue), and cholesterol (white). (b) PANEL plots generated for claudin-5 dimers embedded in symmetric DOPC bilayer and ERc membrane. The ΔPANEL plot shows stabilization (blue) and destabilization (red) of the dimer geometries due to change in lipid environment DOPC to ERc. The stabilized dimer B location is marked (circle, black) on the ΔPANEL plot.
Ijms 20 05583 g008
Figure 9. Effects of palmitoylation on claudin-5 cis dimerization. (a) PANEL plot Figure claudin. and dimer C conformation (green ribbon), (b) PANEL plot for claudin-5P and dimer C conformation (green ribbon) with palmitoylation (pink sticks), and (c) ΔPANEL showing the difference in stability upon palmitoylation. The location of dimers B and C are marked by circles (black) on all plots.
Figure 9. Effects of palmitoylation on claudin-5 cis dimerization. (a) PANEL plot Figure claudin. and dimer C conformation (green ribbon), (b) PANEL plot for claudin-5P and dimer C conformation (green ribbon) with palmitoylation (pink sticks), and (c) ΔPANEL showing the difference in stability upon palmitoylation. The location of dimers B and C are marked by circles (black) on all plots.
Ijms 20 05583 g009
Figure 10. Claudin-5 assembly in two adjacent membranes. A 20 µs snapshot of claudin-5 (purple) trans and cis assembly in a coarse gained MD simulation. The zoomed-in view shows the interacting claudins embedded in the membrane (lipid head groups (blue); lipid tails (light brown). Water and ions in the system are not shown for clarity.
Figure 10. Claudin-5 assembly in two adjacent membranes. A 20 µs snapshot of claudin-5 (purple) trans and cis assembly in a coarse gained MD simulation. The zoomed-in view shows the interacting claudins embedded in the membrane (lipid head groups (blue); lipid tails (light brown). Water and ions in the system are not shown for clarity.
Ijms 20 05583 g010
Figure 11. Model for claudin-5 pore I (left) and pore II (right). (a) Front-view, (b) side-view, and (c) pore diameter along the length of the pore.
Figure 11. Model for claudin-5 pore I (left) and pore II (right). (a) Front-view, (b) side-view, and (c) pore diameter along the length of the pore.
Ijms 20 05583 g011

Share and Cite

MDPI and ACS Style

Rajagopal, N.; Irudayanathan, F.J.; Nangia, S. Computational Nanoscopy of Tight Junctions at the Blood–Brain Barrier Interface. Int. J. Mol. Sci. 2019, 20, 5583. https://doi.org/10.3390/ijms20225583

AMA Style

Rajagopal N, Irudayanathan FJ, Nangia S. Computational Nanoscopy of Tight Junctions at the Blood–Brain Barrier Interface. International Journal of Molecular Sciences. 2019; 20(22):5583. https://doi.org/10.3390/ijms20225583

Chicago/Turabian Style

Rajagopal, Nandhini, Flaviyan Jerome Irudayanathan, and Shikha Nangia. 2019. "Computational Nanoscopy of Tight Junctions at the Blood–Brain Barrier Interface" International Journal of Molecular Sciences 20, no. 22: 5583. https://doi.org/10.3390/ijms20225583

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop