Structure of the human cation-chloride cotransport KCC1 in an outward-open state.

Cation-chloride cotransporters (CCCs) catalyze electroneutral symport of Cl- with Na+ and/or K+ across membranes. CCCs are fundamental in cell volume homeostasis, transepithelia ion movement, maintenance of intracellular Cl- concentration, and neuronal excitability. Here, we present a cryoelectron microscopy structure of human K+-Cl- cotransporter (KCC)1 bound with the VU0463271 inhibitor in an outward-open state. In contrast to many other amino acid-polyamine-organocation transporter cousins, our first outward-open CCC structure reveals that opening the KCC1 extracellular ion permeation path does not involve hinge-bending motions of the transmembrane (TM) 1 and TM6 half-helices. Instead, rocking of TM3 and TM8, together with displacements of TM4, TM9, and a conserved intracellular loop 1 helix, underlie alternate opening and closing of extracellular and cytoplasmic vestibules. We show that KCC1 intriguingly exists in one of two distinct dimeric states via different intersubunit interfaces. Our studies provide a blueprint for understanding the mechanisms of CCCs and their inhibition by small molecule compounds.

Cation Chloride Cotransporter NKCC1 Operates through a Rocking-Bundle Mechanism.

The sodium, potassium, and chloride cotransporter 1 (NKCC1) plays a key role in tightly regulating ion shuttling across cell membranes. Lately, its aberrant expression and function have been linked to numerous neurological disorders and cancers, making it a novel and highly promising pharmacological target for therapeutic interventions. A better understanding of how NKCC1 dynamically operates would therefore have broad implications for ongoing efforts toward its exploitation as a therapeutic target through its modulation. Based on recent structural data on NKCC1, we reveal conformational motions that are key to its function. Using extensive deep-learning-guided atomistic simulations of NKCC1 models embedded into the membrane, we captured complex dynamical transitions between alternate open conformations of the inner and outer vestibules of the cotransporter and demonstrated that NKCC1 has water-permeable states. We found that these previously undefined conformational transitions occur via a rocking-bundle mechanism characterized by the cooperative angular motion of transmembrane helices (TM) 4 and 9, with the contribution of the extracellular tip of TM 10. We found these motions to be critical in modulating ion transportation and in regulating NKCC1's water transporting capabilities. Specifically, we identified interhelical dynamical contacts between TM 10 and TM 6, which we functionally validated through mutagenesis experiments of 4 new targeted NKCC1 mutants. We conclude showing that those 4 residues are highly conserved in most Na+-dependent cation chloride cotransporters (CCCs), which highlights their critical mechanistic implications, opening the way to new strategies for NKCC1's function modulation and thus to potential drug action on selected CCCs.

Correct Nucleotide Selection Is Confined at the Binding Site of Polymerase Enzymes.

DNA polymerases (Pols) add incoming nucleotides (deoxyribonucleoside triphosphate (dNTPs)) to growing DNA strands, a crucial step for DNA synthesis. The insertion of correct (vs incorrect) nucleotides relates to Pols' fidelity, which defines Pols' ability to faithfully replicate DNA strands in a template-dependent manner. We and others have demonstrated that reactant alignment and correct base pairing at the Pols catalytic site are crucial structural features to fidelity. Here, we first used equilibrium molecular simulations to demonstrate that the local dynamics at the protein-DNA interface in the proximity of the catalytic site is different when correct vs incorrect dNTPs are bound to polymerase ß (Pol ß). Formation and dynamic stability of specific interatomic interactions around the incoming nucleotide influence the overall binding site architecture. This explains why certain Pols' mutants can affect the local catalytic environment and influence the selection of correct vs incorrect nucleotides. In particular, this is here demonstrated by analyzing the interaction network formed by the residue R283, whose mutant R283A has an experimentally measured lower capacity of differentiating correct (G:dCTP) vs incorrect (G:dATP) base pairing in Pol ß. We also used alchemical free-energy calculations to quantify the G:dCTP →G:dATP transformation in Pol ß wild-type and mutant R283A. These results correlate well with the experimental trend, thus corroborating our mechanistic insights. Sequence and structural comparisons with other Pols from the same family suggest that these findings may also be valid in similar enzymes.

On the Role of Molecular Conformation of the 8-Oxoguanine Lesion in Damaged DNA Processing by Polymerases.

A common and insidious DNA damage is 8-oxoguanine (8OG), bypassed with low catalytic efficiency and high error frequency by polymerases (Pols) during DNA replication. This is a fundamental process with far-reaching implications in cell function and diseases. However, the molecular determinants of how 8OG exactly affects the catalytic efficiency of Pols remain largely unclear. By examining ternary deoxycytidine triphosphate/DNA/Pol complexes containing the 8OG damage, we found that 8OG consistently adopts different conformations when bound to Pols, compared to when in isolated DNA. Equilibrium molecular dynamics and metadynamics free energy calculations quantified that 8OG is in the lowest energy conformation in isolated DNA. In contrast, 8OG adopts high-energy conformations often characterized by intramolecular steric repulsion when bound to Pols. We show that the 8OG conformation can be regulated by mutating Pol residues interacting with the 8OG phosphate group. These findings propose the 8OG conformation as a factor in Pol-mediated processing of damaged DNA.

Drug Design in the Exascale Era: A Perspective from Massively Parallel QM/MM Simulations.

The initial phases of drug discovery - in silico drug design - could benefit from first principle Quantum Mechanics/Molecular Mechanics (QM/MM) molecular dynamics (MD) simulations in explicit solvent, yet many applications are currently limited by the short time scales that this approach can cover. Developing scalable first principle QM/MM MD interfaces fully exploiting current exascale machines - so far an unmet and crucial goal - will help overcome this problem, opening the way to the study of the thermodynamics and kinetics of ligand binding to protein with first principle accuracy. Here, taking two relevant case studies involving the interactions of ligands with rather large enzymes, we showcase the use of our recently developed massively scalable Multiscale Modeling in Computational Chemistry (MiMiC) QM/MM framework (currently using DFT to describe the QM region) to investigate reactions and ligand binding in enzymes of pharmacological relevance. We also demonstrate for the first time strong scaling of MiMiC-QM/MM MD simulations with parallel efficiency of ∼70% up to >80,000 cores. Thus, among many others, the MiMiC interface represents a promising candidate toward exascale applications by combining machine learning with statistical mechanics based algorithms tailored for exascale supercomputers.

Design, synthesis, docking, and biochemical characterization of non-nucleoside SARS-CoV-2 RdRp inhibitors.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a worldwide pandemic. The identification of effective antiviral drugs remains an urgent medical need. In this context, here we report 17 new 1,4-benzopyrone derivatives, which have been designed, synthesized, and characterized for their ability to block the RNA-dependent RNA polymerase (RdRp) enzyme, a promising target for antiviral drug discovery. This compound series represents a good starting point for developing non-nucleoside inhibitors of RdRp. Compounds 4, 5, and 8 were the most promising drug-like candidates with good potency in inhibiting RdRp, improved in vitro pharmacokinetics compared to the initial hits, and no cytotoxicity effects on normal cell (HEK-293). Compound 8 (ARN25592) stands out as the most promising inhibitor. Our results indicate that this new chemical class of 1,4-benzopyrone derivatives deserves further exploration towards novel and potent antiviral drugs for the treatment of SARS-CoV-2 and potentially other viruses.

Alchemical Free Energy Calculations to Investigate Protein-Protein Interactions: the Case of the CDC42/PAK1 Complex.

Here, we show that alchemical free energy calculations can quantitatively compute the effect of mutations at the protein-protein interface. As a test case, we have used the protein complex formed by the small Rho-GTPase CDC42 and its downstream effector PAK1, a serine/threonine kinase. Notably, the CDC42/PAK1 complex offers a wealth of structural, mutagenesis, and binding affinity data because of its central role in cellular signaling and cancer progression. In this context, we have considered 16 mutations in the CDC42/PAK1 complex and obtained excellent agreement between computed and experimental data on binding affinity. Importantly, we also show that a careful analysis of the side-chain conformations in the mutated amino acids can considerably improve the computed estimates, solving issues related to sampling limitations. Overall, this study demonstrates that alchemical free energy calculations can conveniently be integrated into the design of experimental mutagenesis studies.

Finding the Ion in the RNA-Stack: Can Computational Models Accurately Predict Key Functional Elements in Large Macromolecular Complexes?

This viewpoint discusses the predictive power and impact of computational analyses and simulations to gain prospective, experimentally supported mechanistic insights into complex biological systems. Remarkably, two newly resolved cryoEM structures have confirmed the previous, and independent, prediction of the precise localization and dynamics of key catalytic ions in megadalton-large spliceosomal complexes. This outstanding outcome endorses a prominent synergy of computational and experimental methods in the prospective exploration of such large multicomponent biosystems.

A Gold Nanoparticle Nanonuclease Relying on a Zn(II) Mononuclear Complex.

Similarly to enzymes, functionalized gold nanoparticles efficiently catalyze chemical reactions, hence the term nanozymes. Herein, we present our results showing how surface-passivated gold nanoparticles behave as synthetic nanonucleases, able to cleave pBR322 plasmid DNA with the highest efficiency reported so far for catalysts based on a single metal ion mechanism. Experimental and computational data indicate that we have been successful in creating a catalytic site precisely mimicking that suggested for natural metallonucleases relying on a single metal ion for their activity. It comprises one Zn(II) ion to which a phosphate diester of DNA is coordinated. Importantly, as in nucleic acids-processing enzymes, a positively charged arginine plays a key role by assisting with transition state stabilization and by reducing the pKa of the nucleophilic alcohol of a serine. Our results also show how designing a catalyst for a model substrate (bis-p-nitrophenylphosphate) may provide wrong indications as for its efficiency when it is tested against the real target (plasmid DNA).

Recruiting Mechanism and Functional Role of a Third Metal Ion in the Enzymatic Activity of 5' Structure-Specific Nucleases.

Enzymes of the 5' structure-specific nuclease family are crucial for DNA repair, replication, and recombination. One such enzyme is the human exonuclease 1 (hExo1) metalloenzyme, which cleaves DNA strands, acting primarily as a processive 5'-3' exonuclease and secondarily as a 5'-flap endonuclease. Recently, in crystallo reaction intermediates have elucidated how hExo1 exerts hydrolysis of DNA phosphodiester bonds. These hExo1 structures show a third metal ion intermittently bound close to the two-metal-ion active site, to which recessed ends or 5'-flap substrates bind. Evidence of this third ion has been observed in several nucleic-acid-processing metalloenzymes. However, there is still debate over what triggers the (un)binding of this transient third ion during catalysis and whether this ion has a catalytic function. Using extended molecular dynamics and enhanced sampling free-energy simulations, we observed that the carboxyl side chain of Glu89 (located along the arch motif in hExo1) flips frequently from the reactant state to the product state. The conformational flipping of Glu89 allows one metal ion to be recruited from the bulk and promptly positioned near the catalytic center. This is in line with the structural evidence. Additionally, our simulations show that the third metal ion assists the departure, through the mobile arch, of the nucleotide monophosphate product from the catalytic site. Structural comparisons of nuclease enzymes suggest that this Glu(Asp)-mediated mechanism for third ion recruitment and nucleic acid hydrolysis may be shared by other 5' structure-specific nucleases.

Structure and Dynamics of the Acyl Chains in the Membrane Trafficking and Enzymatic Processing of Lipids.

The regulatory chemical mechanisms of lipid trafficking and degradation are involved in many pathophysiological processes, being implicated in severe pain, inflammation, and cancer. In addition, the processing of lipids is also relevant for industrial and environmental applications. However, there is poor understanding of the chemical features that control lipid membrane trafficking and allow lipid-degrading enzymes to efficiently select and hydrolyze specific fatty acids from a complex cellular milieu of bioactive lipids. This is particularly true for lipid acyl chains, which have diverse structures that can critically affect the many complex reactions needed to elongate, desaturate, or transport fatty acids. Building upon our own contributions in this field, we will discuss how molecular simulations, integrated with experimental evidence, have revealed that the structure and dynamics of the lipid tail are actively involved in modulating membrane trafficking at cellular organelles, and enzymatic reactions at cell membranes. Further evidence comes from recent crystal structures of lipid receptors and remodeling enzymes. Taken together, these recent works have identified those structural features of the lipid acyl chain that are crucial for the regioselectivity and stereospecificity of essential desaturation reactions. In this context, we will first illustrate how atomistic and coarse-grained simulations have elucidated the structure-function relationships between the chemical composition of the lipid's acyl chains and the molecular properties of lipid bilayers. Particular emphasis will be given to the prominent chemical role of the number of double carbon-carbon bonds along the lipid acyl chain, that is, discriminating between saturated, monounsaturated, and polyunsaturated lipids. Different levels of saturation in fatty acid molecules dramatically influence the biophysical properties of lipid assemblies and their interaction with proteins. We will then discuss the processing of lipids by membrane-bound enzymes. Our focus will be on lipids such as anandamide and 2-arachidonoylglycerol. These are the main molecules that act as neurotransmitters in the endocannabinoid system. Specifically, recent findings indicate a crucial interplay between the level of saturation of the lipid tail, its energetically and sterically favored conformations, and the hydrophobic accessory cavities in lipid-degrading enzymes, which help form catalytically active conformations of the selected substrate. This Account will emphasize how the specific chemical structure of acyl chains affects the molecular mechanisms for modulating membrane trafficking and selective hydrolysis. The results examined here show that, by using molecular simulations to investigate lipid plasticity and substrate flexibility, researchers can enrich their interpretation of experimental results about the structure-function relationships of lipids. This could positively impact chemical and biological studies in the field and ultimately support protein engineering studies and structure-based drug discovery to target lipid-processing enzymes.

A Transient and Flexible Cation-π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases.

Metal-dependent DNA and RNA nucleases are enzymes that cleave nucleic acids with great efficiency and precision. These enzyme-mediated hydrolytic reactions are fundamental for the replication, repair, and storage of genetic information within the cell. Here, extensive classical and quantum-based free-energy molecular simulations show that a cation-π interaction is transiently formed in situ at the metal core of Bacteriophage-λ Exonuclease (Exo-λ), during catalysis. This noncovalent interaction (Lys131-Tyr154) triggers nucleophile activation for nucleotide excision. Then, our simulations also show the oscillatory dynamics and swinging of the newly formed cation-π dyad, whose conformational change may favor proton release from the cationic Lys131 to the bulk solution, thus restoring the precatalytic protonation state in Exo-λ. Altogether, we report on the novel mechanistic character of cation-π interactions for catalysis. Structural and bioinformatic analyses support that flexible orientation and transient formation of mobile cation-π interactions may represent a common catalytic strategy to promote nucleic acid hydrolysis in DNA and RNA nucleases.

Smoothed Potential MD Simulations for Dissociation Kinetics of Etoposide To Unravel Isoform Specificity in Targeting Human Topoisomerase II.

Human type II topoisomerases (TopoII) are essential for controlling DNA topology within the cell. For this reason, there are a number of TopoII-targeted anticancer drugs that act by inducing DNA cleavage mediated by both TopoII isoforms (TopoIIα and TopoIIß) in cells. However, recent studies suggest that specific poisoning of TopoIIα may be a safer strategy for treating cancer. This is because poisoning of TopoIIß appears to be linked to the generation of secondary leukemia in patients. We recently reported that enzyme-mediated DNA cleavage complexes (in which TopoII is covalently linked to the cleaved DNA during catalysis) formed in the presence of the anticancer drug etoposide persisted approximately 3-fold longer with TopoIIα than TopoIIß. Notably, enhanced drug-target residence time may reduce the adverse effects of specific TopoIIα poisons. However, it is still not clear how to design drugs that are specific for the α isoform. In this study, we report the results of classical molecular dynamics (MD) simulations to comparatively analyze the molecular interactions formed within the TopoII/DNA/etoposide complex with both isoforms. We also used smoothed potential MD to estimate etoposide dissociation kinetics from the two isoform complexes. These extensive classical and enhanced sampling simulations revealed stabilizing interactions of etoposide with two serine residues (Ser763 and Ser800) in TopoIIα. These interactions are missing in TopoIIß, where both amino acids are alanine residues. This may explain the greater persistence of etoposide-stabilized cleavage complexes formed with Topo TopoIIα. These findings could be useful for the rational design of specific TopoIIα poisons.

Molecular-Dynamics-Simulation-Directed Rational Design of Nanoreceptors with Targeted Affinity.

Here, we demonstrate the possibility of rationally designing nanoparticle receptors with targeted affinity and selectivity for specific small molecules. We used atomistic molecular-dynamics (MD) simulations to gradually mutate and optimize the chemical structure of the molecules forming the coating monolayer of gold nanoparticles (1.7 nm gold-core size). The MD-directed design resulted in nanoreceptors with a 10-fold improvement in affinity for the target analyte (salicylate) and a 100-fold decrease of the detection limit by NMR-chemosensing from the millimolar to the micromolar range. We could define the exact binding mode, which features prolonged contacts and deep penetration of the guest into the monolayer, as well as a distinct shape of the effective binding pockets characterized by exposed interacting points.

Evolutionarily Conserved Allosteric Communication in Protein Tyrosine Phosphatases.

Protein tyrosine phosphatases (PTPs) are an important class of regulatory enzymes that exhibit aberrant activities in a wide range of diseases. A detailed mapping of allosteric communication in these enzymes could, thus, reveal the structural basis of physiologically relevant-and, perhaps, therapeutically informative-perturbations (i.e., mutations, post-translational modifications, or binding events) that influence their catalytic states. This study combines detailed biophysical studies of protein tyrosine phosphatase 1B (PTP1B) with bioinformatic analyses of the PTP family to examine allosteric communication in this class of enzymes. Results of X-ray crystallography, molecular dynamics simulations, and sequence-based statistical analyses indicate that PTP1B possesses a broadly distributed allosteric network that is evolutionarily conserved across the PTP family, and findings from both kinetic studies and mutational analyses show that this network is functionally intact in sequence-diverse PTPs. The allosteric network resolved in this study reveals new sites for targeting allosteric inhibitors of PTPs and helps explain the functional influence of a diverse set of disease-associated mutations.

Abietane-Type Diterpenoids Inhibit Protein Tyrosine Phosphatases by Stabilizing an Inactive Enzyme Conformation.

Protein tyrosine phosphatases (PTPs) contribute to a striking variety of human diseases, yet they remain vexingly difficult to inhibit with uncharged, cell-permeable molecules; no inhibitors of PTPs have been approved for clinical use. This study uses a broad set of biophysical analyses to evaluate the use of abietane-type diterpenoids, a biologically active class of phytometabolites with largely nonpolar structures, for the development of pharmaceutically relevant PTP inhibitors. Results of nuclear magnetic resonance analyses, mutational studies, and molecular dynamics simulations indicate that abietic acid can inhibit protein tyrosine phosphatase 1B, a negative regulator of insulin signaling and an elusive drug target, by binding to its active site in a non-substrate-like manner that stabilizes the catalytically essential WPD loop in an inactive conformation; detailed kinetic studies, in turn, show that minor changes in the structures of abietane-type diterpenoids (e.g., the addition of hydrogens) can improve potency (i.e., lower IC50) by 7-fold. These findings elucidate a previously uncharacterized mechanism of diterpenoid-mediated inhibition and suggest, more broadly, that abietane-type diterpenoids are a promising source of structurally diverse-and, intriguingly, microbially synthesizable-molecules on which to base the design of new PTP-inhibiting therapeutics.

A Strategically Located Arg/Lys Residue Promotes Correct Base Paring During Nucleic Acid Biosynthesis in Polymerases.

Polymerases (Pols) synthesize the double-stranded nucleic acids in the Watson-Crick (W-C) conformation, which is critical for DNA and RNA functioning. Yet, the molecular basis to catalyze the W-C base pairing during Pol-mediated nucleic acids biosynthesis remains unclear. Here, through bioinformatics analyses on a large data set of Pol/DNA structures, we first describe the conserved presence of one positively charged residue (Lys or Arg), which is similarly located near the enzymatic two-metal active site, always interacting directly with the incoming substrate (d)NTP. Incidentally, we noted that some Pol/DNA structures showing the alternative Hoogsteen base pairing were often solved with this specific residue either mutated, displaced, or missing. We then used quantum and classical simulations coupled to free-energy calculations to illustrate how, in human DNA Pol-η, the conserved Arg61 favors W-C base pairing through defined interactions with the incoming nucleotide. Taken together, these structural observations and computational results suggest a structural framework in which this specific residue is critical for stabilizing the incoming (d)NTP nucleotide and base pairing during Pol-mediated nucleic acid biosynthesis. These results may benefit enzyme engineering for nucleic acid processing and encourage new drug discovery strategies to modulate Pols function.

Polyamine-containing etoposide derivatives as poisons of human type II topoisomerases: Differential effects on topoisomerase IIα and IIß.

Etoposide is an anticancer drug that acts by inducing topoisomerase II-mediated DNA cleavage. Despite its wide use, etoposide is associated with some very serious side-effects including the development of treatment-related acute myelogenous leukemias. Etoposide targets both human topoisomerase IIα and IIß. However, the contributions of the two enzyme isoforms to the therapeutic vs. leukemogenic properties of the drug are unclear. In order to develop an etoposide-based drug with specificity for cancer cells that express an active polyamine transport system, the sugar moiety of the drug has been replaced with a polyamine tail. To analyze the effects of this substitution on the specificity of hybrid molecules toward the two enzyme isoforms, we analyzed the activity of a series of etoposide-polyamine hybrids toward human topoisomerase IIα and IIß. All of the compounds displayed an ability to induce enzyme-mediated DNA cleavage that was comparable to or higher than that of etoposide. Relative to the parent drug, the hybrid compounds displayed substantially higher activity toward topoisomerase IIß than IIα. Modeling studies suggest that the enhanced specificity may result from interactions with Gln778 in topoisomerase IIß. The corresponding residue in the α isoform is a methionine.

Cooperative motion of a key positively charged residue and metal ions for DNA replication catalyzed by human DNA Polymerase-η.

Trans-lesion synthesis polymerases, like DNA Polymerase-η (Pol-η), are essential for cell survival. Pol-η bypasses ultraviolet-induced DNA damages via a two-metal-ion mechanism that assures DNA strand elongation, with formation of the leaving group pyrophosphate (PPi). Recent structural and kinetics studies have shown that Pol-η function depends on the highly flexible and conserved Arg61 and, intriguingly, on a transient third ion resolved at the catalytic site, as lately observed in other nucleic acid-processing metalloenzymes. How these conserved structural features facilitate DNA replication, however, is still poorly understood. Through extended molecular dynamics and free energy simulations, we unravel a highly cooperative and dynamic mechanism for DNA elongation and repair, which is here described by an equilibrium ensemble of structures that connect the reactants to the products in Pol-η catalysis. We reveal that specific conformations of Arg61 help facilitate the recruitment of the incoming base and favor the proper formation of a pre-reactive complex in Pol-η for efficient DNA editing. Also, we show that a third transient metal ion, which acts concertedly with Arg61, serves as an exit shuttle for the leaving PPi. Finally, we discuss how this effective and cooperative mechanism for DNA repair may be shared by other DNA-repairing polymerases.

Novel Bacterial Topoisomerase Inhibitors Exploit Asp83 and the Intrinsic Flexibility of the DNA Gyrase Binding Site.

DNA gyrases are enzymes that control the topology of DNA in bacteria cells. This is a vital function for bacteria. For this reason, DNA gyrases are targeted by widely used antibiotics such as quinolones. Recently, structural and biochemical investigations identified a new class of DNA gyrase inhibitors called NBTIs (i.e., novel bacterial topoisomerase inhibitors). NBTIs are particularly promising because they are active against multi-drug resistant bacteria, an alarming clinical issue. Structural data recently demonstrated that these NBTIs bind tightly to a newly identified pocket at the dimer interface of the DNA-protein complex. In the present study, we used molecular dynamics (MD) simulations and docking calculations to shed new light on the binding of NBTIs to this site. Interestingly, our MD simulations demonstrate the intrinsic flexibility of this binding site, which allows the pocket to adapt its conformation and form optimal interactions with the ligand. In particular, we examined two ligands, AM8085 and AM8191, which induced a repositioning of a key aspartate (Asp83B), whose side chain can rotate within the binding site. The conformational rearrangement of Asp83B allows the formation of a newly identified H-bond interaction with an NH on the bound NBTI, which seems important for the binding of NBTIs having such functionality. We validated these findings through docking calculations using an extended set of cognate oxabicyclooctane-linked NBTIs derivatives (~150, in total), screened against multiple target conformations. The newly identified H-bond interaction significantly improves the docking enrichment. These insights could be helpful for future virtual screening campaigns against DNA gyrase.