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.

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.

Influence of Steric Effect on the Pseudo-Multicomponent Synthesis of N-Aroylmethyl-4-Arylimidazoles.

A pseudo-three-component synthesis of N-aroylmethylimidazoles 3 with three new C-N bonds formed regioselectively under microwave conditions was developed. Products were obtained by reacting two equivalents of aroylmethyl bromide (ArCOCH2Br, 1) with the appropriate amidine salt (RCN2H3.HX, 2) and with K2CO3 as a base in acetonitrile. The bicomponent reaction also occurred, giving the expected 4(5)-aryl-1H-imidazoles 4. Notably, the ratio of products 3 and 4 is governed by steric factors of the amidine 2 (i.e., R = H, CH3, Ph). Therefore, a computational study was carried out to understand the reaction course regarding product ratio (3/4), regioselectivity, and the steric effects of the amidine substituent group.

Modeling Kinetics and Thermodynamics of Guest Encapsulation into the [M4L6]12- Supramolecular Organometallic Cage.

The encapsulation of molecular guests into supramolecular hosts is a complex molecular recognition process in which the guest displaces the solvent from the host cavity, while the host deforms to let the guest in. An atomistic description of the association would provide valuable insights on the physicochemical properties that guide it. This understanding may be used to design novel host assemblies with improved properties (i.e., affinities) toward a given class of guests. Molecular simulations may be conveniently used to model the association processes. It is thus of interest to establish efficient protocols to trace the encapsulation process and to predict the associated magnitudes ΔGbind and ΔGbind⧧. Here, we report the calculation of the Gibbs energy barrier and Gibbs binding energy by means of explicit solvent molecular simulations for the [Ga4L6]12- metallocage encapsulating a series of cationic molecules. The ΔGbind⧧ for encapsulation was estimated by means of umbrella sampling simulations. The steps involved were identified, including ion-pair formation and naphthalene rotation (from L ligands of the metallocage) during the guest's entrance. The ΔGbind values were computed using the attach-pull-release method. The results reveal the sensitivity of the estimates on the force field parameters, in particular on atomic charges, showing that higher accuracy is obtained when charges are derived from implicit solvent quantum chemical calculations. Correlation analysis identified some indicators for the binding affinity trends. All computed magnitudes are in very good agreement with experimental observations. This work provides, on one side, a benchmarked way to computationally model a highly charged metallocage encapsulation process. This includes a nonstandard parameterization and charge derivation procedure. On the other hand, it gives specific mechanistic information on the binding processes of [Ga4L6]12- at the molecular level where key motions are depicted. Taken together, the study provides an interesting option for the future design of metal-organic cages.

σ-Silane Platinum(II) Complexes as Intermediates in C-Si Bond-Coupling Processes.

Platinum complexes [Pt(NHC')(NHC)][BArF ] (in which NHC' denotes a cyclometalated N-heterocyclic carbene ligand, NHC) react with primary silanes RSiH3 to afford the cyclometalated platinum(II) silyl complexes [Pt(NHC-SiHR')(NHC)][BArF ] through a process that involves the formation of C-Si and Pt-Si bonds with concomitant extrusion of H2 . Low-temperature NMR studies indicate that the process proceeds through initial formation of the σ-SiH complexes [Pt(NHC')(NHC)(HSiH2 R)][BArF ], which are stable at temperatures below -10 °C. At higher temperatures, activation of one Si-H bond followed by a C-Si coupling reaction generates an agostic SiH platinum hydride derivative [Pt(H)(NHC'-SiH2 R)(NHC)][BArF ], which undergoes a second Si-H bond activation to afford the final products. Computational modeling of the reaction mechanism indicates that the stereochemistry of the silyl/hydride ligands after the first Si-H bond cleavage dictates the nature of the products, favoring the formation of a C-Si bond over a C-H bond, in contrast to previous results obtained for tertiary silanes. Furthermore, the process involves a trans-to-cis isomerization of the NHC ligand before the second Si-H bond cleavage.

Selective Synthesis of Tetrasubstituted Olefins by Copper-Mediated Acetoxythiolation of Internal Alkynes: Scope and Mechanistic Studies.

The Cu-mediated synthesis of tetrasubstituted olefins by the addition of an acetate group and a thiolate to an unactivated internal alkyne is described. The reaction is fully stereoselective, because only the E alkene is obtained. If the alkyne is asymmetric, the reaction also shows a very high degree of regioselectivity. The mechanism of the reaction is elucidated by DFT methods, which show that it takes place through Cu-stabilized radical species. Calculations highlight the crucial role of the dimeric copper(II) diacetate in the process, as it generates the active species in which the sulfur center has an incipient thiyl radical character and accepting, through a series of changes in the oxidation states of the two copper centers, the two electrons released in the addition of two nucleophiles to the alkyne.

Calculation of Reaction Free Energies in Solution: A Comparison of Current Approaches.

The result of the application of different approaches based on the ideal gas/rigid rotor/harmonic oscillator (IGRRHO) model, commonly used in popular software packages, for the calculation of free energies in solution is compared with that of ab initio molecular dynamics for a process involving ligand exchange in palladium complexes. The IGRRHO-based approaches considered differ in most cases in the extent to which the rotational and translational contributions are included in the correction. Our study supports the use the free energy values directly obtained from dispersion-corrected DFT functionals without any correction or with minor corrections at most.

Isolation of a Cationic Platinum(II) σ-Silane Complex.

The platinum complex [Pt(It Bui Pr')(It Bui Pr)][BArF ] interacts with tertiary silanes to form stable (<0 °C) mononuclear PtII σ-SiH complexes [Pt(It Bui Pr')(It Bui Pr)(η1 -HSiR3 )][BArF ]. These compounds have been fully characterized, including X-ray diffraction methods, as the first examples for platinum. DFT calculations (including electronic topological analysis) support the interpretation of the coordination as an unusual η1 -SiH. However, the energies required for achieving a η2 -SiH mode are rather low, and is consistent with the propensity of these derivatives to undergo Si-H cleavage leading to the more stable silyl species [Pt(SiR3 )(It Bui Pr)2 ][BArF ] at room temperature.

Directed Evolution of an Artificial Imine Reductase.

Artificial metalloenzymes, resulting from incorporation of a metal cofactor within a host protein, have received increasing attention in the last decade. The directed evolution is presented of an artificial transfer hydrogenase (ATHase) based on the biotin-streptavidin technology using a straightforward procedure allowing screening in cell-free extracts. Two streptavidin isoforms were yielded with improved catalytic activity and selectivity for the reduction of cyclic imines. The evolved ATHases were stable under biphasic catalytic conditions. The X-ray structure analysis reveals that introducing bulky residues within the active site results in flexibility changes of the cofactor, thus increasing exposure of the metal to the protein surface and leading to a reversal of enantioselectivity. This hypothesis was confirmed by a multiscale approach based mostly on molecular dynamics and protein-ligand dockings.

First-Principles Molecular Dynamics Studies of Organometallic Complexes and Homogeneous Catalytic Processes.

Computational chemistry is a valuable aid to complement experimental studies of organometallic systems and their reactivity. It allows probing mechanistic hypotheses and investigating molecular structures, shedding light on the behavior and properties of molecular assemblies at the atomic scale. When approaching a chemical problem, the computational chemist has to decide on the theoretical approach needed to describe electron/nuclear interactions and the composition of the model used to approximate the actual system. Both factors determine the reliability of the modeling study. The community dedicated much effort to developing and improving the performance and accuracy of theoretical approaches for electronic structure calculations, on which the description of (inter)atomic interactions rely. Here, the importance of the model system used in computational studies is highlighted through examples from our recent research focused on organometallic systems and homogeneous catalytic processes. We show how the inclusion of explicit solvent allows the characterization of molecular events that would otherwise not be accessible in reduced model systems (clusters). These include the stabilization of nascent charged fragments via microscopic solvation (notably, hydrogen bonding), transfer of charge (protons) between distant fragments mediated by solvent molecules, and solvent coordination to unsaturated metal centers. Furthermore, when weak interactions are involved, we show how conformational and solvation properties of organometallic complexes are also affected by the explicit inclusion of solvent molecules. Such extended model systems may be treated under periodic boundary conditions, thus removing the cluster/continuum (or vacuum) boundary, and require a statistical mechanics simulation technique to sample the accessible configurational space. First-principles molecular dynamics, in which atomic forces are computed from electronic structure calculations (namely, density functional theory), is certainly the technique of choice to investigate chemical events in solution. This methodology is well established and thanks to advances in both algorithms and computational resources simulation times required for the modeling of chemical events are nowadays accessible, though the computational requirements use to be high. Specific applications reviewed here include mechanistic studies of the Shilov and Wacker processes, speciation in Pd chemistry, hydrogen bonding to metal centers, and the dynamics of agostic interactions.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases.

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and ß-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

Acid Activation in Phenyliodine Dicarboxylates: Direct Observation, Structures, and Implications.

The use of the hypervalent iodine reagents in oxidative processes has become a staple in modern organic synthesis. Frequently, the reactivity of λ3 iodanes is further enhanced by acids (Lewis or Brønsted). The origin of such activation, however, has remained elusive. Here, we use the common combination of PhI(OAc)2 with BF3·Et2O as a model to fully explore this activation phenomenon. In addition to the spectroscopic assessment of the dynamic acid-base interaction, for the first time the putative PIDA·BF3 complex has been isolated and its structure determined by X-ray diffraction. Consequences of such activation are discussed from a structural and electronic (DFT) points of views, including the origins of the enhanced reactivity.

Rim, Side Arms, and Cavity: Three Sites for the Recognition of Anions by Tetraazolium Resorcinarene Cavitands.

Two tetrabenzoimidazolium-resorcinarene cavitands were prepared and used for the recognition of chloride, bromide, iodide, cyanide, nitrate, perchlorate, hexanoate, benzenesulfonate, and p-toluenesulfonate. Binding affinities of the two cavitands were determined by 1 H NMR titration and computational analysis. The observed spectral changes were related to specific interaction sites, which were supported by the computational studies. In the case of the C2-H tetrabenzoimidazolium-resorcinarene, the recognition region of the inorganic anions and hexanoate was located at the rim of the cavitand, although chloride and bromide also interacted with the aromatic C-H bonds located between adjacent arms of the cavitand. By contrast, the recognition of the two anions with an aromatic ring (benzenesulfonate and p-toluenesulfonate) results from encapsulation of the aromatic part of the anions inside the hydrophobic cavity of the host. In the case of the C2-Me tetrabenzoimizazolium-resorcinarene receptor, the ability of the molecule to bind all inorganic anions and hexanoate was suppressed, but the receptor maintained its ability to strongly bind benzenesulfonate and p-toluenesulfonate. This is interpreted in terms of suppression of the ability of the cavitand to form hydrogen bonds at the rim of the molecule due to replacement of the C2-H proton by a methyl group, while the hydrophobic pocket of the molecule maintains its binding abilities.

NH-Heterocyclic Aryliodonium Salts and their Selective Conversion into N1-Aryl-5-iodoimidazoles.

The synthesis of N-arylimidazoles substituted at the sterically encumbered 5-position is a challenge for modern synthetic approaches. A new family of imidazolyl aryliodonium salts is reported, which serve as a stepping stone on the way to selective formation of N1-aryl-5-iodoimidazoles. Iodine acts as a "universal" placeholder poised for replacement by aryl substituents. These new λ(3) -iodanes are produced by treating the NH-imidazole with ArI(OAc)2 , and are converted to N1-aryl-5-iodoimidazoles by a selective copper-catalyzed aryl migration. The method tolerates a variety of aryl fragments and is also applicable to substituted imidazoles.

A Tetraferrocenyl-Resorcinarene Cavitand as a Redox-Switchable Host of Ammonium Salts.

Tetraannulation of a resorcinarene-octaamino cavitand with ferrocenecarboxaldehyde allows the preparation of a tetrabenzimidazole-resorcinarene cavitand with four ferrocenyl moieties directly linked to the C2 atom of the imidazole units. Oxidation of the four ferrocenyl moieties produces important structural modifications of the molecule, as indicated by DFT calculations performed for the neutral and tetraoxidized forms of the cavitand. By means of (1) H NMR spectroscopic analysis, the encapsulating properties of the new tetraferrocenyl-resorcinarene cavitand toward a series of ammonium salts were evaluated, and a clear cutoff point in binding affinity with respect to size was observed. Cyclic voltammetric studies allowed us to estimate the relative association constants for the neutral and oxidized forms of the cavitand, thus indicating that the guest was bound to the neutral (reduced) state of the cavitand and was released from the oxidized form. These redox-addressable conformational and binding properties of the resorcinarene-tetraferrocenyl cavitand constitute all the necessary features of a redox-switchable molecular gripper. By means of mass-spectrometric analysis, we could unambiguously confirm the molar stoichiometry of the host-guest complex (1:1) and assess the strong guest encapsulation, as indicated by triggering the covalent coupling between host and guest in the gas phase.

Coordination properties of a metal chelator clioquinol to Zn(2+) studied by static DFT and ab initio molecular dynamics.

Several lines of evidence supporting the role of metal ions in amyloid aggregation, one of the hallmarks of Alzheimer's disease (AD), have turned metal ion chelation into a promising therapeutic treatment. The design of efficient chelating ligands requires proper knowledge of the electronic and molecular structure of the complexes formed, including their hydration properties. Among various potential chelators, clioquinol (5-chloro-7-iodo-8-hydroxyquinoline, CQH) has been evaluated with relative success in in vitro experiments and even in phase 2 clinical trials. Clioquinol interacts with Zn(ii) to lead to a binary metal/ligand 1 : 2 stoichiometric complex in which the phenolic group of CQH is deprotonated, resulting in Zn(CQ)2 neutral complexes, to which additional water molecules may coordinate. In the present work, the coordinative properties of clioquinol in aqueous solution have been analyzed by means of static, minimal cluster based DFT calculations and explicit solvent ab initio molecular dynamics simulations. Results from static calculations accounting for solvent effects by means of polarized continuum models suggest that the preferred metal coordination environment is tetrahedral Zn(CQ)2, whereas ab initio molecular dynamics simulations point to quasi degenerate penta Zn(CQ)2(H2O) and hexa Zn(CQ)2(H2O)2 coordinated complexes. The possible reasons for these discrepant results are discussed.

An ionizable active-site tryptophan imparts catalase activity to a peroxidase core.

Catalase peroxidases (KatG's) are bifunctional heme proteins that can disproportionate hydrogen peroxide (catalatic reaction) despite their structural dissimilarity with monofunctional catalases. Using X-ray crystallography and QM/MM calculations, we demonstrate that the catalatic reaction of KatG's involves deprotonation of the active-site Trp, which plays a role similar to that of the distal His in monofunctional catalases. The interaction of a nearby mobile arginine with the distal Met-Tyr-Trp essential adduct (in/out) acts as an electronic switch, triggering deprotonation of the adduct Trp.