Three applications of path integrals: equilibrium and kinetic isotope effects, and the temperature dependence of the rate constant of the [1,5] sigmatropic hydrogen shift in (Z)-1,3-pentadiene

Recent experiments have confirmed the importance of nuclear quantum effects even in large biomolecules at physiological temperature. Here we describe how the path integral formalism can be used to describe rigorously the nuclear quantum effects on equilibrium and kinetic properties of molecules. Specifically, we explain how path integrals can be employed to evaluate the equilibrium (EIE) and kinetic (KIE) isotope effects, and the temperature dependence of the rate constant. The methodology is applied to the [1,5] sigmatropic hydrogen shift in pentadiene. Both the KIE and the temperature dependence of the rate constant confirm the importance of tunneling and other nuclear quantum effects as well as of the anharmonicity of the potential energy surface. Moreover, previous results on the KIE were improved by using a combination of a high level electronic structure calculation within the harmonic approximation with a path integral anharmonicity correction using a lower level method.Comment: 9 pages, 4 figure

Plasmodium falciparum Choline Kinase Inhibition Leads to a Major Decrease in Phosphatidylethanolamine Causing Parasite Death

This work was supported by Agencia Aragonesa para la Investigación y Desarrollo (ARAID), Ministerio de Economía y Competitividad (CTQ2013-44367-C2-2-P to R.H.-G.) and Diputación General de Aragón (DGA; B89 to R.H.-G.) and the EU Seventh Framework Programme (2007–2013) under BioStruct-X (grant agreement 283570 and BIOSTRUCTX 5186, to R.H.-G.). T.K.S. was supported by the Wellcome Trust grant 093228 and European Community’s Seventh Framework Programme under grant agreement No. 602773 (Project KINDRED).Malaria is a life-threatening disease caused by different species of the protozoan parasite Plasmodium, with P. falciparum being the deadliest. Increasing parasitic resistance to existing antimalarials makes the necessity of novel avenues to treat this disease an urgent priority. The enzymes responsible for the synthesis of phosphatidylcholine and phosphatidylethanolamine are attractive drug targets to treat malaria as their selective inhibition leads to an arrest of the parasite’s growth and cures malaria in a mouse model. We present here a detailed study that reveals a mode of action for two P. falciparum choline kinase inhibitors both in vitro and in vivo. The compounds present distinct binding modes to the choline/ethanolamine-binding site of P. falciparum choline kinase, reflecting different types of inhibition. Strikingly, these compounds primarily inhibit the ethanolamine kinase activity of the P. falciparum choline kinase, leading to a severe decrease in the phosphatidylethanolamine levels within P. falciparum, which explains the resulting growth phenotype and the parasites death. These studies provide an understanding of the mode of action, and act as a springboard for continued antimalarial development efforts selectively targeting P. falciparum choline kinase.Publisher PDFPeer reviewe

Atomic-Level Characterization of the Activation Mechanism of SERCA by Calcium

We have performed molecular dynamics (MD) simulations to elucidate, in atomic detail, the mechanism by which the sarcoplasmic reticulum Ca2+-ATPase (SERCA) is activated by Ca2+. Crystal structures suggest that activation of SERCA occurs when the cytoplasmic head-piece, in an open (E1) conformation stabilized by Ca2+, undergoes a large-scale open-to-closed (E1 to E2) transition that is induced by ATP binding. However, spectroscopic measurements in solution suggest that these structural states (E1 and E2) are not tightly coupled to biochemical states (defined by bound ligands); the closed E2 state predominates even in the absence of ATP, in both the presence and absence of Ca2+. How is this loose coupling consistent with the high efficiency of energy transduction in the Ca2+-ATPase? To provide insight into this question, we performed long (500 ns) all-atom MD simulations starting from the open crystal structure, including a lipid bilayer and water. In both the presence and absence of Ca2+, we observed a large-scale open-to-closed conformational transition within 400 ns, supporting the weak coupling between structural and biochemical states. However, upon closer inspection, it is clear that Ca2+ is necessary and sufficient for SERCA to reach the precise geometrical arrangement necessary for activation of ATP hydrolysis. Contrary to suggestions from crystal structures, but in agreement with solution spectroscopy, the presence of ATP is not required for this activating transition. Principal component analysis showed that Ca2+ reshapes the free energy landscape of SERCA to create a path between the open conformation and the activated closed conformation. Thus the malleability of the free energy landscape is essential for SERCA efficiency, ensuring that ATP hydrolysis is tightly coupled to Ca2+ transport. These results demonstrate the importance of real-time dynamics in the formation of catalytically competent conformations of SERCA, with broad implications for understanding enzymatic catalysis in atomic detail

The evolution of multiple active site configurations in a designed enzyme

Developments in computational chemistry, bioinformatics, and laboratory evolution have facilitated the de novo design and catalytic optimization of enzymes. Besides creating useful catalysts, the generation and iterative improvement of designed enzymes can provide valuable insight into the interplay between the many phenomena that have been suggested to contribute to catalysis. In this work, we follow changes in conformational sampling, electrostatic preorganization, and quantum tunneling along the evolutionary trajectory of a designed Kemp eliminase. We observe that in the Kemp Eliminase KE07, instability of the designed active site leads to the emergence of two additional active site configurations. Evolutionary conformational selection then gradually stabilizes the most efficient configuration, leading to an improved enzyme. This work exemplifies the link between conformational plasticity and evolvability and demonstrates that residues remote from the active sites of enzymes play crucial roles in controlling and shaping the active site for efficient catalysis

Unveiling a novel transient druggable pocket in BACE-1 through molecular simulations: conformational analysis and binding mode of multisite inhibitors

The critical role of BACE-1 in the formation of neurotoxic ß-amyloid peptides in the brain makes it an attractive target for an efficacious treatment of Alzheimer’s disease. However, the development of clinically useful BACE-1 inhibitors has proven to be extremely challeng- ing. In this study we examine the binding mode of a novel potent inhibitor (compound 1, with IC50 80 nM) designed by synergistic combination of two fragments—huprine and rhein— that individually are endowed with very low activity against BACE-1. Examination of crystal structures reveals no appropriate binding site large enough to accommodate 1. Therefore we have examined the conformational flexibility of BACE-1 through extended molecular dynamics simulations, paying attention to the highly flexible region shaped by loops 8–14, 154–169 and 307–318. The analysis of the protein dynamics, together with studies of pocket druggability, has allowed us to detect the transient formation of a secondary binding site, which contains Arg307 as a key residue for the interaction with small molecules, at the edge of the catalytic cleft. The formation of this druggable “floppy” pocket would enable the bind- ing of multisite inhibitors targeting both catalytic and secondary sites. Molecular dynamics simulations of BACE-1 bound to huprine-rhein hybrid compounds support the feasibility of this hypothesis. The results provide a basis to explain the high inhibitory potency of the two enantiomeric forms of 1, together with the large dependence on the length of the oligo- methylenic linker. Furthermore, the multisite hypothesis has allowed us to rationalize the inhibitory potency of a series of tacrine-chromene hybrid compounds, specifically regarding the apparent lack of sensitivity of the inhibition constant to the chemical modifications intro- duced in the chromene unit. Overall, these findings pave the way for the exploration of novel functionalities in the design of optimized BACE-1 multisite inhibitors

Biochemical, Structural and Molecular Dynamics Analyses of the Potential Virulence Factor RipA from Yersinia pestis

Human diseases are attributed in part to the ability of pathogens to evade the eukaryotic immune systems. A subset of these pathogens has developed mechanisms to survive in human macrophages. Yersinia pestis, the causative agent of the bubonic plague, is a predominately extracellular pathogen with the ability to survive and replicate intracellularly. A previous study has shown that a novel rip (required for intracellular proliferation) operon (ripA, ripB and ripC) is essential for replication and survival of Y. pestis in postactivated macrophages, by playing a role in lowering macrophage-produced nitric oxide (NO) levels. A bioinformatics analysis indicates that the rip operon is conserved among a distally related subset of macrophage-residing pathogens, including Burkholderia and Salmonella species, and suggests that this previously uncharacterized pathway is also required for intracellular survival of these pathogens. The focus of this study is ripA, which encodes for a protein highly homologous to 4-hydroxybutyrate-CoA transferase; however, biochemical analysis suggests that RipA functions as a butyryl-CoA transferase. The 1.9 Å X-ray crystal structure reveals that RipA belongs to the class of Family I CoA transferases and exhibits a unique tetrameric state. Molecular dynamics simulations are consistent with RipA tetramer formation and suggest a possible gating mechanism for CoA binding mediated by Val227. Together, our structural characterization and molecular dynamic simulations offer insights into acyl-CoA specificity within the active site binding pocket, and support biochemical results that RipA is a butyryl-CoA transferase. We hypothesize that the end product of the rip operon is butyrate, a known anti-inflammatory, which has been shown to lower NO levels in macrophages. Thus, the results of this molecular study of Y. pestis RipA provide a structural platform for rational inhibitor design, which may lead to a greater understanding of the role of RipA in this unique virulence pathway

Biapenem Inactivation by B2 Metallo β-Lactamases: Energy Landscape of the Post-Hydrolysis Reactions

<div><h3>Background</h3><p>The first line of defense by bacteria against <em>β</em>-lactam antibiotics is the expression of β-lactamases, which cleave the amide bond of the β-lactam ring. In the reaction of biapenem inactivation by B2 metallo β-lactamases (MβLs), after the β-lactam ring is opened, the carboxyl group generated by the hydrolytic process and the hydroxyethyl group (common to all carbapenems) rotate around the C5–C6 bond, assuming a new position that allows a proton transfer from the hydroxyethyl group to C2, and a nucleophilic attack on C3 by the oxygen atom of the same side-chain. This process leads to the formation of a bicyclic compound, as originally observed in the X-ray structure of the metallo β-lactamase CphA in complex with product.</p> <h3>Methodology/Principal Findings</h3><p>QM/MM and metadynamics simulations of the post-hydrolysis steps in solution and in the enzyme reveal that while the rotation of the hydroxyethyl group can occur in solution or in the enzyme active site, formation of the bicyclic compound occurs primarily in solution, after which the final product binds back to the enzyme. The calculations also suggest that the rotation and cyclization steps can occur at a rate comparable to that observed experimentally for the enzymatic inactivation of biapenem only if the hydrolysis reaction leaves the N4 nitrogen of the β-lactam ring unprotonated.</p> <h3>Conclusions/Significance</h3><p>The calculations support the existence of a common mechanism (in which ionized N4 is the leaving group) for carbapenems hydrolysis in all MβLs, and suggest a possible revision of mechanisms for B2 MβLs in which the cleavage of the β-lactam ring is associated with or immediately followed by protonation of N4. The study also indicates that the bicyclic derivative of biapenem has significant affinity for B2 MβLs, and that it may be possible to obtain clinically effective inhibitors of these enzymes by modification of this lead compound.</p> </div

A Unified Model of the GABA(A) Receptor Comprising Agonist and Benzodiazepine Binding Sites

We present a full-length α(1)β(2)γ(2) GABA receptor model optimized for agonists and benzodiazepine (BZD) allosteric modulators. We propose binding hypotheses for the agonists GABA, muscimol and THIP and for the allosteric modulator diazepam (DZP). The receptor model is primarily based on the glutamate-gated chloride channel (GluCl) from C. elegans and includes additional structural information from the prokaryotic ligand-gated ion channel ELIC in a few regions. Available mutational data of the binding sites are well explained by the model and the proposed ligand binding poses. We suggest a GABA binding mode similar to the binding mode of glutamate in the GluCl X-ray structure. Key interactions are predicted with residues α(1)R66, β(2)T202, α(1)T129, β(2)E155, β(2)Y205 and the backbone of β(2)S156. Muscimol is predicted to bind similarly, however, with minor differences rationalized with quantum mechanical energy calculations. Muscimol key interactions are predicted to be α(1)R66, β(2)T202, α(1)T129, β(2)E155, β(2)Y205 and β(2)F200. Furthermore, we argue that a water molecule could mediate further interactions between muscimol and the backbone of β(2)S156 and β(2)Y157. DZP is predicted to bind with interactions comparable to those of the agonists in the orthosteric site. The carbonyl group of DZP is predicted to interact with two threonines α(1)T206 and γ(2)T142, similar to the acidic moiety of GABA. The chlorine atom of DZP is placed near the important α(1)H101 and the N-methyl group near α(1)Y159, α(1)T206, and α(1)Y209. We present a binding mode of DZP in which the pending phenyl moiety of DZP is buried in the binding pocket and thus shielded from solvent exposure. Our full length GABA(A) receptor is made available as Model S1