BioMe: biologically relevant metals

Be2+ Causes Hypersensitivity but Mg2+ and Ca2+ Do Not─Favorable Metal Coordination Is the Key for Differential Allosteric Modulation and Binding Affinities

Although the ion selectivity of metalloproteins has been well established, selective metal antigen recognition by immunoproteins remains elusive. One such case is the recognition of the Be2+ ion against its heavier congeners, Mg2+ and Ca2+, by the human leukocyte antigen immunoprotein (HLA-DP2), leading to immunotoxicity. Integrating with our previous mechanistic study on Be2+ toxicity, herein, we have explored the basis of characteristic nontoxicity of Mg2+ and Ca2+ ions despite their in vivo abundance. The ion binding cleft of the HLA-DP2-peptide complex is composed of four acidic residues, p4D and p7E from the peptide and β26E and β69E from the protein. While the tetrahedral coordination site of the smaller Be2+ ion is located deep inside the cavity, hexa- to octa-coordination sites of Mg2+ and Ca2+ ions are located closer to the protein surface. The intrinsic high coordination number of Mg2+/Ca2+ ions induces allosteric modifications on the HLA-DP2_M2 surface, which are atypical for TCR recognition. Furthermore, the lower binding energy of larger Mg2+ and Ca2+ ions with the cavity residues can be correlated to the lower charge density and reduced covalent bonding nature as compared to those of the smaller Be2+ ion. In short, weak binding of Mg2+ and Ca2+ ions and the unfavorable allosteric surface modifications are probably the major determinants for the absence of Mg2+/Ca2+ ion-mediated hypersensitivity in humans.

Molecular Insights into the Calcium Binding in Troponin C through a Molecular Dynamics Study

Calcium-binding proteins play critical roles in various biological processes such as signal transduction, cell growth, and transcription factor regulation. Ion binding and target binding of Ca²⁺-binding proteins are highly related. Therefore, understanding the ion binding mechanism will benefit the relevant inhibitor design toward the Ca²⁺-binding proteins. The EF-hand is the typical ion binding motif in Ca²⁺-binding proteins. Previous studies indicate that the ion binding affinity of the EF-hand increases with the peptide length, but this mechanism has not been fully understood. Herein, using molecular dynamics simulations, thermodynamic integration calculations, and molecular mechanics Poisson-Boltzmann surface area analysis, we systematically investigated four Ca²⁺-binding peptides containing the EF-hand loop in site III of rabbit skeletal troponin C. These four peptides have 13, 21, 26, and 34 residues. Our simulations reproduced the observed trend that the ion binding affinity increases with the peptide length. Our results implied that the E-helix motif preceding the EF-hand loop, likely the Phe99 residue in particular, plays a significant role in this regulation. The E-helix has a significant impact on the backbone and side-chain conformations of the Asp103 residue, rigidifying important hydrogen bonds in the EF-hand and decreasing the solvent exposure of the Ca²⁺ ion, hence leading to more favorable Ca²⁺ binding in longer peptides. The present study provides molecular insights into the ion binding in the EF-hand and establishes an important step toward elucidating the responses of Ca²⁺-binding proteins toward the ion and target availability.

Impacts of targeting different hydration free energy references on the development of ion potentials

Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol^-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H⁺). The choice of ΔG°(H⁺) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K⁺, Mg²⁺, Al³⁺, and Cl^- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H⁺) for the cation and an increase in ΔG°(H⁺) for the anion. The increase in ΔG°(H⁺) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H⁺) value of -259.8 kcal mol^-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H⁺) for Mg²⁺ model is better for a reasonable description of the metal binding network.

MetBP: A Software Tool for Detection of Interaction between Metal Ion-RNA Base Pairs

MOTIVATION

The role of metals in shaping and functioning of RNA is a well established fact and the understanding of that through the analysis of structural data has biological relevance. Often metal ions bind to one or more atoms of the nucleobase of an RNA. This fact becomes more interesting when such bases form a base pair with any other base. Furthermore, when metal ions bind to any residue of an RNA, the secondary structural features of the residue (helix, loop, unpaired etc) are also biologically important. The available metal binding related software tools cannot address such type specific queries.

RESULTS

To fill this limitation, we have designed a software tool, called MetBP, that meets the goal. This tool is a stand-alone command line based tool and has no dependency on the other existing software. It accepts a structure file in mmCIF or PDB format and computes the base pairs and thereafter reports all metals that bind to one or more nucleotides that form pairs with another. It reports binding distance, angles along with base pair stability. It also reports several other important aspects, e.g. secondary structure of the residue in the RNA. MetBP can be used as a generalized metal binding site detection tool for Proteins and DNA as well.

AVAILABILITY

https://github.com/computational-biology/metbp.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Polarizable Cationic Dummy Metal Ion Model

A novel locally polarizable multisite model based on the original cation dummy atom (CDA) model is described for molecular dynamics simulations of ions in condensed phases. Polarization effects are introduced by the electronegativity equalization model (EEM) method where charges on the metal ion and its dummy atoms can fluctuate to respond to the environment. This model includes explicit polarization and ion-induced interactions and can be coupled with nonpolarizable or polarizable water models, making it more transferable to simpler force fields. This approach allows us to enhance the original fixed charge CDA model where the charge distribution cannot adapt to the local solvent structure. To illustrate the new CDApol model, we examined properties of the Zn2+, Al3+, and Zr4+ ions in aqueous solution. The polarizable model and Lennard-Jones parameters were refined for octahedrally coordinated Zn2+, Al3+, and Zr4+ CDAs to reproduce thermodynamic and geometrical properties. Using this locally polarizable model, we were able to obtain the experimental hydration free energy, ion-oxygen distance, and coordination number coupled with the standard 12-6 Lennard-Jones model. This model can be used in myriad additional applications where local polarization and charge transfer is important.

A Comprehensive Review of Computation-Based Metal-Binding Prediction Approaches at the Residue Level

Clear evidence has shown that metal ions strongly connect and delicately tune the dynamic homeostasis in living bodies. They have been proved to be associated with protein structure, stability, regulation, and function. Even small changes in the concentration of metal ions can shift their effects from natural beneficial functions to harmful. This leads to degenerative diseases, malignant tumors, and cancers. Accurate characterizations and predictions of metalloproteins at the residue level promise informative clues to the investigation of intrinsic mechanisms of protein-metal ion interactions. Compared to biophysical or biochemical wet-lab technologies, computational methods provide open web interfaces of high-resolution databases and high-throughput predictors for efficient investigation of metal-binding residues. This review surveys and details 18 public databases of metal-protein binding. We collect a comprehensive set of 44 computation-based methods and classify them into four categories, namely, learning-, docking-, template-, and meta-based methods. We analyze the benchmark datasets, assessment criteria, feature construction, and algorithms. We also compare several methods on two benchmark testing datasets and include a discussion about currently publicly available predictive tools. Finally, we summarize the challenges and underlying limitations of the current studies and propose several prospective directions concerning the future development of the related databases and methods.

Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12-6 Potential

Ions are of central importance in nature, and a variety of potential models was proposed to model ions in different phases for an in-depth exploration of ion-related systems. Here, we developed point charge models of 14 monovalent ions with the traditional 12-6 Lennard-Jones (LJ) potential for use in conjunction with 11 water models of TIP3P, OPC3, SPC/E, SPC/E_b, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models reproduced the real hydration free energy (HFE) of ions and the ion-oxygen distance (IOD) in the first hydration shell accurately and simultaneously, a performance similar to the previously reported 12-6-4 LJ-type ion models (12-6 LJ plus an attractive C₄ term for cations or a repulsive one for anions). This work, along with our previous work on di-, tri-, and tetravalent metal cations (J. Chem. Inf. Model. 2021, 61, 4031-4044; J. Chem. Inf. Model. 2021, 61, 4613-4629), demonstrates the feasibility of the simple 12-6 LJ potential in ion modeling. In order for the 12-6 LJ potential to reproduce both the HFE and IOD, the LJ R parameters need to be close to Shannon's ionic radii for the highly charged cations and to the Stokes's van der Waals (vdW) radii for the monovalent ions. With an additional C₄ term, the R parameters of 12-6-4 LJ ion models agree well with the Stokes's vdW radii and have a more physical meaning. It appears that the C₄ term can be merged into the 12-6 LJ potential by a rational tuning of R and the LJ well depth. Simulations of the osmotic coefficients of alkali chloride solutions and the properties of gaseous and solid alkali halides indicate the necessity of further optimizing ion-ion interactions via, for instance, targeting more properties or using a more physical (polarizable) model.

Determination of van der Waals Parameters Using a Double Exponential Potential for Nonbonded Divalent Metal Cations in TIP3P Solvent

A double exponential (DE) functional form for Lennard-Jones (LJ) interactions, proposed in our previous study, has many advantages over LJ potentials including a natural softcore characteristic for the convenience of the pathway-based free-energy calculations, fast convergence, and flexibility in use. In this work, we put the first step on the application of the DE functional form by identifying a DE potential, coined DE-TIP3P, for molecular simulations using the TIP3P water model. The developed DE-TIP3 potential was better than LJ potential in reproducing the experimental water properties. Afterward, we developed the nonbonded models of 15 divalent metal ions, which frequently appear and play vital roles in biological systems, to be consistent with the DE-TIP3P potential and TIP3P water model. Our nonbonded models were as good as the complicated nonbonded dummy cationic models by Jiang et al. and the nonbonded 12-6-4 LJ models by Li and Merz in reproducing the experimental properties of those ions. Moreover, our nonbonded models achieved a better performance than the compromise (CM) LJ models and 12-6-4 LJ models, developed by Li and Merz, in reproducing the properties of MgCl₂ in aqueous solution.

Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Monovalent ions play significant roles in various biological and material systems. Recently, four new water models (OPC3, OPC, TIP3P-FB, and TIP4P-FB), with significantly improved descriptions of condensed phase water, have been developed. The pairwise interaction between the metal ion and water necessitates the development of ion parameters specifically for these water models. Herein, we parameterized the 12-6 and the 12-6-4 nonbonded models for 12 monovalent ions with the respective four new water models. These monovalent ions contain eight cations including alkali metal ions (Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺), transition-metal ions (Cu⁺ and Ag⁺), and Tl⁺ from the boron family, along with four halide anions (F^-, Cl^-, Br^-, I^-). Our parameters were designed to reproduce the target hydration free energies (the 12-6 hydration free energy (HFE) set), the ion-oxygen distances (the 12-6 ion-oxygen distance (IOD) set), or both of them (the 12-6-4 set). The 12-6-4 parameter set provides highly accurate structural features overcoming the limitations of the routinely used 12-6 nonbonded model for ions. Specifically, we note that the 12-6-4 parameter set is able to reproduce experimental hydration free energies within 1 kcal/mol and experimental ion-oxygen distances within 0.01 Å simultaneously. We further reproduced the experimentally determined activity derivatives for salt solutions, validating the ion parameters for simulations of ion pairs. The improved performance of the present water models over our previous parameter sets for the TIP3P, TIP4P, and SPC/E water models (Li, P. et al J. Chem. Theory Comput. 2015 11 1645 1657) highlights the importance of the choice of water model in conjunction with the metal ion parameter set.

IMPORTANCE OF ZINC CONCENTRATION IN SEMINAL FLUID OF MEN DIAGNOSED WITH INFERTILITY

The prevalence and importance of male infertility make it a serious worldwide medical and social problem. The aim of this study was to determine if there were any significant differences in zinc concentrations in seminal fluid in men diagnosed with infertility compared to control group, and if there were, to determine how it affected the number and percentage of sperm cells with normal motility and morphology. A case-control study was conducted in the Osijek-Baranja County from January 2014 to June of 2015. The study included 276 subjects referred to biochemistry laboratory of the Osijek University Hospital Centre, 100 (36.2%) of them with low sperm count confirmed by spermiogram after 3 months. Semen samples were processed according to the World Health Organization criteria, and zinc concentrations were determined using spectrophotometry and direct colorimetry without deproteinization on a chemistry analyzer (Olympus AU 680, Beckman Coulter, Tokyo, Japan). Study results showed the men with low sperm count to be significantly older (Mann Whitney U test: p=0.013) and to have lower zinc levels (χ²-test: p<0.001). Further analysis included 100 (36.2%) patients with low sperm count. In the group of infertile men, those with low sperm count had higher zinc levels compared to men with azoospermia (Mann Whitney U test: p=0.036), suggesting a connection between lower zinc level in seminal fluid and male fertility. Zinc as a biological marker and an antioxidant affects sperm count, motility and morphology. Zinc supplemental therapy could improve seminal parameters in patients diagnosed with low sperm count. Dietary habits of the subjects differ regionally, and future research could make a solid foundation for complementary approach to treatment of male infertility.

Structure of the Human ACP-ISD11 Heterodimer

In recent years, the mammalian mitochondrial protein complex for iron-sulfur cluster assembly has been the focus of important studies. This is partly because of its high degree of relevance in cell metabolism and because mutations of the involved proteins are the cause of several human diseases. Cysteine desulfurase NFS1 is the key enzyme of the complex. At present, it is well-known that the active form of NFS1 is stabilized by the small protein ISD11. In this work, the structure of the human mitochondrial ACP-ISD11 heterodimer was determined at 2.0 Å resolution. ACP-ISD11 forms a cooperative unit stabilized by several ionic interactions, hydrogen bonds, and apolar interactions. The 4'-phosphopantetheine-acyl chain, which is covalently bound to ACP, interacts with several residues of ISD11, modulating together with ACP the foldability of ISD11. Recombinant human ACP-ISD11 was able to interact with the NFS1 desulfurase, thus yielding an active enzyme, and the NFS1/ACP-ISD11 core complex was activated by frataxin and ISCU proteins. Internal motions of ACP-ISD11 were studied by molecular dynamics simulations, showing the persistence of the interactions between both protein chains. The conformation of the dimer is similar to that found in the context of the (NFS1/ACP-ISD11)₂ supercomplex core, which contains the Escherichia coli ACP instead of the human variant. This fact suggests a sequential mechanism for supercomplex consolidation, in which the ACP-ISD11 complex may fold independently and, after that, the NFS1 dimer would be stabilized.

Crystal structure of Yersinia pestis virulence factor YfeA reveals two polyspecific metal-binding sites

Gram-negative bacteria use siderophores, outer membrane receptors, inner membrane transporters and substrate-binding proteins (SBPs) to transport transition metals through the periplasm. The SBPs share a similar protein fold that has undergone significant structural evolution to communicate with a variety of differentially regulated transporters in the cell. In Yersinia pestis, the causative agent of plague, YfeA (YPO2439, y1897), an SBP, is important for full virulence during mammalian infection. To better understand the role of YfeA in infection, crystal structures were determined under several environmental conditions with respect to transition-metal levels. Energy-dispersive X-ray spectroscopy and anomalous X-ray scattering data show that YfeA is polyspecific and can alter its substrate specificity. In minimal-media experiments, YfeA crystals grown after iron supplementation showed a threefold increase in iron fluorescence emission over the iron fluorescence emission from YfeA crystals grown from nutrient-rich conditions, and YfeA crystals grown after manganese supplementation during overexpression showed a fivefold increase in manganese fluorescence emission over the manganese fluorescence emission from YfeA crystals grown from nutrient-rich conditions. In all experiments, the YfeA crystals produced the strongest fluorescence emission from zinc and could not be manipulated otherwise. Additionally, this report documents the discovery of a novel surface metal-binding site that prefers to chelate zinc but can also bind manganese. Flexibility across YfeA crystal forms in three loops and a helix near the buried metal-binding site suggest that a structural rearrangement is required for metal loading and unloading.

New findings about human dipeptidyl peptidase III based on mutations found in cancer

In this work we investigated the role of two highly conserved residues in the peptidase family M49, whose mutations G313W and R510W were detected in human cancer, using combined experimental and computational approaches.

Minimal Functional Sites in Metalloproteins and Their Usage in Structural Bioinformatics

Metal ions play a functional role in numerous biochemical processes and cellular pathways. Indeed, about 40% of all enzymes of known 3D structure require a metal ion to be able to perform catalysis. The interactions of the metals with the macromolecular framework determine their chemical properties and reactivity. The relevant interactions involve both the coordination sphere of the metal ion and the more distant interactions of the so-called second sphere, i.e., the non-bonded interactions between the macromolecule and the residues coordinating the metal (metal ligands). The metal ligands and the residues in their close spatial proximity define what we call a minimal functional site (MFS). MFSs can be automatically extracted from the 3D structures of metal-binding biological macromolecules deposited in the Protein Data Bank (PDB). They are 3D templates that describe the local environment around a metal ion or metal cofactor and do not depend on the overall macromolecular structure. MFSs provide a different view on metal-binding proteins and nucleic acids, completely focused on the metal. Here we present different protocols and tools based upon the concept of MFS to obtain deeper insight into the structural and functional properties of metal-binding macromolecules. We also show that structure conservation of MFSs in metalloproteins relates to local sequence similarity more strongly than to overall protein similarity.

Sodium and Potassium Interactions with Nucleic Acids

Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.

Fusion RNAs in crystallographic studies of double-stranded RNA from trypanosome RNA editing

Head-to-head fusions of two identical double-stranded fragments of RNA can be designed to self-assemble from a single RNA species and form a double-stranded helix with a twofold rotation axis relating the two strands. These symmetrical RNA molecules are more likely to crystallize without end-on-end statistical packing disorder because the two halves of the molecule are identical. This approach can be used to study many fragments of double-stranded RNA or many isolated helical domains from large single-stranded RNAs that may not yet be amenable to high-resolution studies by crystallography or NMR. We used fusion RNAs to study one (the U-helix) of three functional domains formed when guide RNA binds to its cognate pre-edited mRNA from the trypanosome RNA editing system. The U-helix forms when the 3' oligo(U) tail of the guide RNA (gRNA) binds to the purine-rich, pre-edited mRNA upstream from the current RNA editing site. Fusion RNAs 16-and 32-base pairs in length formed crystals that gave diffraction to 1.37 and 1.05 Å respectively. We provide the composition of a fusion RNA crystallization screen and describe the X-ray data collection, structure determination, and refinement of the crystal structures of fusion RNAs.

Systematic Parameterization of Monovalent Ions Employing the Nonbonded Model

Monovalent ions play fundamental roles in many biological processes in organisms. Modeling these ions in molecular simulations continues to be a challenging problem. The 12-6 Lennard-Jones (LJ) nonbonded model is widely used to model monovalent ions in classical molecular dynamics simulations. A lot of parameterization efforts have been reported for these ions with a number of experimental end points. However, some reported parameter sets do not have a good balance between the two Lennard-Jones parameters (the van der Waals (VDW) radius and potential well depth), which affects their transferability. In the present work, via the use of a noble gas curve we fitted in former work (J. Chem. Theory Comput. 2013, 9, 2733), we reoptimized the 12-6 LJ parameters for 15 monovalent ions (11 positive and 4 negative ions) for three extensively used water models (TIP3P, SPC/E, and TIP4P(EW)). Since the 12-6 LJ nonbonded model performs poorly in some instances for these ions, we have also parameterized the 12-6-4 LJ-type nonbonded model (J. Chem. Theory Comput. 2014, 10, 289) using the same three water models. The three derived parameter sets focused on reproducing the hydration free energies (the HFE set) and the ion-oxygen distance (the IOD set) using the 12-6 LJ nonbonded model and the 12-6-4 LJ-type nonbonded model (the 12-6-4 set) overall give improved results. In particular, the final parameter sets showed better agreement with quantum mechanically calculated VDW radii and improved transferability to ion-pair solutions when compared to previous parameter sets.

Studying allosteric regulation in metal sensor proteins using computational methods

In this chapter, we describe advances made in understanding the mechanism of allosteric regulation of DNA operator binding in the ArsR/SmtB family of metal-sensing proteins using computational methods. The paradigm, zinc-sensing transcriptional repressor Staphylococcus aureus CzrA represents an excellent model system to understand how metal sensor proteins maintain cellular metal homeostasis. Here, we discuss studies that helped to characterize a metal ion-mediated hydrogen-bonding pathway (HBP) that plays a dominant role in the allosteric mechanism of DNA operator binding in these proteins. The chapter discusses computational methods used to provide a molecular basis for the large conformational motions and allosteric coupling free energy (~6kcal/mol) associated with Zn(II) binding in CzrA. We present an accurate and convenient means by which to include metal ions in the nuclear magnetic resonance (NMR) structure determination process using molecular dynamics (MD) constrained by NMR-derived data. The method provides a realistic and physically viable description of the metal-binding site(s) and has potentially broad applicability in the structure determination of metal ion-bound proteins, protein folding, and metal template protein-design studies. Finally, our simulations provide strong support for a proposed HBP that physically connects the metal-binding residue, His97, to the DNA-binding interface through the αR helix that is present only in the Zn(II)-bound state. We find the interprotomer hydrogen bond interaction to be significantly stronger (~8kcal/mol) at functional allosteric metal-binding sites compared to the apo proteins. This interaction works to overcome the considerable disorder at these hydrogen-bonding sites in apo protein and functions as a "switch" to lock in a weak DNA-binding conformation once metal is bound. This interaction is found to be considerably weaker in nonresponsive metal-binding sites. These findings suggest a conserved functional role of metal-mediated second-shell coordination hydrogen bonds at allosterically responsive sites in zinc-sensing transcription regulators.

Analyses of cobalt-ligand and potassium-ligand bond lengths in metalloproteins: trends and patterns

Cobalt and potassium are biologically important metal elements that are present in a large array of proteins. Cobalt is mostly found in vivo associated with a corrin ring, which represents the core of the vitamin B12 molecule. Potassium is the most abundant metal in the cytosol, and it plays a crucial role in maintaining membrane potential as well as correct protein function. Here, we report a thorough analysis of the geometric properties of cobalt and potassium coordination spheres that was performed with high resolution on a representative set of structures from the Protein Data Bank and complemented by quantum mechanical calculations realized at the DFT level of theory (B3LYP/ SDD) on mononuclear model systems. The results allowed us to draw interesting conclusions on the structural characteristics of both Co and K centers, and to evaluate the importance of effects such as their association energies and intrinsic thermodynamic stabilities. Overall, the results obtained provide useful data for enhancing the atomic models normally applied in theoretical and computational studies of Co or K proteins performed at the quantum mechanical level, and for developing molecular mechanical parameters for treating Co or K coordination spheres in molecular mechanics or molecular dynamics studies.

Taking into Account the Ion-induced Dipole Interaction in the Nonbonded Model of Ions

Metal ions exist in almost half of the proteins in the protein databank and they serve as structural, electron-transfer and catalytic elements in the metabolic processes of organisms. Molecular Dynamics (MD) simulation is a powerful tool that provides information about biomolecular systems at the atomic level. Coupled with the growth in computing power, algorithms like the Particle Mesh Ewald (PME) method have become the accepted standard when dealing with long-range interactions in MD simulations. The nonbonded model of metal ions consists of an electrostatic plus 12-6 Lennard Jones (LJ) potential and is used largely because of its speed relative to more accurate models. In previous work we found that ideal parameters do not exist that reproduce several experimental properties for M(II) ions simultaneously using the nonbonded model coupled with the PME method due to the underestimation of metal ion-ligand interactions. Via a consideration of the nature of the nonbonded model, we proposed that the observed error largely arises from overlooking charge-induced dipole interactions. The electrostatic plus 12-6 LJ potential model works reasonably well for neutral systems but does struggle with more highly charged systems. In the present work we designed and parameterized a new nonbonded model for metal ions by adding a 1/r4 term to the 12-6 model. We call it the 12-6-4 LJ-type nonbonded model due to its mathematical construction. Parameters were determined for 16 +2 metal ions for the TIP3P, SPC/E and TIP4PEW water models. The final parameters reproduce the experimental hydration free energies (HFE), ion-oxygen distances (IOD) in the first solvation shell and coordination numbers (CN) accurately for the metal ions investigated. Preliminary tests on MgCl2 at different concentrations in aqueous solution and Mg2+--nucleic acid systems show reasonable results suggesting that the present parameters can work in mixed systems. The 12-6-4 LJ-type nonbonded model is readily adopted into standard force fields like AMBER, CHARMM and OPLS-AA with only a modest computational overhead. The new nonbonded model doesn't consider charge-transfer effects explicitly and, hence, may not suitable for the simulation of systems where charge-transfer effects play a decisive role.

Rational Design of Particle Mesh Ewald Compatible Lennard-Jones Parameters for +2 Metal Cations in Explicit Solvent

Metal ions play significant roles in biological systems. Accurate molecular dynamics (MD) simulations on these systems require a validated set of parameters. Although there are more detailed ways to model metal ions, the nonbonded model, which employs a 12-6 Lennard-Jones (LJ) term plus an electrostatic potential is still widely used in MD simulations today due to its simple form. However, LJ parameters have limited transferability due to different combining rules, various water models and diverse simulation methods. Recently, simulations employing a Particle Mesh Ewald (PME) treatment for long-range electrostatics have become more and more popular owing to their speed and accuracy. In the present work we have systematically designed LJ parameters for 24 +2 metal (M(II)) cations to reproduce different experimental properties appropriate for the Lorentz-Berthelot combining rules and PME simulations. We began by testing the transferability of currently available M(II) ion LJ parameters. The results showed that there are differences between simulations employing Ewald summation with other simulation methods and that it was necessary to design new parameters specific for PME based simulations. Employing the thermodynamic integration (TI) method and performing periodic boundary MD simulations employing PME, allowed for the systematic investigation of the LJ parameter space. Hydration free energies (HFEs), the ion-oxygen distance in the first solvation shell (IOD) and coordination numbers (CNs) were obtained for various combinations of the parameters of the LJ potential for four widely used water models (TIP3P, SPC/E, TIP4P and TIP4PEW). Results showed that the three simulated properties were highly correlated. Meanwhile, M(II) ions with the same parameters in different water models produce remarkably different HFEs but similar structural properties. It is difficult to reproduce various experimental values simultaneously because the nonbonded model underestimates the interaction between the metal ions and water molecules at short range. Moreover, the extent of underestimation increases successively for the TIP3P, SPC/E, TIP4PEW and TIP4P water models. Nonetheless, we fitted a curve to describe the relationship between ε (the well depth) and radius (Rmin/2) from experimental data on noble gases to facilitate the generation of the best possible compromise models. Hence, by targeting different experimental values, we developed three sets of parameters for M(II) cations for three different water models (TIP3P, SPC/E and TIP4PEW). These parameters we feel represent the best possible compromise that can be achieved using the nonbonded model for the ions in combination with simple water models. From a computational uncertainty analysis we estimate that the uncertainty in our computed HFEs is on the order of ±1kcal/mol. Further improvements will require more advanced non-bonded models likely with inclusion of polarization.

Solution NMR refinement of a metal ion bound protein using metal ion inclusive restrained molecular dynamics methods

Correctly calculating the structure of metal coordination sites in a protein during the process of nuclear magnetic resonance (NMR) structure determination and refinement continues to be a challenging task. In this study, we present an accurate and convenient means by which to include metal ions in the NMR structure determination process using molecular dynamics (MD) simulations constrained by NMR-derived data to obtain a realistic and physically viable description of the metal binding site(s). This method provides the framework to accurately portray the metal ions and its binding residues in a pseudo-bond or dummy-cation like approach, and is validated by quantum mechanical/molecular mechanical (QM/MM) MD calculations constrained by NMR-derived data. To illustrate this approach, we refine the zinc coordination complex structure of the zinc sensing transcriptional repressor protein Staphylococcus aureus CzrA, generating over 130 ns of MD and QM/MM MD NMR-data compliant sampling. In addition to refining the first coordination shell structure of the Zn(II) ion, this protocol benefits from being performed in a periodically replicated solvation environment including long-range electrostatics. We determine that unrestrained (not based on NMR data) MD simulations correlated to the NMR data in a time-averaged ensemble. The accurate solution structure ensemble of the metal-bound protein accurately describes the role of conformational sampling in allosteric regulation of DNA binding by zinc and serves to validate our previous unrestrained MD simulations of CzrA. This methodology has potentially broad applicability in the structure determination of metal ion bound proteins, protein folding and metal template protein-design studies.

Cell-free production of integral membrane aspartic acid proteases reveals zinc-dependent methyltransferase activity of the Pseudomonas aeruginosa prepilin peptidase PilD

Integral membrane aspartic acid proteases are receiving growing recognition for their fundamental roles in cellular physiology of eukaryotes and prokaryotes, and may be medically important pharmaceutical targets. The Gram-negative Pseudomonas aeruginosa PilD and the archaeal Methanococcus voltae FlaK were synthesized in the presence of unilamellar liposomes in a cell-free translation system. Cosynthesis of PilD with its full-length substrate, PilA, or of FlaK with its full-length substrate, FlaB2, led to complete cleavage of the substrate signal peptides. Scaled-up synthesis of PilD, followed by solubilization in dodecyl-β-d-maltoside and chromatography, led to a pure enzyme that retained both of its known biochemical activities: cleavage of the PilA signal peptide and S-adenosyl methionine-dependent methylation of the mature pilin. X-ray fluorescence scans show for the first time that PilD is a zinc-binding protein. Zinc is required for the N-terminal methylation of the mature pilin, but not for signal peptide cleavage. Taken together, our work identifies the P. aeruginosa prepilin peptidase PilD as a zinc-dependent N-methyltransferase and provides a new platform for large-scale synthesis of PilD and other integral membrane proteases important for basic microbial physiology and virulence.