Laura L. Perissinotti

Effects of the Donor–Acceptor Distance and Dynamics on Hydride Tunneling in the Dihydrofolate Reductase Catalyzed Reaction

A significant contemporary question in enzymology involves the role of protein dynamics and hydrogen tunneling in enhancing enzyme catalyzed reactions. Here, we report a correlation between the donor-acceptor distance (DAD) distribution and intrinsic kinetic isotope effects (KIEs) for the dihydrofolate reductase (DHFR) catalyzed reaction. This study compares the nature of the hydride-transfer step for a series of active-site mutants, where the size of a side chain that modulates the DAD (I14 in E. coli DHFR) is systematically reduced (I14V, I14A, and I14G). The contributions of the DAD and its dynamics to the hydride-transfer step were examined by the temperature dependence of intrinsic KIEs, hydride-transfer rates, activation parameters, and classical molecular dynamics (MD) simulations. Results are interpreted within the framework of the Marcus-like model where the increase in the temperature dependence of KIEs arises as a direct consequence of the deviation of the DAD from its distribution in the wild type enzyme. Classical MD simulations suggest new populations with larger average DADs, as well as broader distributions, and a reduction in the population of the reactive conformers correlated with the decrease in the size of the hydrophobic residue. The more flexible active site in the mutants required more substantial thermally activated motions for effective H-tunneling, consistent with the hypothesis that the role of the hydrophobic side chain of I14 is to restrict the distribution and dynamics of the DAD and thus assist the hydride-transfer. These studies establish relationships between the distribution of DADs, the hydride-transfer rates, and the DADs rearrangement toward tunneling-ready states. This structure-function correlation shall assist in the interpretation of the temperature dependence of KIEs caused by mutants far from the active site in this and other enzymes, and may apply generally to C-H→C transfer reactions.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.

A Microscopic Study of the Deoxyhemoglobin-Catalyzed Generation of Nitric Oxide from Nitrite Anion†

There is recent evidence suggesting that nitrite anion (NO 2 (-)) represents the major intravascular NO storage molecule whose transduction to NO is facilitated by a reduction mechanism catalyzed by deoxygenated hemoglobin (deoxy-Hb). In this work, we provide a detailed microscopic study of deoxy-Hb nitrite reductase (NIR) activity by combining classical molecular dynamics and hybrid quantum mechanical-molecular mechanical simulations. Our results point out that two alternative mechanisms could be operative and suggest that the most energetic barriers should stem from either reprotonation of the distal histidine or NO dissociation from the ferric heme. In the first proposed mechanism, which is similar to that proposed for bacterial NIRs, nitrite anion or nitrous acid coordinates to the heme through the N atom. This pathway involves HisE7 in a one or two proton transfer process, depending on whether the active species is nitrite anion or nitrous acid, to yield an intermediate Fe(III)NO species which eventually dissociates leading to NO and methemoglobin. In the second mechanism, the nitrite anion coordinates to the heme through the O atom. This pathway requires only one proton transfer from HisE7 and leads directly to the formation of a hydroxo Fe(III) complex and NO.

Ivabradine prolongs phase 3 of cardiac repolarization and blocks the hERG1 (KCNH2) current over a concentration-range overlapping with that required to block HCN4

In Europe, ivabradine has recently been approved to treat patients with angina who have intolerance to beta blockers and/or heart failure. Ivabradine is considered to act specifically on the sinoatrial node by inhibiting the If current (the funny current) to slow automaticity. However, in vitro studies show that ivabradine prolongs phase 3 repolarization in ventricular tissue. No episodes of Torsades de Pointes have been reported in randomized clinical studies. The objective of this study is to assess whether ivabradine blocked the hERG1 current. In the present study we discovered that ivabradine prolongs action potential and blocks the hERG current over a range of concentrations overlapping with those required to block HCN4. Ivabradine produced tonic, rather than use-dependent block. The mutation Y652A significantly suppressed pharmacologic block of hERG by ivabradine. Disruption of C-type inactivation also suppressed block of hERG1 by ivabradine. Molecular docking and molecular dynamics simulations indicate that ivabradine may access the inner cavity of the hERG1 via a lipophilic route and has a well-defined binding site in the closed state of the channel. Structural organization of the binding pockets for ivabradine is discussed. Ivabradine blocks hERG and prolongs action potential duration. Our study is potentially important because it indicates the need for active post marketing surveillance of ivabradine. Importantly, proarrhythmia of a number of other drugs has only been discovered during post marketing surveillance.

Small ligand–globin interactions: Reviewing lessons derived from computer simulation

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

The effect of active-site isoleucine to alanine mutation on the DHFR catalyzed hydride-transfer

Comparison of the nature of hydride transfer in wild-type and active site mutant (I14A) of dihydrofolate reductase suggests that the size of this side chain at position 14 modulates H-tunneling.

NS1643 Interacts around L529 of hERG to Alter Voltage Sensor Movement on the Path to Activation

Activators of hERG1 such as NS1643 are being developed for congenital/acquired long QT syndrome. Previous studies identify the neighborhood of L529 around the voltage-sensor as a putative interacting site for NS1643. With NS1643, the V1/2 of activation of L529I (-34 ± 4 mV) is similar to wild-type (WT) (-37 ± 3 mV; P > 0.05). WT and L529I showed no difference in the slope factor in the absence of NS1643 (8 ± 0 vs. 9 ± 0) but showed a difference in the presence of NS1643 (9 ± 0.3 vs. 22 ± 1; P < 0.01). Voltage-clamp-fluorimetry studies also indicated that in L529I, NS1643 reduces the voltage-sensitivity of S4 movement. To further assess mechanism of NS1643 action, mutations were made in this neighborhood. NS1643 shifts the V1/2 of activation of both K525C and K525C/L529I to hyperpolarized potentials (-131 ± 4 mV for K525C and -120 ± 21 mV for K525C/L529I). Both K525C and K525C/K529I had similar slope factors in the absence of NS1643 (18 ± 2 vs. 34 ± 5, respectively) but with NS1643, the slope factor of K525C/L529I increased from 34 ± 5 to 71 ± 10 (P < 0.01) whereas for K525C the slope factor did not change (18 ± 2 at baseline and 16 ± 2 for NS1643). At baseline, K525R had a slope factor similar to WT (9 vs. 8) but in the presence of NS1643, the slope factor of K525R was increased to 24 ± 4 vs. 9 ± 0 mV for WT (P < 0.01). Molecular modeling indicates that L529I induces a kink in the S4 voltage-sensor helix, altering a salt-bridge involving K525. Moreover, docking studies indicate that NS1643 binds to the kinked structure induced by the mutation with a higher affinity. Combining biophysical, computational, and electrophysiological evidence, a mechanistic principle governing the action of some activators of hERG1 channels is proposed.

Kinetic Model for NS1643 Drug Activation of WT and L529I Variants of Kv11.1 (hERG1) Potassium Channel

Congenital and acquired (drug-induced) forms of the human long-QT syndrome are associated with alterations in Kv11.1 (hERG) channel-controlled repolarizing IKr currents of cardiac action potentials. A mandatory drug screen implemented by many countries led to a discovery of a large group of small molecules that can activate hERG currents and thus may act as potent antiarrhythmic agents. Despite significant progress in identification of channel activators, little is known about their mechanism of action. A combination of electrophysiological studies with molecular and kinetic modeling was used to examine the mechanism of a model activator (NS1643) action on the hERG channel and its L529I mutant. The L529I mutant has gating dynamics similar to that of wild-type while its response to application of NS1643 is markedly different. We propose a mechanism compatible with experiments in which the model activator binds to the closed (C3) and open states (O). We suggest that NS1643 is affecting early gating transitions, probably during movements of the voltage sensor that precede the opening of the activation gate.

Folate binding site of flavin-dependent thymidylate synthase.

The DNA nucleotide thymidylate is synthesized by the enzyme thymidylate synthase, which catalyzes the reductive methylation of deoxyuridylate using the cofactor methylene-tetrahydrofolate (CH2H4folate). Most organisms, including humans, rely on the thyA- or TYMS-encoded classic thymidylate synthase, whereas, certain microorganisms, including all Rickettsia and other pathogens, use an alternative thyX-encoded flavin-dependent thymidylate synthase (FDTS). Although several crystal structures of FDTSs have been reported, the absence of a structure with folates limits understanding of the molecular mechanism and the scope of drug design for these enzymes. Here we present X-ray crystal structures of FDTS with several folate derivatives, which together with mutagenesis, kinetic analysis, and computer modeling shed light on the cofactor binding and function. The unique structural data will likely facilitate further elucidation of FDTSs’ mechanism and the design of structure-based inhibitors as potential leads to new antimicrobial drugs.

A Unique Family of Stable and Water-Soluble S-Nitrosothiol Complexes

In this work, we present a complete and detailed experimental characterization and theoretical study of a variety of coordinated S-nitrosothiols (RSNOs), such as cysteine derivatives, mercaptosuccinic acid, benzyl thiol, and phenyl thiol. Some of them are extremely unstable and sensitive in free form. Strikingly, in contrast with free S-nitrosothiols, we found that, upon coordination to iridium, they become very stable even in aqueous solutions. The study of these coordinated complexes provides further insight on the elucidation of structural aspects dealing with the nature of the S-N bond in RSNOs, a fact which still remains a matter of controversy.