Differential effects of group III metabotropic glutamate receptors on spontaneous inhibitory synaptic currents in spine-innervating double bouquet and parvalbumin-expressing dendrite-targeting GABAergic interneurons in human neocortex.

Diverse neocortical GABAergic neurons specialize in synaptic targeting and their effects are modulated by presynaptic metabotropic glutamate receptors (mGluRs) suppressing neurotransmitter release in rodents, but their effects in human neocortex are unknown. We tested whether activation of group III mGluRs by L-AP4 changes GABAA receptor-mediated spontaneous inhibitory postsynaptic currents (sIPSCs) in 2 distinct dendritic spine-innervating GABAergic interneurons recorded in vitro in human neocortex. Calbindin-positive double bouquet cells (DBCs) had columnar "horsetail" axons descending through layers II-V innervating dendritic spines (48%) and shafts, but not somata of pyramidal and nonpyramidal neurons. Parvalbumin-expressing dendrite-targeting cell (PV-DTC) axons extended in all directions innervating dendritic spines (22%), shafts (65%), and somata (13%). As measured, 20% of GABAergic neuropil synapses innervate spines, hence DBCs, but not PV-DTCs, preferentially select spine targets. Group III mGluR activation paradoxically increased the frequency of sIPSCs in DBCs (to median 137% of baseline) but suppressed it in PV-DTCs (median 92%), leaving the amplitude unchanged. The facilitation of sIPSCs in DBCs may result from their unique GABAergic input being disinhibited via network effect. We conclude that dendritic spines receive specialized, diverse GABAergic inputs, and group III mGluRs differentially regulate GABAergic synaptic transmission to distinct GABAergic cell types in human cortex.

Tonic GABAA Receptor-Mediated Currents of Human Cortical GABAergic Interneurons Vary Amongst Cell Types.

Persistent anion conductances through GABAA receptors (GABAARs) are important modulators of neuronal excitability. However, it is currently unknown how the amplitudes of these currents vary among different cell types in the human neocortex, particularly among diverse GABAergic interneurons. We have recorded 101 interneurons in and near layer 1 from cortical tissue surgically resected from both male and female patients, visualized 84 of them and measured tonic GABAAR currents in 48 cells with an intracellular [Cl-] of 65 mm and in the presence of 5 µm GABA. We compare these tonic currents among five groups of interneurons divided by firing properties and four types of interneuron defined by axonal distributions; rosehip, neurogliaform, stalked-bouton, layer 2-3 innervating and a pool of other cells. Interestingly, the rosehip cell, a type of interneuron only described thus far in human tissue, and layer 2-3 innervating cells exhibit larger tonic currents than other layer 1 interneurons, such as neurogliaform and stalked-bouton cells; the latter two groups showing no difference. The positive allosteric modulators of GABAARs allopregnanolone and DS2 also induced larger current shifts in the rosehip and layer 2-3 innervating cells, consistent with higher expression of the δ subunit of the GABAAR in these neurons. We have also examined how patient parameters, such as age, seizures, type of cancer and anticonvulsant treatment may alter tonic inhibitory currents in human neurons. The cell type-specific differences in tonic inhibitory currents could potentially be used to selectively modulate cortical circuitry.SIGNIFICANCE STATEMENT Tonic currents through GABAA receptors (GABAARs) are a potential therapeutic target for a number of neurologic and psychiatric conditions. Here, we show that these currents in human cerebral cortical GABAergic neurons display cell type-specific differences in their amplitudes which implies differential modulation of their excitability. Additionally, we examine whether the amplitudes of the tonic currents measured in our study show any differences between patient populations, finding some evidence that age, seizures, type of cancer, and anticonvulsant treatment may alter tonic inhibition in human tissue. These results advance our understanding of how pathology affects neuronal excitability and could potentially be used to selectively modulate cortical circuitry.

pDynamo3 Molecular Modeling and Simulation Program.

pDynamo3 is the first formal version of the Dynamo molecular modeling and simulation library that is written in Python 3. It follows from the previous pDynamo versions written in Python 2, the first of which was released in 2007. Both pDynamo and its predecessor, fDynamo, were designed with the aim of providing easy-to-use and flexible frameworks for performing molecular simulations at an atomistic level with a special emphasis on those employing hybrid quantum chemical and molecular mechanical potential methods. Although the use of pDynamo3 is quite similar to that of pDynamo2, it has added significant new capability and also undergone extensive restructuring that will make it much easier to extend with new functionality. The pDynamo3 code is issued under the GNU general public license at https://github.com/pdynamo/pDynamo3 with additional information on the pDynamo website https:www.pdynamo.org.

Density Fitting for QC/MM Interactions.

Density fitting is a standard technique in quantum chemistry as it helps to accelerate certain parts of a calculation, such as the computation of the electron repulsion energy, without significant loss of accuracy. This paper explores the effectiveness of this technique when it is extended to treat interactions with external electrostatic potentials, in particular those that arise from the environment in hybrid quantum chemical/molecular mechanical calculations. It is found that fitted densities are able to well reproduce the energies, forces, and properties obtained using unfitted densities, as long as a suitable operator is employed for the fitting. It is expected that this precision would be improved by the development of basis sets specifically designed to treat these types of interactions and that the use of this approximation could lead to substantial speed-ups in large hybrid potential simulations.

Electrocatalytic Hydrogen Evolution with a Cobalt Complex Bearing Pendant Proton Relays: Acid Strength and Applied Potential Govern Mechanism and Stability.

[Co(bapbpy)Cl]+ (bapbpy: 6,6'-bis(2-aminopyridyl)-2,2'-bipyridine) is a polypyridyl cobalt(II) complex bearing both a redox-active bipyridine ligand and pendant proton relays. This compound catalyzes electro-assisted H2 evolution in DMF with distinct mechanisms depending on the strength of the acid used as the proton source (pKa values ranging from 3.4 to 13.5 in DMF) and the applied potential. Electrochemical studies combining cyclic voltammetry and bulk electrolysis measurements enabled one to bring out four distinct catalytic processes. Where applicable, relevant kinetic information were obtained using either foot-of-the-wave analysis (FOWA) or analytical treatment of bulk electrolysis experiments. Among the different catalytic pathways identified in this study, a clear relationship between the catalyst performances and stability was evidenced. These results draw attention to a number of interesting considerations and may help in the development of future adequately designed catalysts.

Modeling the Hydrolysis of Iron-Sulfur Clusters.

Iron-sulfur (FeS) clusters are essential metal cofactors involved in a wide variety of biological functions. Their catalytic efficiency, biosynthesis, and regulation depend on FeS stability in aqueous solution. Here, molecular modeling is used to investigate the hydrolysis of an oxidized (ferric) mononuclear FeS cluster by bare dissociation and water substitution mechanisms in neutral and acidic solution. First, approximate electronic structure descriptions of FeS reactions by density functional theory are validated against high-level wave function CCSD(T) calculations. Solvation contributions are included by an all-atom model with hybrid quantum chemical/molecular mechanical (QM/MM) potentials and enhanced sampling molecular dynamics simulations. The free energy profile obtained for FeS cluster hydrolysis indicates that the hybrid functional M06 together with an implicit solvent correction capture the most important aspects of FeS cluster reactivity in aqueous solution. Then, 20 reaction channels leading to two consecutive Fe-S bond ruptures were explored with this calibrated model. For all protonation states, nucleophilic substitution with concerted bond breaking and forming to iron is the preferred mechanism, both kinetic and thermodynamically. In neutral solution, proton transfer from water to the sulfur leaving group is also concerted. Dissociative reactions show higher barriers and will not be relevant for FeS reactivity when exposed to solvent. These hydrolysis mechanisms may help to explain the stability and catalytic mechanisms of FeS clusters of multiple sizes and proteins.

Negative Impact of Carbapenem Methylation on the Reactivity of ß-Lactams for Cysteine Acylation as Revealed by Quantum Calculations and Kinetic Analyses.

The Enterococcus faecium l,d-transpeptidase (Ldtfm) mediates resistance to most ß-lactam antibiotics in this bacterium by replacing classical peptidoglycan polymerases. The catalytic Cys of Ldtfm is rapidly acylated by ß-lactams belonging to the carbapenem class but not by penams or cephems. We previously reported quantum calculations and kinetic analyses for Ldtfm and showed that the inactivation profile is not determined by differences in drug binding (KD [equilibrium dissociation constant] values in the 50 to 80 mM range). In this study, we analyzed the reaction of a Cys sulfhydryl with various ß-lactams in the absence of the enzyme environment in order to compare the intrinsic reactivity of drugs belonging to the penam, cephem, and carbapenem classes. For this purpose, we synthesized cyclic Cys-Asn (cCys-Asn) to generate a soluble molecule with a sulfhydryl closely mimicking a cysteine in a polypeptide chain, thereby avoiding free reactive amino and carboxyl groups. Computational studies identified a thermodynamically favored pathway involving a concerted rupture of the ß-lactam amide bond and formation of an amine anion. Energy barriers indicated that the drug reactivity was the highest for nonmethylated carbapenems, intermediate for methylated carbapenems and cephems, and the lowest for penams. Electron-withdrawing groups were key reactivity determinants by enabling delocalization of the negative charge of the amine anion. Acylation rates of cCys-Asn determined by spectrophotometry revealed the same order in the reactivity of ß-lactams. We concluded that the rate of Ldtfm acylation is largely determined by the ß-lactam reactivity with one exception, as the enzyme catalytic pocket fully compensated for the detrimental effect of carbapenem methylation.

Norovirus RNA-dependent RNA polymerase: A computational study of metal-binding preferences.

Norovirus (NV) RNA-dependent RNA polymerase (RdRP) is essential for replicating the genome of the virus, which makes this enzyme a key target for the development of antiviral agents against NV gastroenteritis. In this work, a complex of NV RdRP bound to manganese ions and an RNA primer-template duplex was investigated using X-ray crystallography and hybrid quantum chemical/molecular mechanical simulations. Experimentally, the complex crystallized in a tetragonal crystal form. The nature of the primer/template duplex binding in the resulting structure indicates that the complex is a closed back-tracked state of the enzyme, in which the 3'-end of the primer occupies the position expected for the post-incorporated nucleotide before translocation. Computationally, it is found that the complex can accept a range of divalent metal cations without marked distortions in the active site structure. The highest binding energy is for copper, followed closely by manganese and iron, and then by zinc, nickel, and cobalt. Proteins 2017; 85:1435-1445. © 2017 Wiley Periodicals, Inc.

Structural Determinants of Improved Fluorescence in a Family of Bacteriophytochrome-Based Infrared Fluorescent Proteins: Insights from Continuum Electrostatic Calculations and Molecular Dynamics Simulations.

Using X-ray crystallography, continuum electrostatic calculations, and molecular dynamics simulations, we have studied the structure, protonation behavior, and dynamics of the biliverdin chromophore and its molecular environment in a series of genetically engineered infrared fluorescent proteins (IFPs) based on the chromophore-binding domain of the Deinococcus radiodurans bacteriophytochrome. Our study suggests that the experimentally observed enhancement of fluorescent properties results from the improved rigidity and planarity of the biliverdin chromophore, in particular of the first two pyrrole rings neighboring the covalent linkage to the protein. We propose that the increases in the levels of both motion and bending of the chromophore out of planarity favor the decrease in fluorescence. The chromophore-binding pocket in some of the studied proteins, in particular the weakly fluorescent parent protein, is shown to be readily accessible to water molecules from the solvent. These waters entering the chromophore region form hydrogen bond networks that affect the otherwise planar conformation of the first three rings of the chromophore. On the basis of our simulations, the enhancement of fluorescence in IFPs can be achieved either by reducing the mobility of water molecules in the vicinity of the chromophore or by limiting the interactions of the nearby protein residues with the chromophore. Finally, simulations performed at both low and neutral pH values highlight differences in the dynamics of the chromophore and shed light on the mechanism of fluorescence loss at low pH.

Experimental and Theoretical Insight into Electrocatalytic Hydrogen Evolution with Nickel Bis(aryldithiolene) Complexes as Catalysts.

A series of neutral and monoanionic nickel dithiolene complexes with p-methoxyphenyl-substituted 1,2-dithiolene ligands have been prepared and characterized with physicochemical methods. Two of the complexes, the monoanion of the symmetric [Ni{S2C2(Ph-p-OCH3)2}2] (3(-)) with NBu4(+) as a counterion and the neutral asymmetric [Ni{S2C2(Ph)(Ph-p-OCH3)}2] (2), have been structurally characterized by single-crystal X-ray crystallography. All complexes have been employed as proton-reducing catalysts in N,N-dimethylformamide with trifluoroacetic acid as the proton source. The complexes are active catalysts with good faradaic yields, reaching 83% for 2 but relatively high overpotential requirements (0.91 and 1.55 V measured at the middle of the catalytic wave for two processes observed depending on the different routes of the mechanism). The similarity of the experimental data regardless of whether the neutral or anionic form of the complexes is used indicates that the neutral form acts as a precatalyst. On the basis of detailed density functional theory calculations, the proposed mechanism reveals two different main routes after protonation of the dianion of the catalyst in accordance with the experimental data, indicating the role of the concentration of the acid and the influence of the methoxy groups. Protonation at sulfur seems be more favorable than that at the metal, which is in marked contrast with the catalytic mechanism proposed for analogous cobalt dithiolene complexes.

Multiple proton relay routes in the reaction mechanism of RNAP II: assessing the effect of structural model.

RNA polymerase II catalyzes the nucleotidyl transfer reaction for messenger RNA synthesis in eukaryotes. Two crystal structures of this system have been resolved, each with its own defects in the coordination sphere of Mg(2+) (A) resulting from chemical modifications. We have used both structures and also a novel hybrid of the two that allows a better exploration of the parts of configuration space that reflect substrate-enzyme interactions. MD and QM/MM calculations have been performed, the latter with the semiempirical AM1/d-PhoT method, calibrated against density functional theory. Reaction path scans in 1-D provided insights about the role of Mg(2+) (A) which turns out to be more structural than catalytic. In contrast, 1-D scans of the incorporation of the nucleotidyl group yielded barriers that were much too high, necessitating the use of 2-D reaction coordinates. Three different proton acceptors for the initial reaction step were examined. For those models based on the two PDB structures the 2-D scans continued to yield very high barriers, indicating that the reaction is unlikely to proceed from these configurations. On the other hand, two hybrid models, chosen from the early and late parts of a 12 ns molecular dynamics simulation yielded greatly reduced barriers in the range of ∼ 17 to ∼ 27 kcal/mol for the three proton acceptors, as compared to the experimental estimate of 18 kcal/mol. The final step, release of pyrophosphate, was found to be facile. Our overall mechanism involves only active site residues or water without the need for external reactive agents such as the hydroxide ion previously proposed.

Mechanism of activation of elongation factor Tu by ribosome: catalytic histidine activates GTP by protonation.

Elongation factor Tu (EF-Tu) is central to prokaryotic protein synthesis as it has the role of delivering amino-acylated tRNAs to the ribosome. Release of EF-Tu, after correct binding of the EF-Tu:aa-tRNA complex to the ribosome, is initiated by GTP hydrolysis. This reaction, whose mechanism is uncertain, is catalyzed by EF-Tu, but requires activation by the ribosome. There have been a number of mechanistic proposals, including those spurred by a recent X-ray crystallographic analysis of a ribosome:EF-Tu:aa-tRNA:GTP-analog complex. In this work, we have investigated these and alternative hypotheses, using high-level quantum chemical/molecular mechanical simulations for the wild-type protein and its His85Gln mutant. For both proteins, we find previously unsuggested mechanisms as being preferred, in which residue 85, either His or Gln, directly assists in the reaction. Analysis shows that the RNA has a minor catalytic effect in the wild-type reaction, but plays a significant role in the mutant by greatly stabilizing the reaction's transition state. Given the similarity between EF-Tu and other members of the translational G-protein family, it is likely that these mechanisms of ribosome-activated GTP hydrolysis are pertinent to all of these proteins.

Technical advances in molecular simulation since the 1980s.

This review describes how the theory and practice of molecular simulation have evolved since the beginning of the 1980s when the author started his career in this field. The account is of necessity brief and subjective and highlights the changes that the author considers have had significant impact on his research and mode of working.

Pcetk: A pDynamo-based Toolkit for Protonation State Calculations in Proteins.

Pcetk (a pDynamo-based continuum electrostatic toolkit) is an open-source, object-oriented toolkit for the calculation of proton binding energetics in proteins. The toolkit is a module of the pDynamo software library, combining the versatility of the Python scripting language and the efficiency of the compiled languages, C and Cython. In the toolkit, we have connected pDynamo to the external Poisson-Boltzmann solver, extended-MEAD. Our goal was to provide a modern and extensible environment for the calculation of protonation states, electrostatic energies, titration curves, and other electrostatic-dependent properties of proteins. Pcetk is freely available under the CeCILL license, which is compatible with the GNU General Public License. The toolkit can be found on the Web at the address http://github.com/mfx9/pcetk. The calculation of protonation states in proteins requires a knowledge of pKa values of protonatable groups in aqueous solution. However, for some groups, such as protonatable ligands bound to protein, the pKa aq values are often difficult to obtain from experiment. As a complement to Pcetk, we revisit an earlier computational method for the estimation of pKa aq values that has an accuracy of ± 0.5 pKa-units or better. Finally, we verify the Pcetk module and the method for estimating pKa aq values with different model cases.

Ferric-Thiolate Bond Dissociation Studied with Electronic Structure Calculations.

The stability and reactivity of iron-sulfur clusters are fundamental properties for the biological function of these prosthetic groups. Here, we investigate the ferric-thiolate bond dissociation of model iron-sulfur tetrahedral complexes with high-level ab initio multiconfigurational electronic structure calculations. We find that the reaction mechanism is homolytic with a spin-crossing from the sextet state in the reactant to quartet state in the product. We also compare several density functionals and semiempirical configuration interaction with the high-level ab initio results to find an accurate but computationally more efficient method to describe the reaction. The functionals M06 and those based on the OPTX exchange functional show the best performance and may reasonably describe the various electron correlation effects involved in ferric-thiolate bond dissociation.

Characterization of the divalent metal binding site of bacterial polysaccharide deacetylase using crystallography and quantum chemical calculations.

Peptidoglycan deacetlyase (HP0310, HpPgdA) from the gram-negative pathogen Helicobacter pylori, is the enzyme responsible for a peptidoglycan modification that counteracts the host immune response. In a recent study, we determined the crystallographic structure of the enzyme, which is a homo-tetramer (Shaik et al., PloS One 2011;6:e19207). The metal-binding site, which is essential for the enzyme's catalytic activity, is visible within the structure, but we were unable to identify the nature of the metal itself. In this study, we have obtained a higher-resolution crystal structure of the enzyme, which shows that the ion bound is, in fact, zinc. Analysis of the structure of the four sites, one per monomer, and quantum chemical calculations of models of the site in the presence of different divalent metal ions show an intrinsic preference for zinc, but also significant flexibility of the site so that binding of other ions can eventually occur.

300-Fold increase in production of the Zn2+-dependent dechlorinase TrzN in soluble form via apoenzyme stabilization.

Microbial metalloenzymes constitute a large library of biocatalysts, a number of which have already been shown to catalyze the breakdown of toxic chemicals or industrially relevant chemical transformations. However, while there is considerable interest in harnessing these catalysts for biotechnology, for many of the enzymes, their large-scale production in active, soluble form in recombinant systems is a significant barrier to their use. In this work, we demonstrate that as few as three mutations can result in a 300-fold increase in the expression of soluble TrzN, an enzyme from Arthrobacter aurescens with environmental applications that catalyzes the hydrolysis of triazine herbicides, in Escherichia coli. Using a combination of X-ray crystallography, kinetic analysis, and computational simulation, we show that the majority of the improvement in expression is due to stabilization of the apoenzyme rather than the metal ion-bound holoenzyme. This provides a structural and mechanistic explanation for the observation that many compensatory mutations can increase levels of soluble-protein production without increasing the stability of the final, active form of the enzyme. This study provides a molecular understanding of the importance of the stability of metal ion free states to the accumulation of soluble protein and shows that differences between apoenzyme and holoenzyme structures can result in mutations affecting the stability of either state differently.

Theoretical modeling of low-energy electronic absorption bands in reduced cobaloximes.

The reduced Co(I) states of cobaloximes are powerful nucleophiles that play an important role in the hydrogen-evolving catalytic activity of these species. In this work we analyze the low-energy electronic absorption bands of two cobaloxime systems experimentally and use a variety of density functional theory and molecular orbital ab initio quantum chemical approaches. Overall we find a reasonable qualitative understanding of the electronic excitation spectra of these compounds but show that obtaining quantitative results remains a challenging task.

Cobaloxime-based artificial hydrogenases.

Cobaloximes are popular H2 evolution molecular catalysts but have so far mainly been studied in nonaqueous conditions. We show here that they are also valuable for the design of artificial hydrogenases for application in neutral aqueous solutions and report on the preparation of two well-defined biohybrid species via the binding of two cobaloxime moieties, {Co(dmgH)2} and {Co(dmgBF2)2} (dmgH2 = dimethylglyoxime), to apo Sperm-whale myoglobin (SwMb). All spectroscopic data confirm that the cobaloxime moieties are inserted within the binding pocket of the SwMb protein and are coordinated to a histidine residue in the axial position of the cobalt complex, resulting in thermodynamically stable complexes. Quantum chemical/molecular mechanical docking calculations indicated a coordination preference for His93 over the other histidine residue (His64) present in the vicinity. Interestingly, the redox activity of the cobalt centers is retained in both biohybrids, which provides them with the catalytic activity for H2 evolution in near-neutral aqueous conditions.

QM/MM through the 1990s: The First Twenty Years of Method Development and Applications.

The 2013 Nobel Prize in Chemistry was awarded to the authors of the first two publications utilizing the concept of combined quantum mechanical and molecular mechanical (QM/MM) methods. In celebrating this great event in computational chemistry, we review the early development of combined QM/MM techniques and the associated events that took place through the mid-1990s. We also offer some prospects for the future development of quantum mechanical techniques for macromolecular systems.