Structure and dynamics of the N-terminal domain of the Cu(I) binding protein CusB

Mechanisms of Copper Selectivity and Release by the Metallochaperone CusF: Insights from CO-Binding, Rapid-Freeze-Quench EXAFS, and Unnatural Amino Acid Substitution

Metallochaperones are small proteins that shuttle essential metal ions such as Cu selectively to their cellular targets. CusF has unusual Cu(I) coordination, bound by two methionines, one histidine and a capping tryptophan residue, W44. Here we compare the CO binding reactivity of the wild type (WT) protein and its W44A, F, and M variants. Fourier-transform infrared (FTIR) indicates that W44A provides unhindered access for CO, while W44M is unreactive. WT is also largely unreactive to CO suggesting that the tryptophan cap is effective in shielding the Cu(I) center from exogenous adduct formation, while the Phe variant shows partial reactivity suggestive of an equilibrium between cap-on and cap-off conformers. Rates of metal transfer to the partner CusB are consistent with the π-cation cap providing both selectivity and redox protection. Unnatural amino acid substitutions of the W44 ligand with cyano-Phe and Br-Phe underpin the conclusion that the Phe ligand is a less effective capping residue. Finally, density functional theory (DFT) calculations validate the CO-binding strategy. Overall, the study suggests that CusF uses the tryptophan cap to protect against exogenous ligand (O2) attack while the mechanism of protein-protein complex formation allows the cap to swing out of the way, and thus have minimal effect on the rates of metal transfer.

Structure, Assembly, and Function of Tripartite Efflux and Type 1 Secretion Systems in Gram-Negative Bacteria

Tripartite efflux pumps and the related type 1 secretion systems (T1SSs) in Gram-negative organisms are diverse in function, energization, and structural organization. They form continuous conduits spanning both the inner and the outer membrane and are composed of three principal components-the energized inner membrane transporters (belonging to ABC, RND, and MFS families), the outer membrane factor channel-like proteins, and linking the two, the periplasmic adaptor proteins (PAPs), also known as the membrane fusion proteins (MFPs). In this review we summarize the recent advances in understanding of structural biology, function, and regulation of these systems, highlighting the previously undescribed role of PAPs in providing a common architectural scaffold across diverse families of transporters. Despite being built from a limited number of basic structural domains, these complexes present a staggering variety of architectures. While key insights have been derived from the RND transporter systems, a closer inspection of the operation and structural organization of different tripartite systems reveals unexpected analogies between them, including those formed around MFS- and ATP-driven transporters, suggesting that they operate around basic common principles. Based on that we are proposing a new integrated model of PAP-mediated communication within the conformational cycling of tripartite systems, which could be expanded to other types of assemblies.

Parametrization of Trivalent and Tetravalent Metal Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Commonly seen in rare-earth chemistry and materials science, highly charged metal ions play key roles in many chemical processes. Computer simulations have become an important tool for scientific research nowadays. Meaningful simulations require reliable parameters. In the present work, we parametrized 18 M(III) and 6 M(IV) metal ions for four new water models (OPC3, OPC, TIP3P-FB, TIP4P-FB) in conjunction with each of the 12-6 and 12-6-4 nonbonded models. Similar to what was observed previously, issues with the 12-6 model can be fixed by using the 12-6-4 model. Moreover, the four new water models showed comparable performance or considerable improvement over the previous water models (TIP3P, SPC/E, and TIP4P_EW) in the same category (3-point or 4-point water models, respectively). Finally, we reported a study of a metalloprotein system demonstrating the capability of the 12-6-4 model to model metalloproteins. The reported parameters will facilitate accurate simulations of highly charged metal ions in aqueous solution.

Molecular mechanism of amyloidogenicity and neurotoxicity of a pro-aggregated tau mutant in the presence of histidine tautomerism via replica-exchange simulation

Considering ΔK280 tau mutation, δε isomer with highest sheet content may accelerate aggregation; generating small compounds to inhibit this would help tp prevent tauopathies.

Cu Homeostasis in Bacteria: The Ins and Outs

Copper (Cu) is an essential trace element for all living organisms and used as cofactor in key enzymes of important biological processes, such as aerobic respiration or superoxide dismutation. However, due to its toxicity, cells have developed elaborate mechanisms for Cu homeostasis, which balance Cu supply for cuproprotein biogenesis with the need to remove excess Cu. This review summarizes our current knowledge on bacterial Cu homeostasis with a focus on Gram-negative bacteria and describes the multiple strategies that bacteria use for uptake, storage and export of Cu. We furthermore describe general mechanistic principles that aid the bacterial response to toxic Cu concentrations and illustrate dedicated Cu relay systems that facilitate Cu delivery for cuproenzyme biogenesis. Progress in understanding how bacteria avoid Cu poisoning while maintaining a certain Cu quota for cell proliferation is of particular importance for microbial pathogens because Cu is utilized by the host immune system for attenuating pathogen survival in host cells.

Inhibiting the copper efflux system in microbes as a novel approach for developing antibiotics

Five out of six people receive at least one antibiotic prescription per year. However, the ever-expanding use of antibiotics in medicine, agriculture, and food production has accelerated the evolution of antibiotic-resistant bacteria, which, in turn, made the development of novel antibiotics based on new molecular targets a priority in medicinal chemistry. One way of possibly combatting resistant bacterial infections is by inhibiting the copper transporters in prokaryotic cells. Copper is a key element within all living cells, but it can be toxic in excess. Both eukaryotic and prokaryotic cells have developed distinct copper regulation systems to prevent its toxicity. Therefore, selectively targeting the prokaryotic copper regulation system might be an initial step in developing next-generation antibiotics. One such system is the Gram-negative bacterial CusCFBA efflux system. CusB is a key protein in this system and was previously reported to play an important role in opening the channel for efflux via significant structural changes upon copper binding while also controlling the assembly and disassembly process of the entire channel. In this study, we aimed to develop novel peptide copper channel blockers, designed by in silico calculations based on the structure of CusB. Using a combination of magnetic resonance spectroscopy and various biochemical methods, we found a lead peptide that promotes copper-induced cell toxicity. Targeting copper transport in bacteria has not yet been pursued as an antibiotic mechanism of action. Thus, our study lays the foundation for discovering novel antibiotics.

Catalytic M Center of Copper Monooxygenases Probed by Rational Design. Effects of Selenomethionine and Histidine Substitution on Structure and Reactivity

The M centers of the mononuclear monooxygenases peptidylglycine monooxygenase (PHM) and dopamine β-monooxygenase bind and activate dioxygen en route to substrate hydroxylation. Recently, we reported the rational design of a protein-based model in which the CusF metallochaperone was repurposed via a His to Met mutation to act as a structural and spectroscopic biomimic. The PHM M site exhibits a number of unusual attributes, including a His₂Met ligand set, a fluxional Cu(I)-S(Met) bond, tight binding of exogenous ligands CO and N₃^-, and complete coupling of oxygen reduction to substrate hydroxylation even at extremely low turnover rates. In particular, mutation of the Met ligand to His completely eliminates the catalytic activity despite the propensity of Cu^I-His₃ centers to bind and activate dioxygen in other metalloenzyme systems. Here, we further develop the CusF-based model to explore methionine variants in which Met is replaced by selenomethionine (SeM) and histidine. We examine the effects on coordinate structure and exogenous ligand binding via X-ray absorption spectroscopy and electron paramagnetic resonance and probe the consequences of mutations on redox chemistry via studies of the reduction by ascorbate and oxidation via molecular oxygen. The M-site model is three-coordinate in the Cu(I) state and binds CO to form a four-coordinate carbonyl. In the oxidized forms, the coordination changes to tetragonal five-coordinate with a long axial Met ligand that like the enzymes is undetectable at either the Cu or Se K edges. The EXAFS data at the Se K edge of the SeM variant provide unique information about the nature of the Cu-methionine bond that is likewise weak and fluxional. Kinetic studies document the sluggish reactivity of the Cu(I) complexes with molecular oxygen and rapid rates of reduction of the Cu(II) complexes by ascorbate, indicating a remarkable stability of the Cu(I) state in all three derivatives. The results show little difference between the Met ligand and its SeM and His congeners and suggest that the Met contributes to catalysis in ways that are more complex than simple perturbation of the redox chemistry. Overall, the results stimulate a critical re-examination of the canonical reaction mechanisms of the mononuclear copper monooxygenases.

Exploring the role of the various methionine residues in the Escherichia coli CusB adapter protein

The dissemination of resistant pathogenic microbes has become one of the most challenging problems that modern medicine has faced. Developing novel drugs based on new molecular targets that previously were not targeted, is therefore the highest priority in antibiotics research. One approach that has been recently suggested is to inhibit copper transporters in prokaryotic systems. Copper is required for many biological pathways, but sometimes it can harm the cell. Pathogenic systems have a highly sophisticated copper-regulation network; therefore, a better understanding of how this network operates at the molecular level should assist in developing the next generation of antibiotics. The CusB protein is part of the CusCBA periplasmic Cu(I) efflux system in Gram-negative bacteria, and was recently reported to play a key role in the functioning of the whole CusCBA system, in which conformational changes as well as the assembly/disassembly process control the opening of the transporter. More knowledge of the underlying mechanism is needed to attain a full understanding of CusB functioning, which is associated with targeting specific and crucial residues in CusB. Here, we combine in-vitro structural measurements, which use EPR spectroscopy and UV-Vis measurements, with cell experiments to explore the role of the various methionine residues in CusB. We targeted two methionine residues (M227 and M241) that are essential for the proper functioning of CusB.

Unraveling the Impact of Cysteine-to-Serine Mutations on the Structural and Functional Properties of Cu(I)-Binding Proteins

Appropriate maintenance of Cu(I) homeostasis is an essential requirement for proper cell function because its misregulation induces the onset of major human diseases and mortality. For this reason, several research efforts have been devoted to dissecting the inner working mechanism of Cu(I)-binding proteins and transporters. A commonly adopted strategy relies on mutations of cysteine residues, for which Cu(I) has an exquisite complementarity, to serines. Nevertheless, in spite of the similarity between these two amino acids, the structural and functional impact of serine mutations on Cu(I)-binding biomolecules remains unclear. Here, we applied various biochemical and biophysical methods, together with all-atom simulations, to investigate the effect of these mutations on the stability, structure, and aggregation propensity of Cu(I)-binding proteins, as well as their interaction with specific partner proteins. Among Cu(I)-binding biomolecules, we focused on the eukaryotic Atox1-ATP7B system, and the prokaryotic CueR metalloregulator. Our results reveal that proteins containing cysteine-to-serine mutations can still bind Cu(I) ions; however, this alters their stability and aggregation propensity. These results contribute to deciphering the critical biological principles underlying the regulatory mechanism of the in-cell Cu(I) concentration, and provide a basis for interpreting future studies that will take advantage of cysteine-to-serine mutations in Cu(I)-binding systems.

Trapping intermediates in metal transfer reactions of the CusCBAF export pump of Escherichia coli

Escherichia coli CusCBAF represents an important class of bacterial efflux pump exhibiting selectivity towards Cu(I) and Ag(I). The complex is comprised of three proteins: the CusA transmembrane pump, the CusB soluble adaptor protein, and the CusC outer-membrane pore, and additionally requires the periplasmic metallochaperone CusF. Here we used spectroscopic and kinetic tools to probe the mechanism of copper transfer between CusF and CusB using selenomethionine labeling of the metal-binding Met residues coupled to RFQ-XAS at the Se and Cu edges. The results indicate fast formation of a protein-protein complex followed by slower intra-complex metal transfer. An intermediate coordinated by ligands from each protein forms in 100 ms. Stopped-flow fluorescence of the capping CusF-W44 tryptophan that is quenched by metal transfer also supports this mechanism. The rate constants validate a process in which shared-ligand complex formation assists protein association, providing a driving force that raises the rate into the diffusion-limited regime.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Cadmium(II) inhibition of human uracil-DNA glycosylase by catalytic water supplantation

Toxic metals are known to inhibit DNA repair but the underlying mechanisms of inhibition are still not fully understood. DNA repair enzymes such as human uracil-DNA glycosylase (hUNG) perform the initial step in the base excision repair (BER) pathway. In this work, we showed that cadmium [Cd(II)], a known human carcinogen, inhibited all activity of hUNG at 100 μM. Computational analyses based on 2 μs equilibrium, 1.6 μs steered molecular dynamics (SMD), and QM/MM MD determined that Cd(II) ions entered the enzyme active site and formed close contacts with both D145 and H148, effectively replacing the catalytic water normally found in this position. Geometry refinement by density functional theory (DFT) calculations showed that Cd(II) formed a tetrahedral structure with D145, P146, H148, and one water molecule. This work for the first time reports Cd(II) inhibition of hUNG which was due to replacement of the catalytic water by binding the active site D145 and H148 residues. Comparison of the proposed metal binding site to existing structural data showed that D145:H148 followed a general metal binding motif favored by Cd(II). The identified motif offered structural insights into metal inhibition of other DNA repair enzymes and glycosylases.

Kβ Valence to Core X-ray Emission Studies of Cu(I) Binding Proteins with Mixed Methionine - Histidine Coordination. Relevance to the Reactivity of the M- and H-sites of Peptidylglycine Monooxygenase

Biological systems use copper as a redox center in many metalloproteins, where the role of the metal is to cycle between its +1 and +2 oxidation states. This chemistry requires the redox potential to be in a range that can stabilize both Cu(I) and Cu(II) states and often involves protein-derived ligand sets involving mixed histidine-methionine coordination that balance the preferences of both oxidation states. Transport proteins, on the other hand, utilize copper in the Cu(I) state and often contain sites comprised predominately of the cuprophilic residue methionine. The electronic factors that allow enzymes and transporters to balance their redox requirements are complex and are often elusive due to the dearth of spectroscopic probes of the Cu(I) state. Here we present the novel application of X-ray emission spectroscopy to copper proteins via a study of a series of mixed His-Met copper sites where the ligand set varies in a systematic way between the His3 and Met3 limits. The sites are derived from the wild-type peptidylglycine monooxygenase (PHM), two single-site variants which replicate each of its two copper sites (CuM-site and CuH-site), and the transporters CusF and CusB. Clear differences are observed in the Kβ2,5 region at the Met3 and His3 limits. CusB (Met3) has a distinct peak at 8978.4 eV with a broad shoulder at 8975.6 eV, whereas CuH (His3) has two well-resolved features: a more intense feature at 8974.8 eV and a second at 8977.2 eV. The mixed coordination sphere CusF (Met2His) and the PHM CuM variant (Met1His2) have very similar spectra consisting of two features at 8975.2 and 8977.8 eV. An analysis of DFT calculated spectra indicate that the intensity of the higher energy peak near 8978 eV is mediated by mixing of ligand-based orbitals into the Cu d(10) manifold, with S from Met providing more intensity by facilitating increased Cu p-d mixing. Furthermore, reaction of WT PHM with CO (an oxygen analogue) produced the M site CO complex, which showed a unique XES spectrum that could be computationally reproduced by including interactions between Cu(I) and the CO ligand. The study suggests that the valence-to-core (VtC) region can not only serve as a probe of ligand speciation but also offer insight into the coordination geometry, in a fashion similar to XAS pre-edges, and may be sufficiently sensitive to the coordination of exogenous ligands to be useful in the study of reaction mechanisms.

MCPB.py: A Python Based Metal Center Parameter Builder

MCPB.py, a python based metal center parameter builder, has been developed to build force fields for the simulation of metal complexes employing the bonded model approach. It has an optimized code structure, with far fewer required steps than the previous developed MCPB program. It supports various AMBER force fields and more than 80 metal ions. A series of parametrization schemes to derive force constants and charge parameters are available within the program. We give two examples (one metalloprotein example and one organometallic compound example), indicating the program's ability to build reliable force fields for different metal ion containing complexes. The original version was released with AmberTools15. It is provided via the GNU General Public License v3.0 (GNU_GPL_v3) agreement and is free to download and distribute. MCPB.py provides a bridge between quantum mechanical calculations and molecular dynamics simulation software packages thereby enabling the modeling of metal ion centers. It offers an entry into simulating metal ions in a number of situations by providing an efficient way for researchers to handle the vagaries and difficulties associated with metal ion modeling.

Models for the Metal Transfer Complex of the N-Terminal Region of CusB and CusF

The tripartite CusCFBA pump in Escherichia coli is a very effective heavy metal extrusion system specific for Cu(I) and Ag(I). The N-terminal region of the membrane fusion protein CusB (CusB-NT) is highly disordered, and hence, experimentally characterizing its structure is challenging. In a previous study, this disorder was confirmed with molecular dynamics simulations, although some key structural elements were determined. It was experimentally shown that CusB-NT is fully functional in transferring the metal from the metallochaperone CusF. In this study, we docked these two entities together and formed two representative metal coordination modes, which consist of residues from both proteins. In this way, we created two potential CusB-NT/CusF complexes that share coordination of Cu(I) and thereby represent structural models for the metal transfer process. Each model complex was simulated for 4 μs. The previously observed structural disorder in CusB-NT disappeared upon complexation with CusF. The only differences between the two models occurred in the M21-M36 loop region of CusB-NT and the open flap of CusF: we observed the model with two CusB-NT methionine residues and a CusF methionine as the metal coordination site (termed "MMM") to be more stable than the model with a CusB-NT methionine, a CusF methionine, and a CusF histidine ligating the metal (termed "MMH"). The observed stability of the MMM model was probed for an additional 2 μs, yielding a total simulation time of 6 μs. We hypothesize that both MMM and MMH configurations might take part in the metal exchange process in which the MMH configuration would appear first and would be followed by the MMM configuration. Given the experimental finding of comparable binding affinities of CusB-NT and CusF, the increased stability of the MMM configuration might be a determinant for the transfer from CusF to CusB-NT. The metal would be transferred from the more CusF-dominated metal binding environment (MMH model) to a more CusB-dominated one (MMM model) in which the coordination environment is more stable. From the MMM model, the metal ion would ultimately be coordinated by the CusB methionines only, which would complete the Cu(I) transfer process.

Architecture and roles of periplasmic adaptor proteins in tripartite efflux assemblies

Recent years have seen major advances in the structural understanding of the different components of tripartite efflux assemblies, which encompass the multidrug efflux (MDR) pumps and type I secretion systems. The majority of these investigations have focused on the role played by the inner membrane transporters and the outer membrane factor (OMF), leaving the third component of the system – the Periplasmic Adaptor Proteins (PAPs) – relatively understudied. Here we review the current state of knowledge of these versatile proteins which, far from being passive linkers between the OMF and the transporter, emerge as active architects of tripartite assemblies, and play diverse roles in the transport process. Recognition between the PAPs and OMFs is essential for pump assembly and function, and targeting this interaction may provide a novel avenue for combating multidrug resistance. With the recent advances elucidating the drug efflux and energetics of the tripartite assemblies, the understanding of the interaction between the OMFs and PAPs is the last piece remaining in the complete structure of the tripartite pump assembly puzzle.

EPR spectroscopy identifies Met and Lys residues that are essential for the interaction between the CusB N-terminal domain and metallochaperone CusF

Copper plays a key role in all living organisms by serving as a cofactor for a large variety of proteins and enzymes involved in electron transfer, oxidase and oxygenase activities, and the detoxification of oxygen radicals. Due to its toxicity, a conserved homeostasis mechanism is required. In E. coli, the CusCFBA efflux system is a copper-regulating system and is responsible for transferring Cu(I) and Ag(I) out of the periplasm domain into the extracellular domain. Two of the components of this efflux system, the CusF metallochaperone and the N-terminal domain of CusB, have been thought to play significant roles in the function of this efflux system. Resolving the metal ion transport mechanism through this efflux system is vital for understanding metal- and multidrug-resistant microorganisms. This work explores one aspect of the E. coli resistance mechanism by observing the interaction between the N-terminal domain of CusB and the CusF protein, using electron paramagnetic resonance (EPR) spectroscopy, circular dichroism (CD), and chemical cross-linking. The data summarized here show that M36 and M38 of CusB are important residues for both the Cu(I) coordination to the CusB N-terminal domain and the interaction with CusF, and K32 is essential for the interaction with CusF. In contrast, the K29 residue is less consequential for the interaction with CusF, whereas M21 is mostly important for the proper interaction with CusF.

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.

Role of substrate dynamics in protein prenylation reactions

CONSPECTUS: The role dynamics plays in proteins is of intense contemporary interest. Fundamental insights into how dynamics affects reactivity and product distributions will facilitate the design of novel catalysts that can produce high quality compounds that can be employed, for example, as fuels and life saving drugs. We have used molecular dynamics (MD) methods and combined quantum mechanical/molecular mechanical (QM/MM) methods to study a series of proteins either whose substrates are too far away from the catalytic center or whose experimentally resolved substrate binding modes cannot explain the observed product distribution. In particular, we describe studies of farnesyl transferase (FTase) where the farnesyl pyrophosphate (FPP) substrate is ∼8 Å from the zinc-bound peptide in the active site of FTase. Using MD and QM/MM studies, we explain how the FPP substrate spans the gulf between it and the active site, and we have elucidated the nature of the transition state (TS) and offered an alternate explanation of experimentally observed kinetic isotope effects (KIEs). Our second story focuses on the nature of substrate dynamics in the aromatic prenyltransferase (APTase) protein NphB and how substrate dynamics affects the observed product distribution. Through the examples chosen we show the power of MD and QM/MM methods to provide unique insights into how protein substrate dynamics affects catalytic efficiency. We also illustrate how complex these reactions are and highlight the challenges faced when attempting to design de novo catalysts. While the methods used in our previous studies provided useful insights, several clear challenges still remain. In particular, we have utilized a semiempirical QM model (self-consistent charge density functional tight binding, SCC-DFTB) in our QM/MM studies since the problems we were addressing required extensive sampling. For the problems illustrated, this approach performed admirably (we estimate for these systems an uncertainty of ∼2 kcal/mol), but it is still a semiempirical model, and studies of this type would benefit greatly from more accurate ab initio or DFT models. However, the challenge with these methods is to reach the level of sampling needed to study systems where large conformational changes happen in the many nanoseconds to microsecond time regimes. Hence, how to couple expensive and accurate QM methods with sophisticated sampling algorithms is an important future challenge especially when large-scale studies of catalyst design become of interest. The use of MD and QM/MM models to elucidate enzyme catalytic pathways and to design novel catalytic agents is in its infancy but shows tremendous promise. While this Account summarizes where we have been, we also discuss briefly future directions that improve our fundamental ability to understand enzyme catalysis.

Parameterization of highly charged metal ions using the 12-6-4 LJ-type nonbonded model in explicit water

Highly charged metal ions act as catalytic centers and structural elements in a broad range of chemical complexes. The nonbonded model for metal ions is extensively used in molecular simulations due to its simple form, computational speed, and transferability. We have proposed and parametrized a 12-6-4 LJ (Lennard-Jones)-type nonbonded model for divalent metal ions in previous work, which showed a marked improvement over the 12-6 LJ nonbonded model. In the present study, by treating the experimental hydration free energies and ion–oxygen distances of the first solvation shell as targets for our parametrization, we evaluated 12-6 LJ parameters for 18 M(III) and 6 M(IV) metal ions for three widely used water models (TIP3P, SPC/E, and TIP4P_EW). As expected, the interaction energy underestimation of the 12-6 LJ nonbonded model increases dramatically for the highly charged metal ions. We then parametrized the 12-6-4 LJ-type nonbonded model for these metal ions with the three water models. The final parameters reproduced the target values with good accuracy, which is consistent with our previous experience using this potential. Finally, tests were performed on a protein system, and the obtained results validate the transferability of these nonbonded model parameters.

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.