Effects of the RNA-Polymerase Inhibitors Remdesivir and Favipiravir on the Structure of Lipid Bilayers-An MD Study.

The structure and dynamics of membranes are crucial to ensure the proper functioning of cells. There are some compounds used in therapeutics that show nonspecific interactions with membranes in addition to their specific molecular target. Among them, two compounds recently used in therapeutics against COVID-19, remdesivir and favipiravir, were subjected to molecular dynamics simulation assays. In these, we demonstrated that the compounds can spontaneously bind to model lipid membranes in the presence or absence of cholesterol. These findings correlate with the corresponding experimental results recently reported by our group. In conclusion, insertion of the compounds into the membrane is observed, with a mean position close to the phospholipid head groups.

Class III Peroxidases PRX01, PRX44, and PRX73 Control Root Hair Growth in Arabidopsis thaliana.

Root hair cells are important sensors of soil conditions. They grow towards and absorb water-soluble nutrients. This fast and oscillatory growth is mediated by continuous remodeling of the cell wall. Root hair cell walls contain polysaccharides and hydroxyproline-rich glycoproteins, including extensins (EXTs). Class-III peroxidases (PRXs) are secreted into the apoplastic space and are thought to trigger either cell wall loosening or polymerization of cell wall components, such as Tyr-mediated assembly of EXT networks (EXT-PRXs). The precise role of these EXT-PRXs is unknown. Using genetic, biochemical, and modeling approaches, we identified and characterized three root-hair-specific putative EXT-PRXs, PRX01, PRX44, and PRX73. prx01,44,73 triple mutation and PRX44 and PRX73 overexpression had opposite effects on root hair growth, peroxidase activity, and ROS production, with a clear impact on cell wall thickness. We use an EXT fluorescent reporter with contrasting levels of cell wall insolubilization in prx01,44,73 and PRX44-overexpressing background plants. In this study, we propose that PRX01, PRX44, and PRX73 control EXT-mediated cell wall properties during polar expansion of root hair cells.

Depletion of Mannose Receptor-Positive Tumor-associated Macrophages via a Peptide-targeted Star-shaped Polyglutamate Inhibits Breast Cancer Progression in Mice.

Although many studies have explored the depletion of tumor-associated macrophages (TAM) as a therapeutic strategy for solid tumors, currently available compounds suffer from poor efficacy and dose-limiting side effects. Here, we developed a novel TAM-depleting agent ("OximUNO") that specifically targets CD206+ TAMs and demonstrated efficacy in a triple-negative breast cancer (TNBC) mouse model. OximUNO comprises a star-shaped polyglutamate (St-PGA) decorated with the CD206-targeting peptide mUNO that carries the chemotherapeutic drug doxorubicin (DOX). In the TNBC model, a fluorescently labeled mUNO-decorated St-PGA homed to CD206+ TAMs within primary lesions and metastases. OximUNO exhibited no acute liver or kidney toxicity in vivo. Treatment with OximUNO reduced the progression of primary tumor lesions and pulmonary metastases, significantly diminished the number of CD206+ TAMs and increased the CD8/FOXP3 expression ratio (indicating immunomodulation). Our findings suggest the potential benefit of OximUNO as a TAM-depleting agent for TNBC treatment. Importantly, our studies also represent a novel design of a peptide-targeted St-PGA as a targeted therapeutic nanoconjugate. Significance: A peptide-targeted nanoformulation of DOX exclusively eliminates mannose receptor+ TAMs in breast cancer models, generating response without off-target effects (a drawback of many TAM-depleting agents under clinical study).

Functional asymmetry and chemical reactivity of CsoR family persulfide sensors.

CstR is a persulfide-sensing member of the functionally diverse copper-sensitive operon repressor (CsoR) superfamily. While CstR regulates the bacterial response to hydrogen sulfide (H2S) and more oxidized reactive sulfur species (RSS) in Gram-positive pathogens, other dithiol-containing CsoR proteins respond to host derived Cu(I) toxicity, sometimes in the same bacterial cytoplasm, but without regulatory crosstalk in cells. It is not clear what prevents this crosstalk, nor the extent to which RSS sensors exhibit specificity over other oxidants. Here, we report a sequence similarity network (SSN) analysis of the entire CsoR superfamily, which together with the first crystallographic structure of a CstR and comprehensive mass spectrometry-based kinetic profiling experiments, reveal new insights into the molecular basis of RSS specificity in CstRs. We find that the more N-terminal cysteine is the attacking Cys in CstR and is far more nucleophilic than in a CsoR. Moreover, our CstR crystal structure is markedly asymmetric and chemical reactivity experiments reveal the functional impact of this asymmetry. Substitution of the Asn wedge between the resolving and the attacking thiol with Ala significantly decreases asymmetry in the crystal structure and markedly impacts the distribution of species, despite adopting the same global structure as the parent repressor. Companion NMR, SAXS and molecular dynamics simulations reveal that the structural and functional asymmetry can be traced to fast internal dynamics of the tetramer. Furthermore, this asymmetry is preserved in all CstRs and with all oxidants tested, giving rise to markedly distinct distributions of crosslinked products. Our exploration of the sequence, structural, and kinetic features that determine oxidant-specificity suggest that the product distribution upon RSS exposure is determined by internal flexibility.

A Highly Conserved Iron-Sulfur Cluster Assembly Machinery between Humans and Amoeba Dictyostelium discoideum: The Characterization of Frataxin.

Several biological activities depend on iron-sulfur clusters ([Fe-S]). Even though they are well-known in several organisms their function and metabolic pathway were poorly understood in the majority of the organisms. We propose to use the amoeba Dictyostelium discoideum, as a biological model to study the biosynthesis of [Fe-S] at the molecular, cellular and organism levels. First, we have explored the D. discoideum genome looking for genes corresponding to the subunits that constitute the molecular machinery for Fe-S cluster assembly and, based on the structure of the mammalian supercomplex and amino acid conservation profiles, we inferred the full functionality of the amoeba machinery. After that, we expressed the recombinant mature form of D. discoideum frataxin protein (DdFXN), the kinetic activator of this pathway. We characterized the protein and its conformational stability. DdFXN is monomeric and compact. The analysis of the secondary structure content, calculated using the far-UV CD spectra, was compatible with the data expected for the FXN fold, and near-UV CD spectra were compatible with the data corresponding to a folded protein. In addition, Tryptophan fluorescence indicated that the emission occurs from an apolar environment. However, the conformation of DdFXN is significantly less stable than that of the human FXN, (4.0 vs. 9.0 kcal mol-1, respectively). Based on a sequence analysis and structural models of DdFXN, we investigated key residues involved in the interaction of DdFXN with the supercomplex and the effect of point mutations on the energetics of the DdFXN tertiary structure. More than 10 residues involved in Friedreich's Ataxia are conserved between the human and DdFXN forms, and a good correlation between mutational effect on the energetics of both proteins were found, suggesting the existence of similar sequence/function/stability relationships. Finally, we integrated this information in an evolutionary context which highlights particular variation patterns between amoeba and humans that may reflect a functional importance of specific protein positions. Moreover, the complete pathway obtained forms a piece of evidence in favor of the hypothesis of a shared and highly conserved [Fe-S] assembly machinery between Human and D. discoideum.

Relationship between activity and stability: Design and characterization of stable variants of human frataxin.

The relationships between conformational dynamics, stability and protein function are not obvious. Frataxin (FXN) is an essential protein that forms part of a supercomplex dedicated to the iron-sulfur (Fe-S) cluster assembly within the mitochondrial matrix. In humans, the loss of FXN expression or a decrease in its functionality results in Friedreich's Ataxia, a cardio-neurodegenerative disease. Recently, the way in which FXN interacts with the rest of the subunits of the supercomplex was uncovered. This opens a window to explore relationships between structural dynamics and function. In this study, we prepared a set of FXN variants spanning a broad range of conformational stabilities. Variants S160I, S160M and A204R were more stable than the wild-type and showed similar biological activity. Additionally, we prepared SILCAR, a variant that combines S160I, L203C and A204R mutations. SILCAR was 2.4 kcal mol-1 more stable and equally active. Some of the variants were significantly more resistant to proteolysis than the wild-type FXN. SILCAR showed the highest resistance, suggesting a more rigid structure. It was corroborated by means of molecular dynamics simulations. Relaxation dispersion NMR experiments comparing SILCAR and wild-type variants suggested similar internal motions in the microsecond to millisecond timescale. Instead, variant S157I showed higher denaturation resistance but a significant lower function, similarly to that observed for the FRDA variant N146K. We concluded that the contribution of particular side chains to the conformational stability of FXN might be highly subordinated to their impact on both the protein function and the stability of the functional supercomplex.

Hemeproteins as Targets for Sulfide Species.

Significance: Sulfides are endogenous and ubiquitous signaling species that share the hemeproteins as biochemical targets with O2, nitric oxide, and carbon monoxide. The description of the binding mechanisms is mandatory to anticipate the biochemical relevance of the interaction. Recent Advances: The binding of sulfide to ferric hemeproteins has been described in more than 40 systems, including native proteins, mutants, and model systems. Mechanisms of sulfide binding to ferric hemeproteins have been examined by a combination of kinetic and computational experiments. The distal control of the association process, dissected into the migration of the ligand to the active site and the binding event, reveals that neutral hydrogen sulfide (H2S) reaches the active site and is the predominant binding ligand, while the HS- is excluded by the protein matrix. Experiments with model compounds, devoid of a protein scaffold, reveal that both H2S and HS- can bind the ferric heme if accessing the site. A critical role of the proximal ligand in the prevention of the metal-centered reduction has been experimentally assessed. For metmyoglobin and methemoglobin, the coordination of sulfide leads to noncanonical functions: sulfide storage and its oxidative detoxification have been evidenced under physiological and excess sulfide concentrations, respectively. Critical Issues: The bound species is suggested to predominate in the monoprotonated form, although spectroscopic evidence is pending. Future Directions: A description of the role of hemeproteins as biochemical targets for inorganic sulfide requires understanding the reactivity of bound sulfide, for example: the metal-centered reduction, the reaction with excess sulfide, oxidants, or other gasotransmitters, among other biomolecules.

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations.

Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O2 and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.

Mechanism of Sulfide Binding by Ferric Hemeproteins.

The reaction of hydrogen sulfide (H2S) with hemeproteins is a key physiological reaction; still, its mechanism and implications are not completely understood. In this work, we propose a combination of experimental and theoretical tools to shed light on the reaction in model system microperoxidase 11 (MP11-FeIII) and myoglobin (Mb-FeIII), from the estimation of the intrinsic binding constants of the species H2S and hydrosulfide (HS-), and the computational description of the overall binding process. Our results show that H2S and HS- are the main reactive species in Mb-FeIII and MP11-FeIII, respectively, and that the magnitude of their intrinsic binding constants are similar to most of the binding constants reported so far for hemeproteins systems and model compounds. However, while the binding of HS- to Mb-FeIII was negligible, the binding of H2S to MP11-FeIII was significant, providing a frame for a discriminated analysis of both species and revealing differential mechanistic aspects. A joint inspection of the kinetic data and the free energy profiles of the binding processes suggests that a dissociative mechanism with the release of a coordinated water molecule as rate limiting step is operative in the binding of H2S to Mb-FeIII and that the binding of HS- is prevented in the access to the protein matrix. For the MP11-FeIII case, where no access restrictions for the ligands are present, an associative component in the mechanism seems to be operative. Overall, the results suggest that if accessing the active site then both H2S and HS- are capable of binding a ferric heme moiety.

Filling the Gaps to Solve the Extensin Puzzle.

Extensins (EXTs) are highly repetitive plant O-glycoproteins that require several post-translational modifications (PTMs) to become functional in plant cell walls. First, they are hydroxylated on contiguous proline residues; then they are O-glycosylated on hydroxyproline and serine. After secretion into the apoplast, O-glycosylated EXTs form a tridimensional network organized by inter- and intra-Tyr linkages. Recent studies have made significant progress in the identification of the enzymatic machinery required to process EXTs, which includes prolyl 4-hydroxylases, glycosyltransferases, papain-type cysteine endopeptidases, and peroxidases. EXTs are abundant in plant tissues and are particularly important in rapidly expanding root hairs and pollen tubes, which grow in a polar manner. Small changes in EXT PTMs affect fast-growing cells, although the molecular mechanisms underlying this regulation are unknown. In this review, we highlight recent advances in our understanding of EXT modifications throughout the secretory pathway, EXT assembly in cell walls, and possible sensing mechanisms involving the Catharanthus roseus cell surface sensor receptor-like kinases located at the interface between the apoplast and the cytoplasmic side of the plasma membrane.

Tertiary and quaternary structural basis of oxygen affinity in human hemoglobin as revealed by multiscale simulations.

Human hemoglobin (Hb) is a benchmark protein of structural biology that shaped our view of allosterism over 60 years ago, with the introduction of the MWC model based on Perutz structures of the oxy(R) and deoxy(T) states and the more recent Tertiary Two-State model that proposed the existence of individual subunit states -"r" and "t"-, whose structure is yet unknown. Cooperative oxygen binding is essential for Hb function, and despite decades of research there are still open questions related to how tertiary and quaternary changes regulate oxygen affinity. In the present work, we have determined the free energy profiles of oxygen migration and for HisE7 gate opening, with QM/MM calculations of the oxygen binding energy in order to address the influence of tertiary differences in the control of oxygen affinity. Our results show that in the α subunit the low to high affinity transition is achieved by a proximal effect that mostly affects oxygen dissociation and is the driving force of the allosteric transition, while in the ß subunit the affinity change results from a complex interplay of proximal and distal effects, including an increase in the HE7 gate opening, that as shown by free energy profiles promotes oxygen uptake.

Theoretical investigation of the mechanism of nitroxyl decomposition in aqueous solution.

Nitroxyl (HNO) is a species that has been proposed recently to play different roles in nitrosative stress processes. HNO decomposition in aqueous solution leading to N2O is a fast reaction that competes with many biochemical reactions in which HNO may be involved. Since molecular determinants of this reaction are still not fully understood, we present in this work an exhaustive analysis of the mechanism in terms of electronic-structure calculations as well as state of the art hybrid quantum mechanics/molecular mechanics molecular dynamics simulations. We characterized the reaction mechanism and computed free energy profiles for the reaction steps using an umbrella sampling procedure. We propose a first dimerization step followed by an acid-base equilibria. Afterwards, the product is formed from two main pathways involving cis-hyponitrous acid (cis-HONNOH) and its conjugate basis as intermediate. Our calculations show preference for the anionic pathway under physiological conditions and allow us to rationalize the results in terms of a molecular description of specific interactions with the solvent. These interactions turn out to be determinant in the stabilization of transition states and, thereby, modifying the free energy barriers. We predict a strong pH-dependence of the overall kinetics of N2O formation, related with the fraction of reactive species available in solution. Finally, we suggest experimental procedures which could validate this mechanism.

Coarse-Grained Simulations of Heme Proteins: Validation and Study of Large Conformational Transitions.

Heme proteins are ubiquitous in nature and perform many diverse functions in all kingdoms of life. Many of these functions are related to large-scale conformational transitions and allosteric processes. Sampling of these large conformational changes is computationally very challenging. In this context, coarse-grain simulations emerge as an efficient approach to explore the conformational landscape. In this work, we present a coarse-grained model of the heme group and thoroughly validate this model in different benchmark examples, which include the monomeric heme proteins myoglobin and neuroglobin and the tetrameric human hemoglobin where we evaluated the method's ability to explore conformational changes (as the formation of hexacoordinated species) and allosteric transitions (as the well-known R → T transition). The obtained results are compared with atomistic molecular dynamics simulations. Overall, the results indicate that this approach conserves the essential dynamical information on different allosteric processes.