C-N Coupling of DNA-Conjugated (Hetero)aryl Bromides and Chlorides for DNA-Encoded Chemical Library Synthesis.

DNA-encoded chemical library (DECL) screens are a rapid and economical tool to identify chemical starting points for drug discovery. As a robust transformation for drug discovery, palladium-catalyzed C-N coupling is a valuable synthetic method for the construction of DECL chemical matter; however, currently disclosed methods have only been demonstrated on DNA-attached (hetero)aromatic iodide and bromide electrophiles. We developed conditions utilizing an N-heterocyclic carbene-palladium catalyst that extends this reaction to the coupling of DNA-conjugated (hetero)aromatic chlorides with (hetero)aromatic and select aliphatic amine nucleophiles. In addition, we evaluated steric and electronic effects within this catalyst series, carried out a large substrate scope study on two representative (hetero)aryl bromides, and applied this newly developed method within the construction of a 63 million-membered DECL.

Molecular Dynamics Study of Twister Ribozyme: Role of Mg(2+) Ions and the Hydrogen-Bonding Network in the Active Site.

The recently discovered twister ribozyme is thought to utilize general acid-base catalysis in its self-cleavage mechanism, but the roles of nucleobases and metal ions in the mechanism are unclear. Herein, molecular dynamics simulations of the env22 twister ribozyme are performed to elucidate the structural and equilibrium dynamical properties, as well as to examine the role of Mg(2+) ions and possible candidates for the general base and acid in the self-cleavage mechanism. The active site region and the ends of the pseudoknots were found to be less mobile than other regions of the ribozyme, most likely providing structural stability and possibly facilitating catalysis. A purported catalytic Mg(2+) ion and the closest neighboring Mg(2+) ion remained chelated and relatively immobile throughout the microsecond trajectories, although removal of these Mg(2+) ions did not lead to any significant changes in the structure or equilibrium motions of the ribozyme on the microsecond time scale. In addition, a third metal ion, a Na(+) ion remained close to A1(O5'), the leaving group atom, during the majority of the microsecond trajectories, suggesting that it might stabilize the negative charge on A1(O5') during self-cleavage. The locations of these cations and their interactions with key nucleotides in the active site suggest that they may be catalytically relevant. The P1 stem is partially melted at its top and bottom in the crystal structure and further unwinds in the trajectories. The simulations also revealed an interconnected network comprised of hydrogen-bonding and π-stacking interactions that create a relatively rigid network around the self-cleavage site. The nucleotides involved in this network are among the highly conserved nucleotides in twister ribozymes, suggesting that this interaction network may be important to structure and function.

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.

Relative Binding Free Energies of Adenine and Guanine to Damaged and Undamaged DNA in Human DNA Polymerase η: Clues for Fidelity and Overall Efficiency.

Human DNA polymerase η (Pol Î·) plays an essential protective role against skin cancer caused by cyclobutane thymine-thymine dimers (TTDs), a frequent form of DNA damage arising from exposure to the sun. This enzyme rescues stalled replication forks at the TTDs by inserting bases opposite these DNA defects. Herein we calculate binding free energies for a free deoxyribose nucleotide triphosphate, dATP or dGTP, to Pol Î· complexed with undamaged or damaged DNA. The calculations indicate that the binding of dATP to the enzyme-DNA complex is thermodynamically favored for TTD-containing DNA over undamaged DNA, most likely because of more extensive hydrogen-bonding interactions between the TTD and the enzyme that hold the TTD more rigidly in place. The calculations also illustrate that dATP binding is thermodynamically favored over dGTP binding at both thymine positions of the TTD, most likely due to more persistent and stable hydrogen-bonding interactions between the TTD and dATP than between the TTD and dGTP. This free energy difference is slightly greater for binding at the 5' thymine position than at the 3' thymine position, presumably because of stabilization arising from the A:T base pair formed at the 3' position of the TTD in the previous step of Pol Î· function. All of these trends in binding free energies are consistent with experimental measurements of binding strength, fidelity, processivity, and overall efficiency. The insights gained from this analysis have implications for drug design efforts aimed at modifying the binding properties of this enzyme for improving cancer chemotherapy treatments.

Comparative Molecular Dynamics Studies of Human DNA Polymerase η.

High-energy ultraviolet radiation damages DNA through the formation of cyclobutane pyrimidine dimers, which stall replication. When the lesion is a thymine-thymine dimer (TTD), human DNA polymerase η (Pol η) assists in resuming the replication process by inserting nucleotides opposite the damaged site. We performed extensive molecular dynamics (MD) simulations to investigate the structural and dynamical effects of four different Pol η complexes with or without a TTD and with either dATP or dGTP as the incoming base. No major differences in the overall structures and equilibrium dynamics were detected among the four systems, suggesting that the specificity of this enzyme is due predominantly to differences in local interactions in the binding regions. Analysis of the hydrogen-bonding interactions between the enzyme and the DNA and dNTP provided molecular-level insights. Specifically, the TTD was observed to engage in more hydrogen-bonding interactions with the enzyme than its undamaged counterpart of two normal thymines. The resulting greater rigidity and specific orientation of the TTD are consistent with the experimental observation of higher processivity and overall efficiency at TTD sites than at analogous sites with two normal thymines. The similarities between the systems containing dATP and dGTP are consistent with the experimental observation of relatively low fidelity with respect to the incoming base. Moreover, Q38 and R61, two strictly conserved amino acids across the Pol η family, were found to exhibit persistent hydrogen-bonding interactions with the TTD and cation-π interactions with the free base, respectively. Thus, these simulations provide molecular level insights into the basis for the selectivity and efficiency of this enzyme, as well as the roles of the two most strictly conserved residues.

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

CusCFBA is one of the metal efflux systems in Escherichia coli that is highly specific for its substrates, Cu(I) and Ag(I). It serves to protect the bacteria in environments that have lethal concentrations of these metals. The membrane fusion protein CusB is the periplasmic piece of CusCFBA, which has not been fully characterized by crystallography because of its extremely disordered N-terminal region. This region has both structural and functional importance because it has been experimentally proven to transfer the metal by itself from the metallochaperone CusF and to induce a structural change in the rest of CusB to increase Cu(I)/Ag(I) resistance. Understanding metal uptake from the periplasm is critical to gain insight into the mechanism of the whole CusCFBA pump, which makes resolving a structure for the N-terminal region necessary because it contains the metal binding site. We ran extensive molecular dynamics simulations to reveal the structural and dynamic properties of both the apo and Cu(I)-bound versions of the CusB N-terminal region. In contrast to its functional companion CusF, Cu(I) binding to the N-terminus of CusB causes only a slight, local stabilization around the metal site. The trajectories were analyzed in detail, revealing extensive structural disorder in both the apo and holo forms of the protein. CusB was further analyzed by breaking the protein up into three subdomains according to the extent of the observed disorder: the N- and C-terminal tails, the central beta strand motif, and the M21-M36 loop connecting the two metal-coordinating methionine residues. Most of the observed disorder was traced back to the tail regions, leading us to hypothesize that the latter two subdomains (residues 13-45) may form a functionally competent metal-binding domain because the tail regions appear to play no role in metal binding.

Insights into the mechanistic dichotomy of the protein farnesyltransferase peptide substrates CVIM and CVLS.

Protein farnesyltransferase (FTase) catalyzes farnesylation of a variety of peptide substrates. (3)H α-secondary kinetic isotope effect (α-SKIE) measurements of two peptide substrates, CVIM and CVLS, are significantly different and have been proposed to reflect a rate-limiting S(N)2-like transition state with dissociative characteristics for CVIM, while, due to the absence of an isotope effect, CVLS was proposed to have a rate-limiting peptide conformational change. Potential of mean force quantum mechanical/molecular mechanical studies coupled with umbrella sampling techniques were performed to further probe this mechanistic dichotomy. We observe the experimentally proposed transition state (TS) for CVIM but find that CVLS has a symmetric S(N)2 TS, which is also consistent with the absence of a (3)H α-SKIE. These calculations demonstrate facile substrate-dependent alterations in the transition state structure catalyzed by FTase.

Prediction of trypsin/molecular fragment binding affinities by free energy decomposition and empirical scores.

Two families of binding affinity estimation methodologies are described which were utilized in the SAMPL3 trypsin/fragment binding affinity challenge. The first is a free energy decomposition scheme based on a thermodynamic cycle, which included separate contributions from enthalpy and entropy of binding as well as a solvent contribution. Enthalpic contributions were estimated with PM6-DH2 semiempirical quantum mechanical interaction energies, which were modified with a statistical error correction procedure. Entropic contributions were estimated with the rigid-rotor harmonic approximation, and solvent contributions to the free energy were estimated with several different methods. The second general methodology is the empirical score LISA, which contains several physics-based terms trained with the large PDBBind database of protein/ligand complexes. Here we also introduce LISA+, an updated version of LISA which, prior to scoring, classifies systems into one of four classes based on a ligand's hydrophobicity and molecular weight. Each version of the two methodologies (a total of 11 methods) was trained against a compiled set of known trypsin binders available in the Protein Data Bank to yield scaling parameters for linear regression models. Both raw and scaled scores were submitted to SAMPL3. Variants of LISA showed relatively low absolute errors but also low correlation with experiment, while the free energy decomposition methods had modest success when scaling factors were included. Nonetheless, re-scaled LISA yielded the best predictions in the challenge in terms of RMS error, and six of these models placed in the top ten best predictions by RMS error. This work highlights some of the difficulties of predicting binding affinities of small molecular fragments to protein receptors as well as the benefit of using training data.

Insight into the cation-π interaction at the metal binding site of the copper metallochaperone CusF.

The periplasmic Cu(+)/Ag(+) chaperone CusF features a novel cation-π interaction between a Cu(+)/Ag(+) ion and Trp44 at the metal binding site. The nature and strength of the Cu(+)/Ag(+)-Trp44 interactions were investigated using computational methodologies. Quantum-mechanical (QM) calculations showed that the Cu(+) and Ag(+) interactions with Trp44 are of similar strength (~14 kcal/mol) and bond order. Quantum-mechanical/molecular-mechanical (QM/MM) calculations showed that Cu(+) binds in a distorted tetrahedral coordination environment in the Trp44Met mutant, which lacks the cation-π interaction. Molecular dynamics (MD) simulations of CusF in the apo and Cu(+)-bound states emphasized the importance of the Cu(+)-Trp44 interaction in protecting Cu(+) from water oxidation. The protein structure does not change over the time scale of hundreds of nanoseconds in the metal-bound state. The metal recognition site exhibits small motions in the apo state but remains largely preorganized toward metal binding. Trp44 remains oriented to form the cation-π interaction in the apo state and faces an energetic penalty to move away from the metal ion. Cu(+) binding quenches the protein's internal motions in regions linked to binding CusB, suggesting that protein motions play an essential role in Cu(+) transfer to CusB.

Pairwise additivity of energy components in protein-ligand binding: the HIV II protease-Indinavir case.

An energy expansion (binding energy decomposition into n-body interaction terms for n ≥ 2) to express the receptor-ligand binding energy for the fragmented HIV II protease-Indinavir system is described to address the role of cooperativity in ligand binding. The outcome of this energy expansion is compared to the total receptor-ligand binding energy at the Hartree-Fock, density functional theory, and semiempirical levels of theory. We find that the sum of the pairwise interaction energies approximates the total binding energy to ∼82% for HF and to >95% for both the M06-L density functional and PM6-DH2 semiempirical method. The contribution of the three-body interactions amounts to 18.7%, 3.8%, and 1.4% for HF, M06-L, and PM6-DH2, respectively. We find that the expansion can be safely truncated after n=3. That is, the contribution of the interactions involving more than three parties to the total binding energy of Indinavir to the HIV II protease receptor is negligible. Overall, we find that the two-body terms represent a good approximation to the total binding energy of the system, which points to pairwise additivity in the present case. This basic principle of pairwise additivity is utilized in fragment-based drug design approaches and our results support its continued use. The present results can also aid in the validation of non-bonded terms contained within common force fields and in the correction of systematic errors in physics-based score functions.

Identifying Oxacillinase-48 Carbapenemase Inhibitors Using DNA-Encoded Chemical Libraries.

Bacterial resistance to ß-lactam antibiotics is largely mediated by ß-lactamases, which catalyze the hydrolysis of these drugs and continue to emerge in response to antibiotic use. ß-Lactamases that hydrolyze the last resort carbapenem class of ß-lactam antibiotics (carbapenemases) are a growing global health threat. Inhibitors have been developed to prevent ß-lactamase-mediated hydrolysis and restore the efficacy of these antibiotics. However, there are few inhibitors available for problematic carbapenemases such as oxacillinase-48 (OXA-48). A DNA-encoded chemical library approach was used to rapidly screen for compounds that bind and potentially inhibit OXA-48. Using this approach, a hit compound, CDD-97, was identified with submicromolar potency (Ki = 0.53 ± 0.08 µM) against OXA-48. X-ray crystallography showed that CDD-97 binds noncovalently in the active site of OXA-48. Synthesis and testing of derivatives of CDD-97 revealed structure-activity relationships and informed the design of a compound with a 2-fold increase in potency. CDD-97, however, synergizes poorly with ß-lactam antibiotics to inhibit the growth of bacteria expressing OXA-48 due to poor accumulation into E. coli. Despite the low in vivo activity, CDD-97 provides new insights into OXA-48 inhibition and demonstrates the potential of using DNA-encoded chemistry technology to rapidly identify ß-lactamase binders and to study ß-lactamase inhibition, leading to clinically useful inhibitors.

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.

Effects of Active Site Mutations on Specificity of Nucleobase Binding in Human DNA Polymerase η.

Human DNA polymerase η (Pol η) plays a vital role in protection against skin cancer caused by damage from ultraviolet light. This enzyme rescues stalled replication forks at cyclobutane thymine-thymine dimers (TTDs) by inserting nucleotides opposite these DNA lesions. Residue R61 is conserved in the Pol η enzymes across species, but the corresponding residue, as well as its neighbor S62, is different in other Y-family polymerases, Pol ι and Pol κ. Herein, R61 and S62 are mutated to their Pol ι and Pol κ counterparts. Relative binding free energies of dATP to mutant Pol η•DNA complexes with and without a TTD were calculated using thermodynamic integration. The binding free energies of dATP to the Pol η•DNA complex with and without a TTD are more similar for all of these mutants than for wild-type Pol η, suggesting that these mutations decrease the ability of this enzyme to distinguish between a TTD lesion and undamaged DNA. Molecular dynamics simulations of the mutant systems provide insights into the molecular level basis for the changes in relative binding free energies. The simulations identified differences in hydrogen-bonding, cation-π, and π-π interactions of the side chains with the dATP and the TTD or thymine-thymine (TT) motif. The simulations also revealed that R61 and Q38 act as a clamp to position the dATP and the TTD or TT and that the mutations impact the balance among the interactions related to this clamp. Overall, these calculations suggest that R61 and S62 play key roles in the specificity and effectiveness of Pol η for bypassing TTD lesions during DNA replication. Understanding the basis for this specificity is important for designing drugs aimed at cancer treatment.

Bringing Clarity to the Prediction of Protein-Ligand Binding Free Energies via "Blurring"

We present a method to evaluate the free energies of ligand binding utilizing a Monte Carlo estimation of the configuration integrals concomitant with uncertainty quantification. Ensembles for integration are built through systematically perturbing an initial ligand conformation in a rigid binding pocket, which is optimized separately prior to incorporation of the ligand. We call the procedure producing the ensembles "blurring", and it is carried out using an in-house developed code. The Boltzmann factor contribution of each pose to the configuration integral is computed and from there the free energy is obtained. Potential function uncertainties are estimated using a fragment-based error propagation method. This method has been applied to a set of small aromatic ligands complexed with T4 Lysozyme L99A mutant. Microstate energies have been determined with the force fields ff99SB and ff94, and the semiempirical method PM6DH2 in conjunction with continuum solvation models including Generalized Born (GB), the Conductor-like Screening Model (COSMO), and SMD. Of the methods studied, PM6DH2-based scoring gave binding free energy estimates, which yielded a good correlation to the experimental binding affinities (R2 = 0.7). All methods overestimated the calculated binding affinities. We trace this to insufficient sampling, the single static protein structure, and inaccuracies in the solvent models we have used in this study.

Molecular Dynamics Study of Helicobacter pylori Urease.

Helicobacter pylori have been implicated in an array of gastrointestinal disorders including, but not limited to, gastric and duodenal ulcers and adenocarcinoma. This bacterium utilizes an enzyme, urease, to produce copious amounts of ammonia through urea hydrolysis in order to survive the harsh acidic conditions of the stomach. Molecular dynamics (MD) studies on the H. pylori urease enzyme have been employed in order to study structural features of this enzyme that may shed light on the hydrolysis mechanism. A total of 400 ns of MD simulation time were collected and analyzed in this study. A wide-open flap state previously observed in MD simulations on Klebsiella aerogenes [Roberts et al. J. Am. Chem. Soc.2012, 134, 9934] urease has been identified in the H. pylori enzyme that has yet to be experimentally observed. Critical distances between residues on the flap, contact points in the closed state, and the separation between the active site Ni2+ ions and the critical histidine α322 residue were used to characterize flap motion. An additional flap in the active site was elaborated upon that we postulate may serve as an exit conduit for hydrolysis products. Finally we discuss the internal hollow cavity and present analysis of the distribution of sodium ions over the course of the simulation.

The Movable Type Method Applied to Protein-Ligand Binding.

Accurately computing the free energy for biological processes like protein folding or protein-ligand association remains a challenging problem. Both describing the complex intermolecular forces involved and sampling the requisite configuration space make understanding these processes innately difficult. Herein, we address the sampling problem using a novel methodology we term "movable type". Conceptually it can be understood by analogy with the evolution of printing and, hence, the name movable type. For example, a common approach to the study of protein-ligand complexation involves taking a database of intact drug-like molecules and exhaustively docking them into a binding pocket. This is reminiscent of early woodblock printing where each page had to be laboriously created prior to printing a book. However, printing evolved to an approach where a database of symbols (letters, numerals, etc.) was created and then assembled using a movable type system, which allowed for the creation of all possible combinations of symbols on a given page, thereby, revolutionizing the dissemination of knowledge. Our movable type (MT) method involves the identification of all atom pairs seen in protein-ligand complexes and then creating two databases: one with their associated pairwise distant dependent energies and another associated with the probability of how these pairs can combine in terms of bonds, angles, dihedrals and non-bonded interactions. Combining these two databases coupled with the principles of statistical mechanics allows us to accurately estimate binding free energies as well as the pose of a ligand in a receptor. This method, by its mathematical construction, samples all of configuration space of a selected region (the protein active site here) in one shot without resorting to brute force sampling schemes involving Monte Carlo, genetic algorithms or molecular dynamics simulations making the methodology extremely efficient. Importantly, this method explores the free energy surface eliminating the need to estimate the enthalpy and entropy components individually. Finally, low free energy structures can be obtained via a free energy minimization procedure yielding all low free energy poses on a given free energy surface. Besides revolutionizing the protein-ligand docking and scoring problem this approach can be utilized in a wide range of applications in computational biology which involve the computation of free energies for systems with extensive phase spaces including protein folding, protein-protein docking and protein design.