Structural and biochemical characterization of Arabidopsis alcohol dehydrogenases reveals distinct functional properties but similar redox sensitivity.

Alcohol dehydrogenases (ADHs) are a group of zinc-binding enzymes belonging to the medium-length dehydrogenase/reductase (MDR) protein superfamily. In plants, these enzymes fulfill important functions involving the reduction of toxic aldehydes to the corresponding alcohols (as well as catalyzing the reverse reaction, i.e., alcohol oxidation; ADH1) and the reduction of nitrosoglutathione (GSNO; ADH2/GSNOR). We investigated and compared the structural and biochemical properties of ADH1 and GSNOR from Arabidopsis thaliana. We expressed and purified ADH1 and GSNOR and determined two new structures, NADH-ADH1 and apo-GSNOR, thus completing the structural landscape of Arabidopsis ADHs in both apo- and holo-forms. A structural comparison of these Arabidopsis ADHs revealed a high sequence conservation (59% identity) and a similar fold. In contrast, a striking dissimilarity was observed in the catalytic cavity supporting substrate specificity and accommodation. Consistently, ADH1 and GSNOR showed strict specificity for their substrates (ethanol and GSNO, respectively), although both enzymes had the ability to oxidize long-chain alcohols, with ADH1 performing better than GSNOR. Both enzymes contain a high number of cysteines (12 and 15 out of 379 residues for ADH1 and GSNOR, respectively) and showed a significant and similar responsivity to thiol-oxidizing agents, indicating that redox modifications may constitute a mechanism for controlling enzyme activity under both optimal growth and stress conditions.

Variants in Human ATP Synthase Mitochondrial Genes: Biochemical Dysfunctions, Associated Diseases, and Therapies.

Mitochondrial ATP synthase (Complex V) catalyzes the last step of oxidative phosphorylation and provides most of the energy (ATP) required by human cells. The mitochondrial genes MT-ATP6 and MT-ATP8 encode two subunits of the multi-subunit Complex V. Since the discovery of the first MT-ATP6 variant in the year 1990 as the cause of Neuropathy, Ataxia, and Retinitis Pigmentosa (NARP) syndrome, a large and continuously increasing number of inborn variants in the MT-ATP6 and MT-ATP8 genes have been identified as pathogenic. Variants in these genes correlate with various clinical phenotypes, which include several neurodegenerative and multisystemic disorders. In the present review, we report the pathogenic variants in mitochondrial ATP synthase genes and highlight the molecular mechanisms underlying ATP synthase deficiency that promote biochemical dysfunctions. We discuss the possible structural changes induced by the most common variants found in patients by considering the recent cryo-electron microscopy structure of human ATP synthase. Finally, we provide the state-of-the-art of all therapeutic proposals reported in the literature, including drug interventions targeting mitochondrial dysfunctions, allotopic gene expression- and nuclease-based strategies, and discuss their potential translation into clinical trials.

Conformational response to ligand binding of TMPRSS2, a protease involved in SARS-CoV-2 infection: Insights through computational modeling.

Thanks to the considerable research which has been undertaken in the last few years to improve our understanding of the biology and mechanism of action of SARS-CoV-2, we know how the virus uses its surface spike protein to infect host cells. The transmembrane prosthesis, serine 2 (TMPRSS2) protein, located on the surface of human cells, recognizes the cleavage site in the spike protein, leading to the release of the fusion peptide and entry of the virus into the host cells. Because of its role, TMPRSS2 has been proposed as a drug target to prevent infection by the virus. In this study, we aim to increase our understanding of TMPRSS2 using long scale microsecond atomistic molecular dynamics simulations, focusing on the conformational changes over time. The comparison between simulations conducted on the protein in the native (apo) and inhibited form (holo), has shown that in the holo form the inhibitor stabilizes the catalytic site and induces rearrangements in the extracellular domain of the protein. In turn, it leads to the formation of a new cavity in the vicinity of the ligand binding pocket that is stable in the microsecond time scale. Given the low specificity of known protease inhibitors, these findings suggest a new potential drug target site that can be used to improve TMPRSS2 specific recognition by newly designed inhibitors.

Design, Pharmacological Characterization, and Molecular Docking of Minimalist Peptidomimetic Antagonists of α4ß1 Integrin.

Integrin receptors mediate cell-cell interactions via the recognition of cell-adhesion glycoproteins, as well as via the interactions of cells with proteins of the extracellular matrix, and upon activation they transduce signals bi-directionally across the cell membrane. In the case of injury, infection, or inflammation, integrins of ß2 and α4 families participate in the recruitment of leukocytes, a multi-step process initiated by the capturing of rolling leukocytes and terminated by their extravasation. In particular, α4ß1 integrin is deeply involved in leukocyte firm adhesion preceding extravasation. Besides its well-known role in inflammatory diseases, α4ß1 integrin is also involved in cancer, being expressed in various tumors and showing an important role in cancer formation and spreading. Hence, targeting this integrin represents an opportunity for the treatment of inflammatory disorders, some autoimmune diseases, and cancer. In this context, taking inspiration from the recognition motives of α4ß1 integrin with its natural ligands FN and VCAM-1, we designed minimalist α/ß hybrid peptide ligands, with our approach being associated with a retro strategy. These modifications are expected to improve the compounds' stability and bioavailability. As it turned out, some of the ligands were found to be antagonists, being able to inhibit the adhesion of integrin-expressing cells to plates coated with the natural ligands without inducing any conformational switch and any activation of intracellular signaling pathways. An original model structure of the receptor was generated using protein-protein docking to evaluate the bioactive conformations of the antagonists via molecular docking. Since the experimental structure of α4ß1 integrin is still unknown, the simulations might also shed light on the interactions between the receptor and its native protein ligands.

Determination and Kinetic Characterization of a New Potential Inhibitor for AmsI Protein Tyrosine Phosphatase from the Apple Pathogen Erwinia amylovora.

Erwinia amylovora is a Gram-negative bacterium, responsible for the fire blight disease in Rosaceae plants. Its virulence is correlated with the production of an exopolysaccharide (EPS) called amylovoran, which protects the bacterium from the surrounding environment and helps its diffusion inside the host. Amylovoran biosynthesis relies on the expression of twelve genes clustered in the ams operon. One of these genes, amsI, encodes for a Low Molecular Weight Protein Tyrosine Phosphatase (LMW-PTP) called EaAmsI, which plays a key role in the regulation of the EPS production pathway. For this reason, EaAmsI was chosen in this work as a target for the development of new antibacterial agents against E. amylovora. To achieve this aim, a set of programs (DOCK6, OpenEye FRED) was selected to perform a virtual screening using a database of ca. 700 molecules. The six best-scoring compounds identified were tested in in vitro assays. A complete inhibition kinetic characterization carried out on the most promising molecule (n-Heptyl ß-D-glucopyranoside, N7G) showed an inhibition constant of 7.8 ± 0.6 µM. This study represents an initial step towards the development of new EaAmsI inhibitors able to act as antibacterial agents against E. amylovora infections.

DNMT1 mutations leading to neurodegeneration paradoxically reflect on mitochondrial metabolism.

ADCA-DN and HSN-IE are rare neurodegenerative syndromes caused by dominant mutations in the replication foci targeting sequence (RFTS) of the DNA methyltransferase 1 (DNMT1) gene. Both phenotypes resemble mitochondrial disorders, and mitochondrial dysfunction was first observed in ADCA-DN. To explore mitochondrial involvement, we studied the effects of DNMT1 mutations in fibroblasts from four ADCA-DN and two HSN-IE patients. We documented impaired activity of purified DNMT1 mutant proteins, which in fibroblasts results in increased DNMT1 amount. We demonstrated that DNMT1 is not localized within mitochondria, but it is associated with the mitochondrial outer membrane. Concordantly, mitochondrial DNA failed to show meaningful CpG methylation. Strikingly, we found activated mitobiogenesis and OXPHOS with significant increase of H2O2, sharply contrasting with a reduced ATP content. Metabolomics profiling of mutant cells highlighted purine, arginine/urea cycle and glutamate metabolisms as the most consistently altered pathways, similar to primary mitochondrial diseases. The most severe mutations showed activation of energy shortage AMPK-dependent sensing, leading to mTORC1 inhibition. We propose that DNMT1 RFTS mutations deregulate metabolism lowering ATP levels, as a result of increased purine catabolism and urea cycle pathways. This is associated with a paradoxical mitochondrial hyper-function and increased oxidative stress, possibly resulting in neurodegeneration in non-dividing cells.

Computational Study of Helicase from SARS-CoV-2 in RNA-Free and Engaged Form.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain-domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains.

PacDOCK: A Web Server for Positional Distance-Based and Interaction-Based Analysis of Docking Results.

Molecular docking is a key method for structure-based drug design used to predict the conformations assumed by small drug-like ligands when bound to their target. However, the evaluation of molecular docking studies can be hampered by the lack of a free and easy to use platform for the complete analysis of results obtained by the principal docking programs. To this aim, we developed PacDOCK, a freely available and user-friendly web server that comprises a collection of tools for positional distance-based and interaction-based analysis of docking results, which can be provided in several file formats. PacDOCK allows a complete analysis of molecular docking results through root mean square deviation (RMSD) calculation, molecular visualization, and cluster analysis of docked poses. The RMSD calculation compares docked structures with a reference structure, also when atoms are randomly labelled, and their conformational and positional differences can be visualised. In addition, it is possible to visualise a ligand into the target binding pocket and investigate the key receptor-ligand interactions. Moreover, PacDOCK enables the clustering of docking results by identifying a restrained number of clusters from many docked poses. We believe that PacDOCK will contribute to facilitating the analysis of docking results to improve the efficiency of computer-aided drug design.

New Insights on Rotenone Resistance of Complex I Induced by the m.11778G>A/MT-ND4 Mutation Associated with Leber's Hereditary Optic Neuropathy.

The finding that the most common mitochondrial DNA mutation m.11778G>A/MT-ND4 (p.R340H) associated with Leber's hereditary optic neuropathy (LHON) induces rotenone resistance has produced a long-standing debate, because it contrasts structural evidence showing that the ND4 subunit is far away from the quinone-reaction site in complex I, where rotenone acts. However, recent cryo-electron microscopy data revealed that rotenone also binds to the ND4 subunit. We investigated the possible structural modifications induced by the LHON mutation and found that its amino acid replacement would disrupt a possible hydrogen bond between native R340 and Q139 in ND4, thereby destabilizing rotenone binding. Our analysis thus explains rotenone resistance in LHON patients as a biochemical signature of its pathogenic effect on complex I.

Definition of the Binding Architecture to a Target Promoter of HP1043, the Essential Master Regulator of Helicobacter pylori.

HP1043 is an essential orphan response regulator of Helicobacter pylori orchestrating multiple crucial cellular processes. Classified as a member of the OmpR/PhoB family of two-component systems, HP1043 exhibits a highly degenerate receiver domain and evolved to function independently of phosphorylation. Here, we investigated the HP1043 binding mode to a target sequence in the hp1227 promoter (Php1227). Scanning mutagenesis of HP1043 DNA-binding domain and consensus sequence led to the identification of residues relevant for the interaction of the protein with a target DNA. These determinants were used as restraints to guide a data-driven protein-DNA docking. Results suggested that, differently from most other response regulators of the same family, HP1043 binds in a head-to-head conformation to the Php1227 target promoter. HP1043 interacts with DNA largely through charged residues and contacts with both major and minor grooves of the DNA are required for a stable binding. Computational alanine scanning on molecular dynamics trajectory was performed to corroborate our findings. Additionally, in vitro transcription assays confirmed that HP1043 positively stimulates the activity of RNA polymerase.

Combining Different Docking Engines and Consensus Strategies to Design and Validate Optimized Virtual Screening Protocols for the SARS-CoV-2 3CL Protease.

The 3CL-Protease appears to be a very promising medicinal target to develop anti-SARS-CoV-2 agents. The availability of resolved structures allows structure-based computational approaches to be carried out even though the lack of known inhibitors prevents a proper validation of the performed simulations. The innovative idea of the study is to exploit known inhibitors of SARS-CoV 3CL-Pro as a training set to perform and validate multiple virtual screening campaigns. Docking simulations using four different programs (Fred, Glide, LiGen, and PLANTS) were performed investigating the role of both multiple binding modes (by binding space) and multiple isomers/states (by developing the corresponding isomeric space). The computed docking scores were used to develop consensus models, which allow an in-depth comparison of the resulting performances. On average, the reached performances revealed the different sensitivity to isomeric differences and multiple binding modes between the four docking engines. In detail, Glide and LiGen are the tools that best benefit from isomeric and binding space, respectively, while Fred is the most insensitive program. The obtained results emphasize the fruitful role of combining various docking tools to optimize the predictive performances. Taken together, the performed simulations allowed the rational development of highly performing virtual screening workflows, which could be further optimized by considering different 3CL-Pro structures and, more importantly, by including true SARS-CoV-2 3CL-Pro inhibitors (as learning set) when available.

Inhibition of Urease, a Ni-Enzyme: The Reactivity of a Key Thiol With Mono- and Di-Substituted Catechols Elucidated by Kinetic, Structural, and Theoretical Studies.

The inhibition of urease from Sporosarcina pasteurii (SPU) and Canavalia ensiformis (jack bean, JBU) by a class of six aromatic poly-hydroxylated molecules, namely mono- and dimethyl-substituted catechols, was investigated on the basis of the inhibitory efficiency of the catechol scaffold. The aim was to probe the key step of a mechanism proposed for the inhibition of SPU by catechol, namely the sulfanyl radical attack on the aromatic ring, as well as to obtain critical information on the effect of substituents of the catechol aromatic ring on the inhibition efficacy of its derivatives. The crystal structures of all six SPU-inhibitors complexes, determined at high resolution, as well as kinetic data obtained on JBU and theoretical studies of the reaction mechanism using quantum mechanical calculations, revealed the occurrence of an irreversible inactivation of urease by means of a radical-based autocatalytic multistep mechanism, and indicate that, among all tested catechols, the mono-substituted 3-methyl-catechol is the most efficient inhibitor for urease.

The structure-based reaction mechanism of urease, a nickel dependent enzyme: tale of a long debate.

This review is an attempt to retrace the chronicle that starts from the discovery of the role of nickel as the essential metal ion in urease for the enzymatic catalysis of urea, a key step in the biogeochemical cycle of nitrogen on Earth, to the most recent progress in understanding the chemistry of this historical enzyme. Data and facts are presented through the magnifying lenses of the authors, using their best judgment to filter and elaborate on the many facets of the research carried out on this metalloenzyme over the years. The tale is divided in chapters that discuss and describe the results obtained in the subsequent leaps in the knowledge that led from the discovery of a biological role for Ni to the most recent advancements in the comprehension of the relationship between the structure and function of urease. This review is intended not only to focus on the bioinorganic chemistry of this beautiful metal-based catalysis, but also, and maybe primarily, to evoke inspiration and motivation to further explore the realm of bio-based coordination chemistry.

The model structure of the copper-dependent ammonia monooxygenase.

Ammonia monooxygenase is a copper-dependent membrane-bound enzyme that catalyzes the first step of nitrification in ammonia-oxidizing bacteria to convert ammonia to hydroxylamine, through the reductive insertion of a dioxygen-derived O atom in an N-H bond. This reaction is analogous to that carried out by particulate methane monooxygenase, which catalyzes the conversion of methane to methanol. The enzymatic activity of ammonia monooxygenase must be modulated to reduce the release of nitrogen-based soil nutrients for crop production into the atmosphere or underground waters, a phenomenon known to significantly decrease the efficiency of primary production as well as increase air and water pollution. The structure of ammonia monooxygenase is not available, rendering the rational design of enzyme inhibitors impossible. This study describes a successful attempt to build a structural model of ammonia monooxygenase, and its accessory proteins AmoD and AmoE, from Nitrosomonas europaea, taking advantage of the high sequence similarity with particulate methane monooxygenase and the homologous PmoD protein, for which crystal structures are instead available. The results obtained not only provide the structural details of the proteins ternary and quaternary structures, but also suggest a location for the copper-containing active site for both ammonia and methane monooxygenases, as well as support a proposed structure of a CuA-analogue dinuclear copper site in AmoD and PmoD.

Targeting the Protein Tunnels of the Urease Accessory Complex: A Theoretical Investigation.

Urease is a nickel-containing enzyme that is essential for the survival of several and often deadly pathogenic bacterial strains, including Helicobacter pylori. Notwithstanding several attempts, the development of direct urease inhibitors without side effects for the human host remains, to date, elusive. The recently solved X-ray structure of the HpUreDFG accessory complex involved in the activation of urease opens new perspectives for structure-based drug discovery. In particular, the quaternary assembly and the presence of internal tunnels for nickel translocation offer an intriguing possibility to target the HpUreDFG complex in the search of indirect urease inhibitors. In this work, we adopted a theoretical framework to investigate such a hypothesis. Specifically, we searched for putative binding sites located at the protein-protein interfaces on the HpUreDFG complex, and we challenged their druggability through structure-based virtual screening. We show that, by virtue of the presence of tunnels, some protein-protein interfaces on the HpUreDFG complex are intrinsically well suited for hosting small molecules, and, as such, they possess good potential for future drug design endeavors.

A Solvent-Exposed Cysteine Forms a Peculiar NiII -Binding Site in the Metallochaperone CooT from Rhodospirillum rubrum.

In Rhodospirillum rubrum, the maturation of carbon monoxide dehydrogenase (CODH) requires three nickel chaperones, namely RrCooC, RrCooT and RrCooJ. Recently, the biophysical characterisation of the RrCooT homodimer and the X-ray structure of its apo form revealed the existence of a solvent-exposed NiII -binding site at the dimer interface, involving the strictly conserved Cys2. Here, a multifaceted approach that used NMR and X-ray absorption spectroscopies, complemented with structural bio-modelling methodologies, was used to characterise the binding mode of NiII in RrCooT. This study suggests that NiII adopts a square-planar geometry through a N2 S2 coordinating environment that comprises the two thiolate and amidate groups of both Cys2 residues at the dimer interface. The existence of a diamagnetic mononuclear NiII centre with bis-amidate/bis-thiolate ligands, coordinated by a single-cysteine motif, is unprecedented in biology and raises the question of its role in the activation of CODH at the molecular level.

Predicting ligand binding poses for low-resolution membrane protein models: Perspectives from multiscale simulations.

Membrane receptors constitute major targets for pharmaceutical intervention. Drug design efforts rely on the identification of ligand binding poses. However, the limited experimental structural information available may make this extremely challenging, especially when only low-resolution homology models are accessible. In these cases, the predictions may be improved by molecular dynamics simulation approaches. Here we review recent developments of multiscale, hybrid molecular mechanics/coarse-grained (MM/CG) methods applied to membrane proteins. In particular, we focus on our in-house MM/CG approach. It is especially tailored for G-protein coupled receptors, the largest membrane receptor family in humans. We show that our MM/CG approach is able to capture the atomistic details of the receptor/ligand binding interactions, while keeping the computational cost low by representing the protein frame and the membrane environment in a highly simplified manner. We close this review by discussing ongoing improvements and challenges of the current implementation of our MM/CG code.

Structure and dynamics of Helicobacter pylori nickel-chaperone HypA: an integrated approach using NMR spectroscopy, functional assays and computational tools.

Helicobacter pylori HypA (HpHypA) is a metallochaperone necessary for maturation of [Ni,Fe]-hydrogenase and urease, the enzymes required for colonization and survival of H. pylori in the gastric mucosa. HpHypA contains a structural Zn(II) site and a unique Ni(II) binding site at the N-terminus. X-ray absorption spectra suggested that the Zn(II) coordination depends on pH and on the presence of Ni(II). This study was performed to investigate the structural properties of HpHypA as a function of pH and Ni(II) binding, using NMR spectroscopy combined with DFT and molecular dynamics calculations. The solution structure of apo,Zn-HpHypA, containing Zn(II) but devoid of Ni(II), was determined using 2D, 3D and 4D NMR spectroscopy. The structure suggests that a Ni-binding and a Zn-binding domain, joined through a short linker, could undergo mutual reorientation. This flexibility has no physiological effect on acid viability or urease maturation in H. pylori. Atomistic molecular dynamics simulations suggest that Ni(II) binding is important for the conformational stability of the N-terminal helix. NMR chemical shift perturbation analysis indicates that no structural changes occur in the Zn-binding domain upon addition of Ni(II) in the pH 6.3-7.2 range. The structure of the Ni(II) binding site was probed using 1H NMR spectroscopy experiments tailored to reveal hyperfine-shifted signals around the paramagnetic metal ion. On this basis, two possible models were derived using quantum-mechanical DFT calculations. The results provide a comprehensive picture of the Ni(II) mode to HpHypA, important to rationalize, at the molecular level, the functional interactions of this chaperone with its protein partners.

Polymer-Mediated Delivery of siRNAs to Hepatocellular Carcinoma: Variables Affecting Specificity and Effectiveness.

Despite the advances in anticancer therapies, their effectiveness for many human tumors is still far from being optimal. Significant improvements in treatment efficacy can come from the enhancement of drug specificity. This goal may be achieved by combining the use of therapeutic molecules with tumor specific effects and delivery carriers with tumor targeting ability. In this regard, nucleic acid-based drug (NABD) and particularly small interfering RNAs (siRNAs), are attractive molecules due to the possibility to be engineered to target specific tumor genes. On the other hand, polymeric-based delivery systems are emerging as versatile carriers to generate tumor-targeted delivery systems. Here we will focus on the most recent findings in the selection of siRNA/polymeric targeted delivery systems for hepatocellular carcinoma (HCC), a human tumor for which currently available therapeutic approaches are poorly effective. In addition, we will discuss the most attracting and, in our opinion, promising siRNA-polymer combinations for HCC in relation to the biological features of HCC tissue. Attention will be also put on the mathematical description of the mechanisms ruling siRNA-carrier delivery, this being an important aspect to improve effectiveness reducing the experimental work.

Urease Inhibition in the Presence of N-(n-Butyl)thiophosphoric Triamide, a Suicide Substrate: Structure and Kinetics.

The nickel-dependent enzyme urease is a virulence factor for a large number of pathogenic and antibiotic-resistant bacteria, as well as a negative factor for the efficiency of soil nitrogen fertilization for crop production. The use of urease inhibitors to offset these effects requires knowledge, at a molecular level, of their mode of action. The 1.28 Å resolution structure of the enzyme-inhibitor complex obtained upon incubation of Sporosarcina pasteurii urease with N-(n-butyl)thiophosphoric triamide (NBPT), a molecule largely utilized in agriculture, reveals the presence of the monoamidothiophosphoric acid (MATP) moiety, obtained upon enzymatic hydrolysis of the diamide derivative of NBPT (NBPD) to yield n-butyl amine. MATP is bound to the two Ni(II) ions in the active site of urease using a µ2-bridging O atom and terminally bound O and NH2 groups, with the S atom of the thiophosphoric amide pointing away from the metal center. The mobile flap modulating the size of the active site cavity is found in the closed conformation. Docking calculations suggest that the interaction between urease in the open flap conformation and NBPD involves a role for the conserved αArg339 in capturing and orienting the inhibitor prior to flap closure. Calorimetric and spectrophotometric determinations of the kinetic parameters of this inhibition indicate the occurrence of a reversible slow inhibition mode of action, characterized, for both bacterial and plant ureases, by a very small value of the dissociation constant of the urease-MATP complex. No need to convert NBPT to its oxo derivative NBPTO, as previously proposed, is necessary for urease inhibition.