Dynamical allosterism in the mechanism of action of DNA mismatch repair protein MutS

Reconnaissance of Allostery via the Restoration of Native p53 DNA-Binding Domain Dynamics in Y220C Mutant p53 Tumor Suppressor Protein

Allosteric regulation of protein dynamics infers a long-range deliberate propagation of information via micro- and macroscale interactions. The Y220C structural mutant is one of the most frequent cancerous p53 mutants. The mutation is distally located from the DNA-binding site of the p53 DNA-binding domain yet causes changes in DNA recognition. This system presents a unique opportunity to examine the allosteric control of mutated proteins under a drug design paradigm. We focus on the key case study of p53 Y220C mutation restoration by a series of new compounds suggested to have Y220C reactivation properties in comparison to our previous findings on the restorative potential of PK11000, a compound studied extensively for reactivation in vitro and in vivo. Previously, we implemented all-atom molecular dynamics (MD) simulations and our lab's techniques of MD-Sectors and MD-Markov state models on the wild type, the Y220C mutant, and Y220C with PK11000 to characterize the effector's restorative properties in terms of conformational dynamics and hydrogen bonding. In this study, we turn to probing the effects made by docking the battery of a new but less well-tested set of aminobenzothiazole derivative compounds reported by Baud et al., which show promise of Y220C rescue. We find that while complete and precise reconstitution of p53 WT molecular dynamics may not be observed as was the case with PK11000, dispersed local reconstitution of loop dynamics provides evidence of rescuing effects by aminobenzothiazole derivative N,2-dihydroxy-3,5-diiodo-4-(1H-pyrrol-1-yl)benzamide, Effector 22, like what we observed for PK11000. Generalizable insights into the mutation and allosteric reactivation of p53 by various effectors by reconstitution of WT dynamics observed in statistical conformational ensemble analysis and network inference are discussed, considering the development of allosteric drug design rooted in first principles.

Mutagenic Activation of Glutathione Peroxidase-4: Approaches toward Rational Design of Allosteric Drugs

Glutathione peroxidase 4 (GPX4) reduces lipid hydroperoxides in lipid membranes, effectively inhibiting iron-dependent cell death or ferroptosis. The upregulation of the enzyme by the mutations at residues D21 and D23 has been suggested to be associated with higher protein activity, which confers more protection against neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases. Therefore, it has become an attractive target for treating and preventing neurodegenerative diseases. However, identifying means of mimicking the beneficial effects of these mutations distant from the active site constitutes a formidable challenge in moving toward therapeutics. In this study, we explore using molecular dynamics simulations to computationally map the conformational and energetic landscape of the wild-type GPX4 protein and three mutant variants to identify the allosteric networks of the enzyme. We present the conformational dynamic profile providing the desired signature behavior of the enzyme. We also discuss the implications of these findings for drug design efforts.

Insights into Rational Design of a New Class of Allosteric Effectors with Molecular Dynamics Markov State Models and Network Theory

The development of drugs to restore protein function has been a major advance facilitated by molecular medicine. Allosteric regulation, a phenomenon widely observed in nature, in which a molecule binds to control a distance active site, holds great promise for regulating proteins, yet how to rationally design such a molecule remains a mystery. Over the past few years, we and others have developed several techniques based on molecular dynamics (MD) simulations: MD-Markov state models to capture global conformational substates, and network theory approach utilizing the interaction energy within the protein to confer local allosteric control. We focus on the key case study of the p53 Y220C mutation restoration by PK11000, a compound experimentally shown to reactivate p53 native function in Y220C mutant present tumors. We gain insights into the mutation and allosteric reactivation of the protein, which we anticipate will be applicable to de novo design to engineer new compounds not only for this mutation, but in other macromolecular systems as well.

A computational study of Tat-CDK9-Cyclin binding dynamics and its implication in transcription-dependent HIV latency

HIV is a virus that attacks the T cells. HIV may either actively replicate or become latent within host cells for years. Since HIV uses its own protein Tat to hijack the host CDK9-Cyclin complex for transcription, Tat is implicated in transcription-dependent HIV latency. To quantify the impact of Tat binding, we propose a computational framework to probe the dynamics of the CDK9-Cyclin interface and the ATP pocket reorganization upon binding by different Tat mutants. Specifically, we focus on mutations at three Tat residues P10, W11, and N12 that are known to interact directly with CDK9 based on the crystal structure of the Tat-CDK9-Cyclin complex. Our molecular dynamics simulations show that the CDK9-Cyclin interface becomes slightly weaker for P10S and W11R mutants but tighter for the K12N mutant. Furthermore, the side chain orientation of residue K48 in the ATP pocket of CDK9 is similar to the inactive state in P10S and W11R simulations, but similar to the active state in K12N simulations. These are consistent with some existing but puzzling observations of latency for these mutants. This framework may hence help gain a better understanding of the role of Tat in the transcription-dependent HIV latency establishment.

Spectral analysis of molecular dynamics simulations on PDZ: MD sectors

The idea of protein "sectors" posits that sparse subsets of amino acid residues form cooperative networks that are key elements of protein stability, ligand binding, and allosterism. To date, protein sectors have been calculated by the statistical coupling analysis (SCA) method of Ranganathan and co-workers via the spectral analysis of conservation-weighted evolutionary covariance matrices obtained from a multiple sequence alignments of homologous families of proteins. SCA sectors, a knowledge-based protocol, have been indentified with functional properties and allosterism for a number of systems. In this study, we investigate the utility of the sector idea for the analysis of physics-based molecular dynamics (MD) trajectories of proteins. Our test case for this procedure is PSD95- PDZ3, one of the smallest proteins for which allosterism has been observed. It has served previously as a model system for a number of prediction algorithms, and is well characterized by X-ray crystallography, NMR spectroscopy and site specific mutagenisis. All-atom MD simulations were performed for a total of 500 nanoseconds using AMBER, and MD-calculated covariance matrices for the fluctuations of residue displacements and non-bonded interaction energies were subjected to spectral analysis in a manner analogous to that of SCA. The composition of MD sectors was compared with results from SCA, site specific mutagenesis, and allosterism. The concordance indicates that MD sectors are a viable protocol for analyzing MD trajectories and provide insight into the physical origin of the phenomenon.Communicated by Ramaswamy H. Sarma.

Coordinated protein and DNA conformational changes govern mismatch repair initiation by MutS

MutS homologs identify base-pairing errors made in DNA during replication and initiate their repair. In the presence of adenosine triphosphate, MutS induces DNA bending upon mismatch recognition and subsequently undergoes conformational transitions that promote its interaction with MutL to signal repair. In the absence of MutL, these transitions lead to formation of a MutS mobile clamp that can move along the DNA. Previous single-molecule FRET (smFRET) studies characterized the dynamics of MutS DNA-binding domains during these transitions. Here, we use protein–DNA and DNA–DNA smFRET to monitor DNA conformational changes, and we use kinetic analyses to correlate DNA and protein conformational changes to one another and to the steps on the pathway to mobile clamp formation. The results reveal multiple sequential structural changes in both MutS and DNA, and they suggest that DNA dynamics play a critical role in the formation of the MutS mobile clamp. Taking these findings together with data from our previous studies, we propose a unified model of coordinated MutS and DNA conformational changes wherein initiation of mismatch repair is governed by a balance of DNA bending/unbending energetics and MutS conformational changes coupled to its nucleotide binding properties.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Applications of Ruthenium Complex in Tumor Diagnosis and Therapy

Ruthenium complexes are a new generation of metal antitumor drugs that are currently of great interest in multidisciplinary research. In this review article, we introduce the applications of ruthenium complexes in the diagnosis and therapy of tumors. We focus on the actions of ruthenium complexes on DNA, mitochondria, and endoplasmic reticulum of cells, as well as signaling pathways that induce tumor cell apoptosis, autophagy, and inhibition of angiogenesis. Furthermore, we highlight the use of ruthenium complexes as specific tumor cell probes to dynamically monitor the active biological component of the microenvironment and as excellent photosensitizer, catalyst, and bioimaging agents for phototherapies that significantly enhance the diagnosis and therapeutic effect on tumors. Finally, the combinational use of ruthenium complexes with existing clinical antitumor drugs to synergistically treat tumors is discussed.

Dependence of prevalence of contiguous pathways in proteins on structural complexity

Allostery is a regulatory mechanism in proteins where an effector molecule binds distal from an active site to modulate its activity. Allosteric signaling may occur via a continuous path of residues linking the active and allosteric sites, which has been suggested by large conformational changes evident in crystal structures. An alternate possibility is that the signal occurs in the realm of ensemble dynamics via an energy landscape change. While the latter was first proposed on theoretical grounds, increasing evidence suggests that such a control mechanism is plausible. A major difficulty for testing the two methods is the ability to definitively determine that a residue is directly involved in allosteric signal transduction. Statistical Coupling Analysis (SCA) is a method that has been successful at predicting pathways, and experimental tests involving mutagenesis or domain substitution provide the best available evidence of signaling pathways. However, ascertaining energetic pathways which need not be contiguous is far more difficult. To date, simple estimates of the statistical significance of a pathway in a protein remain to be established. The focus of this work is to estimate such benchmarks for the statistical significance of contiguous pathways for the null model of selecting residues at random. We found that when 20% of residues in proteins are randomly selected, contiguous pathways at the 6 Å cutoff level were found with success rates of 51% in PDZ, 30% in p53, and 3% in MutS. The results suggest that the significance of pathways may have system specific factors involved. Furthermore, the possible existence of false positives for contiguous pathways implies that signaling could be occurring via alternate routes including those consistent with the energetic landscape model.

Evolutionary Covariance Combined with Molecular Dynamics Predicts a Framework for Allostery in the MutS DNA Mismatch Repair Protein

Mismatch repair (MMR) is an essential, evolutionarily conserved pathway that maintains genome stability by correcting base-pairing errors in DNA. Here we examine the sequence and structure of MutS MMR protein to decipher the amino acid framework underlying its two key activities—recognizing mismatches in DNA and using ATP to initiate repair. Statistical coupling analysis (SCA) identified a network (sector) of coevolved amino acids in the MutS protein family. The potential functional significance of this SCA sector was assessed by performing molecular dynamics (MD) simulations for alanine mutants of the top 5% of 160 residues in the distribution, and control nonsector residues. The effects on three independent metrics were monitored: (i) MutS domain conformational dynamics, (ii) hydrogen bonding between MutS and DNA/ATP, and (iii) relative ATP binding free energy. Each measure revealed that sector residues contribute more substantively to MutS structure–function than nonsector residues. Notably, sector mutations disrupted MutS contacts with DNA and/or ATP from a distance via contiguous pathways and correlated motions, supporting the idea that SCA can identify amino acid networks underlying allosteric communication. The combined SCA/MD approach yielded novel, experimentally testable hypotheses for unknown roles of many residues distributed across MutS, including some implicated in Lynch cancer syndrome.

Use of restrained molecular dynamics to predict the conformations of phosphorylated receiver domains in two-component signaling systems

Two‐component signaling (TCS) is the primary means by which bacteria, as well as certain plants and fungi, respond to external stimuli. Signal transduction involves stimulus‐dependent autophosphorylation of a sensor histidine kinase and phosphoryl transfer to the receiver domain of a downstream response regulator. Phosphorylation acts as an allosteric switch, inducing structural and functional changes in the pathway's components. Due to their transient nature, phosphorylated receiver domains are challenging to characterize structurally. In this work, we provide a methodology for simulating receiver domain phosphorylation to predict conformations that are nearly identical to experimental structures. Using restrained molecular dynamics, phosphorylated conformations of receiver domains can be reliably sampled on nanosecond timescales. These simulations also provide data on conformational dynamics that can be used to identify regions of functional significance related to phosphorylation. We first validated this approach on several well‐characterized receiver domains and then used it to compare the upstream and downstream components of the fungal Sln1 phosphorelay. Our results demonstrate that this technique provides structural insight, obtained in the absence of crystallographic or NMR information, regarding phosphorylation‐induced conformational changes in receiver domains that regulate the output of their associated signaling pathway. To our knowledge, this is the first time such a protocol has been described that can be broadly applied to TCS proteins for predictive purposes. Proteins 2016; 85:155–176. © 2016 Wiley Periodicals, Inc.

Long-Range Signaling in MutS and MSH Homologs via Switching of Dynamic Communication Pathways

Allostery is conformation regulation by propagating a signal from one site to another distal site. This study focuses on the long-range communication in DNA mismatch repair proteins MutS and its homologs where intramolecular signaling has to travel over 70 Å to couple lesion detection to ATPase activity and eventual downstream repair. Using dynamic network analysis based on extensive molecular dynamics simulations, multiple preserved communication pathways were identified that would allow such long-range signaling. The pathways appear to depend on the nucleotides bound to the ATPase domain as well as the type of DNA substrate consistent with previously proposed functional cycles of mismatch recognition and repair initiation by MutS and homologs. A mechanism is proposed where pathways are switched without major conformational rearrangements allowing for efficient long-range signaling and allostery.

Mechanism for verification of mismatched and homoduplex DNAs by nucleotides-bound MutS analyzed by molecular dynamics simulations

In order to understand how MutS recognizes mismatched DNA and induces the reaction of DNA repair using ATP, the dynamics of the complexes of MutS (bound to the ADP and ATP nucleotides, or not) and DNA (with mismatched and matched base-pairs) were investigated using molecular dynamics simulations. As for DNA, the structure of the base-pairs of the homoduplex DNA which interacted with the DNA recognition site of MutS was intermittently disturbed, indicating that the homoduplex DNA was unstable. As for MutS, the disordered loops in the ATPase domains, which are considered to be necessary for the induction of DNA repair, were close to (away from) the nucleotide-binding sites in the ATPase domains when the nucleotides were (not) bound to MutS. This indicates that the ATPase domains changed their structural stability upon ATP binding using the disordered loop. Conformational analysis by principal component analysis showed that the nucleotide binding changed modes which have structurally solid ATPase domains and the large bending motion of the DNA from higher to lower frequencies. In the MutS-mismatched DNA complex bound to two nucleotides, the bending motion of the DNA at low frequency modes may play a role in triggering the formation of the sliding clamp for the following DNA-repair reaction step. Moreover, MM-PBSA/GBSA showed that the MutS-homoduplex DNA complex bound to two nucleotides was unstable because of the unfavorable interactions between MutS and DNA. This would trigger the ATP hydrolysis or separation of MutS and DNA to continue searching for mismatch base-pairs. Proteins 2016; 84:1287-1303. © 2016 Wiley Periodicals, Inc.

Recapturing the Correlated Motions of Protein Using Coarse- Grained Models

Long-range interactions and allostery are important for many biological processes. Increasing numbers of studies, both experimental and computational, show that internal dynamics may play an important role in such behaviors. Investigating the dynamical effects of proteins, how- ever, is a challenging problem using all-atom molecular dynamics because of the length-scales and timescales involved. As a result, coarse-grained models are often implemented. Herein, we use three well-defined coarse-grained models: Go, Martini and Cafemol, and a small model protein Eglin C, which is readily studied via all-atom molecular dynamics, to examine if these coarse grained models can explore the dynamics of Eglin C accurately as well as to see how these models respond to mutations. We found that all three models can recapture the dynamics of Eglin C to a significant extent - where we focus on root-mean square fluctuations and correlated motions as dynamical measures - but that the Cafemol and Go models are superior. The best agreement with all-atom simulations is for structured regions of Eglin C.

Mismatch binding, ADP-ATP exchange and intramolecular signaling during mismatch repair

The focus of this article is on the DNA binding and ATPase activities of the mismatch repair (MMR) protein, MutS-our current understanding of how this protein uses ATP to fuel its actions on DNA and initiate repair via interactions with MutL, the next protein in the pathway. Structure-function and kinetic studies have yielded detailed views of the MutS mechanism of action in MMR. How MutS and MutL work together after mismatch recognition to enable strand-specific nicking, which leads to strand excision and synthesis, is less clear and remains an active area of investigation.

Molecular Dynamics Simulations of Ternary Complexes: Comparisons of LEAFY Protein Binding to Different DNA Motifs

LEAFY (LFY) is a plant-specific transcription factor, with a variety of roles in different species. LFY contains a conserved DNA-binding domain (DBD) that determines its DNA-binding specificity. Recently, the structures of the dimeric LFY-DBD bound to different DNA motifs were successively solved by X-ray crystallography. In this article, molecular dynamics (MD) simulations are employed to study two crystal structures of DNA-bound LFY protein from angiosperms and the moss Physcomitrella patens, respectively. The comparison of stabilities of the two systems is consistent with the experimental data of binding affinities. The calculation of hydrogen bonds showed that position 312 in LFY determines the difference of DNA-binding specificity. By using principal component analysis (PCA) and free energy landscape (FEL) methods, the open-close conformational change of the dimerization interface was found to be important for the system stability. At the dimerization interface, the protein-protein interaction has multiple influences on the cooperative DNA binding of LFY. The following analysis of DNA structural parameters further revealed that the protein-protein interaction contributes varying roles according to the specific DNA-binding efficiency. We propose that the protein-protein interaction serves a dual function as a connector between LFY monomers and a regulator of DNA-binding specificity. It will improve the robustness and adaptivity of the LFY-DNA ternary structure. This study provides some new insights into the understanding of the dynamics and interaction mechanism of dimeric LFY-DBD bound to DNA at the atomic level.

ATP binding and hydrolysis by Saccharomyces cerevisiae Msh2-Msh3 are differentially modulated by mismatch and double-strand break repair DNA substrates

In Saccharomyces cerevisiae, Msh2-Msh3-mediated mismatch repair (MMR) recognizes and targets insertion/deletion loops for repair. Msh2-Msh3 is also required for 3' non-homologous tail removal (3'NHTR) in double-strand break repair. In both pathways, Msh2-Msh3 binds double-strand/single-strand junctions and initiates repair in an ATP-dependent manner. However, we recently demonstrated that the two pathways have distinct requirements with respect to Msh2-Msh3 activities. We identified a set of aromatic residues in the nucleotide binding pocket (FLY motif) of Msh3 that, when mutated, disrupted MMR, but left 3'NHTR largely intact. One of these mutations, msh3Y942A, was predicted to disrupt the nucleotide sandwich and allow altered positioning of ATP within the pocket. To develop a mechanistic understanding of the differential requirements for ATP binding and/or hydrolysis in the two pathways, we characterized Msh2-Msh3 and Msh2-msh3Y942A ATP binding and hydrolysis activities in the presence of MMR and 3'NHTR DNA substrates. We observed distinct, substrate-dependent ATP hydrolysis and nucleotide turnover by Msh2-Msh3, indicating that the MMR and 3'NHTR DNA substrates differentially modify the ATP binding/hydrolysis activities of Msh2-Msh3. Msh2-msh3Y942A retained the ability to bind DNA and ATP but exhibited altered ATP hydrolysis and nucleotide turnover. We propose that both ATP and structure-specific repair substrates cooperate to direct Msh2-Msh3-mediated repair and suggest an explanation for the msh3Y942A separation-of-function phenotype.

Thumb inhibitor binding eliminates functionally important dynamics in the hepatitis C virus RNA polymerase

Hepatitis C virus (HCV) has infected almost 200 million people worldwide, typically causing chronic liver damage and severe complications such as liver failure. Currently, there are few approved treatments for viral infection. Thus, the HCV RNA-dependent RNA polymerase (gene product NS5B) has emerged as an important target for small molecule therapeutics. Potential therapeutic agents include allosteric inhibitors that bind distal to the enzyme active site. While their mechanism of action is not conclusively known, it has been suggested that certain inhibitors prevent a conformational change in NS5B that is crucial for RNA replication. To gain insight into the molecular origin of long-range allosteric inhibition of NS5B, we employed molecular dynamics simulations of the enzyme with and without an inhibitor bound to the thumb domain. These studies indicate that the presence of an inhibitor in the thumb domain alters both the structure and internal motions of NS5B. Principal components analysis identified motions that are severely attenuated by inhibitor binding. These motions may have functional relevance by facilitating interactions between NS5B and RNA template or nascent RNA duplex, with presence of the ligand leading to enzyme conformations with narrower and thus less accessible RNA binding channels. This study provides the first evidence for a mechanistic basis of allosteric inhibition in NS5B. Moreover, we present evidence that allosteric inhibition of NS5B results from intrinsic features of the enzyme free energy landscape, suggesting a common mechanism for the action of diverse allosteric ligands.