Kv7.2 (KCNQ2) is the principal molecular component of the slow voltage gated M-channel, which strongly influences neuronal excitability. Calmodulin (CaM) binds to two intracellular C-terminal segments of Kv7.2 channels, helices A and B, and it is required for exit from the endoplasmic reticulum. However, the molecular mechanisms by which CaM controls channel trafficking are currently unknown. Here we used two complementary approaches to explore the molecular events underlying the association between CaM and Kv7.2 and their regulation by Ca(2+). First, we performed a fluorometric assay using dansylated calmodulin (D-CaM) to characterize the interaction of its individual lobes to the Kv7.2 CaM binding site (Q2AB). Second, we explored the association of Q2AB with CaM by NMR spectroscopy, using (15)N-labeled CaM as a reporter. The combined data highlight the interdependency of the N- and C-lobes of CaM in the interaction with Q2AB, suggesting that when CaM binds Ca(2+) the binding interface pivots between the N-lobe whose interactions are dominated by helix B and the C-lobe where the predominant interaction is with helix A. In addition, Ca(2+) makes CaM binding to Q2AB more difficult and, reciprocally, the channel weakens the association of CaM with Ca(2+).

The accepted model for the interaction of α and β integrins in the transmembrane (TM) domain is based on the pair αIIbβ3. This involves the so-called outer and inner membrane association clasps (OMC and IMC, respectively). In the α chain, the OMC involves a GxxxG-like motif, whereas in the IMC a conserved juxtamembrane GFFKR motif experiences a backbone reversal that partially fills the void generated by TM separation towards the cytoplasmic half. However, the GFFKR motif of several α integrin cytoplasmic tails in non-bicelle environments has been shown to adopt an α-helical structure that is not membrane-embedded and which was shown to bind a variety of cytoplasmic proteins. Thus it is not known if a membrane-embedded backbone reversal is a conserved structural feature in α integrins. We have studied the system αLβ2 because of its importance in leukocytes, where integrin deactivation is particularly important. Herein we show that the backbone reversal feature is not only present in αIIb but also in αL-TM when reconstituted in bicelles. Additionally, titration with β2 TM showed eight residues clustering along one side of αL-TM, forming a plausible interacting face with β2. The latter orientation is consistent with a previously predicted reported polar interaction between αL Ser-1071 and β2 Thr-686.

Voltage-gated Na⁺-channels are transmembrane proteins that are responsible for the fast depolarizing phase of the action potential in nerve and muscular cells. Selective permeability of Na⁺ over Ca²⁺ or K⁺ ions is essential for the biological function of Na⁺-channels. After the emergence of the first high-resolution structure of a Na⁺-channel, an anionic coordination site was proposed to confer Na⁺ selectivity through partial dehydration of Na⁺ via its direct interaction with conserved glutamate side chains. By combining molecular dynamics simulations and free-energy calculations, a low-energy permeation pathway for Na⁺ ion translocation through the selectivity filter of the recently determined crystal structure of a prokaryotic sodium channel from Arcobacter butzleri is characterised. The picture that emerges is that of a pore preferentially occupied by two ions, which can switch between different configurations by crossing low free-energy barriers. In contrast to K⁺-channels, the movements of the ions appear to be weakly coupled in Na⁺-channels. When the free-energy maps for Na⁺ and K⁺ ions are compared, a selective site is characterised in the narrowest region of the filter, where a hydrated Na⁺ ion, and not a hydrated K⁺ ion, is energetically stable.

Conditionally disordered proteins are either ordered or disordered depending on the environmental context. The substrates of the mitochondrial intermembrane space (IMS) oxidoreductase Mia40 are synthesized on cytosolic ribosomes and diffuse as intrinsically disordered proteins to the IMS, where they fold into their functional conformations; behaving thus as conditionally disordered proteins. It is not clear how the sequences of these polypeptides encode at the same time for their ability to adopt a folded structure and to remain unfolded. Here we characterize the disorder-to-order transition of a Mia40 substrate, the human small copper chaperone Cox17. Using an integrated real-time approach, including chromatography, fluorescence, CD, FTIR, SAXS, NMR, and MS analysis, we demonstrate that in this mitochondrial protein, the conformational switch between disordered and folded states is controlled by the formation of a single disulfide bond, both in the presence and in the absence of Mia40. We provide molecular details on how the folding of a conditionally disordered protein is tightly regulated in time and space, in such a way that the same sequence is competent for protein translocation and activity.

The interaction of human galectin-1 with a variety of oligosaccharides, from di-(N-acetyllactosamine) to tetra-saccharides (blood B type-II antigen) has been scrutinized by using a combined approach of different NMR experiments, molecular dynamics (MD) simulations, and isothermal titration calorimetry. Ligand- and receptor-based NMR experiments assisted by computational methods allowed proposing three-dimensional structures for the different complexes, which explained the lack of enthalpy gain when increasing the chemical complexity of the glycan. Interestingly, and independently of the glycan ligand, the entropy term does not oppose the binding event, a rather unusual feature for protein-sugar interactions. CLEANEX-PM and relaxation dispersion experiments revealed that sugar binding affected residues far from the binding site and described significant changes in the dynamics of the protein. In particular, motions in the microsecond-millisecond timescale in residues at the protein dimer interface were identified in the presence of high affinity ligands. The dynamic process was further explored by extensive MD simulations, which provided additional support for the existence of allostery in glycan recognition by human galectin-1.

Congenital erythropoietic porphyria (CEP) is an autosomal recessive disorder of the heme biosynthetic pathway that is characterized by uroporphyrinogen III synthase (UROS) deficiency and the accumulation of non-physiological isomer I porphyrins. These phototoxic metabolites predominantly produced by the erythron result in ineffective erythropoiesis, chronic hemolysis and splenomegaly, but they also disseminate in tissues causing bullous photosensitivity to UV light and skin fragility that may progress to scarring with photo mutilation. Therapeutic management is currently limited to supportive care and bone marrow transplantation is reserved for the most severe cases. We describe here a 26-year-old women previously diagnosed with CEP harbouring two novel UROS gene mutations whose pathogenic mechanism was investigated by extensive molecular analysis. Clinical features included disabling hypertrichosis and skin photosensitivity without hemolysis. The first and rate-limiting 5-aminolevulinate synthase 2 (ALAS2) enzyme controls heme synthesis and porphyrin production in erythroid cells, while iron availability modulates its expression through a post-transcriptional mechanism. We performed iterative phlebotomies over 26 months to induce iron depletion in the patient and investigated the effectiveness and tolerance of this cost-effective approach. We observed a progressive decrease in plasma ferritin and urinary porphyrins upon treatment without inducing anemia. The patient reported improved quality of life and photosensitivity. Our data confirm recent reports highlighting the benefit of iron restriction on the disease phenotype through a reduction in porphyrin accumulation. This new strategy may represent an efficient and well-tolerated treatment for CEP patients with skin involvement and limited hematological component if iron restriction is carefully monitored.

Improved methods are required for investigating the systemic metabolic effects of SARS-CoV-2 infection and patient stratification for precision treatment. We aimed to develop an effective method using lipid profiles for discriminating between SARS-CoV-2 infection, healthy controls, and non-SARS-CoV-2 respiratory infections. Targeted liquid chromatography-mass spectrometry lipid profiling was performed on discovery (20 SARS-CoV-2-positive; 37 healthy controls; 22 COVID-19 symptoms but SARS-CoV-2negative) and validation (312 SARS-CoV-2-positive; 100 healthy controls) cohorts. Orthogonal projection to latent structure-discriminant analysis (OPLS-DA) and Kruskal-Wallis tests were applied to establish discriminant lipids, significance, and effect size, followed by logistic regression to evaluate classification performance. OPLS-DA reported separation of SARS-CoV-2 infection from healthy controls in the discovery cohort, with an area under the curve (AUC) of 1.000. A refined panel of discriminant features consisted of six lipids from different subclasses (PE, PC, LPC, HCER, CER, and DCER). Logistic regression in the discovery cohort returned a training ROC AUC of 1.000 (sensitivity = 1.000, specificity = 1.000) and a test ROC AUC of 1.000. The validation cohort produced a training ROC AUC of 0.977 (sensitivity = 0.855, specificity = 0.948) and a test ROC AUC of 0.978 (sensitivity = 0.948, specificity = 0.922). The lipid panel was also able to differentiate SARS-CoV-2-positive individuals from SARS-CoV-2-negative individuals with COVID-19-like symptoms (specificity = 0.818). Lipid profiling and multivariate modelling revealed a signature offering mechanistic insights into SARS-CoV-2, with strong predictive power, and the potential to facilitate effective diagnosis and clinical management.

Fast C-type inactivation confers distinctive functional properties to the hERG potassium channel, and its association to inherited and acquired cardiac arrythmias makes the study of the inactivation mechanism of hERG at the atomic detail of paramount importance. At present, two models have been proposed to describe C-type inactivation in K+-channels. Experimental data and computational work on the bacterial KcsA channel support the hypothesis that C-type inactivation results from a closure of the selectivity filter that sterically impedes ion conduction. Alternatively, recent experimental structures of a mutated Shaker channel revealed a widening of the extracellular portion of the selectivity filter, which might diminish conductance by interfering with the mechanism of ion permeation. Here, we performed molecular dynamics simulations of the wild-type hERG, a non-inactivating mutant (hERG-N629D), and a mutant that inactivates faster than the wild-type channel (hERG-F627Y) to find out which and if any of the two reported C-type inactivation mechanisms applies to hERG. Closure events of the selectivity filter were not observed in any of the simulated trajectories but instead, the extracellular section of the selectivity filter deviated from the canonical conductive structure of potassium channels. The degree of widening of the potassium binding sites at the extracellular entrance of the channel was directly related to the degree of inactivation with hERG-F627Y > wild-type hERG > hERG-N629D. These findings support the hypothesis that C-type inactivation in hERG entails a widening of the extracellular entrance of the channel rather than a closure of the selectivity filter.

The pathogenesis of liver fibrosis requires activation of hepatic stellate cells (HSCs); once activated, HSCs lose intracellular fatty acids but the role of fatty acid oxidation and carnitine palmitoyltransferase 1A (CPT1A) in this process remains largely unexplored.

Transient receptor potential (TRP) ion channels are important pharmacological targets because of their role in the perception of pain, and so, understanding their chemical regulation is essential for the development of analgesic drugs. Among the currently known TRP channel chemical agonists, capsaicin, the active compound of chili pepper, is probably the most exhaustively studied. The availability of the three-dimensional structure of the vanilloid receptor 1 (TRPV1) has fueled computational studies revealing the molecular details of capsaicin binding modes. Although this is a significant step, a comprehensible binding mechanism or pathway is invaluable for targeting TRP channels in modern pharmacology. In the present work, free-energy and enhanced sampling techniques have been used to explore a possible membrane-mediated pathway for capsaicin to enter the TRPV1 binding pocket where capsaicin accesses the protein starting at the extracellular milieu through the outer leaflet and into its binding site in the protein. The main states visited along this route have been characterized and include (i) a bound state in agreement with the binding mode "head-down, tail-up" and (ii) an alternative state corresponding to a "head-up, tail-down" binding mode. In agreement with previous reports, binding is mediated by both hydrogen bonds and van der Waals interactions, and residue Y511 is crucial for stabilizing the bound state and during the binding process. Together, these results provide a foundation to further understand TRPV channels, and they could be used to guide therapeutic design of selective inhibitors potentially leading to novel avenues for pharmacological applications targeting the TRPV1 channel.

Potassium channels are a diverse family of integral membrane proteins through which K(+) can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K(+) channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K(+) concentration is thought to be representative of a K(+) channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K(+) KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K(+) KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K(+) structure of KcsA and KirBac. The low-K(+) KcsA structure is not permeable by Na(+), K(+), or Rb(+), and the selectivity of its binding sites is different from that of the high-K(+) structure.

The lactose repressor protein may bind DNA in two possible configurations: a specific one, if the DNA sequence corresponds to a binding site, and a non-specific one otherwise. To find its target sequences, the lactose repressor first binds non-specifically to DNA, and subsequently, it rapidly searches for a binding site. Atomic structures of non-specific and specific complexes are available from crystallographic and nuclear magnetic resonance experiments. However, what remains unknown is a detailed description of the steps that transform the non-specific complex into the specific one. Here, how the protein first recognizes its binding site has been studied using molecular dynamics simulations. The picture that emerges is that of a protein that is as mobile when interacting with non-specific DNA sequences as when free in solution. This high degree of mobility allows the protein to rapidly sample different DNA sequences. In contrast, when the protein encounters a binding site, the configuration ensemble collapses, and the protein sliding movements along the DNA sequence become scarce. The binding energies in the specific and non-specific complexes were analysed using the Molecular Mechanics Poisson Boltzmann Surface Area approach. These results represent a first step towards a throughout characterization of the DNA-recognition process.

Interactions of Na(+), K(+), Rb(+), and Cs(+) ions within the selectivity filter of a potassium channel have been investigated via multiple molecular dynamics simulations (total simulation time, 48 ns) based on the high resolution structure of KcsA, embedded in a phospholipid bilayer. As in simulations based on a lower resolution structure of KcsA, concerted motions of ions and water within the filter are seen. Despite the use of a higher resolution structure and the inclusion of four buried water molecules thought to stabilize the filter, this region exhibits a significant degree of flexibility. In particular, pronounced distortion of filter occurs if no ions are present within it. The two most readily permeant ions, K(+) and Rb(+), are similar in their interactions with the selectivity filter. In contrast, Na(+) ions tend to distort the filter by binding to a ring of four carbonyl oxygens. The larger Cs(+) ions result in a small degree of expansion of the filter relative to the x-ray structure. Cs(+) ions also appear to interact differently with the gate region of the channel, showing some tendency to bind within a predominantly hydrophobic pocket. The four water molecules buried between the back of the selectivity filter and the remainder of the protein show comparable mobility to the surrounding protein and do not exchange with water molecules within the filter or the central cavity. A preliminary comparison of the use of particle mesh Ewald versus cutoff protocols for the treatment of long-range electrostatics suggests some difference in the kinetics of ion translocation within the filter.

We have designed, synthesised and characterised the effects of a number of novel anthraquinone derivatives and assessed their effects on large conductance, Ca(2+) activated K(+) (BK) channels recorded from rabbit bladder smooth muscle cells using the excised, inside/out configuration of the patch clamp technique. These compounds are members of the GoSlo-SR family of compounds, which potently open BK channels and shift the voltage required for half maximal activation (V1/2) negatively. The efficacy of the anilinoanthraquinone derivatives was enhanced when the size of ring D was increased, since the cyclopentane and cyclohexane derivatives shifted the V1/2, by -24 ± 6 mV and -54 ± 8 mV, respectively, whereas the cycloheptane and cyclooctane derivatives shifted the V1/2 by -61 ± 6 mV and -106 ± 6 mV. To examine if a combination of hydrophobicity and steric bulking of this region further enhanced their ability to open BK channels, we synthesised a number of naphthalene and tetrahydro-naphthalene derivatives. The tetrahydro-2-naphthalene derivative GoSlo-SR-5-69 was the most potent and efficacious of the series since it was able to shift the activation V1/2 by greater than -100 mV when applied at a concentration of 1 μM and had an EC50 of 251 nM, making it one of the most potent and efficacious BK channel openers synthesised to date.

The response to hypoxia in animals involves the expression of multiple genes regulated by the αβ-hypoxia-inducible transcription factors (HIFs). The hypoxia-sensing mechanism involves oxygen limited hydroxylation of prolyl residues in the N- and C-terminal oxygen-dependent degradation domains (NODD and CODD) of HIFα isoforms, as catalysed by prolyl hydroxylases (PHD 1-3). Prolyl hydroxylation promotes binding of HIFα to the von Hippel-Lindau protein (VHL)-elongin B/C complex, thus signalling for proteosomal degradation of HIFα. We reveal that certain PHD2 variants linked to familial erythrocytosis and cancer are highly selective for CODD or NODD. Crystalline and solution state studies coupled to kinetic and cellular analyses reveal how wild-type and variant PHDs achieve ODD selectivity via different dynamic interactions involving loop and C-terminal regions. The results inform on how HIF target gene selectivity is achieved and will be of use in developing selective PHD inhibitors.

Congenital erythropoietic porphyria (CEP), also known as Günther's disease, results from a deficient activity in the fourth enzyme, uroporphyrinogen III synthase (UROIIIS), of the heme pathway. Ciclopirox (CPX) is an off-label drug, topically prescribed as an antifungal. It has been recently shown that it also acts as a pharmacological chaperone in CEP, presenting a specific activity in deleterious mutations in UROIIIS. Despite CPX is active at subtoxic concentrations, acute gastrointestinal (GI) toxicity was found due to the precipitation in the stomach of the active compound and subsequent accumulation in the intestine. To increase its systemic availability, we carried out pharmacokinetic (PK) and pharmacodynamic (PD) studies using alternative formulations for CPX. Such strategy effectively suppressed GI toxicity in WT mice and in a mouse model of the CEP disease (UROIIISP248Q/P248Q). In terms of activity, phosphorylation of CPX yielded good results in CEP cellular models but showed limited activity when administered to the CEP mouse model. These results highlight the need of a proper formulation for pharmacological chaperones used in the treatment of rare diseases.

Liver fibrosis must be evaluated in patients with hepatitis C virus (HCV) after liver transplantation because its severity affects their prognosis and the recurrence of HCV. Since invasive biopsy is still the gold standard to identify patients at risk of graft loss from rapid fibrosis progression, it becomes crucial the development of new accurate, non-invasive methods that allow repetitive examination of the patients. Therefore, we have developed a non-invasive, accurate model to distinguish those patients with different liver fibrosis stages. Two hundred and three patients with HCV were histologically classified (METAVIR) into five categories of fibrosis one year after liver transplantation. In this cross-sectional study, patients at fibrosis stages F0-F1 (n = 134) were categorised as "slow fibrosers" and F2-F4 (n = 69) as "rapid fibrosers". Chloroform/methanol serum extracts were analysed by reverse ultra-high performance liquid chromatography coupled to mass spectrometry. A diagnostic model was built through linear discriminant analyses. An algorithm consisting of two sphingomyelins and two phosphatidylcholines accurately classifies rapid and slow fibrosers after transplantation. The proposed model yielded an AUROC of 0.92, 71% sensitivity, 85% specificity, and 84% accuracy. Moreover, specific bile acids and sphingomyelins increased notably along with liver fibrosis severity, differentiating between rapid and slow fibrosers.

Fluorine (19 F) incorporation into glycan-binding proteins (lectins) has been achieved and exploited to monitor the binding to carbohydrate ligands by nuclear magnetic resonance (NMR) spectroscopy. Galectins are a family of lectins that bind carbohydrates, generally with weak affinities, through a combination of intermolecular interactions including a key CH-π stacking involving a conserved tryptophan residue. Herein, Galectin-3 (Gal3) and Galectin-8 (Gal8) with one and two carbohydrate recognition domains (CRDs), respectively, were selected. Gal3 contains one Trp, whereas Gal8 contains three, one at each binding site and a third one not involved in sugar binding; these were substituted by the corresponding F-Trp analogues. The presence of fluorine did not significantly modify the affinity for glycan binding, which was in slow exchange on the 19 F NMR chemical-shift timescale, even for weak ligands, and allowed binding events taking place at two different binding sites within the same lectin to be individualized.

Mesoamerican nephropathy (MeN) is a form of chronic kidney disease found predominantly in young men in Mesoamerica. Strenuous agricultural labor is a consistent risk factor for MeN, but the pathophysiologic mechanism leading to disease is poorly understood. We compared the urine metabolome among men in Nicaragua engaged in sugarcane harvest and seed cutting (n = 117), a group at high risk for MeN, against three referents: Nicaraguans working less strenuous jobs at the same sugarcane plantations (n = 78); Nicaraguans performing non-agricultural work (n = 102); and agricultural workers in Spain (n = 78). Using proton nuclear magnetic resonance, we identified 136 metabolites among participants. Our non-hypothesis-based approach identified distinguishing urine metabolic features in the high-risk group, revealing increased levels of hippurate and other gut-derived metabolites and decreased metabolites related to central energy metabolism when compared to referent groups. Our complementary hypothesis-based approach, focused on nicotinamide adenine dinucleotide (NAD+) related metabolites, and revealed a higher kynurenate/tryptophan ratio in the high-risk group (p = 0.001), consistent with a heightened inflammatory state. Workers in high-risk occupations are distinguishable by urinary metabolic features that suggest increased gut permeability, inflammation, and altered energy metabolism. Further study is needed to explore the pathophysiologic implications of these findings.

Neuronal KV7 channels, important regulators of cell excitability, are among the most sensitive proteins to reactive oxygen species. The S2S3 linker of the voltage sensor was reported as a site-mediating redox modulation of the channels. Recent structural insights reveal potential interactions between this linker and the Ca2+-binding loop of the third EF-hand of calmodulin (CaM), which embraces an antiparallel fork formed by the C-terminal helices A and B, constituting the calcium responsive domain (CRD). We found that precluding Ca2+ binding to the EF3 hand, but not to EF1, EF2, or EF4 hands, abolishes oxidation-induced enhancement of KV7.4 currents. Monitoring FRET (Fluorescence Resonance Energy Transfer) between helices A and B using purified CRDs tagged with fluorescent proteins, we observed that S2S3 peptides cause a reversal of the signal in the presence of Ca2+ but have no effect in the absence of this cation or if the peptide is oxidized. The capacity of loading EF3 with Ca2+ is essential for this reversal of the FRET signal, whereas the consequences of obliterating Ca2+ binding to EF1, EF2, or EF4 are negligible. Furthermore, we show that EF3 is critical for translating Ca2+ signals to reorient the AB fork. Our data are consistent with the proposal that oxidation of cysteine residues in the S2S3 loop relieves KV7 channels from a constitutive inhibition imposed by interactions between the EF3 hand of CaM which is crucial for this signaling.