Influence of solution ionic strength on the stabilities of M20 loop conformations in apo E. coli dihydrofolate reductase.

Interactions among ions and their specific interactions with macromolecular solutes are known to play a central role in biomolecular stability. However, similar effects in the conformational stability of protein loops that play functional roles, such as binding ligands, proteins, and DNA/RNA molecules, remain relatively unexplored. A well-characterized enzyme that has such a functional loop is Escherichia coli dihydrofolate reductase (ecDHFR), whose so-called M20 loop has been observed in three ordered conformations in crystal structures. To explore how solution ionic strengths may affect the M20 loop conformation, we proposed a reaction coordinate that could quantitatively describe the loop conformation and used it to classify the loop conformations in representative ecDHFR x-ray structures crystallized in varying ionic strengths. The Protein Data Bank survey indicates that at ionic strengths (I) below the intracellular ion concentration-derived ionic strength in E. coli (I ≤ 0.237M), the ecDHFR M20 loop tends to adopt open/closed conformations, and rarely an occluded loop state, but when I is >0.237M, the loop tends to adopt closed/occluded conformations. Distance-dependent electrostatic potentials around the most mobile M20 loop region from molecular dynamics simulations of ecDHFR in equilibrated CaCl2 solutions of varying ionic strengths show that high ionic strengths (I = 0.75/1.5M) can preferentially stabilize the loop in closed/occluded conformations. These results nicely correlate with conformations derived from ecDHFR structures crystallized in varying ionic strengths. Altogether, our results suggest caution in linking M20 loop conformations derived from crystal structures solved at ionic strengths beyond that tolerated by E. coli to the ecDHFR function.

The Zinc Linchpin Motif in the DNA Repair Glycosylase MUTYH: Identifying the Zn2+ Ligands and Roles in Damage Recognition and Repair.

The DNA base excision repair (BER) glycosylase MUTYH prevents DNA mutations by catalyzing adenine (A) excision from inappropriately formed 8-oxoguanine (8-oxoG):A mismatches. The importance of this mutation suppression activity in tumor suppressor genes is underscored by the association of inherited variants of MUTYH with colorectal polyposis in a hereditary colorectal cancer syndrome known as MUTYH-associated polyposis, or MAP. Many of the MAP variants encompass amino acid changes that occur at positions surrounding the two-metal cofactor-binding sites of MUTYH. One of these cofactors, found in nearly all MUTYH orthologs, is a [4Fe-4S]2+ cluster coordinated by four Cys residues located in the N-terminal catalytic domain. We recently uncovered a second functionally relevant metal cofactor site present only in higher eukaryotic MUTYH orthologs: a Zn2+ ion coordinated by three Cys residues located within the extended interdomain connector (IDC) region of MUTYH that connects the N-terminal adenine excision and C-terminal 8-oxoG recognition domains. In this work, we identified a candidate for the fourth Zn2+ coordinating ligand using a combination of bioinformatics and computational modeling. In addition, using in vitro enzyme activity assays, fluorescence polarization DNA binding assays, circular dichroism spectroscopy, and cell-based rifampicin resistance assays, the functional impact of reduced Zn2+ chelation was evaluated. Taken together, these results illustrate the critical role that the "Zn2+ linchpin motif" plays in MUTYH repair activity by providing for proper engagement of the functional domains on the 8-oxoG:A mismatch required for base excision catalysis. The functional importance of the Zn2+ linchpin also suggests that adjacent MAP variants or exposure to environmental chemicals may compromise Zn2+ coordination, and ability of MUTYH to prevent disease.

Modeling Zn²âº release from metallothionein.

Mammalian metallothioneins (MTs) comprise a Zn3Cys9 cluster in the ß domain and a Zn4Cys11 cluster in the α domain. They play a crucial role in storing and donating Zn(2+) ions to target metalloproteins and have been implicated in several diseases, thus understanding how MTs release Zn(2+) is of widespread interest. In this work, we present a strategy to compute the free energy for releasing Zn(2+) from MTs using a combination of classical molecular dynamics (MD) simulations, quantum-mechanics/molecular-mechanics (QM/MM) minimizations, and continuum dielectric calculations. The methodology is shown to reproduce the experimental observations that (1) the Zn-binding sites do not have equal Zn(2+) affinity and (2) the isolated ß domain is thermodynamically less stable and releases Zn(2+) faster with oxidizing agents than the isolated α domain. It was used to compute the free energies for Zn(2+) release from the metal cluster in the absence and presence of the protein matrix (protein architecture and coupled protein-water interactions) to yield the respective disulfide-bonded product. The results show the importance of the protein matrix as well as protein dynamics and coupled conformational changes in accounting for the differential Zn(2+)-releasing propensity of the two domains with oxidizing agents.

Solution Ionic Strength Can Modulate Functional Loop Conformations in E. coli Dihydrofolate Reductase.

The observation of multiple conformations of a functional loop (termed M20) in the Escherichia coli dihydrofolate reductase (ecDHFR) enzyme triggered the proposition that large-scale motions of protein structural elements contribute to enzyme catalysis. The transition of the M20 loop from a closed conformation to an occluded conformation was thought to aid the rate-limiting release of the products. However, the influence of charged species in the solution environment on the observed M20 loop conformations, independent of charged ligands bound to the enzyme, had not been considered. Molecular dynamics simulations of ecDHFR in model CaCl2 solutions of varying molar ionic strengths IM reveal a substantial free energy barrier between occluded and closed M20 loop states at IM exceeding the E. coli threshold (∼0.24 M). This barrier may facilitate crystallization of ecDHFR in the occluded state, consistent with ecDHFR structures obtained at IM exceeding 0.3 M. At lower IM (≤0.15 M), the M20 loop can explore the occluded state, but prefers an open/partially closed conformation, again consistent with ecDHFR structures. Our findings caution against using ecDHFR structures obtained at nonphysiological ionic strengths in interpreting catalytic events or in structure-based drug design.

Differential role of the protein matrix on the binding of a catalytic aspartate to Mg2+ vs Ca2+: application to ribonuclease H.

Divalent metal cations are essential cofactors for many enzyme functions. Although Mg(2+) is the native cofactor in many enzymes such as ribonuclease H, its competitor Ca(2+) may also bind to the enzyme but inhibit catalysis. Thus, the competition between Mg(2+) and Ca(2+) for a given metal-binding site in an enzyme and their effects on enzyme activity are of great interest. Most studies have focused on the interactions between Mg(2+) or Ca(2+) and the metal ligands in the first and sometimes second coordination shell. However, no study (to our knowledge) has examined the role of the protein architecture and surrounding aqueous environment on the binding of Mg(2+) vs Ca(2+) to a given protein metal-binding site. In this work, the free energy barriers for the binding of a catalytically essential aspartate to Mg(2+) or Ca(2+) in ribonuclease H from two organisms were computed using umbrella sampling with a classical force field ("classical" model). The corresponding free energy barriers in water were computed using the "classical" model as well as density functional theory combined with a self-consistent reaction field. The results reveal that, relative to water, the protein architecture and coupled protein-water interactions raise the free energy barrier for binding of the catalytically essential aspartate to the native Mg(2+) cofactor more than the respective binding to Ca(2+). They also reveal the physical basis for the different observed binding modes of Mg(2+) and Ca(2+) and highlight limitations of simulations with classical force fields that do not explicitly account for charge transfer and polarization effects.

Protein/solvent medium effects on Mg(2+)-carboxylate interactions in metalloenzymes.

We employed umbrella sampling molecular dynamics simulations in explicit water to study the binding of the Mg(2+) cofactor to ribonuclease H (RNase H) from three different organisms. We show that the enzyme can differentiate between different Mg(2+)-binding modes that are nearly equally stable by creating a free-energy barrier between a water-rich mode and a water-depleted mode. Through a comparison with the corresponding free-energy barrier in water, this effect is shown to emanate from the enzymes's three-dimensional architecture and its associated environment. Implications of these protein medium effects in RNase H function and in structure-based drug/molecular design are discussed.

Sensitivity of Functional Loop Conformations on Long-Range Electrostatics: Implications for M20 Loop Dynamics in E. coli Dihydrofolate Reductase.

In E. coli dihydrofolate reductase, unusual conformational motions of a functional M20 loop that interacts with substrate and coenzyme have been construed as evidence for dynamical effects in enzyme catalysis. By computing this loop's conformational free energies in the apoenzyme, we show that it is sensitive to the treatment of long-range electrostatic interactions and the solvation box size in modeling/simulations. These results provide important guidelines for computing reaction/binding free energy profiles of proteins with functional loops.

Structural Basis for the Magnesium-Dependent Activation and Hexamerization of the Lon AAA+ Protease.

The Lon AAA+ protease (LonA) plays important roles in protein homeostasis and regulation of diverse biological processes. LonA behaves as a homomeric hexamer in the presence of magnesium (Mg(2+)) and performs ATP-dependent proteolysis. However, it is also found that LonA can carry out Mg(2+)-dependent degradation of unfolded protein substrate in an ATP-independent manner. Here we show that in the presence of Mg(2+) LonA forms a non-secluded hexameric barrel with prominent openings, which explains why Mg(2+)-activated LonA can operate as a diffusion-based chambered protease to degrade unstructured protein and peptide substrates efficiently in the absence of ATP. A 1.85 Å crystal structure of Mg(2+)-activated protease domain reveals Mg(2+)-dependent remodeling of a substrate-binding loop and a potential metal-binding site near the Ser-Lys catalytic dyad, supported by biophysical binding assays and molecular dynamics simulations. Together, these findings reveal the specific roles of Mg(2+) in the molecular assembly and activation of LonA.

Conserved structural motif for recognizing nicotinamide adenine dinucleotide in poly(ADP-ribose) polymerases and ADP-ribosylating toxins: implications for structure-based drug design.

The ADP-ribosylating toxins (ADPRTs) and poly(ADP-ribose) polymerases (PARPs) are two important drug-target protein families. Although the Y-X(10)-Y motif for the diphtheria toxin group and the STS motif for the other ADPRTs have been found to recognize the NAD(+) substrate, it is not known (i) if these two different motifs share any structural similarity, (ii) the key forces/residues contributing to NAD(+) binding, and (iii) if they recognize the same or different NAD(+) conformations. Here, we show that even though the different toxin groups and PARPs share insignificant sequence identity, they share a similar 3D structure shaped like a scorpion (the "scorpion" motif) whose first three and last residues interact mainly with the NAD(+) nicotinamide ring via van der Waals forces. This locally conserved structure binds the nicotinamide mononucleotide moiety in a structurally conserved ringlike conformation. The biological implications/applications of locally conserved structures for toxins/PARPs and the nicotinamide mononucleotide are discussed.

Empirical force fields for biologically active divalent metal cations in water.

We have presented a strategy for deriving ion-water van der Waals (vdW) parameters that implicitly include the microscopic solvent molecular effects around the ion. The strategy can be used to obtain vdW parameters for metal cations of the same formal charge and known experimental hydration free energies. In this work, it was applied to derive the vdW parameters for 24 divalent metal ions with measured hydration free energies ranging from -300 to -572 kcal/mol, coordination numbers (CNs) from 4 to 15, and ion-O (water) distances from 1.67 to 2.90 angstroms. The strategy used to derive the vdW parameters employs (1) a numerical procedure that links the coupling parameter used in free energy simulations with the experimental hydration free energies and (2) the first-shell CNs and structure for the entire series of divalent cations. One of the parameter sets obtained (referred to as MWc) simultaneously reproduces the observed (i) relative hydration free energies, (ii) first-shell CNs, and (iii) average ion-water distances of all the dications studied. In particular, the MWc parameters reproduce the observed (i) decrease in the CN from 6 for Cu2+ to 4 for Be2+, (ii) no change in the CN of 6 for dications with hydration free energies between those of Cu2+ and Cd2+, and (iii) an expansion of the CN from 6 for Cd2+ to 9.5 for Ba2+. The ion-water parameters derived herein represent a first step in the simulations of metalloproteins, which will also require potential energy functions incorporating polarizability, charge transfer, and other electronic effects to accurately model the protein-metal interactions in aqueous solution.

A combined experimental and theoretical study of divalent metal ion selectivity and function in proteins: application to E. coli ribonuclease H1.

Structural and thermodynamic aspects of alkaline earth metal dication (Mg(2+), Ca(2+), Sr(2+), Ba(2+)) binding to E. coli ribonuclease H1 (RNase H1) have been investigated using both experimental and theoretical methods. The various metal-binding modes of the enzyme were explored using classical molecular dynamics simulations, and relative binding free energies were subsequently evaluated by free energy simulations. The trends in the free energies of model systems based on the simulation structures were subsequently verified using a combination of density functional theory and continuum dielectric methods. The calculations provide a physical basis for the experimental results and suggest plausible role(s) for the metal cation and the catalytically important acidic residues in protein function. Magnesium ion indirectly activates water attack of the phosphorus atom by freeing one of the active site carboxylate residues, D70, to act as a general base through its four first-shell water molecules, which prevent D70 from binding directly to Mg(2+). Calcium ion, on the other hand, inhibits enzyme activity by preventing D70 from acting as a general base through bidentate interactions with both carboxylate oxygen atoms of D70. These additional interactions to D70, in addition to the D10 and E48 monodentate interactions found for Mg(2+), enable Ca(2+) to bind tighter than the other divalent ions. However, a bare Mg(2+) ion with two or less water molecules in the first shell could bind directly to the three active-site carboxylates, in particular D70, thus inhibiting enzymatic activity. The present analyses and results could be generalized to other members of the RNase H family that possess the same structural fold and show similar metal-binding site and Mg(2+)-dependent activity.