Energetics of zinc-mediated interactions in the allosteric pathways of metal sensor proteins

AsIII Selectively Induces a Disorder-to-Order Transition in the Metalloid Binding Region of the AfArsR Protein

Arsenic is highly toxic and a significant threat to human health, but certain bacteria have developed defense mechanisms initiated by AsIII binding to AsIII-sensing proteins of the ArsR family. The transcriptional regulator AfArsR responds to AsIII and SbIII by coordinating the metalloids with three cysteines, located in a short sequence of the same monomer chain. Here, we characterize the binding of AsIII and HgII to a model peptide encompassing this fragment of the protein via solution equilibrium and spectroscopic/spectrometric techniques (pH potentiometry, UV, CD, NMR, PAC, EXAFS, and ESI-MS) combined with DFT calculations and MD simulations. Coordination of AsIII changes the peptide structure from a random-coil to a well-defined structure of the complex. A trigonal pyramidal AsS3 binding site is formed with almost exactly the same structure as observed in the crystal structure of the native protein, implying that the peptide possesses all of the features required to mimic the AsIII recognition and response selectivity of AfArsR. Contrary to this, binding of HgII to the peptide does not lead to a well-defined structure of the peptide, and the atoms near the metal binding site are displaced and reoriented in the HgII model. Our model study suggests that structural organization of the metal site by the inducer ion is a key element in the mechanism of the metalloid-selective recognition of this protein.

Natural Polymorphisms D60E and I62V Stabilize a Closed Conformation in HIV-1 Protease in the Absence of an Inhibitor or Substrate

HIV infection remains a global health issue plagued by drug resistance and virological failure. Natural polymorphisms (NPs) contained within several African and Brazilian protease (PR) variants have been shown to induce a conformational landscape of more closed conformations compared to the sequence of subtype B prevalent in North America and Western Europe. Here we demonstrate through experimental pulsed EPR distance measurements and molecular dynamic (MD) simulations that the two common NPs D60E and I62V found within subtypes F and H can induce a closed conformation when introduced into HIV-1PR subtype B. Specifically, D60E alters the conformation in subtype B through the formation of a salt bridge with residue K43 contained within the nexus between the flap and hinge region of the HIV-1 PR fold. On the other hand, I62V modulates the packing of the hydrophobic cluster of the cantilever and fulcrum, also resulting in a more closed conformation.

Calorimetric analysis of AdcR and its interactions with zinc(II) and DNA

Zinc(II) ions play critical roles in all known life as structurally important stabilizing ions in proteins, catalytically active metals in enzymes, and signaling agents impacting physiological changes. To maintain homeostasis, the intracellular concentration of zinc(II) is strictly controlled by a family of metal-regulatory proteins in both prokaryotic and eukaryotic organisms. In S. pneumoniae, there are two proteins that share responsibility for Zn2+ homeostasis, one of them is the Adhesin Competence Repressor (AdcR) and it binds to a specific double-stranded DNA binding domain (dsDNA). AdcR has been structurally characterized containing two zinc(II) metal centers per monomeric unit. Here we report data collected from differential scanning calorimetry (DSC) experiments aimed to measure the structural stability of AdcR, the fully complimented Zn2AdcR complex, and the protein/DNA complex Zn2AdcR/dsDNA. Thermograms collected from DSC experiments yielded endothermic unfolding events for AdcR, Zn2AdcR, and Zn2AdcR/dsDNA complex at 55.6, 70.2, and 56.6 °C, respectively. A non-two state unfolding model best fits the data, giving ΔH terms associated with these thermal unfolding events of 5.1, 7.1, and 4.9 kcal/mol. These data allow for the development of a thermodynamic cycle connecting both zinc(II) and DNA binding to AdcR. Furthermore, pairing this newly reported data with known association constants for zinc(II) and DNA binding allowed for the generation of thermodynamic profiles for both zinc(II) binding to AdcR and Zn2AdcR binding to DNA, which show both are decisively entropy-driven processes.

Direct C–H amination reactions of arenes with N-hydroxyphthalimides catalyzed by cuprous bromide

An efficient Cu-catalyzed strategy for the direct C–H amination of arenes in high yields using N-hydroxyphthalimide as the amidyl radical precursor under air is reported. A possible mechanism is proposed that proceeds via a radical reaction in the presence of CuBr and triethyl phosphite.

Recent advances in computational methods for biosensor design

Biosensors are analytical tools with a great application in healthcare, food quality control, and environmental monitoring. They are of considerable interest to be designed by using cost-effective and efficient approaches. Designing biosensors with improved functionality or application in new target detection has been converted to a fast-growing field of biomedicine and biotechnology branches. Experimental efforts have led to valuable successes in the field of biosensor design; however, some deficiencies restrict their utilization for this purpose. Computational design of biosensors is introduced as a promising key to eliminate the gap. A set of reliable structure prediction of the biosensor segments, their stability, and accurate descriptors of molecular interactions are required to computationally design biosensors. In this review, we provide a comprehensive insight into the progress of computational methods to guide the design and development of biosensors, including molecular dynamics simulation, quantum mechanics calculations, molecular docking, virtual screening, and a combination of them as the hybrid methodologies. By relying on the recent advances in the computational methods, an opportunity emerged for them to be complementary or an alternative to the experimental methods in the field of biosensor design.

NMR and computational methods for molecular resolution of allosteric pathways in enzyme complexes

Allostery is a ubiquitous biological mechanism in which a distant binding site is coupled to and drastically alters the function of a catalytic site in a protein. Allostery provides a high level of spatial and temporal control of the integrity and activity of biomolecular assembles composed of proteins, nucleic acids, or small molecules. Understanding the physical forces that drive allosteric coupling is critical to harnessing this process for use in bioengineering, de novo protein design, and drug discovery. Current microscopic models of allostery highlight the importance of energetics, structural rearrangements, and conformational fluctuations, and in this review, we discuss the synergistic use of solution NMR spectroscopy and computational methods to probe these phenomena in allosteric systems, particularly protein-nucleic acid complexes. This combination of experimental and theoretical techniques facilitates an unparalleled detection of subtle changes to structural and dynamic equilibria in biomolecules with atomic resolution, and we provide a detailed discussion of specialized NMR experiments as well as the complementary methods that provide valuable insight into allosteric pathways in silico. Lastly, we highlight two case studies to demonstrate the adaptability of this approach to enzymes of varying size and mechanistic complexity.

Tuning site-specific dynamics to drive allosteric activation in a pneumococcal zinc uptake regulator

MarR (multiple antibiotic resistance repressor) family proteins are bacterial repressors that regulate transcription in response to a wide range of chemical signals. Although specific features of MarR family function have been described, the role of atomic motions in MarRs remains unexplored thus limiting insights into the evolution of allostery in this ubiquitous family of repressors. Here, we provide the first experimental evidence that internal dynamics play a crucial functional role in MarR proteins. Streptococcus pneumoniae AdcR (adhesin-competence repressor) regulates Zn^II homeostasis and Zn^II functions as an allosteric activator of DNA binding. Zn^II coordination triggers a transition from somewhat independent domains to a more compact structure. We identify residues that impact allosteric activation on the basis of Zn^II-induced perturbations of atomic motions over a wide range of timescales. These findings appear to reconcile the distinct allosteric mechanisms proposed for other MarRs and highlight the importance of conformational dynamics in biological regulation.

Functional Role of Solvent Entropy and Conformational Entropy of Metal Binding in a Dynamically Driven Allosteric System

Allostery is a regulatory phenomenon whereby ligand binding to one site influences the binding of the same or a different ligand to another site on a macromolecule. The physical origins of allosteric regulation remain under intense investigation. In general terms, ligand-induced structural changes, perturbations of residue-specific dynamics, and surrounding solvent molecules all potentially contribute to the global energetics of allostery. While the role of solvent is generally well understood in regulatory events associated with major protein structural rearrangements, the degree to which protein dynamics impact solvent degrees of freedom is unclear, particularly in cases of dynamically driven allostery. With the aid of new crystal structures, extensive calorimetric and residue-specific dynamics studies over a range of time scales and temperatures, we dissect for the first time the relative degree to which changes in solvent entropy and residue-specific dynamics impact dynamically driven, allosteric inhibition of DNA binding by Zn in the zinc efflux repressor, CzrA (chromosomal zinc-regulated repressor). We show that non-native residue-specific dynamics in allosterically impaired CzrA mutants are accompanied by significant perturbations in solvent entropy that cannot be predicted from crystal structures. We conclude that functional dynamics are not necessarily restricted to protein residues but involve surface water molecules that may be responding to ligand (Zn)-mediated perturbations in protein internal motions that define the conformational ensemble, rather than major structural rearrangements.

Negative regulation of daptomycin production by DepR2, an ArsR-family transcriptional factor

Daptomycin, a lipopeptide antibiotic potently active against Gram-positive bacterial pathogens, is produced by Streptomyces roseosporus, but the transcriptional regulation on its biosynthesis is not fully understood. Here, we report that DepR2, an ArsR-family transcriptional regulator isolated previously by DNA-affinity purification, interacts directly with dptEp, the major promoter of the daptomycin gene cluster. DepR2 binds to an imperfect palindromic sequence at the very upstream of dptEp. Meanwhile, higher dptEp activities were consistently observed in the ΔdepR2 mutant, correlating with a nearly 2.5-fold increased production of daptomycin and three structurally related secondary metabolites A21978C1–3. Thus, our data suggest that the ArsR-family transcriptional regulator DepR2 negatively regulates production of daptomycin by directly repressing the expression of its gene cluster in S. roseosporus. To the best of our knowledge, this is the first report to show the involvement of an ArsR-family regulator in the direct regulation of secondary metabolite biosynthesis in Streptomyces.

Metallochaperones and metalloregulation in bacteria

Bacterial transition metal homoeostasis or simply 'metallostasis' describes the process by which cells control the intracellular availability of functionally required metal cofactors, from manganese (Mn) to zinc (Zn), avoiding both metal deprivation and toxicity. Metallostasis is an emerging aspect of the vertebrate host-pathogen interface that is defined by a 'tug-of-war' for biologically essential metals and provides the motivation for much recent work in this area. The host employs a number of strategies to starve the microbial pathogen of essential metals, while for others attempts to limit bacterial infections by leveraging highly competitive metals. Bacteria must be capable of adapting to these efforts to remodel the transition metal landscape and employ highly specialized metal sensing transcriptional regulators, termed metalloregulatory proteins,and metallochaperones, that allocate metals to specific destinations, to mediate this adaptive response. In this essay, we discuss recent progress in our understanding of the structural mechanisms and metal specificity of this adaptive response, focusing on energy-requiring metallochaperones that play roles in the metallocofactor active site assembly in metalloenzymes and metallosensors, which govern the systems-level response to metal limitation and intoxication.

Solution NMR Spectroscopy for the Study of Enzyme Allostery

Allostery is a ubiquitous biological regulatory process in which distant binding sites within a protein or enzyme are functionally and thermodynamically coupled. Allosteric interactions play essential roles in many enzymological mechanisms, often facilitating formation of enzyme-substrate complexes and/or product release. Thus, elucidating the forces that drive allostery is critical to understanding the complex transformations of biomolecules. Currently, a number of models exist to describe allosteric behavior, taking into account energetics as well as conformational rearrangements and fluctuations. In the following Review, we discuss the use of solution NMR techniques designed to probe allosteric mechanisms in enzymes. NMR spectroscopy is unequaled in its ability to detect structural and dynamical changes in biomolecules, and the case studies presented herein demonstrate the range of insights to be gained from this valuable method. We also provide a detailed technical discussion of several specialized NMR experiments that are ideally suited for the study of enzymatic allostery.

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.

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.

Characterizing solution surface loop conformational flexibility of the GM2 activator protein

GM2AP has a β-cup topology with numerous X-ray structures showing multiple conformations for some of the surface loops, revealing conformational flexibility that may be related to function, where function is defined as either membrane binding associated with ligand binding and extraction or interaction with other proteins. Here, site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy and molecular dynamic (MD) simulations are used to characterize the mobility and conformational flexibility of various structural regions of GM2AP. A series of 10 single cysteine amino acid substitutions were generated, and the constructs were chemically modified with the methanethiosulfonate spin label. Continuous wave (CW) EPR line shapes were obtained and subsequently simulated using the microscopic order macroscopic disorder (MOMD) program. Line shapes for sites that have multiple conformations in the X-ray structures required two spectral components, whereas spectra of the remaining sites were adequately fit with single-component parameters. For spin labeled sites L126C and I66C, spectra were acquired as a function of temperature, and simulations provided for the determination of thermodynamic parameters associated with conformational change. Binding to GM2 ligand did not alter the conformational flexibility of the loops, as evaluated by EPR and NMR spectroscopies. These results confirm that the conformational flexibility observed in the surface loops of GM2AP crystals is present in solution and that the exchange is slow on the EPR time scale (>ns). Furthermore, MD simulation results are presented and agree well with the conformational heterogeneity revealed by SDSL.

Recent developments in copper and zinc homeostasis in bacterial pathogens

Copper and zinc homeostasis systems in pathogenic bacteria are required to resist host efforts to manipulate the availability and toxicity of these metal ions. Central to this microbial adaptive response is the involvement of metal-trafficking and metal-sensing proteins that ultimately exercise control of metal speciation in the cell. Cu-specific and Zn-specific metalloregulatory proteins regulate the transcription of metal-responsive genes while metallochaperones and related proteins ensure that these metals are appropriately buffered by the intracellular milieu and delivered to correct intracellular targets. In this review, we summarize recent findings on how bacterial pathogens mount a metal-specific response to derail host efforts to win the 'fight over metals.'

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.