Salt Bridge Network Research Articles

Molecular Dynamics (MD) Simulations Provide Insights into the Activation Mechanisms of 5-HT2A Receptors

Recent breakthroughs in the determination of atomic resolution 3-D cryo-electron microscopy structures of membrane proteins present an unprecedented opportunity for drug discovery. Structure-based drug discovery utilizing in silico methods enables the study of dynamic connectivity of stable conformations induced by the drug in achieving its effect. With the ever-expanding computational power, simulations of this type reveal protein dynamics in the nano-, micro-, and even millisecond time scales. In the present study, aiming to characterize the protein dynamics of the 5HT2A receptor stimulated by ligands (agonist/antagonist), we performed 1 µs MD simulations on 5HT2A/DOI (agonist), 5HT2A/GSK215083 (antagonist), and 5HT2A (APO, no ligand) systems. The crystal structure of 5HT2A/zotepine (antagonist) (PDB: 6A94) was used to set up the simulation systems in a lipid bilayer environment. We found the monitoring of the ionic lock residue pair (R3.50-E6.30) of 5HT2A in MD simulations to be a good approximation of the effects of agonists (ionic lock breakage) or antagonists (ionic lock formation) on receptor activation. We further performed analyses of the MD trajectories, including Principal Component Analysis (PCA), hydrogen bond, salt bridge, and hydrophobic interaction network analyses, and correlation between residues to identify key elements of receptor activation. Our results suggest that in order to trigger receptor activation, DOI must interact with 5HT2A through residues V5.39, G5.42, S5.43, and S5.46 on TM5, inducing significant conformational changes in the backbone angles of G5.42 and S5.43. DOI also interacted with residues W6.48 (toggle switch) and F6.51 on TM6, causing major conformational shifts in the backbone angles of F6.44 and V6.45. These structural changes were transmitted to the intracellular ends of TM5, TM6, and ICL3, resulting in the breaking of the ionic lock and subsequent G protein activation. The studies could be helpful in future design of selective agonists/antagonists for various serotonin receptors (5HT1A, 5HT2A, 5HT2B, 5HT2C, and 5HT7) involved in detrimental disorders, such as addiction and schizophrenia.

Sodium Ion-Induced Structural Transition on the Surface of a DNA-Interacting Protein.

Protein surfaces have pivotal roles in interactions between proteins and other biological molecules. However, the structural dynamics of protein surfaces have rarely been explored and are poorly understood. Here, the surface of a single-stranded DNA (ssDNA) binding protein (SSB) with four DNA binding domains that bind ssDNA in binding site sizes of 35, 56, and 65 nucleotides per tetramer is investigated. Using oligonucleotides as probes to sense the charged surface, NaCl induces a two-state structural transition on the SSB surface even at moderate concentrations. Chelation of sodium ions with charged amino acids alters the network of hydrogen bonds and/or salt bridges on the surface. Such changes are associated with changes in the electrostatic potential landscape and interaction mode. These findings advance the understanding of the molecular mechanism underlying the enigmatic salt-induced transitions between different DNA binding site sizes of SSBs. This work demonstrates that monovalent salt is a key regulator of biomolecular interactions that not only play roles in non-specific electrostatic screening effects as usually assumed but also may configure the surface of proteins to contribute to the effective regulation of biomolecular recognition and other downstream events.

ELUCIDATING THE DIFFERENTIAL IMPACTS OF EQUIVALENT GATING-CHARGE MUTATIONS IN VOLTAGE-GATED SODIUM CHANNELS.

Voltage-gated sodium (Nav) channels are pivotal for cellular signaling and mutations in Nav channels can lead to excitability disorders in cardiac, muscular, and neural tissues. A major cluster of pathological mutations localizes in the voltage-sensing domains (VSDs), resulting in either gain-of-function (GoF), loss-of-function (LoF) effects, or both. However, the mechanism behind this functional divergence of mutations at equivalent positions remains elusive. Through hotspot analysis, we identified three gating charges (R1, R2, and R3) as major mutational hotspots in VSDs. The same amino-acid substitutions at equivalent gating-charge positions in VSDI and VSDII of the cardiac sodium channel Nav1.5 show differential gating-property impacts in electrophysiology measurements. We conducted 120 μs molecular dynamics (MD) simulations on wild-type and six mutants to elucidate the structural basis of their differential impacts. Our μs-scale MD simulations with applied external electric fields captured VSD state transitions and revealed the differential structural dynamics between equivalent R-to-Q mutants. Notably, we observed transient leaky conformations in some mutants during structural transitions, offering a detailed structural explanation for gating-pore currents. Our salt-bridge network analysis uncovered VSD-specific and state-dependent interactions among gating charges, countercharges, and lipids. This detailed analysis elucidated how mutations disrupt critical electrostatic interactions, thereby altering VSD permeability and modulating gating properties. By demonstrating the crucial importance of considering the specific structural context of each mutation, our study represents a significant leap forward in understanding structure-function relationships in Nav channels. Our work establishes a robust framework for future investigations into the molecular basis of ion channel-related disorders.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1.

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 in apo. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganizations in key salt bridge and hydrogen bond networks contributing to the CB1 activation/inactivation. For instance, D213-Y224 hydrogen bond and D184-K192 salt bridge showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Insights into the role of glycerophospholipids on the iron export function of SLC40A1 and the molecular mechanisms of ferroportin disease.

SLC40A1 is the sole iron export protein reported in mammals. In humans, its dysfunction is responsible for ferroportin disease, an inborn error of iron metabolism transmitted as an autosomal dominant trait and observed in different ethnic groups. As a member of the major facilitator superfamily, SLC40A1 requires a series of conformational changes to enable iron translocation across the plasma membrane. The influence of lipids on protein stability and its conformational changes has been little investigated to date. Here, we combine molecular dynamics simulations of SLC40A1 embedded in membrane bilayers with experimental alanine scanning mutagenesis to analyze the specific role of glycerophospholipids. We identify four basic residues (Lys90, Arg365, Lys366, and Arg371) that are located at the membrane-cytosol interface and consistently interact with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) molecules. These residues surround a network of salt bridges and hydrogens bonds that play a critical role in stabilizing SLC40A1 in its basal outward-facing conformation. More deeply embedded in the plasma membrane, we identify Arg179 as a charged amino acid residue also tightly interacting with lipid polar heads. This results in a local deformation of the lipid bilayer. Interestingly, Arg179 is adjacent to Arg178, which forms a functionally important salt-bridge with Asp473 and is a recurrently associated with ferroportin disease when mutated to glutamine. We demonstrate that the two p.Arg178Gln and p.Arg179Thr missense variants have similar functional behaviors. These observations provide insights into the role of phospholipids in the formation/disruption of the SLC40A1 inner gate, and give a better understanding of the diversity of molecular mechanisms of ferroportin disease.

Targeting the protein folding transition state by mutation: Large scale (un)folding rate accelerations without altering native stability.

Proteins are constantly undergoing folding and unfolding transitions, with rates that determine their homeostasis in vivo and modulate their biological function. The ability to optimize these rates without affecting overall native stability is hence highly desirable for protein engineering and design. The great challenge is, however, that mutations generally affect folding and unfolding rates with inversely complementary fractions of the net free energy change they inflict on the native state. Here we address this challenge by targeting the folding transition state (FTS) of chymotrypsin inhibitor 2 (CI2), a very slow and stable two-state folding protein with an FTS known to be refractory to change by mutation. We first discovered that the CI2's FTS is energetically taxed by the desolvation of several, highly conserved, charges that form a buried salt bridge network in the native structure. Based on these findings, we designed a CI2 variant that bears just four mutations and aims to selectively stabilize the FTS. This variant has >250-fold faster rates in both directions and hence identical native stability, demonstrating the success of our FTS-centric design strategy. With an optimized FTS, CI2 also becomes 250-fold more sensitive to proteolytic degradation by its natural substrate chymotrypsin, and completely loses its activity as inhibitor. These results indicate that CI2 has been selected through evolution to have a very unstable FTS in order to attain the kinetic stability needed to effectively function as protease inhibitor. Moreover, the CI2 case showcases that protein (un)folding rates can critically pivot around a few key residues-interactions, which can strongly modify the general effects of known structural factors such as domain size and fold topology. From a practical standpoint, our results suggest that future efforts should perhaps focus on identifying such critical residues-interactions in proteins as best strategy to significantly improve our ability to predict and engineer protein (un)folding rates.

Effect of the Lys62Ala Mutation on the Thermal Stability of BstHPr Protein by Molecular Dynamics.

We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 μs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3-Lys62-Glu36 triad, made up of the Glu3-Lys62 and Glu36-Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.

Mobile barrier mechanisms for Na+-coupled symport in an MFS sugar transporter

While many 3D structures of cation-coupled transporters have been determined, the mechanistic details governing the obligatory coupling and functional regulations still remain elusive. The bacterial melibiose transporter (MelB) is a prototype of major facilitator superfamily transporters. With a conformation-selective nanobody, we determined a low-sugar affinity inward-facing Na+-bound cryoEM structure. The available outward-facing sugar-bound structures showed that the N- and C-terminal residues of the inner barrier contribute to the sugar selectivity. The inward-open conformation shows that the sugar selectivity pocket is also broken when the inner barrier is broken. Isothermal titration calorimetry measurements revealed that this inward-facing conformation trapped by this nanobody exhibited a greatly decreased sugar-binding affinity, suggesting the mechanisms for substrate intracellular release and accumulation. While the inner/outer barrier shift directly regulates the sugar-binding affinity, it has little or no effect on the cation binding, which is supported by molecular dynamics simulations. Furthermore, the hydron/deuterium exchange mass spectrometry analyses allowed us to identify dynamic regions; some regions are involved in the functionally important inner barrier-specific salt-bridge network, which indicates their critical roles in the barrier switching mechanisms for transport. These complementary results provided structural and dynamic insights into the mobile barrier mechanism for cation-coupled symport.

Mobile barrier mechanisms for Na+-coupled symport in an MFS sugar transporter.

While many 3D structures of cation-coupled transporters have been determined, the mechanistic details governing the obligatory coupling and functional regulations still remain elusive. The bacterial melibiose transporter (MelB) is a prototype of major facilitator superfamily transporters. With a conformation-selective nanobody, we determined a low-sugar affinity inward-facing Na+-bound cryoEM structure. The available outward-facing sugar-bound structures showed that the N- and C-terminal residues of the inner barrier contribute to the sugar selectivity. The inward-open conformation shows that the sugar selectivity pocket is also broken when the inner barrier is broken. Isothermal titration calorimetry measurements revealed that this inward-facing conformation trapped by this nanobody exhibited a greatly decreased sugar-binding affinity, suggesting the mechanisms for substrate intracellular release and accumulation. While the inner/outer barrier shift directly regulates the sugar-binding affinity, it has little or no effect on the cation binding, which is supported by molecular dynamics simulations. Furthermore, the hydron/deuterium exchange mass spectrometry analyses allowed us to identify dynamic regions; some regions are involved in the functionally important inner barrier-specific salt-bridge network, which indicates their critical roles in the barrier switching mechanisms for transport. These complementary results provided structural and dynamic insights into the mobile barrier mechanism for cation-coupled symport.

The unique salt bridge network in GlacPETase: a key to its stability.

The extensive accumulation of polyethylene terephthalate (PET) has become a critical environmental issue. PET hydrolases can break down PET into its building blocks. Recently, we identified a glacial PET hydrolase GlacPETase sharing less than 31% amino acid identity with any known PET hydrolases. In this study, the crystal structure of GlacPETase was determined at 1.8 Å resolution, revealing unique structural features including a distinctive N-terminal disulfide bond and a specific salt bridge network. Site-directed mutagenesis demonstrated that the disruption of the N-terminal disulfide bond did not reduce GlacPETase's thermostability or its catalytic activity on PET. However, mutations in the salt bridges resulted in changes in melting temperature ranging from -8°C to +2°C and the activity on PET ranging from 17.5% to 145.5% compared to the wild type. Molecular dynamics simulations revealed that these salt bridges stabilized the GlacPETase's structure by maintaining their surrounding structure. Phylogenetic analysis indicated that GlacPETase represented a distinct branch within PET hydrolases-like proteins, with the salt bridges and disulfide bonds in this branch being relatively conserved. This research contributed to the improvement of our comprehension of the structural mechanisms that dictate the thermostability of PET hydrolases, highlighting the diverse characteristics and adaptability observed within PET hydrolases.IMPORTANCEThe pervasive problem of polyethylene terephthalate (PET) pollution in various terrestrial and marine environments is widely acknowledged and continues to escalate. PET hydrolases, such as GlacPETase in this study, offered a solution for breaking down PET. Its unique origin and less than 31% identity with any known PET hydrolases have driven us to resolve its structure. Here, we report the correlation between its unique structure and biochemical properties, focusing on an N-terminal disulfide bond and specific salt bridges. Through site-directed mutagenesis experiments and molecular dynamics simulations, the roles of the N-terminal disulfide bond and salt bridges were elucidated in GlacPETase. This research enhanced our understanding of the role of salt bridges in the thermostability of PET hydrolases, providing a valuable reference for the future engineering of PET hydrolases.

Salt-bridge mediated conformational dynamics in the figure-of-eight knotted ketol acid reductoisomerase (KARI).

The utility of knotted proteins in biological activities has been ambiguous since their discovery. From their evolutionary significance to their functionality in stabilizing the native protein structure, a unilateral conclusion hasn't been achieved yet. While most studies have been performed to understand the stabilizing effect of the knotted fold on the protein chain, more ideas are yet to emerge regarding the interactions in stabilizing the knot. Using classical molecular dynamics (MD) simulations, we have explored the dynamics of the figure-of-eight knotted domain present in ketol acid reductoisomerase (KARI). Our main focus was on the presence of a salt bridge network evident within the knotted region and its role in shaping the conformational dynamics of the knotted chain. Through the potential of mean forces (PMFs) calculation, we have also marked the specific salt bridges that are pivotal in stabilizing the knotted structure. The correlated motions have been further monitored with the help of principal component analysis (PCA) and dynamic cross-correlation maps (DCCM). Furthermore, mutation of the specific salt bridges led to a change in their conformational stability, vindicating their importance.

A cytosolic mutp53(E285K) variant confers chemoresistance of malignant melanoma

Malignant melanoma (MM) is known to be intrinsically chemoresistant, even though only ~20% of MM carry mutations of the tumor suppressor p53. Despite improvement of systemic therapy the mortality rate of patients suffering from metastatic MM is still ~70%, highlighting the need for alternative treatment options or for the re-establishment of conventional therapeutic approaches, including chemotherapy. Screening the p53 mutation status in a cohort of 19 patient-derived melanoma samples, we identified one rarely described missense mutation of p53 leading to E285K amino acid exchange (mutp53(E285K)). Employing structural and computational analysis we revealed a major role of E285 residue in maintaining stable conformation of wild-type p53 (wtp53). E285K mutation was predicted to cause interruption of a salt-bridge network affecting the conformation of the C-terminal helix of the DNA-binding domain (DBD) thereby preventing DNA interaction. In this context, a cluster of frequently mutated amino acid residues in cancer was identified to putatively lead to similar structural effects as E285K substitution (E285 cluster). Functional analysis, including knockdown of endogenous p53 and reconstitution with diverse p53 missense mutants confirmed mutp53(E285K) to have lost transcriptional activity, to be localized in the cytosol of cancer cells, by both means conferring chemoresistance. Re-sensitization to cisplatin-induced cell death was achieved using clinically approved compounds aiming to restore p53 wild-type function (PRIMA1-Met), or inhibition of AKT-driven MAPK survival pathways (afuresertib), in both cases being partially due to ferroptosis induction. Consequently, active ferroptosis induction using the GPX4 inhibitor RSL3 proved superior in tumorselectively fighting MM cells. Due to high prevalence of the E285-cluster mutations in MM as well as in a variety of other tumor types, we conclude this cluster to serve an important function in tumor development and therapy and suggest new implications for ferroptosis induction in therapeutic applications fighting MM in particular and cancer in general.

Exploring the molecular mechanism of cold-adaption of an alkaline protease mutant by molecular dynamics simulations and residue interaction network.

Psychrophilic proteases have attracted enormous attention in past decades, due to their high catalytic activity at low temperatures in a wide range of industrial processes, especially in the detergent and leather industries. Among them, H5 is an alkaline protease mutant, which featuring psychrophilic-like behavior, but the reasons that H5 with higher activity at low temperatures are still poorly understood. Herein, the molecular dynamics (MD) simulations combined with residue interaction network (RIN) were utilized to investigate the mechanisms of the cold-adaption of mutant H5. The results demonstrated that two loops involved in the substrate binding G100-S104 and S125-S129 in H5 had higher mobility, and the distance enlargement between the two loops modulated the substrate's accessibility compared with wild type counterpart. Besides, H5 enhanced conformational flexibility by weakening salt bridges and increasing interaction with the solvent. In particular, the absence of Lys251-Asp197-Arg247 salt bridge network may contribute to the structural mobility. Based on the free energy landscape and molecular mechanics Poisson-Boltzmann surface area of the wild type and H5, it was elucidated that H5 possessed a large population of interconvertible conformations, resulting in the weaker substrate binding free energy. The calculated RIN topology parameters such as the average degree, average cluster coefficient, and average path length further verified that the mutant H5 attenuated residue-to-residue interactions. The investigation of the mechanisms by which how the residue mutation affects the stability and activity of enzymes provides a theoretical basis for the development of cold-adapted protease.

The role of salt bridge networks in the stability of the yeast hexose transporter 1

BackgroundThe yeast S. cerevisiae preferably metabolizes glucose through aerobic glycolysis. Glucose transport is facilitated by multiple hexose transporters (Hxts), and their expression and activity are tightly regulated by multiple mechanisms. However, detailed structural and functional analyses of Hxts remain limited, largely due to the lack of crystal structure. MethodsHomology modeling was used to build a 3D structural model for the yeast glucose transporter Hxt1 and investigate the effects of site directed mutations on Hxt1 stability and glucose transport activity. ResultsThe conserved salt bridge-forming residues observed in the human Glut4 and the yeast glucose receptor Rgt2 were identified within and between the two 6-transmembrane spanning segments of Hxt1. Most of the RGT2 mutations that disrupt the salt bridge networks were known to cause constitutive signal generation, whereas the corresponding substitutions in HXT1 were shown to decrease Hxt1 stability. While substitutions of the two residues in the salt bridge 2 in Glut4—E329Q and E393D—were reported to abolish glucose transport, the equivalent substitutions in Hxt1 (D382Q and E454D) did not affect Hxt1 glucose transport activity. ConclusionsSubstitutions of equivalent salt bridge-forming residues in Hxt1, Rgt2, and Glut4 are predicted to lock them in an inward-facing conformation but lead to different functional consequences. General significanceThe salt bridge networks in yeast and human glucose transporters and yeast glucose receptors may play different roles in maintaining their structural and functional integrity.

Computation-Based Design of Salt Bridges in PETase for Enhanced Thermostability and Performance for PET Degradation.

Polyethylene terephthalate (PET) is one of the most widely used plastics, and the accumulation of PET poses a great threat to the environment. IsPETase can degrade PET rapidly at moderate temperatures, but its application is greatly limited by the low stability. Herein, a molecular dynamics (MD) simulations combined with a sequence alignment strategy was adopted to introduce salt bridges into the flexible region of IsPETase to improve its thermal stability. In the designed variants, the Tm values of IsPETaseI168R/S188D and IsPETaseI168R/S188E were 7.4 and 8.7 °C higher than that of the wild type, respectively. The release of products degraded by IsPETaseI168R/S188E was 4.3 times that of the wild type. Tertiary structure characterization demonstrated that the structure of the variants IsPETaseI168R/S188D and IsPETaseI168R/S188E became more compact. Extensive MD simulations verified that a stable salt bridge was formed between the residue R168 and D186 in IsPETaseI168R/S188D, while in IsPETaseI168R/S188E an R168-D186-E188 salt bridge network was observed. These results confirmed that the proposed computation-based salt bridge design strategy could efficiently generate variants with enhanced thermal stability for the long-term degradation of PET, which would be helpful for the design of enzymes with improved stability.

Molecular Dynamics Simulations of HPr Proteins from a Thermophilic and a Mesophilic Organism: A Comparative Thermal Study.

The histidine-containing phosphocarrier (HPr) is a monomeric protein conserved in Gram-positive bacteria, which may be of mesophilic or thermophilic nature. In particular, the HPr protein from the thermophilic organism B. stearothermophilus is a good model system for thermostability studies, since experimental data, such as crystal structure and thermal stability curves, are available. However, its unfolding mechanism at higher temperatures is yet unclear at a molecular level. Therefore, in this work, we researched the thermal stability of this protein using molecular dynamics simulations, subjecting it to five different temperatures during a time span of 1 μs. The analyses of the structural parameters and molecular interactions were compared with those of the mesophilic homologue HPr protein from B. subtilis. Each simulation was run in triplicate using identical conditions for both proteins. The results showed that the two proteins lose stability as the temperature increases, but the mesophilic structure is more affected. We found that the salt bridge network formed by the triad of Glu3-Lys62-Glu36 residues and the salt bridge made up of Asp79-Lys83 ion pair are key factors to keep stable the thermophilic protein, maintaining the hydrophobic core protected and the structure packed. In addition, these molecular interactions neutralize the negative surface charge, acting as "natural molecular staples".

Classification of GTP-dependent K-Ras4B active and inactive conformational states.

Classifying reliably active and inactive molecular conformations of wildtype (WT) and mutated oncogenic proteins is a key, ongoing challenge in molecular cancer studies. Here, we probe the GTP-bound K-Ras4B conformational dynamics using long-time atomistic molecular dynamics (MD) simulations. We extract and analyze the detailed underlying free energy landscape of WT K-Ras4B. We use two key reaction coordinates, labeled d1 and d2 (i.e., distances coordinating the Pβ atom of the GTP ligand with two key residues, T35 and G60), shown to correlate closely with activities of WT and mutated K-Ras4B. However, our new K-Ras4B conformational kinetics study reveals a more complex network of equilibrium Markovian states. We show that a new reaction coordinate is required to account for the orientation of acidic K-Ras4B sidechains such as D38 with respect to the interface with binding effector RAF1 and rationalize the activation/inactivation propensities and the corresponding molecular binding mechanisms. We use this understanding to unveil how a relatively conservative mutation (i.e., D33E, in the switch I region) can lead to significantly different activation propensities compared with WT K-Ras4B. Our study sheds new light on the ability of residues near the K-Ras4B-RAF1 interface to modulate the network of salt bridges at the binding interface with the RAF1 downstream effector and, thus, to influence the underlying GTP-dependent activation/inactivation mechanism. Altogether, our hybrid MD-docking modeling approach enables the development of new in silico methods for quantitative assessment of activation propensity changes (e.g., due to mutations or local binding environment). It also unveils the underlying molecular mechanisms and facilitates the rational design of new cancer drugs.

Structure of the planar cell polarity cadherins Fat4 and Dachsous1

The atypical cadherins Fat and Dachsous are key regulators of cell growth and animal development. In contrast to classical cadherins, which form homophilic interactions to segregate cells, Fat and Dachsous cadherins form heterophilic interactions to induce cell polarity within tissues. Here, we determine the co-crystal structure of the human homologs Fat4 and Dachsous1 (Dchs1) to establish the molecular basis for Fat-Dachsous interactions. The binding domains of Fat4 and Dchs1 form an extended interface along extracellular cadherin (EC) domains 1-4 of each protein. Biophysical measurements indicate that Fat4-Dchs1 affinity is among the highest reported for cadherin superfamily members, which is attributed to an extensive network of salt bridges not present in structurally similar protocadherin homodimers. Furthermore, modeling suggests that unusual extracellular phosphorylation modifications directly modulate Fat-Dachsous binding by introducing charged contacts across the interface. Collectively, our analyses reveal how the molecular architecture of Fat4-Dchs1 enables them to form long-range, high-affinity interactions to maintain planar cell polarity.

Conversion of Helix 1 into a Loop in Prion Protein Misfolding.

Cellular prion protein PrPC consists of three α-helices, one β-sheet, and an unstructured N-terminal domain. Misfolding of this protein into the scrapie form (PrPSc) increases dramatically the β-sheet content. H1 is the most stable helix on PrPC and contains an unusual number of hydrophilic amino acids. Its fate in PrPSc is not clear. We performed replica exchange molecular dynamics simulations on H1 alone, H1 together with an N-terminally flanking H1B1 loop and H1 in complex with other hydrophilic regions of the prion protein. In the presence of the H99SQWNKPSKPKTNMK113 sequence, H1 is almost completely converted to a loop structure stabilized by a network of salt bridges. On the other hand, H1 retains its helical structure alone or together with the other sequences considered in this study. We carried out an additional simulation by restraining the distance between the two ends of H1, mimicking a possible geometric restriction by the rest of the protein. Even though the loop was the major conformation, a significant amount of helical structure was also observed. This suggests that the interaction with H99SQWNKPSKPKTNMK113 is necessary for complete helix-to-loop conversion.

Elucidating the differential mutational effects of gating charges in voltage-gated sodium channels.

Recent advances in genetics and pathophysiology have identified a large number of mutations in voltage-gated sodium (NaV) channels that cause various heart, muscle, and brain excitability disorders. A significant feature of mutation distribution is that the disease-associated mutations are enriched in gating charges in voltage-sensing domains (VSDs). Intriguingly, the same mutation at the equivalent positions in different VSDs brings distinctive phenotypes.To address this mechanistic mystery, we studied the mutational effects on structural transitions of two gating charges (R2 and R3) in both VSD1 and VSD2 using molecular dynamics simulation and conformational free energy calculation. The application of an external electric field allows us to simulate the VSD transitions between “up” and “down” conformations at a microsecond scale. Our computational study shows that the structural impact of the same mutation (R to Q substitution) not only alters between different gating charges within the same VSD but also differs at the equivalent position between VSD1 and VSD2. The differential impacts are determined by a salt-bridge network. This network formed by gating charges and countercharges differs from one VSD to another. Our salt-bridge analysis reveals that the perturbations induced by mutations on this VSD-specific network bring unique impacts to the formation and stability of the intermediate state. The salt bridges in the intermediate play an unexpected role in determining the transition rate of VSDs, which elucidates a general principle behind the diverse mutational effects and will help us to predict the effects of more disease-associated mutations in voltage-gated ion channels.