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.

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.

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.

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

Persistent anion conductances through GABA receptors (GABARs) 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 GABAR currents in 48 cells with an intracellular [Cl] of 65 mm and in the presence of 5 μm GABA.

Physiological role for GABA receptor desensitization in the induction of long-term potentiation at inhibitory synapses.

GABA receptors (GABARs) are pentameric ligand-gated ion channels distributed throughout the brain where they mediate synaptic and tonic inhibition. Following activation, these receptors undergo desensitization which involves entry into long-lived agonist-bound closed states. Although the kinetic effects of this state are recognised and its structural basis has been uncovered, the physiological impact of desensitization on inhibitory neurotransmission remains unknown.

Theoretical characterization of the shikimate 5-dehydrogenase reaction from Mycobacterium tuberculosis by hybrid QC/MM simulations and quantum chemical descriptors.

In this study, we have investigated the enzyme shikimate 5-dehydrogenase from the causative agent of tuberculosis, Mycobacterium tuberculosis. We have employed a mixture of computational techniques, including molecular dynamics, hybrid quantum chemical/molecular mechanical potentials, relaxed surface scans, quantum chemical descriptors and free-energy simulations, to elucidate the enzyme's reaction pathway. Overall, we find a two-step mechanism, with a single transition state, that proceeds by an energetically uphill hydride transfer, followed by an energetically downhill proton transfer.

Hydrogen evolution reaction mediated by an all-sulfur trinuclear nickel complex.

We report the synthesis and the characterization of a trinuclear nickel complex. Solid state and solution studies using X-ray diffraction, NMR and UV-vis spectroscopy highlight the square planar geometries around the metal centers and an all-sulfur coordination sphere. It exhibits significant electrocatalytic activity for hydrogen evolution in DMF using EtNHCl as the proton source.

Photoswitching mechanism of a fluorescent protein revealed by time-resolved crystallography and transient absorption spectroscopy.

Reversibly switchable fluorescent proteins (RSFPs) serve as markers in advanced fluorescence imaging. Photoswitching from a non-fluorescent off-state to a fluorescent on-state involves trans-to-cis chromophore isomerization and proton transfer. Whereas excited-state events on the ps timescale have been structurally characterized, conformational changes on slower timescales remain elusive.

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.

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

[Co(bapbpy)Cl] (: 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 H evolution in DMF with distinct mechanisms depending on the strength of the acid used as the proton source (p values ranging from 3.4 to 13.

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

The l,d-transpeptidase (Ldt) mediates resistance to most β-lactam antibiotics in this bacterium by replacing classical peptidoglycan polymerases. The catalytic Cys of Ldt 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 Ldt and showed that the inactivation profile is not determined by differences in drug binding ( [equilibrium dissociation constant] values in the 50 to 80 mM range).

Chromophore twisting in the excited state of a photoswitchable fluorescent protein captured by time-resolved serial femtosecond crystallography.

Chromophores absorb light in photosensitive proteins and thereby initiate fundamental biological processes such as photosynthesis, vision and biofluorescence. An important goal in their understanding is the provision of detailed structural descriptions of the ultrafast photochemical events that they undergo, in particular of the excited states that connect chemistry to biological function. Here we report on the structures of two excited states in the reversibly photoswitchable fluorescent protein rsEGFP2.

Crystal structures of a GABA-receptor chimera reveal new endogenous neurosteroid-binding sites.

γ-Aminobutyric acid receptors (GABARs) are vital for controlling excitability in the brain. This is emphasized by the numerous neuropsychiatric disorders that result from receptor dysfunction. A critical component of most native GABARs is the α subunit.

Proton reduction reaction catalyzed by homoleptic nickel bis-1,2-dithiolate complexes: Experimental and theoretical mechanistic investigations.

A series of homoleptic monoanionic nickel dithiolene complexes [Ni(bdt)](NBu), [Ni(tdt)](NBu), and [Ni(mnt)](NBu) containing the ligands benzene-1,2-dithiolate (bdt), toluene-3,4-dithiolate (tdt) and maleonitriledithiolate (mnt), respectively, have been employed as electrocatalysts in the hydrogen evolution reaction with trifluoroacetic acid as proton source in acetonitrile. All complexes were active catalysts with TONs reaching 113, 158 and 6 for [Ni(bdt)](NBu), [Ni(tdt)](NBu), and [Ni(mnt)](NBu), respectively. Faradaic yield for hydrogen evolution reaction reaches 88 % for , which also displays the minimal overpotential requirement value (467 mV) within the series.

An Algorithm for Adaptive QC/MM Simulations.

An algorithm is proposed for the simulation of molecular systems with hybrid quantum chemical (QC) and molecular mechanical (MM) potentials that permits the adaptive partitioning of the atoms in the system between QC and MM regions. In contrast to existing methods, the algorithm requires only a single QC calculation of the QC/MM system per energy calculation and yet has consistent energy and forces, which makes it suitable for geometry optimizations and molecular dynamics calculations within the microcanonical ensemble, in addition to other types of simulation. This article describes the algorithm and its implementation, presents some simple test cases with both semiempirical and density functional theory QC/MM potentials, and discusses perspectives for future work.

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.

Catalytic Mechanism of Peptidoglycan Deacetylase: A Computational Study.

Bacterial peptidoglycan deacetylase enzymes are potentially important targets for the design of new drugs. In pathogenic bacteria, they modify cell-wall peptidoglycan by removing the acetyl group, which makes the bacteria more resistant to the host's immune response and other forms of attack, such as degradation by lysozyme. In this study, we have investigated the mechanism of reaction of acetyl removal from a model substrate, the N-acetylglucosamine/N-acetylmuramic acid dimer, by peptidogylcan deacetylase from Helicobacter pylori.

Artificial Hydrogenases Based on Cobaloximes and Heme Oxygenase.

The insertion of cobaloxime catalysts in the heme-binding pocket of heme oxygenase (HO) yields artificial hydrogenases active for H evolution in neutral aqueous solutions. These novel biohybrids have been purified and characterized by using UV/visible and EPR spectroscopy. These analyses revealed the presence of two distinct binding conformations, thereby providing the cobaloxime with hydrophobic and hydrophilic environments, respectively.

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.

Hybrid Potential Simulation of the Acylation of Enterococcus faecium l,d-Transpeptidase by Carbapenems.

The l,d-transpeptidases, Ldts, catalyze peptidoglycan cross-linking in β-lactam-resistant mutant strains of several bacteria, including Enterococcus faecium and Mycobacterium tuberculosis. Although unrelated to the essential d,d-transpeptidases, which are inactivated by the β-lactam antibiotics, they are nevertheless inhibited by the carbapenem antibiotics, making them potentially useful targets in the treatment of some important diseases. In this work, we have investigated the acylation mechanism of the Ldt from E.

Aminoacetylation Reaction Catalyzed by Leucyl-tRNA Synthetase Operates via a Self-Assisted Mechanism Using a Conserved Residue and the Aminoacyl Substrate.

Leucyl-tRNA synthetase catalyzes attachment of leucine amino acid to its cognate tRNA. During the second, aminoacetylation, step of the reaction, the leucyl moiety is transferred from leucyl-adenylate to the terminal A76 adenosine of tRNA. In this work, we have investigated the aminoacetylation step catalyzed by leucyl-tRNA synthase, using ab initio quantum chemical/molecular mechanical hybrid potentials in conjunction with reaction-path-location algorithms and molecular dynamics free energy simulations.

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.

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.

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.

Force-induced chemical reactions on the metal centre in a single metalloprotein molecule.

Metalloproteins play indispensable roles in biology owing to the versatile chemical reactivity of metal centres. However, studying their reactivity in many metalloproteins is challenging, as protein three-dimensional structure encloses labile metal centres, thus limiting their access to reactants and impeding direct measurements. Here we demonstrate the use of single-molecule atomic force microscopy to induce partial unfolding to expose metal centres in metalloproteins to aqueous solution, thus allowing for studying their chemical reactivity in aqueous solution for the first time.