Recent experiments suggest that cellular crowding facilitates the target search dynamics of proteins on DNA, the mechanism of which is not yet known. By using large scale computer simulations, we show that two competing factors, namely the width of the depletion layer that separates the crowder cloud from the DNA molecule and the degree of protein-crowder crosstalk, act in harmony to affect the target search dynamics of proteins. The impacts vary from nonspecific to specific target search regime. During a nonspecific search, dynamics of a protein is only minimally affected, whereas, a significantly different behaviour is observed when the protein starts forming a specific protein-DNA complex. We also find that the severity of impacts largely depends upon physiological crowder concentration and deviation from it leads to attenuation in the binding kinetics. Based on extensive kinetic study and binding energy landscape analysis, we further present a comprehensive molecular description of the search process that allows us to interpret the experimental findings.

Abstract The glucose transporter from Staphylococcus epidermidis, GlcPSe, is a homolog of the human GLUT sugar transporters of the major facilitator superfamily. Together with the xylose transporter from Escherichia coli, XylEEc, the other prominent prokaryotic GLUT homolog, GlcPSe, is equipped with a conserved proton-binding site arguing for an electrogenic transport mode. However, the electrophysiological analysis of GlcPSe presented here reveals important differences between the two GLUT homologs. GlcPSe, unlike XylEEc, does not perform steady-state electrogenic transport at symmetrical pH conditions. Furthermore, when a pH gradient is applied, partially uncoupled transport modes can be generated. In contrast to other bacterial sugar transporters analyzed so far, in GlcPSe sugar binding, translocation and release are also accomplished by the deprotonated transporter. Based on these experimental results, we conclude that coupling of sugar and H+ transport is incomplete in GlcPSe. To verify the viability of the observed partially coupled GlcPSe transport modes, we propose a universal eight-state kinetic model in which any degree of coupling is realized and H+/sugar symport represents only a specific instance. Furthermore, using sequence comparison with strictly coupled XylEEc and similar sugar transporters, we identify an additional charged residue that may be essential for effective H+/sugar symport.

Abstract The transport channel of nuclear pore complexes (NPCs) contains a high density of intrinsically disordered proteins that are rich in phenylalanine-glycine (FG)-repeat motifs (FG Nups). The FG Nups interact promiscuously with various nuclear transport receptors (NTRs), such as karyopherins (Kaps), that mediate the trafficking of nucleocytoplasmic cargoes while also generating a selectively permeable barrier against other macromolecules. Although the binding of NTRs to FG Nups increases molecular crowding in the NPC transport channel, it is unclear how this impacts FG Nup barrier function or the movement of other molecules, such as the Ran importer NTF2. Here, we use surface plasmon resonance to evaluate FG Nup conformation, binding equilibria, and interaction kinetics associated with the multivalent binding of NTF2 and karyopherinβ1 (Kapβ1) to Nsp1p molecular brushes. NTF2 and Kapβ1 show different long- and short-lived binding characteristics that emerge from varying degrees of molecular retention and FG repeat binding avidity within the Nsp1p brush. Physiological concentrations of NTF2 produce a collapse of Nsp1p brushes, whereas Kapβ1 binding generates brush extension. However, the presence of prebound Kapβ1 inhibits Nsp1p brush collapse during NTF2 binding, which is dominated by weak, short-lived interactions that derive from steric hindrance and diminished avidity with Nsp1p. This suggests that binding promiscuity confers kinetic advantages to NTF2 by expediting its facilitated diffusion and reinforces the proposal that Kapβ1 contributes to the integral barrier function of the NPC.

Abstract The stimulator-of-interferon-genes (STING) protein is involved in innate immunity. It has recently been shown that modulation of STING can lead to an aggressive antitumor response. DMXAA is an antitumor agent that had shown great promise in murine models but failed in human clinical trials. The molecular target of DMXAA was subsequently shown to be murine STING (mSTING); however, human STING (hSTING) is insensitive to DMXAA. Molecular dynamics simulations were employed to investigate the differences between hSTING and mSTING that could influence DMXAA binding. An initial set of simulations was performed to investigate a single lid region mutation G230I in hSTING (corresponding residue in mSTING is an Ile), which rendered the protein sensitive to DMXAA. The simulations found that an Ile side chain was enough to form a steric barrier that prevents exit of DMXAA, whereas in WT hSTING, the Gly residue that lacks a side chain formed a porous lid region that allowed DMXAA to exit. A second set of molecular dynamics simulations compared the tendency of STING to be in an open-inactive conformation or a closed-active conformation. The results show that hSTING prefers to be in an open-inactive conformation even with cGAMP, the native ligand, bound. On the other hand, mSTING prefers a closed-active conformation even without a ligand bound. These results highlight the challenges in translating a mouse active STING compound into a human active compound, while also providing avenues to pursue for designing a small-molecule drug targeting human STING.

Abstract The structure and functions of the M2 protein from Influenza A are sensitive to pH, cholesterol, and the antiinfluenza drug Amantadine. This is a tetrameric membrane protein of 97 amino-acid residues that has multiple functions, among them as a proton-selective channel and facilitator of viral budding, replacing the need for the ESCRT proteins that other viruses utilize. Here, various amino-acid-specific-labeled samples of the full-length protein were prepared and mixed, so that only interresidue 13C-13C cross peaks between two differently labeled proteins representing interhelical interactions are observed. This channel is activated at slightly acidic pH values in the endosome when the His37 residues in the middle of the transmembrane domain take on a +2 or +3 charged state. Changes observed here in interhelical distances in the N-terminus can be accounted for by modest structural changes, and no significant changes in structure were detected in the C-terminal portion of the channel upon activation of the channel. Amantadine, which blocks proton conductance by binding in the aqueous pore near the N-terminus, however, significantly modifies the tetrameric structure on the opposite side of the membrane. The interactions between the juxtamembrane amphipathic helix of one monomer and its neighboring monomer observed in the absence of drug are disrupted in its presence. However, the addition of cholesterol prevents this structural disruption. In fact, strong interactions are observed between cholesterol and residues in the amphipathic helix, accounting for cholesterol binding adjacent to a native palmitoylation site and near to an interhelix crevice that is typical of cholesterol binding sites. The resultant stabilization of the amphipathic helix deep in the bilayer interface facilitates the bilayer curvature that is essential for viral budding.

Abstract Icosahedral viral capsids are made of a large number of symmetrically organized protein subunits whose local movements can be essential for infection. In the capsid of the minute virus of mice, events required for infection that involve translocation of peptides through capsid pores are associated with a subtle conformational change. In vitro, this change can be reversibly induced by overcoming the energy barrier through mild heating of the capsid, but little is known about the capsid regions involved in the process. Here, we use hydrogen-deuterium exchange coupled to mass spectrometry to analyze the dynamics of the minute virus of mice capsid at increasing temperatures. Our results indicate that the transition associated with peptide translocation involves the structural rearrangement of regions distant from the capsid pores. These alterations are reflected in an increased dynamics of some secondary-structure elements in the capsid shell from which spikes protrude, and a decreased dynamics in the long intertwined loops that form the large capsid spikes. Thus, the translocation events through capsid pores involve a global conformational rearrangement of the capsid and a complex alteration of its equilibrium dynamics. This study additionally demonstrates the potential of hydrogen-deuterium exchange coupled to mass spectrometry to explore in detail temperature-dependent structural dynamics in large macromolecular protein assemblies. Most importantly, it paves the way for undertaking novel studies of the relationship between structure, dynamics, and biological function in virus particles and other large protein cages.

Abstract Interactions of the hydrophobic anticancer drug tamoxifen (TAM) with lipid model membranes were studied using calcein-encapsulated vesicle leakage, attenuated total reflection Fourier transform infrared (FTIR) spectroscopy, small-angle neutron scattering (SANS), atomic force microscopy (AFM) based force spectroscopy, and all-atom molecular dynamics (MD) simulations. The addition of TAM enhances membrane permeability, inducing calcein to translocate from the interior to the exterior of lipid vesicles. A large decrease in the FTIR absorption band’s magnitude was observed in the hydrocarbon chain region, suggesting suppressed bond vibrational dynamics. Bilayer thickening was determined from SANS data. Force spectroscopy measurements indicate that the lipid bilayer area compressibility modulus KA is increased by a large amount after the incorporation of TAM. MD simulations show that TAM decreases the lipid area and increases chain order parameters. Moreover, orientational and positional analyses show that TAM exhibits a highly dynamic conformation within the lipid bilayer. Our detailed experimental and computational studies of TAM interacting with model lipid membranes shed new light on membrane modulation by TAM.

Abstract Glucosylceramide (GlcCer), one of the simplest glycosphingolipids, plays key roles in physiology and pathophysiology. It has been suggested that GlcCer modulates cellular events by forming specialized domains. In this study, we investigated the interplay between GlcCer and cholesterol (Chol), an important lipid involved in the formation of liquid-ordered (lo) phases. Using fluorescence microscopy and spectroscopy, and dynamic and electrophoretic light scattering, we characterized the interaction between these lipids in different pH environments. A quantitative description of the phase behavior of the ternary unsaturated phospholipid/Chol/GlcCer mixture is presented. The results demonstrate coexistence between lo and liquid-disordered (ld) phases. However, the extent of lo/ld phase separation is sparse, mainly due to the ability of GlcCer to segregate into tightly packed gel domains. As a result, the phase diagram of these mixtures is characterized by an extensive three-phase coexistence region of fluid (ld-phospholipid enriched)/lo (Chol enriched)/gel (GlcCer enriched). Moreover, the results show that upon acidification, GlcCer solubility in the lo phase is increased, leading to a larger lo/ld coexistence region. Quantitative analyses allowed us to determine the differences in the composition of the phases at neutral and acidic pH. These results predict the impact of GlcCer on domain formation and membrane organization in complex biological membranes, and provide a background for unraveling the relationship between the biophysical properties of GlcCer and its biological action.

Abstract We studied the viscoelastic properties of homogeneous and inhomogeneous levan-DNA mixtures using optical tweezers and a rotational rheometer. Levan and DNA are important components of the extracellular matrix of bacterial biofilms. Their viscoelastic properties influence the mechanical as well as molecular-transport properties of biofilm. Both macro- and microrheology measurements in homogeneous levan-DNA mixtures revealed pseudoplastic behavior. When the concentration of DNA reached a critical value, levan started to aggregate, forming clusters of a few microns in size. Microrheology using optical tweezers enabled us to measure local viscoelastic properties within the clusters as well as in the DNA phase surrounding the levan aggregates. In phase-separated levan-DNA mixtures, the results of macro- and microrheology differed significantly. The local viscosity and elasticity of levan increased, whereas the local viscosity of DNA decreased. On the other hand, the results of bulk viscosity measurements suggest that levan clusters do not interact strongly with DNA. Upon treatment with DNase, levan aggregates dispersed. These results demonstrate the advantages of microrheological measurements compared to bulk viscoelastic measurements when the materials under investigation are complex and inhomogeneous, as is often the case in biological samples.

Abstract PDEδ is a small protein that binds and controls the trafficking of RAS subfamily proteins. Its inhibition protects initiation of RAS signaling, and it is one of the common targets considered for oncological drug development. In this study, we used solved x-ray structures of inhibitor-bound PDEδ targets to investigate mechanisms of action of six independent all-atom MD simulations. An analysis of atomic simulations combined with the molecular mechanic-Poisson-Boltzmann solvent accessible surface area/generalized Born solvent accessible surface area calculations led to the identification of action mechanisms for a panel of novel PDEδ inhibitors. To the best of our knowledge, this study is one of the first in silico investigations on co-crystallized PDEδ protein. A detailed atomic-scale understanding of the molecular mechanism of PDEδ inhibition may assist in the design of novel PDEδ inhibitors. One of the most common side effects for diverse small molecules/kinase inhibitors is their off-target interactions with cardiac ion channels and human-ether-a-go-go channel specifically. Thus, all of the studied PDEδ inhibitors are also screened in silico at the central cavities of hERG1 potassium channels.

Recent experiments suggest that cellular crowding facilitates the target search dynamics of proteins on DNA, the mechanism of which is not yet known. By using large scale computer simulations, we show that two competing factors, namely the width of the depletion layer that separates the crowder cloud from the DNA molecule and the degree of protein-crowder crosstalk, act in harmony to affect the target search dynamics of proteins. The impacts vary from nonspecific to specific target search regime. During a nonspecific search, dynamics of a protein is only minimally affected, whereas, a significantly different behaviour is observed when the protein starts forming a specific protein-DNA complex. We also find that the severity of impacts largely depends upon physiological crowder concentration and deviation from it leads to attenuation in the binding kinetics. Based on extensive kinetic study and binding energy landscape analysis, we further present a comprehensive molecular description of the search process that allows us to interpret the experimental findings.

Abstract The glucose transporter from Staphylococcus epidermidis, GlcPSe, is a homolog of the human GLUT sugar transporters of the major facilitator superfamily. Together with the xylose transporter from Escherichia coli, XylEEc, the other prominent prokaryotic GLUT homolog, GlcPSe, is equipped with a conserved proton-binding site arguing for an electrogenic transport mode. However, the electrophysiological analysis of GlcPSe presented here reveals important differences between the two GLUT homologs. GlcPSe, unlike XylEEc, does not perform steady-state electrogenic transport at symmetrical pH conditions. Furthermore, when a pH gradient is applied, partially uncoupled transport modes can be generated. In contrast to other bacterial sugar transporters analyzed so far, in GlcPSe sugar binding, translocation and release are also accomplished by the deprotonated transporter. Based on these experimental results, we conclude that coupling of sugar and H+ transport is incomplete in GlcPSe. To verify the viability of the observed partially coupled GlcPSe transport modes, we propose a universal eight-state kinetic model in which any degree of coupling is realized and H+/sugar symport represents only a specific instance. Furthermore, using sequence comparison with strictly coupled XylEEc and similar sugar transporters, we identify an additional charged residue that may be essential for effective H+/sugar symport.

Abstract The transport channel of nuclear pore complexes (NPCs) contains a high density of intrinsically disordered proteins that are rich in phenylalanine-glycine (FG)-repeat motifs (FG Nups). The FG Nups interact promiscuously with various nuclear transport receptors (NTRs), such as karyopherins (Kaps), that mediate the trafficking of nucleocytoplasmic cargoes while also generating a selectively permeable barrier against other macromolecules. Although the binding of NTRs to FG Nups increases molecular crowding in the NPC transport channel, it is unclear how this impacts FG Nup barrier function or the movement of other molecules, such as the Ran importer NTF2. Here, we use surface plasmon resonance to evaluate FG Nup conformation, binding equilibria, and interaction kinetics associated with the multivalent binding of NTF2 and karyopherinβ1 (Kapβ1) to Nsp1p molecular brushes. NTF2 and Kapβ1 show different long- and short-lived binding characteristics that emerge from varying degrees of molecular retention and FG repeat binding avidity within the Nsp1p brush. Physiological concentrations of NTF2 produce a collapse of Nsp1p brushes, whereas Kapβ1 binding generates brush extension. However, the presence of prebound Kapβ1 inhibits Nsp1p brush collapse during NTF2 binding, which is dominated by weak, short-lived interactions that derive from steric hindrance and diminished avidity with Nsp1p. This suggests that binding promiscuity confers kinetic advantages to NTF2 by expediting its facilitated diffusion and reinforces the proposal that Kapβ1 contributes to the integral barrier function of the NPC.

Abstract The stimulator-of-interferon-genes (STING) protein is involved in innate immunity. It has recently been shown that modulation of STING can lead to an aggressive antitumor response. DMXAA is an antitumor agent that had shown great promise in murine models but failed in human clinical trials. The molecular target of DMXAA was subsequently shown to be murine STING (mSTING); however, human STING (hSTING) is insensitive to DMXAA. Molecular dynamics simulations were employed to investigate the differences between hSTING and mSTING that could influence DMXAA binding. An initial set of simulations was performed to investigate a single lid region mutation G230I in hSTING (corresponding residue in mSTING is an Ile), which rendered the protein sensitive to DMXAA. The simulations found that an Ile side chain was enough to form a steric barrier that prevents exit of DMXAA, whereas in WT hSTING, the Gly residue that lacks a side chain formed a porous lid region that allowed DMXAA to exit. A second set of molecular dynamics simulations compared the tendency of STING to be in an open-inactive conformation or a closed-active conformation. The results show that hSTING prefers to be in an open-inactive conformation even with cGAMP, the native ligand, bound. On the other hand, mSTING prefers a closed-active conformation even without a ligand bound. These results highlight the challenges in translating a mouse active STING compound into a human active compound, while also providing avenues to pursue for designing a small-molecule drug targeting human STING.

Abstract The structure and functions of the M2 protein from Influenza A are sensitive to pH, cholesterol, and the antiinfluenza drug Amantadine. This is a tetrameric membrane protein of 97 amino-acid residues that has multiple functions, among them as a proton-selective channel and facilitator of viral budding, replacing the need for the ESCRT proteins that other viruses utilize. Here, various amino-acid-specific-labeled samples of the full-length protein were prepared and mixed, so that only interresidue 13C-13C cross peaks between two differently labeled proteins representing interhelical interactions are observed. This channel is activated at slightly acidic pH values in the endosome when the His37 residues in the middle of the transmembrane domain take on a +2 or +3 charged state. Changes observed here in interhelical distances in the N-terminus can be accounted for by modest structural changes, and no significant changes in structure were detected in the C-terminal portion of the channel upon activation of the channel. Amantadine, which blocks proton conductance by binding in the aqueous pore near the N-terminus, however, significantly modifies the tetrameric structure on the opposite side of the membrane. The interactions between the juxtamembrane amphipathic helix of one monomer and its neighboring monomer observed in the absence of drug are disrupted in its presence. However, the addition of cholesterol prevents this structural disruption. In fact, strong interactions are observed between cholesterol and residues in the amphipathic helix, accounting for cholesterol binding adjacent to a native palmitoylation site and near to an interhelix crevice that is typical of cholesterol binding sites. The resultant stabilization of the amphipathic helix deep in the bilayer interface facilitates the bilayer curvature that is essential for viral budding.

Abstract Icosahedral viral capsids are made of a large number of symmetrically organized protein subunits whose local movements can be essential for infection. In the capsid of the minute virus of mice, events required for infection that involve translocation of peptides through capsid pores are associated with a subtle conformational change. In vitro, this change can be reversibly induced by overcoming the energy barrier through mild heating of the capsid, but little is known about the capsid regions involved in the process. Here, we use hydrogen-deuterium exchange coupled to mass spectrometry to analyze the dynamics of the minute virus of mice capsid at increasing temperatures. Our results indicate that the transition associated with peptide translocation involves the structural rearrangement of regions distant from the capsid pores. These alterations are reflected in an increased dynamics of some secondary-structure elements in the capsid shell from which spikes protrude, and a decreased dynamics in the long intertwined loops that form the large capsid spikes. Thus, the translocation events through capsid pores involve a global conformational rearrangement of the capsid and a complex alteration of its equilibrium dynamics. This study additionally demonstrates the potential of hydrogen-deuterium exchange coupled to mass spectrometry to explore in detail temperature-dependent structural dynamics in large macromolecular protein assemblies. Most importantly, it paves the way for undertaking novel studies of the relationship between structure, dynamics, and biological function in virus particles and other large protein cages.

Abstract Interactions of the hydrophobic anticancer drug tamoxifen (TAM) with lipid model membranes were studied using calcein-encapsulated vesicle leakage, attenuated total reflection Fourier transform infrared (FTIR) spectroscopy, small-angle neutron scattering (SANS), atomic force microscopy (AFM) based force spectroscopy, and all-atom molecular dynamics (MD) simulations. The addition of TAM enhances membrane permeability, inducing calcein to translocate from the interior to the exterior of lipid vesicles. A large decrease in the FTIR absorption band’s magnitude was observed in the hydrocarbon chain region, suggesting suppressed bond vibrational dynamics. Bilayer thickening was determined from SANS data. Force spectroscopy measurements indicate that the lipid bilayer area compressibility modulus KA is increased by a large amount after the incorporation of TAM. MD simulations show that TAM decreases the lipid area and increases chain order parameters. Moreover, orientational and positional analyses show that TAM exhibits a highly dynamic conformation within the lipid bilayer. Our detailed experimental and computational studies of TAM interacting with model lipid membranes shed new light on membrane modulation by TAM.

Abstract Glucosylceramide (GlcCer), one of the simplest glycosphingolipids, plays key roles in physiology and pathophysiology. It has been suggested that GlcCer modulates cellular events by forming specialized domains. In this study, we investigated the interplay between GlcCer and cholesterol (Chol), an important lipid involved in the formation of liquid-ordered (lo) phases. Using fluorescence microscopy and spectroscopy, and dynamic and electrophoretic light scattering, we characterized the interaction between these lipids in different pH environments. A quantitative description of the phase behavior of the ternary unsaturated phospholipid/Chol/GlcCer mixture is presented. The results demonstrate coexistence between lo and liquid-disordered (ld) phases. However, the extent of lo/ld phase separation is sparse, mainly due to the ability of GlcCer to segregate into tightly packed gel domains. As a result, the phase diagram of these mixtures is characterized by an extensive three-phase coexistence region of fluid (ld-phospholipid enriched)/lo (Chol enriched)/gel (GlcCer enriched). Moreover, the results show that upon acidification, GlcCer solubility in the lo phase is increased, leading to a larger lo/ld coexistence region. Quantitative analyses allowed us to determine the differences in the composition of the phases at neutral and acidic pH. These results predict the impact of GlcCer on domain formation and membrane organization in complex biological membranes, and provide a background for unraveling the relationship between the biophysical properties of GlcCer and its biological action.

Abstract We studied the viscoelastic properties of homogeneous and inhomogeneous levan-DNA mixtures using optical tweezers and a rotational rheometer. Levan and DNA are important components of the extracellular matrix of bacterial biofilms. Their viscoelastic properties influence the mechanical as well as molecular-transport properties of biofilm. Both macro- and microrheology measurements in homogeneous levan-DNA mixtures revealed pseudoplastic behavior. When the concentration of DNA reached a critical value, levan started to aggregate, forming clusters of a few microns in size. Microrheology using optical tweezers enabled us to measure local viscoelastic properties within the clusters as well as in the DNA phase surrounding the levan aggregates. In phase-separated levan-DNA mixtures, the results of macro- and microrheology differed significantly. The local viscosity and elasticity of levan increased, whereas the local viscosity of DNA decreased. On the other hand, the results of bulk viscosity measurements suggest that levan clusters do not interact strongly with DNA. Upon treatment with DNase, levan aggregates dispersed. These results demonstrate the advantages of microrheological measurements compared to bulk viscoelastic measurements when the materials under investigation are complex and inhomogeneous, as is often the case in biological samples.

Abstract PDEδ is a small protein that binds and controls the trafficking of RAS subfamily proteins. Its inhibition protects initiation of RAS signaling, and it is one of the common targets considered for oncological drug development. In this study, we used solved x-ray structures of inhibitor-bound PDEδ targets to investigate mechanisms of action of six independent all-atom MD simulations. An analysis of atomic simulations combined with the molecular mechanic-Poisson-Boltzmann solvent accessible surface area/generalized Born solvent accessible surface area calculations led to the identification of action mechanisms for a panel of novel PDEδ inhibitors. To the best of our knowledge, this study is one of the first in silico investigations on co-crystallized PDEδ protein. A detailed atomic-scale understanding of the molecular mechanism of PDEδ inhibition may assist in the design of novel PDEδ inhibitors. One of the most common side effects for diverse small molecules/kinase inhibitors is their off-target interactions with cardiac ion channels and human-ether-a-go-go channel specifically. Thus, all of the studied PDEδ inhibitors are also screened in silico at the central cavities of hERG1 potassium channels.