Tautomer-Specific Deacylation and Ω-Loop Flexibility Explain the Carbapenem-Hydrolyzing Broad-Spectrum Activity of the KPC-2 ß-Lactamase.

KPC-2 (Klebsiella pneumoniae carbapenemase-2) is a globally disseminated serine-ß-lactamase (SBL) responsible for extensive ß-lactam antibiotic resistance in Gram-negative pathogens. SBLs inactivate ß-lactams via a mechanism involving a hydrolytically labile covalent acyl-enzyme intermediate. Carbapenems, the most potent ß-lactams, evade the activity of many SBLs by forming long-lived inhibitory acyl-enzymes; however, carbapenemases such as KPC-2 efficiently deacylate carbapenem acyl-enzymes. We present high-resolution (1.25-1.4 Å) crystal structures of KPC-2 acyl-enzymes with representative penicillins (ampicillin), cephalosporins (cefalothin), and carbapenems (imipenem, meropenem, and ertapenem) obtained utilizing an isosteric deacylation-deficient mutant (E166Q). The mobility of the Ω-loop (residues 165-170) negatively correlates with antibiotic turnover rates (kcat), highlighting the role of this region in positioning catalytic residues for efficient hydrolysis of different ß-lactams. Carbapenem-derived acyl-enzyme structures reveal the predominance of the Δ1-(2R) imine rather than the Δ2 enamine tautomer. Quantum mechanics/molecular mechanics molecular dynamics simulations of KPC-2:meropenem acyl-enzyme deacylation used an adaptive string method to differentiate the reactivity of the two isomers. These identify the Δ1-(2R) isomer as having a significantly (7 kcal/mol) higher barrier than the Δ2 tautomer for the (rate-determining) formation of the tetrahedral deacylation intermediate. Deacylation is therefore likely to proceed predominantly from the Δ2, rather than the Δ1-(2R) acyl-enzyme, facilitated by tautomer-specific differences in hydrogen-bonding networks involving the carbapenem C-3 carboxylate and the deacylating water and stabilization by protonated N-4, accumulating a negative charge on the Δ2 enamine-derived oxyanion. Taken together, our data show how the flexible Ω-loop helps confer broad-spectrum activity upon KPC-2, while carbapenemase activity stems from efficient deacylation of the Δ2-enamine acyl-enzyme tautomer.

Multiscale Workflow for Modeling Ligand Complexes of Zinc Metalloproteins.

Zinc metalloproteins are ubiquitous, with protein zinc centers of structural and functional importance, involved in interactions with ligands and substrates and often of pharmacological interest. Biomolecular simulations are increasingly prominent in investigations of protein structure, dynamics, ligand interactions, and catalysis, but zinc poses a particular challenge, in part because of its versatile, flexible coordination. A computational workflow generating reliable models of ligand complexes of biological zinc centers would find broad application. Here, we evaluate the ability of alternative treatments, using (nonbonded) molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) at semiempirical (DFTB3) and density functional theory (DFT) levels of theory, to describe the zinc centers of ligand complexes of six metalloenzyme systems differing in coordination geometries, zinc stoichiometries (mono- and dinuclear), and the nature of interacting groups (specifically the presence of zinc-sulfur interactions). MM molecular dynamics (MD) simulations can overfavor octahedral geometries, introducing additional water molecules to the zinc coordination shell, but this can be rectified by subsequent semiempirical (DFTB3) QM/MM MD simulations. B3LYP/MM geometry optimization further improved the accuracy of the description of coordination distances, with the overall effectiveness of the approach depending upon factors, including the presence of zinc-sulfur interactions that are less well described by semiempirical methods. We describe a workflow comprising QM/MM MD using DFTB3 followed by QM/MM geometry optimization using DFT (e.g., B3LYP) that well describes our set of zinc metalloenzyme complexes and is likely to be suitable for creating accurate models of zinc protein complexes when structural information is more limited.

Molecular Simulations suggest Vitamins, Retinoids and Steroids as Ligands of the Free Fatty Acid Pocket of the SARS-CoV-2 Spike Protein*.

We investigate binding of linoleate and other potential ligands to the recently discovered fatty acid binding site in the SARS-CoV-2 spike protein, using docking and molecular dynamics simulations. Simulations suggest that linoleate and dexamethasone stabilize the locked spike conformation, thus reducing the opportunity for ACE2 interaction. In contrast, cholesterol may expose the receptor-binding domain by destabilizing the closed structure, preferentially binding to a different site in the hinge region of the open structure. We docked a library of FDA-approved drugs to the fatty acid site using an approach that reproduces the structure of the linoleate complex. Docking identifies steroids (including dexamethasone and vitamin D); retinoids (some known to be active in vitro, and vitamin A); and vitamin K as potential ligands that may stabilize the closed conformation. The SARS-CoV-2 spike fatty acid site may bind a diverse array of ligands, including dietary components, and therefore provides a promising target for therapeutics or prophylaxis.

A transient receptor potential ion channel in Chlamydomonas shares key features with sensory transduction-associated TRP channels in mammals.

Sensory modalities are essential for navigating through an ever-changing environment. From insects to mammals, transient receptor potential (TRP) channels are known mediators for cellular sensing. Chlamydomonas reinhardtii is a motile single-celled freshwater green alga that is guided by photosensory, mechanosensory, and chemosensory cues. In this type of alga, sensory input is first detected by membrane receptors located in the cell body and then transduced to the beating cilia by membrane depolarization. Although TRP channels seem to be absent in plants, C. reinhardtii possesses genomic sequences encoding TRP proteins. Here, we describe the cloning and characterization of a C. reinhardtii version of a TRP channel sharing key features present in mammalian TRP channels associated with sensory transduction. In silico sequence-structure analysis unveiled the modular design of TRP channels, and electrophysiological experiments conducted on Human Embryonic Kidney-293T cells expressing the Cr-TRP1 clone showed that many of the core functional features of metazoan TRP channels are present in Cr-TRP1, suggesting that basic TRP channel gating characteristics evolved early in the history of eukaryotes.

Voltage sensor gating charge transfer in a hERG potassium channel model.

Relaxation of a hERG K(+) channel model during molecular-dynamics simulation in a hydrated POPC bilayer was accompanied by transitions of an arginine gating charge across a charge transfer center in two voltage sensor domains. Inspection of the passage of arginine side chains across the charge transfer center suggests that the unique hydration properties of the arginine guanidine cation facilitates charge transfer during voltage sensor responses to changes in membrane potential, and underlies the preference of Arg over Lys as a mobile charge carrier in voltage-sensitive ion channels.

Assessing hERG pore models as templates for drug docking using published experimental constraints: the inactivated state in the context of drug block.

Many structurally and therapeutically diverse drugs interact with the human heart K+ channel hERG by binding within the K+ permeation pathway of the open channel, leading to drug-induced 'long QT syndrome'. Drug binding to hERG is often stabilized by inactivation gating. In the absence of a crystal structure, hERG pore homology models have been used to characterize drug interactions. Here we assess potentially inactivated states of the bacterial K+ channel, KcsA, as templates for inactivated state hERG pore models in the context of drug binding using computational docking. Although Flexidock and GOLD docking produced low energy score poses in the models tested, each method selected a MthK K+ channel-based model over models based on the putative inactivated state KcsA structures for each of the 9 drugs tested. The variety of docking poses found indicates that an optimal arrangement for drug binding of aromatic side chains in the hERG pore can be achieved in several different configurations. This plasticity of the drug "binding site" is likely to be a feature of the hERG inactivated state. The results demonstrate that experimental data on specific drug interactions can be used as structural constraints to assess and refine hERG homology models.

Interactions between voltage sensor and pore domains in a hERG K+ channel model from molecular simulations and the effects of a voltage sensor mutation.

The hERG K(+) channel is important for establishing normal electrical activity in the human heart. The channel's unique gating response to membrane potential changes indicates specific interactions between voltage sensor and pore domains that are poorly understood. In the absence of a crystal structure we constructed a homology model of the full hERG membrane domain and performed 0.5 µs molecular dynamics (MD) simulations in a hydrated membrane. The simulations identify potential interactions involving residues at the extracellular surface of S1 in the voltage sensor and at the N-terminal end of the pore helix in the hERG model. In addition, a diffuse interface involving hydrophobic residues on S4 (voltage sensor) and pore domain S5 of an adjacent subunit was stable during 0.5 µs of simulation. To assess the ability of the model to give insight into the effects of channel mutation we simulated a hERG mutant that contains a Leu to Pro substitution in the voltage sensor S4 helical segment (hERG L532P). Consistent with the retention of gated K(+) conductance, the L532P mutation was accommodated in the S4 helix with little disruption of helical structure. The mutation reduced the extent of interaction across the S4-S5 interface, suggesting a structural basis for the greatly enhanced deactivation rate in hERG L532P. The study indicates that pairwise comparison of wild-type and mutated channel models is a useful approach to interpreting functional data where uncertainty in model structures exist.

The hERG K(+) channel S4 domain L532P mutation: characterization at 37°C.

hERG (human Ether-à-go-go Related Gene) is responsible for ion channels mediating rapid delayed rectifier potassium current, I(Kr), which is key to cardiac action potential repolarization. Gain-of-function hERG mutations give rise to the SQT1 variant of the Short QT Syndrome (SQTS). Reggae mutant zebrafish, with a S4 zERG mutation (Leucine499Proline; L499P), display arrhythmic features analogous to those seen in the SQTS. The affected S4 domain ERG residue is highly conserved. This study was executed to determine how the homologous hERG mutation (L532P) influences channel function at 37°C. Whole-cell measurements of current (I(hERG)) were made from HEK 293 cells expressing WT or L532P hERG. The half maximal activation voltage (V(0.5)) of L532P I(hERG) was positively shifted by ~+36mV compared to WT I(hERG); however at negative voltages a pronounced L532P I(hERG) was observed. Both activation and deactivation time-courses were accelerated for L532P I(hERG). The inactivation V(0.5) for L532P I(hERG) was shifted by ~+32mV. Under action potential (AP) voltage-clamp, L532P I(hERG) exhibited a dome-shaped current peaking at ~+16mV, compared to ~-31mV for WT-I(hERG). The L532P mutation produced an ~5-fold increase in the IC(50) for dronedarone inhibition of I(hERG). Homology modeling indicated that the L532 residue within the S4 helix lies closely apposed to the S5 region of an adjacent hERG subunit. Alterations to the S4 domain structure and, potentially, to interactions between adjacent hERG subunits are likely to account for the functional effects of this mutation.

Sequence and structural conservation reveal fingerprint residues in TRP channels.

Transient receptor potential (TRP) proteins are a large family of cation-selective channels, surpassed in variety only by voltage-gated potassium channels. Detailed molecular mechanisms governing how membrane voltage, ligand binding, or temperature can induce conformational changes promoting the open state in TRP channels are still a matter of debate. Aiming to unveil distinctive structural features common to the transmembrane domains within the TRP family, we performed phylogenetic reconstruction, sequence statistics, and structural analysis over a large set of TRP channel genes. Here, we report an exceptionally conserved set of residues. This fingerprint is composed of twelve residues localized at equivalent three-dimensional positions in TRP channels from the different subtypes. Moreover, these amino acids are arranged in three groups, connected by a set of aromatics located at the core of the transmembrane structure. We hypothesize that differences in the connectivity between these different groups of residues harbor the apparent differences in coupling strategies used by TRP subgroups.

Evolutionary analyses reveal independent origins of gene repertoires and structural motifs associated to fast inactivation in calcium-selective TRPV channels.

Essential for calcium homeostasis, TRPV5 and TRPV6 are calcium-selective channels belonging to the transient receptor potential (TRP) gene family. In this study, we investigated the evolutionary history of these channels to add an evolutionary context to the already available physiological information. Phylogenetic analyses revealed that paralogs found in mammals, sauropsids, amphibians, and chondrichthyes, are the product of independent duplication events in the ancestor of each group. Within amniotes, we identified a traceable signature of three amino acids located at the amino-terminal intracellular region. The signature correlates with both the duplication events and the phenotype of fast inactivation observed in mammalian TRPV6 channels. Electrophysiological recordings and mutagenesis revealed that the signature sequence modulates the phenotype of fast inactivation in all clades of vertebrates but reptiles. A transcriptome analysis showed a change in tissue expression from gills, in marine vertebrates, to kidneys in terrestrial vertebrates. Our results highlight a cytoplasmatic structural triad composed by the Helix-Loop-Helix domain, the S2-S3 linker, and the TRP domain helix that is important on modulating the activity of calcium-selective TRPV channels.

ß-Lactamases and ß-Lactamase Inhibitors in the 21st Century.

The ß-lactams retain a central place in the antibacterial armamentarium. In Gram-negative bacteria, ß-lactamase enzymes that hydrolyze the amide bond of the four-membered ß-lactam ring are the primary resistance mechanism, with multiple enzymes disseminating on mobile genetic elements across opportunistic pathogens such as Enterobacteriaceae (e.g., Escherichia coli) and non-fermenting organisms (e.g., Pseudomonas aeruginosa). ß-Lactamases divide into four classes; the active-site serine ß-lactamases (classes A, C and D) and the zinc-dependent or metallo-ß-lactamases (MBLs; class B). Here we review recent advances in mechanistic understanding of each class, focusing upon how growing numbers of crystal structures, in particular for ß-lactam complexes, and methods such as neutron diffraction and molecular simulations, have improved understanding of the biochemistry of ß-lactam breakdown. A second focus is ß-lactamase interactions with carbapenems, as carbapenem-resistant bacteria are of grave clinical concern and carbapenem-hydrolyzing enzymes such as KPC (class A) NDM (class B) and OXA-48 (class D) are proliferating worldwide. An overview is provided of the changing landscape of ß-lactamase inhibitors, exemplified by the introduction to the clinic of combinations of ß-lactams with diazabicyclooctanone and cyclic boronate serine ß-lactamase inhibitors, and of progress and strategies toward clinically useful MBL inhibitors. Despite the long history of ß-lactamase research, we contend that issues including continuing unresolved questions around mechanism; opportunities afforded by new technologies such as serial femtosecond crystallography; the need for new inhibitors, particularly for MBLs; the likely impact of new ß-lactam:inhibitor combinations and the continuing clinical importance of ß-lactams mean that this remains a rewarding research area.

Interactions between amiodarone and the hERG potassium channel pore determined with mutagenesis and in silico docking.

The antiarrhythmic drug amiodarone delays cardiac repolarisation through inhibition of hERG-encoded potassium channels responsible for the rapid delayed rectifier potassium current (IKr). This study aimed to elucidate molecular determinants of amiodarone binding to the hERG channel. Whole-cell patch-clamp recordings were made at 37°C of ionic current (IhERG) carried by wild-type (WT) or mutant hERG channels expressed in HEK293 cells. Alanine mutagenesis and ligand docking were used to investigate the roles of pore cavity amino-acid residues in amiodarone binding. Amiodarone inhibited WT outward IhERG tails with a half-maximal inhibitory concentration (IC50) of ∼45nM, whilst inward IhERG tails in a high K(+) external solution ([K(+)]e) of 94mM were blocked with an IC50 of 117.8nM. Amiodarone's inhibitory action was contingent upon channel gating. Alanine-mutagenesis identified multiple residues directly or indirectly involved in amiodarone binding. The IC50 for the S6 aromatic Y652A mutation was increased to ∼20-fold that of WT IhERG, similar to the pore helical mutant S624A (∼22-fold WT control). The IC50 for F656A mutant IhERG was ∼17-fold its corresponding WT control. Computational docking using a MthK-based hERG model differentiated residues likely to interact directly with drug and those whose Ala mutation may affect drug block allosterically. The requirements for amiodarone block of aromatic residues F656 and Y652 within the hERG pore cavity are smaller than for other high affinity IhERG inhibitors, with relative importance to amiodarone binding of the residues investigated being S624A∼Y652A>F656A>V659A>G648A>T623A.

A TRP conductance modulates repolarization after sensory-dependent depolarization in Chlamydomonas reinhardtii.

Sensory integration is vital for motile organisms constantly exposed to changing surroundings. Chlamydomonas reinhardtii is a single-celled green alga found swimming in freshwater. In this type of alga, sensory input is first detected by membrane receptors located in the cell body, and then transduced to the beating cilia by membrane depolarization. Many components of the machinery associated with sensory integration in C. reinhardtii, such as chemoreceptors and repolarization-associated channels, are yet uncharacterized. TRP channels are known mediators for cellular sensing in animal cells and it has been suggested that the C. reinhardtii genome encodes for a set of TRP proteins. Here, by combining behavioral studies with electrophysiological experiments conducted on both population and single alga, we test whether TRP channel blockers affect algal swimming behavior. Our results suggest that a TRP conductance is associated to the repolarization that follows a depolarizing receptor potential, highlighting a primitive function of TRP proteins.

Cloning, identification and functional characterization of bovine free fatty acid receptor-1 (FFAR1/GPR40) in neutrophils.

Long chain fatty acids (LCFAs), which are ligands for the G-protein coupled receptor FFAR1 (GPR40), are increased in cow plasma after parturition, a period in which they are highly susceptible to infectious diseases. This study identified and analyzed the functional role of the FFAR1 receptor in bovine neutrophils, the first line of host defense against infectious agents. We cloned the putative FFAR1 receptor from bovine neutrophils and analyzed the sequence to construct a homology model. Our results revealed that the sequence of bovine FFAR1 shares 84% identity with human FFAR1 and 31% with human FFAR3/GPR41. Therefore, we constructed a homology model of bovine FFAR1 using human as the template. Expression of the bovine FFAR1 receptor in Chinese hamster ovary (CHO)-K1 cells increased the levels of intracellular calcium induced by the LCFAs, oleic acid (OA) and linoleic acid (LA); no increase in calcium mobilization was observed in the presence of the short chain fatty acid propionic acid. Additionally, the synthetic agonist GW9508 increased intracellular calcium in CHO-K1/bFFAR1 cells. OA and LA increased intracellular calcium in bovine neutrophils. Furthermore, GW1100 (antagonist of FFAR1) and U73122 (phospholipase C (PLC) inhibitor) reduced FFAR1 ligand-induced intracellular calcium in CHO-K1/bFFAR1 cells and neutrophils. Additionally, inhibition of FFAR1, PLC and PKC reduced the FFAR1 ligand-induced release of matrix metalloproteinase (MMP)-9 granules and reactive oxygen species (ROS) production. Thus, we identified the bovine FFAR1 receptor and demonstrate a functional role for this receptor in neutrophils activated with oleic or linoleic acid.