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 H₂ evolution in DMF with distinct mechanisms depending on the strength of the acid used as the proton source (pK_a 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.

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)2](NBu4), [Ni(tdt)2](NBu4), and [Ni(mnt)2](NBu4) containing the ligands benzene-1,2-dithiolate (bdt2-), toluene-3,4-dithiolate (tdt2-) and maleonitriledithiolate (mnt2-), 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)2](NBu4), [Ni(tdt)2](NBu4), and [Ni(mnt)2](NBu4), respectively. Faradaic yield for hydrogen evolution reaction reaches 88 % for 2- , which also displays the minimal overpotential requirement value (467 mV) within the series. Two pathways for H2 evolution can be hypothesized that differ on on the sequence of protonation and reduction steps. DFT calculations are in agreement with experimental data and indicate that protonation at sulfur follows reduction to the dianion. Hydrogen evolves from the direduced-diprotonated form via a highly distorted nickel hydride intermediate. The effects of acid strength and concentration in the hydrogen-evolving mechanism are also discussed.

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.

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.

Simple Link Atom Saccharide Hybrid (SLASH) Treatment for Glycosidic Bonds at the QM/MM Boundary

We investigated link atom approaches for treating the polar C-O bond with particular reference to the glycosidic bond found in complex carbohydrates. We show that cutting this bond after the oxygen in the QM region and saturating the QM system with a hydrogen link atom leads to greater conformational and configurational accuracy at the boundary compared with cutting the bond before oxygen and saturating the QM system with a halogen link atom to represent the oxygen. Furthermore, we find that balancing the MM atom charges and redistributing the boundary atom charges at the QM/MM boundary minimizes the effect of the link atom, both energetically and structurally. This is illustrated via a series of calculations on a set of carbohydrate and carbohydrate-like model compounds. Finally, we confirm the validity of our model by performing molecular dynamics simulations for a typical disaccharide model compound in water. Our postsimulation conformational and configurational analyses show that the oxygen-to-water hydrogen pair distribution functions and the Φ,Ψ distributions at the glycosidic boundary between the quantum and classical regions compare favorably with results obtained from complete QM and complete MM treatments of the saccharide.

The pDynamo Program for Molecular Simulations using Hybrid Quantum Chemical and Molecular Mechanical Potentials

The pDynamo program has been developed for the simulation of molecular systems using hybrid quantum chemical (QC) and molecular mechanical (MM) potentials. pDynamo is written in a mixture of the computer languages Python and C and is a successor to the previous version of Dynamo, now denoted fDynamo, that was written in Fortran 90 (J. Comput. Chem. 2000, 21, 1088). The current version of Dynamo has a similar range of functionality to the older one but extends it in some significant ways, including the addition of a density functional theory QC capability. This paper gives a general description of pDynamo and outlines some of the advantages and disadvantages that have been encountered in switching computer languages. Some technical aspects of the implementation of pDynamo's algorithms are also discussed and illustrated with the results of example calculations. pDynamo is available on the Web at the address http://www.pdynamo.org and is released under the CeCILL license which is equivalent to the GNU general public license but conforms to the principles of French law.

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. As a proof-of-principle, we demonstrate two chemical reactions for the FeS₄ centre in rubredoxin: electrophilic protonation and nucleophilic ligand substitution. Our results show that protonation and ligand substitution result in mechanical destabilization of the FeS₄ centre. Quantum chemical calculations corroborated experimental results and revealed detailed reaction mechanisms. We anticipate that this novel approach will provide insights into chemical reactivity of metal centres in metalloproteins under biologically more relevant conditions.

The investigation of the chemical reactivity of metal centres in metalloproteins in aqueous solution is challenging. Here, the authors demonstrate the use of single molecule force spectroscopy to study the chemical reactivity of the iron-sulfur centre in rubredoxin in aqueous solution.

Dye-sensitized PS-b-P2VP-templated nickel oxide films for photoelectrochemical applications

Moving from homogeneous water-splitting photocatalytic systems to photoelectrochemical devices requires the preparation and evaluation of novel p-type transparent conductive photoelectrode substrates. We report here on the sensitization of polystyrene-block-poly-(2-vinylpyridine) (PS-b-P2VP) diblock copolymer-templated NiO films with an organic push-pull dye. The potential of these new templated NiO film preparations for photoelectrochemical applications is compared with NiO material templated by F108 triblock copolymers. We conclude that NiO films are promising materials for the construction of dye-sensitized photocathodes to be inserted into photoelectrochemical (PEC) cells. However, a combined effort at the interface between materials science and molecular chemistry, ideally funded within a Global Artificial Photosynthesis Project, is still needed to improve the overall performance of the photoelectrodes and progress towards economically viable PEC devices.

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.

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.

All-atom empirical potential for molecular modeling and dynamics studies of proteins

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent-solvent, solvent-solute, and solute-solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volumes, and crystal pressures and structures were used in the optimization. The resulting protein parameters were tested by applying them to noncyclic tripeptide crystals, cyclic peptide crystals, and the proteins crambin, bovine pancreatic trypsin inhibitor, and carbonmonoxy myoglobin in vacuo and in crystals. A detailed analysis of the relationship between the alanine dipeptide potential energy surface and calculated protein φ, χ angles was made and used in optimizing the peptide group torsional parameters. The results demonstrate that use of ab initio structural and energetic data by themselves are not sufficient to obtain an adequate backbone representation for peptides and proteins in solution and in crystals. Extensive comparisons between molecular dynamics simulations and experimental data for polypeptides and proteins were performed for both structural and dynamic properties. Energy minimization and dynamics simulations for crystals demonstrate that the latter are needed to obtain meaningful comparisons with experimental crystal structures. The presented parameters, in combination with the previously published CHARMM all-atom parameters for nucleic acids and lipids, provide a consistent set for condensed-phase simulations of a wide variety of molecules of biological interest.

Hydrolysis of DFP and the nerve agent (S)-sarin by DFPase proceeds along two different reaction pathways: implications for engineering bioscavengers

Organophosphorus (OP) nerve agents such as (S)-sarin are among the most highly toxic compounds that have been synthesized. Engineering enzymes that catalyze the hydrolysis of nerve agents (“bioscavengers”) is an emerging prophylactic approach to diminish their toxic effects. Although its native function is not known, diisopropyl fluorophosphatase (DFPase) from Loligo vulgaris catalyzes the hydrolysis of OP compounds. Here, we investigate the mechanisms of diisopropylfluorophosphate (DFP) and (S)-sarin hydrolysis by DFPase with quantum mechanical/molecular mechanical umbrella sampling simulations. We find that the mechanism for hydrolysis of DFP involves nucleophilic attack by Asp229 on phosphorus to form a pentavalent intermediate. P–F bond dissociation then yields a phosphoacyl enzyme intermediate in the rate-limiting step. The simulations suggest that a water molecule, coordinated to the catalytic Ca²⁺, donates a proton to Asp121 and then attacks the tetrahedral phosphoacyl intermediate to liberate the diisopropylphosphate product. In contrast, the calculated free energy barrier for hydrolysis of (S)-sarin by the same mechanism is highly unfavorable, primarily because of the instability of the pentavalent phosphoenzyme species. Instead, simulations suggest that hydrolysis of (S)-sarin proceeds by a mechanism in which Asp229 could activate an intervening water molecule for nucleophilic attack on the substrate. These findings may lead to improved strategies for engineering DFPase and related six-bladed β-propeller folds for more efficient degradation of OP compounds.

GTKDynamo: a PyMOL plug-in for QC/MM hybrid potential simulations

Hybrid quantum chemical/molecular mechanical (QCMM) potentials are very powerful tools for molecular simulation. They are especially useful for studying processes in condensed phase systems, such as chemical reactions that involve a relatively localized change in electronic structure and where the surrounding environment contributes to these changes but can be represented with more computationally efficient functional forms. Despite their utility, however, these potentials are not always straightforward to apply since the extent of significant electronic structure changes occurring in the condensed phase process may not be intuitively obvious. To facilitate their use, we have developed an open-source graphical plug-in, GTKDynamo that links the PyMOL visualization program and the pDynamo QC/MM simulation library. This article describes the implementation of GTKDynamo and its capabilities and illustrates its application to QC/MM simulations.

Pentacoordinated phosphorus revisited by high-level QM/MM calculations

Enzymes catalyzing phosphoryl transfer reactions are extremely efficient and are involved in crucial biochemical processes. The mechanisms of these enzymes are complex due to the diversity of substrates that are involved. The reaction can proceed through a pentacoordinated phosphorus species that is either a stable intermediate or a transition state (TS). Because of this, the first X-ray structure of a pentacoordinated phosphorus intermediate in the beta-phosphoglucomutase enzyme aroused great interest but also much controversy. To provide new insights into the nature of that structure, we have determined the reaction path of the phosphorylation step using high-level QM/MM calculations, and have also calculated the geometry of a complex with a transition state analogue (TSA) that has been suggested to be the actual species in the crystal. The protein crystalline environment has been modeled so as to mimic the experimental conditions. We conclude that the pentacoordinated phosphorus formed in this enzyme is not a stable species but a TS, which gives an activation energy for phosphorylation in agreement with kinetic results. We also show that the TSA is a good mimic of the true TS. We have performed a new crystallographic refinement of the original diffraction map of the pentacoordinated phosphorus structure with the MgF(3)(-) TSA. The new fit improves significantly with respect to the original one, which strongly supports that Allen and coworkers wrongly assigned the X-ray structure to a pentavalent phosphorane.

Editing mechanism of aminoacyl-tRNA synthetases operates by a hybrid ribozyme/protein catalyst

Aminoacyl-tRNA synthetases (aaRSs) are critical for the translational process, catalyzing the attachment of specific amino acids to their cognate tRNAs. To ensure formation of the correct aminoacyl-tRNA, and thereby enhance the reliability of translation, several aaRSs have an editing capability that hinders formation of misaminoacylated tRNAs. We investigated theoretically the mechanism of the editing reaction for a class I enzyme, leucyl-tRNA synthetase (LeuRS), complexed with a misaminoacylated tRNA(Leu), employing ab initio hybrid quantum mechanical/molecular mechanical potentials in conjunction with molecular dynamics simulations. It is shown that the water molecule that acts as the nucleophile in the editing reaction is activated by a 3'-hydroxyl group at the 3'-end of tRNA(Leu) and that the O2' atom of the leaving group of the substrate is capped by one of the water's hydrogen atoms. Thus, it is shown that editing is a self-cleavage reaction of the tRNA and so it is the tRNA, and not the protein, that drives the reaction. The protein does, however, have an important stabilizing effect on some high-energy intermediates along the reaction path, which is more efficient than the ribozyme would be alone. This indicates that editing is achieved by a novel "hybrid ribozyme/protein catalyst". Analysis of existing experimental data and additional modeling shows that this ribozymal mechanism appears to be widespread, occurring in the ribosome as well as in other aaRSs. It also suggests transitional forms that could have played an important role in the RNA world hypothesis for the origin of life.

Cyclopentadienyl ruthenium-nickel catalysts for biomimetic hydrogen evolution: electrocatalytic properties and mechanistic DFT studies

The new dinuclear nickel-ruthenium complexes [Ni(xbsms)RuCp(L)][PF(6)] (H(2)xbsms = 1,2-bis(4-mercapto-3,3-dimethyl-2-thiabutyl)benzene; Cp(-) = cyclopentadienyl; L = DMSO, CO, PPh(3), and PCy(3)) are reported and are bioinspired mimics of NiFe hydrogenases. These compounds were characterized by X-ray diffraction techniques and display novel structural motifs. Interestingly, [Ni(xbsms)RuCpCO][PF(6)] is stereochemically nonrigid in solution and an isomerization mechanism was derived with the help of density functional theory (DFT) calculations. Because of an increased electron density on the metal centers [Eur. J. Inorg. Chem. 2007, 18, 2613-2626] with respect to the previously described [Ni(xbsms)Ru(CO)(2)Cl(2)] and [Ni(xbsms)Ru(p-cymene)Cl](+) complexes, [Ni(xbsms)RuCp(dmso)][PF(6)] catalyzes hydrogen evolution from Et(3)NH(+) in DMF with an overpotential reduced by 180 mV and thus represents the most efficient NiFe hydrogenase functional mimic. DFT calculations were carried out with several methods to investigate the catalytic cycle and, coupled with electrochemical measurements, allowed a mechanism to be proposed. A terminal or bridging hydride derivative was identified as the active intermediate, with the structure of the bridging form similar to that of the Ni-C active state of NiFe hydrogenases.

Photoconversion of the fluorescent protein EosFP: a hybrid potential simulation study reveals intersystem crossings

Fluorescent proteins undergoing green to red photoconversion have proved to be essential tools in cell biology, notably in superlocalization nanoscopy. However, the exact mechanism governing photoconversion, which overall involves irreversible cleavage of the protein backbone and elongation of the chromophore pi-conjugation, remains unclear. In this paper we present a theoretical investigation of the photoconversion reaction in the fluorescent protein EosFP, using excited-state hybrid quantum chemical and molecular mechanical potentials, in conjunction with reaction-path-finding techniques. Our results reveal a mechanism in which the hydroxybenzylidene moiety of the chromophore remains protonated and there is an excited state proton transfer from His62 to Phe61 that promotes peptide bond cleavage. Excitation of the neutral green form of EosFP to the first singlet excited state is followed by two intersystem crossing events, first to a triplet state and then back to the ground state singlet surface. From there, a number of rearrangements occur in the ground state and lead to the red form. Analyses of the structures and energies of the intermediates along the reaction path enable us to identify the critical role of the chromophore environment in promoting photoinduced backbone cleavage. Possible ways in which photoconvertible fluorescent proteins can be engineered to facilitate photoconversion are considered.

High glucose induced endothelial cell growth inhibition is associated with an increase in TGFbeta1 secretion and inhibition of Ras prenylation via suppression of the mevalonate pathway

OBJECTIVE

Ras proteins are known to affect cellular growth and function. The influence of the prenylation status of Ras on the observed changes in endothelial cell growth under high glucose conditions has not previously been examined.

METHODS

Human umbilical vein endothelial cells were exposed to normal or high glucose conditions for 72 h. They were then examined for proliferative and hypertrophic effects, transforming growth factor beta(1) (TGFbeta(1)) release, and phosphorylated p38 expression. The importance of prenylation was explored by the addition of mevalonate, isoprenoids or farnesyltransferase inhibitors to control the high glucose media and by measuring changes induced by high glucose and exogenous TGFbeta(1) in Ras prenylation and farnesyltransferase activity. Kidneys from diabetic rats treated with atorvastatin were also compared to specimens from untreated animals and the expression of the Ras effector p-Akt examined.

RESULTS

High glucose conditions caused a reduction in cell number. This was reversed in the presence of mevalonate or farnesylpyrophosphate (FPP), suggesting that the cell growth abnormalities observed are due to high glucose induced inhibition of the mevalonate pathway and subsequent prenylation of proteins. Endothelial cells exposed to high glucose increased their secretion of TGFbeta(1) and the phosphorylation of p38 both of which were reversed by concurrent exposure to FPP. A reduction in farnesyltransferase activity was observed after exposure to both high glucose and TGFbeta(1). Exposure to a farnesyltransferase inhibitor in control conditions mimicked the growth response observed with high glucose exposure and prenylated Ras was reduced by exposure to both high glucose and TGFbeta(1). Finally, interruption of the mevalonate pathway with a statin reduced the expression of p-Akt in diabetic rat kidneys.

CONCLUSION

This study demonstrates that high glucose induced significant alterations in endothelial cell growth by inhibition of the mevalonate pathway, which subsequently mediates the increase in TGFbeta(1) and inhibition of Ras prenylation.

Binding of calcium and other metal ions to the EF-hand loops of calmodulin studied by quantum chemical calculations and molecular dynamics simulations

Calcium ion binding by the four EF-hand motifs of the protein calmodulin (CaM) is a central event in Ca2+-based cellular signaling. To understand molecular details of this complex process, isolated Ca2+-binding loops can be studied, by use of both experiments and calculations. In this work, we explore the metal specificity of the four Ca2+-binding loops of CaM using density functional theory (DFT) quantum chemical calculations and molecular dynamics simulations. We study CaM complexes with the physiologically important ions of calcium (Ca2+) and magnesium (Mg2+) and also with two other ions, strontium (Sr2+) and lanthanum (La3+). The former is of interest in the area of radioactive waste bioremediation, whereas the latter is often used as a probe of Ca2+-binding sites. We obtain intrinsic metal ion-loop binding energies as well as their components: vacuum, charge-transfer, solvation, entropy, and deformation terms. A detailed analysis of the results reveals that the total binding energy depends on a delicate balance among these energy components. They, in turn, are determined by the cation's charge and size as well as the amino acid composition and flexibility of the loops and the identity of the metal-chelating residues.

Quantification of the rat spinal microglial response to peripheral nerve injury as revealed by immunohistochemical image analysis and flow cytometry

Microgliosis is implicated in the pathophysiology of several neurological disorders, including neuropathic pain. Consequently, perturbation of microgliosis is a mechanistic and drug development target in neuropathic pain, which highlights the requirement for specific, sensitive and reproducible methods of microgliosis measurement. In this study, we used the spinal microgliosis associated with L5 spinal nerve transection and minocycline-induced attenuation thereof to: (1) evaluate novel software based semi-quantitative image analysis paradigms for the assessment of immunohistochemical images. Microgliosis was revealed by immunoreactivity to OX42. Several image analysis paradigms were assessed and compared to a previously validated subjective categorical rating scale. This comparison revealed that grey scale measurement of the proportion of a defined area of spinal cord occupied by OX42 immunoreactive cells is a robust image analysis paradigm. (2) Develop and validate a flow cytometric approach for quantification of spinal microgliosis. The flow cytometric technique reliably quantified microgliosis in spinal cord cell suspensions, using OX42 and ED9 immunoreactivity to identify microglia.

The results suggest that image analysis of immunohistochemical revelation of microgliosis reliably detects the spinal microgliosis in response to peripheral nerve injury and pharmacological attenuation thereof. In addition, flow cytometry provides an alternative approach for quantitative analysis of spinal microgliosis elicited by nerve injury.

A Transition Path Sampling Study of the Reaction Catalyzed by the Enzyme Chorismate Mutase

The study of the chemical steps in enzyme-catalyzed reactions represents a challenge for molecular simulation techniques. One concern is how to calculate paths for the reaction. Common techniques include the definition of a reaction coordinate in terms of a small set of (normally) geometrical variables or the determination of minimum energy paths on the potential energy surface of the reacting system. Both have disadvantages, the former because it presupposes knowledge of which variables are likely to be important for reaction and the latter because it provides a static picture and dynamical effects are ignored. In this paper, we employ the transition path sampling method developed by Chandler and co-workers, which overcomes some of these limitations. The reaction that we have chosen is the chorismate-mutase-catalyzed conversion of chorismate into prephenate, which has become something of a test case for simulation studies of enzyme mechanisms. We generated an ensemble of approximately 1000 independent transition paths for the reaction in the enzyme and another approximately 500 for the corresponding reaction in solution. A large variety of analyses of these paths was performed, but we have concentrated on characterizing the transition state ensemble, particularly the flexibility of its structures with respect to other ligands of the enzyme and the time evolution of various geometrical and energetic properties as the reaction proceeds. We have also devised an approximate technique for locating transition state structures along the paths.

New insights into the reaction mechanism catalyzed by the glutamate racemase enzyme: pH titration curves and classical molecular dynamics simulations

The mechanism of the reactions catalyzed by the pyridoxal-phosphate-independent amino acid racemases and epimerases faces the difficult task of deprotonating a relatively low acidicity proton, the amino acid's alpha-hydrogen, with a relatively poor base, a cysteine. In this work, we propose a mechanism for one of these enzymes, glutamate racemase (MurI), about which many controversies exist, and the roles that its active site residues may play. The titration curves and the pK1/2 values of all of the ionizable residues for different structures leading from reactants to products have been analyzed. From these results a concerted mechanism has been proposed in which the Cys70 residue would deprotonate the alpha-hydrogen of the substrate while, at the same time, being deprotonated by the Asp7 residue. To study the consistency of this mechanism classical molecular dynamics (MD) simulations have been carried out along with pK1/2 calculations on the MD-generated structures.