Repository logo
 
Loading...
Profile Picture

Search Results

Now showing 1 - 10 of 26
  • Chemiosmotic misunderstandings
    Publication . Silva, Pedro J.
    Recent publications have questioned the appropriateness of the chemiosmotic theory, a key tenet of modern bioenergetics originally described by Mitchell and since widely improved upon and applied. In one of them, application of Gauss’ law to a model charge distribution in mitochondria was argued to refute the possibility of ATP generation through H+ movement in the absence of a counterion, whereas a different author advocated, for other reasons, the impossibility of chemiosmosis and proposed that a novel energy-generation scheme (referred to as “murburn”) relying on superoxide-catalyzed (or superoxide-promoted) ADP phosphorylation would operate instead. In this letter, those proposals are critically examined and found to be inconsistent with established experimental data and new theoretical calculations.
  • BBr3-assisted cleavage of most ethers does not follow the commonly assumed mechanism
    Publication . Silva, Carla Sousa e; Silva, Pedro J.
    Density‐functional computations were used to probe the reaction mechanism of BBr3‐assisted ether cleavage. After the initial formation of an ether–BBr3 adduct, secondary and tertiary alkyl ethers are cleaved through Br– transfer from the activated BBr3 to the alkyl moiety, as postulated in the literature. In contrast, all other ethers studied react through a novel pathway involving two ether–BBr3 adducts, one of which acts as Br– donor, and the second as the reaction substrate. The identification of the novel bimolecular mechanism for this classical reaction has further applications, because it implies that BBr3‐assisted ether cleavage may become impossible if the ether is surrounded by bulky portions of the molecule that prevent the approach of the attacking BBr3 adduct. Our data also allow the construction of an order of reactivity of alkyl ether deprotection: isopropyl, benzyl, tertiary alkyl, allyl, isobutyl and ethyl can be removed sequentially as their bromo derivatives; phenyl, cyanomethyl and chloromethyl groups can be sequentially removed as their corresponding alcohols.
  • Computational exploration of the reaction mechanism of the Cu+- catalysed synthesis of indoles from N-aryl enaminones
    Publication . Bernardo, Carlos E. P.; Silva, Pedro J.
    We have studied the role of Cu+-phenantroline as a catalyst in the cyclization of N-aryl-enaminones using density-functional theory computations. The catalyst was found to bind the substrate upon deprotonation of its eneaminone, and to dramatically increase the acidity of the carbon adjacent to the ketone functionality. The deprotonation of this carbon atom yields a carbanion which attacks the aryl moiety, thereby closing the heterocycle in the rate-determining step. This C–C bond forming reaction was found to proceed much more rapidly when preceded by re-protonation of the substrate N-atom (which had lost H+ in the initial step). Hydride transfer to the catalyst then completes the indole synthesis, in a very fast step. The influence of Li+ and K+ on the regioselectivity of the cyclization of bromo-substituted analogues could not, however, be reproduced by our model. Alternative pathways involving either single-electron transfer from the catalyst to the substrate or ring cyclization without previous carbon α-deprotonation were found to be kinetically or thermodynamically inaccessible.
  • Evaluation of density functional methods on the geometric and energetic descriptions of species involved in Cu+-promoted catalysis
    Publication . Bernardo, Carlos E. P.; Bauman, Nicholas P.; Piecuch, Piotr; Silva, Pedro J.
    We have evaluated the performance of 15 density functionals of diverse complexity on the geometry optimization and energetic evaluation of model reaction steps present in the proposed reaction mechanisms of Cu(I)-catalyzed indole synthesis and click chemistry of iodoalkynes and azides. The relative effect of the Cu(+) ligand on the relative strength of Cu(+)-alkyne interactions, and the strong preference for a π-bonding mode is captured by all functionals. The best energetic correlations with MP2 are obtained with PBE0, M06-L, and PBE1PW91, which also provide good quality geometries. Furthermore, PBE0 and PBE1PW91 afford the best agreement with the high-level CCSD(T) computations of the deprotonation energies of Cu(+)-coordinated eneamines, where MP2 strongly disagrees with CCSD(T) and the examined DFT functionals. PBE0 also emerged as the most suitable functional for the study of the energetics and geometries of Cu(+) hydrides, while at the same time correctly capturing the influence of the Cu(+) ligands on the metal reactivity.
  • Mechanistic pathways of mercury removal from the organomercurial lyase active site
    Publication . Silva, Pedro J.; Rodrigues, Viviana
    Bacterial populations present in Hg-rich environments have evolved biological mechanisms to detoxify methylmercury and other organometallic mercury compounds. The most common resistance mechanism relies on the H(+)-assisted cleavage of the Hg-C bond of methylmercury by the organomercurial lyase MerB. Although the initial reaction steps which lead to the loss of methane from methylmercury have already been studied experimentally and computationally, the reaction steps leading to the removal of Hg(2+) from MerB and regeneration of the active site for a new round of catalysis have not yet been elucidated. In this paper, we have studied the final steps of the reaction catalyzed by MerB through quantum chemical computations at the combined MP2/CBS//B3PW91/6-31G(d) level of theory. While conceptually simple, these reaction steps occur in a complex potential energy surface where several distinct pathways are accessible and may operate concurrently. The only pathway which clearly emerges as forbidden in our analysis is the one arising from the sequential addition of two thiolates to the metal atom, due to the accumulation of negative charges in the active site. The addition of two thiols, in contrast, leads to two feasible mechanistic possibilities. The most straightforward pathway proceeds through proton transfer from the attacking thiol to Cys159 , leading to its removal from the mercury coordination sphere, followed by a slower attack of a second thiol, which removes Cys96. The other pathway involves Asp99 in an accessory role similar to the one observed earlier for the initial stages of the reaction and affords a lower activation enthalpy, around 14 kcal mol(-1), determined solely by the cysteine removal step rather than by the thiol ligation step. Addition of one thiolate to the intermediates arising from either thiol attack occurs without a barrier and produces an intermediate bound to one active site cysteine and from which Hg(SCH3)2 may be removed only after protonation by solvent-provided H3O(+). Thiolate addition to the active site (prior to any attack by thiols) leads to pathways where the removal of the first cysteine becomes the rate-determining step, irrespective of whether Cys159 or Cys96 leaves first. Comparisons with the recently computed mechanism of the related enzyme MerA further underline the important role of Asp99 in the energetics of the MerB reaction. Kinetic simulation of the mechanism derived from our computations strongly suggests that in vivo the thiolate-only pathway is operative, and the Asp-assisted pathway (as well as the conversion of intermediates of the thiolate pathway into intermediates of the Cys-assisted pathway) is prevented by steric factors absent from our model and related to the precise geometry of the organomercurial binding-pocket.
  • Refining the reaction mechanism of O2 towards its co-substrate in cofactor-free dioxygenases
    Publication . Silva, Pedro J.
    Cofactor-less oxygenases perform challenging catalytic reactions between singlet co-substrates and triplet oxygen, in spite of apparently violating the spin-conservation rule. In 1-H-3-hydroxy-4-oxoquinaldine-2,4-dioxygenase, the active site has been suggested by quantum chemical computations to fine tune triplet oxygen reactivity, allowing it to interact rapidly with its singlet substrate without the need for spin inversion, and in urate oxidase the reaction is thought to proceed through electron transfer from the deprotonated substrate to an aminoacid sidechain, which then feeds the electron to the oxygen molecule. In this work, we perform additional quantum chemical computations on these two systems to elucidate several intriguing features unaddressed by previous workers. These computations establish that in both enzymes the reaction proceeds through direct electron transfer from co-substrate to O2 followed by radical recombination, instead of minimum-energy crossing points between singlet and triplet potential energy surfaces without formal electron transfer. The active site does not affect the reactivity of oxygen directly but is crucial for the generation of the deprotonated form of the co-substrates, which have redox potentials far below those of their protonated forms and therefore may transfer electrons to oxygen without sizeable thermodynamic barriers. This mechanism seems to be shared by most cofactor-less oxidases studied so far.
  • Response to “molecular-level understanding of biological energy coupling and transduction: response to “chemiosmotic misunderstandings”
    Publication . Silva, Pedro J.
    The most recent contribution by Sunil Nath in these pages is, mostly, a repetition of his previous claims regarding failures of the chemiosmotic hypotheses, supplemented with some fresh misunderstandings of the points I had sought to clarify in my previous critique⁠. Considerable portions rehash 50-60 years-old controversies, with no apparent understanding that the current chemiosmotic hypothesis, while birthed by Mitchell, differs from Mitchell's details in many respects. As such, Nath has devoted much time dealing with a few errors (or wrong hypotheses) by Mitchell (in a few places I would almost venture to say "typographical mistakes in typesetting") and presents the ensuing conclusions as "refutations" of the chemiosmotic paradigm, completely neglecting that such details (such as the precise H+/ATP or H+:O ratios) are completely irrelevant to the reality (or not) of an electron-transport chain that uses the free energy liberated by electron-transfer to remove H+ from a compartment, to which it returns through and ATP synthase which uses the energy in that spontaneous return to drive ATP synthesis. The thermodynamical mistakes and misunderstandings of the relevant literature present in Nath's new contribution are so numerous, though, that I feel forced to call the attention of the readers of "Biophysical Chemistry" to them.
  • Computational development of rubromycin-based lead compounds for HIV-1 reverse transcriptase inhibition
    Publication . Bernardo, Carlos E. P.; Silva, Pedro J.
    The binding of several rubromycin-based ligands to HIV1-reverse transcriptase was analyzed using molecular docking and molecular dynamics simulations. MM-PBSA analysis and examination of the trajectories allowed the identification of several promising compounds with predicted high affinity towards reverse transcriptase mutants which have proven resistant to current drugs. Important insights on the complex interplay of factors determining the ability of ligands to selectively target each mutant have been obtained.
  • Influence of alkyne and azide substituents on the choice of the reaction mechanism of the Cu+-catalyzed addition of azides to iodoalkynes
    Publication . Silva, Pedro J.; Bernardo, Carlos E. P.
    The cycloaddition of azides to iodoalkynes is strongly enhanced by some Cu+-complexes. We have studied computationally six reaction pathways for the cycloaddition of 24 combinations of azide and iodoalkyne to identify the dominant pathways and the influence of reactant structure on the evolution of the reaction. Two pathways were found to be operating for distinct sets of reactants. In the first pathway, initial complexation of iodoalkyne by Cu+ is followed by the binding of the azide to the metal through its substituted nitrogen atom, followed by attack of the nonhalogenated alkyne carbon by the terminal nitrogen atom. This pathway is generally followed by aromatic or electron-deficient azides, unless the iodoalkyne bears an electron-withdrawing group. The second pathway is a single-step mechanism similar (apart from the alkyne bond weakening caused by complexation) to that observed in the absence of catalyst. Electron-deficient iodoalkynes and methyl azides strongly prefer this mechanism, regardless of the identity of the reaction partners. The catalytic gain obtained through the use of Cu+ depends only partially on its direct effect on the energy of the transition state (relative to that of the infinitely separated reactants) and may be lost if the iodoalkyne itself strongly interacts with the catalyst through the formation of too strong a π-complex.
  • Will 1,2-dihydro-1,2-azaborine-based drugs resist metabolism by cytochrome P450 compound I?
    Publication . Silva, Pedro J.
    1,2-dihydro-1,2-azaborine is a structural and electronic analogue of benzene which is able to occupy benzene-binding pockets in T4 lysozyme and has been proposed as suitable arene-mimicking group for biological and pharmaceutical applications. Its applicability in a biological context requires it to be able to resist modification by xenobiotic-degrading enzymes like the P450 cytochromes. Quantum chemical computations described in this work show that 1,2-dihydro-1,2-azaborine is much more prone to modification by these enzymes than benzene, unless steric crowding of the ring prevents it from reaching the active site, or otherwise only allows reaction at the less reactive C4-position. This novel heterocyclic compound is therefore expected to be of limited usefulness as an aryl bioisostere.