Search

The accumulation of amyloid-beta ( Aβ ) oligomers is recognized as a potential culprit in Alzheimer's disease (AD). Experimental studies show that melatonin, a hormone that mainly regulates circadian rhythm and sleep, can interact with Aβ peptides and disrupt the formation of oligomers. However, how melatonin inhibits the oligomerization of soluble Aβ is unclear. Here, by computational simulations, we investigate the effect of different levels of melatonin on the conformation of the Aβ 42 dimer. We find that the conformation of the Aβ 42 dimer is dependent on melatonin levels. When melatonin is absent, the dimer mainly forms a parallel β -sheet in the CHC region. When one melatonin molecule is present, the overall conformation of the dimer does not change much, but the N-terminal of the dimer tends to adopt antiparallel β -sheets. When two melatoinin molecules are present, the Aβ 42 dimer exhibits significant structural change, especially in its central region, resulting in a more compact conformation, and forms parallel β -sheets in the C-terminal. This conformational difference induced by different levels of melatoinin can shed light on the protective role of melatonin.

Thrombomodulin (TM), a transmembrane receptor integral to the anticoagulant pathway, governs thrombin's substrate specificity via interaction with thrombin's anion-binding exosite I. Despite its established role, the precise mechanisms underlying this regulatory function are yet to be fully unraveled. In this study, we deepen the understanding of these mechanisms through eight independent 1 μs all-atom simulations, analyzing thrombin both in its free form and when bound to TM fragments TM456 and TM56. Our investigations revealed distinct and significant conformational changes in thrombin mediated by the binding of TM56 and TM456. While TM56 predominantly influences motions within exosite I, TM456 orchestrates coordinated alterations across various loop regions, thereby unveiling a multifaceted modulatory role that extends beyond that of TM56. A highlight of our study is the identification of critical hydrogen bonds that undergo transformations during TM56 and TM456 binding, shedding light on the pivotal allosteric influence exerted by TM4 on thrombin's structural dynamics. This work offers a nuanced appreciation of TM's regulatory role in blood coagulation, paving the way for innovative approaches in the development of anticoagulant therapies and expanding the horizons in oncology therapeutics through a deeper understanding of molecular interactions in the coagulation pathway., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

Hydrogen bonds play a critical role in the folding and stability of proteins, such as proteins and nucleic acids, by providing strong and directional interactions. They help to maintain the secondary and 3D structure of proteins, and structural changes in these molecules often result from the formation or breaking of hydrogen bonds. To gain insights into these hydrogen bonding networks, we applied two machine learning models - a logistic regression model and a decision tree model - to study four variants of thrombin: wild-type, ΔK9, E8K, and R4A. Our results showed that both models have their unique advantages. The logistic regression model highlighted potential key residues (GLU295) in thrombin's allosteric pathways, while the decision tree model identified important hydrogen bonding motifs. This information can aid in understanding the mechanisms of folding in proteins and has potential applications in drug design and other therapies. The use of these two models highlights their usefulness in studying hydrogen bonding networks in proteins.

The proteolytic cleavage of C99 by γ-secretase is the last step in the production of amyloid-β (Aβ) peptides. Previous studies have shown that membrane lipid composition, cholesterol concentration, and mutation in the transmembrane helix modified the structures and fluctuations of C99. In this study, we performed atomistic molecular dynamics simulations of the homodimer of the 55-residue congener of the C-terminal domain of the amyloid protein precursor, C99(1-55), in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine-cholesterol lipid bilayer and compared the conformational ensemble of wild-type (WT) sequence to those of the A2T and D23N variants. These mutations are particularly interesting as the protective Alzheimer's disease (AD) A2T mutation is known to decrease Aβ production, whereas the early onset AD D23N mutation does not affect Aβ production. We found noticeable differences in the structural ensembles of the three sequences. In particular, A2T varies from both WT and D23N by having long-range effects on the population of the extracellular juxtamembrane helix, the interface between the G29xxx-G33xxx-G37 motifs, and the fluctuations of the transmembrane helical topologies.

Thrombin is a Na[Formula: see text]-activated serine protease existing in two forms targeted to procoagulant and anticoagulant activities, respectively. There is one Na[Formula: see text]-binding site that has been the focus of the study of the thrombin. However, molecular dynamics (MD) simulations suggest that there might be actually two Na[Formula: see text]-binding sites in thrombin and that Na[Formula: see text] ions can even bind to two sites simultaneously. In this study, we performed 12 independent 2-µs all-atom MD simulations for the wild-type (WT) thrombin and we studied the effects of the different Na[Formula: see text] binding modes on thrombin. From the root-mean-square fluctuations (RMSF) for the [Formula: see text]-carbons, we see that the atomic fluctuations mainly change in the 60s, 170s, and 220s loops, and the connection (residue 167 to 170). The correlation matrices for different binding modes suggest regions that may play an important role in thrombin's allosteric response and provide us a possible allosteric pathway for the sodium binding. Amorim-Hennig (AH) clustering tells us how the structure of the regions of interest changes on sodium binding. Principal component analysis (PCA) shows us how the different regions of thrombin change conformation together with sodium binding. Solvent-accessible surface area (SASA) exposes the conformational change in exosite I and catalytic triad. Finally, we argue that the double binding mode might be an inactive mode and that the kinetic scheme for the Na[Formula: see text] binding to thrombin might be a multiple-step mechanism rather than a 2-step mechanism., (© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Mutations in the TREX1 3' → 5' exonuclease are associated with a spectrum of autoimmune disease phenotypes in humans and mice. Failure to degrade DNA activates the cGAS-STING DNA-sensing pathway signaling a type-I interferon (IFN) response that ultimately drives immune system activation. TREX1 and the cGAS-STING DNA-sensing pathway have also been implicated in the tumor microenvironment, where TREX1 is proposed to degrade tumor-derived DNA that would otherwise activate cGAS-STING. If tumor-derived DNA were not degraded, the cGAS-STING pathway would be activated to promote IFN-dependent antitumor immunity. Thus, we hypothesize TREX1 exonuclease inhibition as a novel immunotherapeutic strategy. We present data demonstrating antitumor immunity in the TREX1 D18N mouse model and discuss theory surrounding the best strategy for TREX1 inhibition. Potential complications of TREX1 inhibition as a therapeutic strategy are also discussed., Competing Interests: WH, TH, and FP declare the filing of U.S. Provisional Application No. 62/706,167 Trex1 Inhibitors and Uses Thereof. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Hemphill, Simpson, Liu, Salsbury, Hollis, Grayson and Perrino.)

Thrombin plays an important role in the process of hemostasis and blood coagulation. Studies in thrombin can help us find ways to treat cancer because thrombin is able to reduce the characteristic hypercoagulability of cancer. Thrombin is composed of two chains, the light chain and the heavy chain. The function of the heavy chain has been largely explored, while the function of the light chain was obscured until several disease-associated mutations in the light chain come to light. In this study, we want to explore the dynamic and conformation effects of mutations on the light chain further to determine possible associations between mutation, conformational changes, and disease. The study, which is a follow-up for our studies on apo thrombin and the mutant, ΔK9, mainly focuses on the mutants E8K and R4A. E8K is a disease-associated mutation, and R4A is used to study the role of Arg4, which is suggested experimentally to play a critical role for thrombin's catalytic activities. We performed five all-atom one microsecond-scale molecular dynamics (MD) simulations for both E8K and R4A, and quantified the changes in the conformational ensemble of the mutants. From the root-mean-square fluctuations (RMSF) for the α-carbons, we find that the atomic fluctuations change in the mutants in the 60s loop and γ loop. The correlation coefficients for the α-carbons indicate that the correlation relation for atom-pairs in the protein is also impacted. The clustering analysis and the principal component analysis (PCA) consistently tell us that the catalytic pocket and the regulatory loops are destabilized by the mutations. We also find that there are two binding modes for Na + by clustering the vector difference between the Na + ions and the 220s loop. After further analysis, we find that there is a relation between the Na + binding and the rigidification of the γ loop, which may shed light on the mysterious role of the γ loop in thrombin.

The omega-3 and omega-6 polyunsaturated fatty acids are two important components of cell membranes in human brains. When incorporated into phospholipids, omega-3 slows the progression of Alzheimer's disease (AD), whereas omega-6 is linked to increased risk of AD. Little is known on the amyloid-β (Aβ) conformations in membranes rich in omega-3 and omega-6 phospholipids. Herein, the structural properties of the Aβ 29-42 dimer embedded in both fatty acid membranes were comparatively studied to a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine (POPC) bilayer using all-atom molecular dynamics (MD) simulations. Starting from α-helix, both omega-6 and omega-3 membranes promote new orientations and conformations of the dimer, in agreement with the observed dependence of Aβ production upon addition of these two fatty acids. This conformational result is corroborated by atomistic MD simulations of the dimer of the 99 amino acid C-terminal fragment of amyloid precursor protein spanning the residues 15-55. Starting from β-sheet, omega-6 membrane promotes helical and disordered structures of Aβ 29-42 dimer, whereas omega-3 membrane preserves the β-sheet structures differing however from those observed in POPC. Remarkably, the mixture of the two fatty acids and POPC depicts another conformational ensemble of the Aβ 29-42 dimer. This finding demonstrates that variation in the abundance of the molecular phospholipids, which changes with age, modulates membrane-embedded Aβ oligomerization.

Thrombin is a key component for chemotherapeutic and antithrombotic therapy development. As the physiologic and pathologic roles of the light chain still remain vague, here, we continue previous efforts to understand the impacts of the disease-associated single deletion of LYS9 in the light chain. By combining supervised and unsupervised machine learning methodologies and more traditional structural analyses on data from 10 μs molecular dynamics simulations, we show that the conformational ensemble of the ΔK9 mutant is significantly perturbed. Our analyses consistently indicate that LYS9 deletion destabilizes both the catalytic cleft and regulatory functional regions and result in some conformational changes that occur in tens to hundreds of nanosecond scaled motions. We also reveal that the two forms of thrombin each prefer a distinct binding mode of a Na + ion. We expand our understanding of previous experimental observations and shed light on the mechanisms of the LYS9 deletion associated bleeding disorder by providing consistent but more quantitative and detailed structural analyses than early studies in literature. With a novel application of supervised learning, i.e. the decision tree learning on the hydrogen bonding features in the wild-type and ΔK9 mutant forms of thrombin, we predict that seven pairs of critical hydrogen bonding interactions are significant for establishing distinct behaviors of wild-type thrombin and its ΔK9 mutant form. Our calculations indicate the LYS9 in the light chain has both localized and long-range allosteric effects on thrombin, supporting the opinion that light chain has an important role as an allosteric effector.

The fidelity of protein synthesis is largely dominated by the accurate recognition of transfer RNAs (tRNAs) by their cognate aminoacyl-tRNA synthetases. Aminoacylation of each tRNA with its cognate amino acid is necessary to maintain the accuracy of genetic code input. Aminoacylated tRNA Met functions in both initiation and elongation steps during protein synthesis. As a precursor to the investigation of a methionyl-tRNA synthetase-tRNA Met complex, presented here are the results of molecular dynamics (MD) for single nucleotide substitutions in the D-loop of tRNA Met (G15A, G18A, and G19A) probing structure/function relationships. The core of tRNA Met likely mediates an effective communication between the tRNA anticodon and acceptor ends, contributing an acceptor stem rearrangement to fit into the enzyme-active site. Simulations of Escherichia coli tRNA Met were performed for 1 μs four times each. The MD simulations showed changes in tRNA flexibility and long-range communication most prominently in the G18A variant. The results indicate that the overall tertiary structure of tRNA Met remains unchanged with these substitutions; yet, there are perturbations to the secondary structure. Network-based analysis of the hydrogen bond structure and correlated motion indicates that the secondary structure elements of the tRNA are highly intraconnected, but loosely interconnected. Specific nucleotides, including U8 and G22, stabilize the mutated structures and are candidates for substitution in future studies., Competing Interests: The authors declare no competing financial interest.

An important regulatory domain of NF-[Formula: see text]B Essential Modulator (NEMO) is a ubiquitin-binding zinc finger, with a tetrahedral CYS3HIS1 zinc-coordinating binding site. Two variations of NEMO's zinc finger are implicated in various disease states including ectodermal dysplasia and adult-onset glaucoma. To discern structural and dynamical differences between these disease states, we present results of 48-[Formula: see text]s of molecular dynamics simulations for three zinc finger systems each in two states, with and without zinc-bound and correspondingly appropriate cysteine thiol/thiolate configurations. The wild-type protein, often studied for its role in cancer, maintains the most rigid and conformationally stable zinc-bound configuration compared with the diseased counterparts. The glaucoma-related protein has persistent loss of secondary structure except within the dominant conformation. Conformational overlap between wild-type and glaucoma isoforms indicate a competitive binding mechanism may be substantial in the malfunctioning configuration, while the alpha-helical disruption of the ectodermal dysplasia suggests a loss of binding selectivity is responsible for aberrant function.

Aims: Peroxiredoxins (Prxs) are ubiquitous cysteine-based peroxidases involved in oxidant defense and signal transduction. Despite much study, the precise roles of conserved residues remain poorly defined. In this study, we carried out extensive functional and structural characterization of 10 variants of such residues in a model decameric bacterial Prx., Results: Three active site proximal mutations of Salmonella typhimurium AhpC, T43V, R119A, and E49Q, lowered catalytic efficiency with hydrogen peroxide by 4-5 orders of magnitude, but did not affect reactivity toward their reductant, AhpF. pK a values of the peroxidatic cysteine were also shifted up by 1-1.3 pH units for these and a decamer disruption mutant, T77I. Except for the decamer-stabilizing T77V, all mutations destabilized decamers in the reduced form. In the oxidized form, three mutants-T77V, T43A, and T43S-exhibited stabilized decamers and were more efficiently reduced by AhpF than wild-type AhpC. Crystal structures of most mutants were solved and many showed alterations in stability of the fully folded active site loop., Innovation: This is the first study of Prx mutants to comprehensively assess the effects of mutations on catalytic activities, the active site cysteine pK a , and the protein structure and oligomeric status., Conclusion: The Arg119 side chain must be properly situated for efficient catalysis, but for other debilitating variants, the functional defects could be explained by structural perturbations and/or associated decamer destabilization rather than direct effects. This underscores the importance of our comprehensive approach. A remarkable new finding was the preference of the reductant for decamers. Antioxid. Redox Signal. 28, 521-536.

Alzheimer's disease is linked to various types of aggregates of amyloid-β (Aβ) peptide and their interactions with protein receptors and neuronal cell membranes. Little is known on the impact of the electric field on membrane-embedded Aβ. Here we use atomistic molecular dynamics simulations to study the effects of a constant electric field on the conformations of Aβ 29-42 dimer inside a membrane, where the electric field has a strength of 20 mV/nm which exists across the membrane of a human neuron. Starting from α-helix peptides, the transmembrane electric field (TMEF) accelerates the conversion from the Gly-out substate to the Gly-side and Gly-in substates. Starting from β-sheet peptides, TMEF induces changes of the kink and tilt angles at Gly33 and Gly37. Overall, in the simulations totaling 10 μs, TMEF establishes new ground states for the dimer, similar to induced-fit in ligand binding. Our findings indicate that TMEF can stabilize rare conformations of amyloid peptides, and this could influence the cleavage of the amyloid precursor protein and the formation of β-sheet oligomers in membrane bilayers.

Correlated motion analysis provides a method for understanding communication between and dynamic similarities of biopolymer residues and domains. The typical equal-time correlation matrices-frequently visualized with pseudo-colorings or heat maps-quickly convey large regions of highly correlated motion but hide more subtle similarities of motion. Here we propose a complementary method for visualizing correlations within proteins (or general biopolymers) that quickly conveys intuition about which residues have a similar dynamic behavior. For grouping residues, we use the recently developed non-parametric clustering algorithm HDBSCAN. Although the method we propose here can be used to group residues using correlation as a similarity matrix-the most straightforward and intuitive method-it can also be used to more generally determine groups of residues which have similar dynamic properties. We term these latter groups "Dynamic Domains", as they are based not on spatial closeness but rather closeness in the column space of a correlation matrix. We provide examples of this method across three human proteins of varying size and function-the Nf-Kappa-Beta essential modulator, the clotting promoter Thrombin and the mismatch repair protein (dimer) complex MutS-alpha. Although the examples presented here are from all-atom molecular dynamics simulations, this visualization technique can also be used on correlations matrices built from any ensembles of conformations from experiment or computation., (© 2017 The Protein Society.)

Given their increasingly frequent usage, understanding the chemical and structural properties which allow therapeutic nucleic acids to promote the death of cancer cells is critical for medical advancement. One molecule of interest is a 10-mer of FdUMP (5-fluoro-2'-deoxyuridine-5'-O-monophosphate) also called F10. To investigate causes of structural stability, we have computationally restored the 2' oxygen on each ribose sugar of the phosphodiester backbone, creating FUMP[10]. Microsecond time-scale, all-atom, simulations of FUMP[10] in the presence of 150 mM MgCl 2 predict that the strand has a 45% probability of folding into a stable hairpin-like secondary structure. Analysis of 16 μs of data reveals phosphate interactions as likely contributors to the stability of this folded state. Comparison with polydT and polyU simulations predicts that FUMP[10]'s lowest order structures last for one to 2 orders of magnitude longer than similar nucleic acid strands. Here we provide a brief structural and conformational analysis of the predicted structures of FUMP[10], and suggest insights into its stability via comparison to F10, polydT, and polyU.

The formation of senile plaques in central neural system resulting from the aggregation of the amyloid β (Aβ) of 40 and 42 residues is one of the two hallmarks of Alzheimer's disease. Numerous experiments and computational studies have shown that the aggregation of Aβ peptides in vitro is very complex and depends on many factors such as pH, agitation, temperature, and peptide concentration. The impact of a static electric field (EF) on amyloid peptide aggregation has been much less studied, although EFs may have some applications to treat Parkinson's disease symptoms. Here, we study the influence of an EF strength of 20 mV/nm, present in the human brains, on the conformation of the Aβ 29-42 dimer. Our 7 μs non-equilibrium atomistic simulations in aqueous solution show that this field-strength promotes substantially the formation of β-hairpins, believed to be a very important intermediate state during aggregation. This work also suggests that structural biology experiments conducted under appropriate EF strengths may help reduce the conformational heterogeneity of Aβ 1-40 /Aβ 1-42 dimers and provide significant insights into their structures that may be disease-causing.

MutS α is a key component in the mismatch repair (MMR) pathway. This protein is responsible for initiating the signaling pathways for DNA repair or cell death. Herein we investigate this heterodimer's post-recognition, post-binding response to three types of DNA damage involving cytotoxic, anti-cancer agents-carboplatin, cisplatin, and FdU. Through a combination of supervised and unsupervised machine learning techniques along with more traditional structural and kinetic analysis applied to all-atom molecular dynamics (MD) calculations, we predict that MutS α has a distinct response to each of the three damage types. Via a binary classification tree (a supervised machine learning technique), we identify key hydrogen bond motifs unique to each type of damage and suggest residues for experimental mutation studies. Through a combination of a recently developed clustering (unsupervised learning) algorithm, RMSF calculations, PCA, and correlated motions we predict that each type of damage causes MutS α to explore a specific region of conformation space. Detailed analysis suggests a short range effect for carboplatin-primarily altering the structures and kinetics of residues within 10 angstroms of the damaged DNA-and distinct longer-range effects for cisplatin and FdU. In our simulations, we also observe that a key phenylalanine residue-known to stack with a mismatched or unmatched bases in MMR-stacks with the base complementary to the damaged base in 88.61% of MD frames containing carboplatinated DNA. Similarly, this Phe71 stacks with the base complementary to damage in 91.73% of frames with cisplatinated DNA. This residue, however, stacks with the damaged base itself in 62.18% of trajectory frames with FdU-substituted DNA and has no stacking interaction at all in 30.72% of these frames. Each drug investigated here induces a unique perturbation in the MutS α complex, indicating the possibility of a distinct signaling event and specific repair or death pathway (or set of pathways) for a given type of damage., Competing Interests: Conflict of Interest Statement: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Zinc-finger proteins are regulators of critical signaling pathways for various cellular functions, including apoptosis and oncogenesis. Here, we investigate how binding site protonation states and zinc coordination influence protein structure, dynamics, and ultimately function, as these pivotal regulatory proteins are increasingly important for protein engineering and therapeutic discovery. To better understand the thermodynamics and dynamics of the zinc finger of NEMO (NF-κB essential modulator), as well as the role of zinc, we present results of 20 μs molecular dynamics trajectories, 5 μs for each of four active site configurations. Consistent with experimental evidence, the zinc ion is essential for mechanical stabilization of the functional, folded conformation. Hydrogen bond motifs are unique for deprotonated configurations yet overlap in protonated cases. Correlated motions and principal component analysis corroborate the similarity of the protonated configurations and highlight unique relationships of the zinc-bound configuration. We hypothesize a potential mechanism for zinc binding from results of the thiol configurations. The deprotonated, zinc-bound configuration alone predominantly maintains its tertiary structure throughout all 5 μs and alludes rare conformations potentially important for (im)proper zinc-finger-related protein-protein or protein-DNA interactions.

As the length of molecular dynamics (MD) trajectories grows with increasing computational power, so does the importance of clustering methods for partitioning trajectories into conformational bins. Of the methods available, the vast majority require users to either have some a priori knowledge about the system to be clustered or to tune clustering parameters through trial and error. Here we present non-parametric uses of two modern clustering techniques suitable for first-pass investigation of an MD trajectory. Being non-parametric, these methods require neither prior knowledge nor parameter tuning. The first method, HDBSCAN, is fast-relative to other popular clustering methods-and is able to group unstructured or intrinsically disordered systems (such as intrinsically disordered proteins, or IDPs) into bins that represent global conformational shifts. HDBSCAN is also useful for determining the overall stability of a system-as it tends to group stable systems into one or two bins-and identifying transition events between metastable states. The second method, iMWK-Means, with explicit rescaling followed by K-Means, while slower than HDBSCAN, performs well with stable, structured systems such as folded proteins and is able to identify higher resolution details such as changes in relative position of secondary structural elements. Used in conjunction, these clustering methods allow a user to discern quickly and without prior knowledge the stability of a simulated system and identify both local and global conformational changes.

Advancing the use of therapeutic nucleic acids requires understanding the chemical and structural properties that allow these polymers to promote the death of malignant cells. Here we explore Zn 2+ complexation by the fluoropyrimidine polymer F10, which has strong activities in multiple preclinical models of cancer. Delivery of fluoropyrimidine FdUMP in the 10-residue polymer F10 rather than the nucleobase (5-fluorouracil) allows consideration of metal ion binding effects on drug delivery. The differences in metal ion interactions with fluoropyrimidine compared to normal DNA results in conformation changes that affect protein binding, cell uptake, and codelivery of metals such as zinc, and the cytoxicity thereof. Microsecond-time-scale, all-atom simulations of F10 predict that zinc selectively stabilizes the polymer via interactions with backbone phosphate groups and suggest a mechanism of complexation for the zinc-base interactions shown in previous experimental work. The positive zinc ions are attracted to the negatively charged phosphate groups. Once the Zn 2+ ions are near F10, they cause the base's N3 nitrogen to deprotonate. Subsequently, magnesium atoms displace zinc from their interactions with phosphate, freeing the zinc ions to interact with the FdU bases by forming weak interactions with the O4 oxygen and the fluorine attached to C5. These interactions of magnesium with phosphate groups and zinc with nucleobases agree with previous experimental results and are seen in MD simulations only when magnesium is introduced after N3 deprotonation, indicating a specific order of metal binding events. Additionally, we predict interactions between zinc and F10's O2 atoms, which were not previously observed. By comparison to 10mers of polyU and polydT, we also predict that the presence of fluorine increases the binding affinity of zinc to F10 relative to analogous strands of RNA and DNA consisting of only native nucleotides.

Displaying a single representative conformation of a biopolymer rather than an ensemble of states mistakenly conveys a static nature rather than the actual dynamic personality of biopolymers. However, there are few apparent options due to the fixed nature of print media. Here we suggest a standardized methodology for visually indicating the distribution width, standard deviation and uncertainty of ensembles of states with little loss of the visual simplicity of displaying a single representative conformation. Of particular note is that the visualization method employed clearly distinguishes between isotropic and anisotropic motion of polymer subunits. We also apply this method to ligand binding, suggesting a way to indicate the expected error in many high throughput docking programs when visualizing the structural spread of the output. We provide several examples in the context of nucleic acids and proteins with particular insights gained via this method. Such examples include investigating a therapeutic polymer of FdUMP (5-fluoro-2-deoxyuridine-5-O-monophosphate) - a topoisomerase-1 (Top1), apoptosis-inducing poison - and nucleotide-binding proteins responsible for ATP hydrolysis from Bacillus subtilis. We also discuss how these methods can be extended to any macromolecular data set with an underlying distribution, including experimental data such as NMR structures., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Molecular dynamics (MD) simulation methods have seen significant improvement since their inception in the late 1950s. Constraints of simulation size and duration that once impeded the field have lessened with the advent of better algorithms, faster processors, and parallel computing. With newer techniques and hardware available, MD simulations of more biologically relevant timescales can now sample a broader range of conformational and dynamical changes including rare events. One concern in the literature has been under which circumstances it is sufficient to perform many shorter timescale simulations and under which circumstances fewer longer simulations are necessary. Herein, our simulations of the zinc finger NEMO (2JVX) using multiple simulations of length 15, 30, 1000, and 3000 ns are analyzed to provide clarity on this point.

ERG, an ETS-family transcription factor, acts as a regulator of differentiation of early hematopoietic cells. It contains an autoinhibitory domain, which negatively regulates DNA-binding. The mechanism of autoinhibitory is still illusive. To understand the mechanism, we study the dynamical properties of ERG protein by molecular dynamics simulations. These simulations suggest that DNA binding autoinhibition associates with the internal dynamics of ERG. Specifically, we find that (1), The N-C terminal correlation in the inhibited ERG is larger than that in uninhibited ERG that contributes to the autoinhibition of DNA-binding. (2), DNA-binding changes the property of the N-C terminal correlation from being anti-correlated to correlated, that is, changing the relative direction of the correlated motions and (3), For the Ets-domain specifically, the inhibited and uninhibited forms exhibit essentially the same dynamics, but the binding of the DNA decreases the fluctuation of the Ets-domain. We also find from PCA analysis that the three systems, even with quite different dynamics, do have highly similar free energy surfaces, indicating that they share similar conformations.

Doxorubicin (Dox) is widely used for breast cancer treatment but causes serious side effects including cardiotoxicity that may adversely impact patient lifespan even if treatment is successful. Herein, we describe selective conjugation of Dox to a single site in a DNA hairpin resulting in a highly stable complex that enables Dox to be used more effectively. Selective conjugation of Dox to G15 in the hairpin loop was verified using site-specific labeling with [2-(15)N]-2'-deoxyguanosine in conjunction with [(1)H-(15)N] 2D NMR, while 1:1 stoichiometry for the conjugate was validated by ESI-QTOF mass spectrometry and UV spectroscopy. Molecular modeling indicated covalently bound Dox also intercalated into the stem of the hairpin and stability studies demonstrated the resulting Dox-conjugated hairpin (DCH) complex had a half-life >30 h, considerably longer than alternative covalent and noncovalent complexes. Secondary conjugation of DCH with folic acid (FA) resulted in increased internalization into breast cancer cells. The dual conjugate, DCH-FA, can be used for safer and more effective chemotherapy with Dox and this conjugation strategy can be expanded to include additional anticancer drugs.

MutSα is the most abundant mismatch-binding factor of human DNA mismatch repair (MMR) proteins. MMR maintains genetic stability by recognizing and repairing DNA defects. Failure to accomplish their function may lead to cancer. In addition, MutSα recognizes at least some types of DNA damage making it a target for anticancer agents. Here, complementing scarce experimental data, we report unique hydrogen-bonding motifs associated with the recognition of the carboplatin induced DNA damage by MutSα. These data predict that carboplatin and cisplatin induced damaging DNA adducts are recognized by MutSα in a similar manner. Our simulations also indicate that loss of base pairing at the damage site results in (1) non-specific binding and (2) changes in the atomic flexibility at the lesion site and beyond. To further quantify alterations at MutSα-DNA interface in response to damage recognition, non-bonding interactions and salt bridges were investigated. These data indicate (1) possible different packing and (2) disruption of the salt bridges at the MutSα-DNA interface in the damaged complex. These findings (1) underscore the general observation of disruptions at the MutSα-DNA interface and (2) highlight the nature of the anticancer effect of the carboplatin agent. The analysis was carried out from atomistic simulations.

DNA mismatch repair (MMR) proteins maintain genetic integrity in all organisms by recognizing and repairing DNA errors. Such alteration of hereditary information can lead to various diseases, including cancer. Besides their role in DNA repair, MMR proteins detect and initiate cellular responses to certain type of DNA damage. Its response to the damaged DNA has made the human MMR pathway a useful target for anticancer agents such as carboplatin. This study indicates that strong, specific interactions at the interface of MutSα in response to the mismatched DNA recognition are replaced by weak, non-specific interactions in response to the damaged DNA recognition. Data suggest a severe impairment of the dimerization of MutSα in response to the damaged DNA recognition. While the core of MutSα is preserved in response to the damaged DNA recognition, the loss of contact surface and the rearrangement of contacts at the protein interface suggest a different packing in response to the damaged DNA recognition. Coupled in response to the mismatched DNA recognition, interaction energies, hydrogen bonds, salt bridges, and solvent accessible surface areas at the interface of MutSα and within the subunits are uncoupled or asynchronously coupled in response to the damaged DNA recognition. These pieces of evidence suggest that the loss of a synchronous mode of response in the MutSα's surveillance for DNA errors would possibly be one of the mechanism(s) of signaling the MMR-dependent programed cell death much wanted in anticancer therapies. The analysis was drawn from dynamics simulations.

We, and others, have previously shown that mismatch repair proteins, in addition to their repair function, contribute to cell death initiation. In response to some drugs, this cell death activity is independent of the repair function of the proteins. Rescinnamine, a derivative of the indole alkaloid reserpine, a drug used to treat hypertension several decades ago, was shown to target the cell death-initiating activity of mismatch repair proteins. When used in animals, the hypotensive action of this drug prevents applying appropriate concentrations for statistically significant tumor reduction. Using a combination of computational modeling, chemical synthesis and cell assays, we determine how rescinnamine can be structurally modified and what effect these modifications have on cell survival. These results inform further computational modeling to suggest new synthetic lead molecules to move toward further biological testing.

The simultaneous binding of netropsin in the minor groove and Zn(2+) in the major groove of a DNA hairpin that includes 10 consecutive FdU nucleotides at the 3'-terminus (3'FdU) was demonstrated based upon NMR spectroscopy, circular dichroism (CD), and computational modeling studies. The resulting Zn(2+)/netropsin: 3'FdU complex had very high thermal stability with aspects of the complex intact at 85 °C, conditions that result in complete dissociation of Mg(2+) complexes. CD and (19)F NMR spectroscopy were consistent with Zn(2+) binding in the major groove of the DNA duplex and utilizing F5 and O4 of consecutive FdU nucleotides as ligands with FdU nucleotides hemi-deprotonated in the complex. Netropsin is bound in the minor groove of the DNA duplex based upon 2D NOESY data demonstrating contacts between AH2 (1)H and netropsin (1)H resonances. The Zn(2+)/netropsin: 3'FdU complex displayed increased cytotoxicity towards PC3 prostate cancer (PCa) cells relative to the constituent components or separate complexes (e.g. Zn(2+):3'FdU) indicating that this new structural motif may be therapeutically useful for PCa treatment. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:32.

One of the most popular and simple models for the calculation of pK(a) s from a protein structure is the semi-macroscopic electrostatic model MEAD. This model requires empirical parameters for each residue to calculate pK(a) s. Analysis of current, widely used empirical parameters for cysteine residues showed that they did not reproduce expected cysteine pK(a) s; thus, we set out to identify parameters consistent with the CHARMM27 force field that capture both the behavior of typical cysteines in proteins and the behavior of cysteines which have perturbed pK(a) s. The new parameters were validated in three ways: (1) calculation across a large set of typical cysteines in proteins (where the calculations are expected to reproduce expected ensemble behavior); (2) calculation across a set of perturbed cysteines in proteins (where the calculations are expected to reproduce the shifted ensemble behavior); and (3) comparison to experimentally determined pK(a) values (where the calculation should reproduce the pK(a) within experimental error). Both the general behavior of cysteines in proteins and the perturbed pK(a) in some proteins can be predicted reasonably well using the newly determined empirical parameters within the MEAD model for protein electrostatics. This study provides the first general analysis of the electrostatics of cysteines in proteins, with specific attention paid to capturing both the behavior of typical cysteines in a protein and the behavior of cysteines whose pK(a) should be shifted, and validation of force field parameters for cysteine residues., (Copyright © 2012 Wiley Periodicals, Inc.)

The peroxiredoxins (Prx) are ubiquitous peroxidases involved in important biological processes; however, details of their enzymatic mechanism remain elusive. To probe potential dynamics-function relationships, molecular dynamics simulations and electrostatic calculations were performed on the atypical 2-cysteine thiol peroxidase (Tpx) from Streptococcus pneumoniae and results compared to a previous study of a typical 2-cysteine Prx from Trypanosoma cruzi. The analyses indicate a commonality between both typical and atypical Prx: dynamic asymmetry. Asymmetry is observed in structure, fluctuations, and active site electrostatics. Key residues, including Glu150 and Phe153, play roles in the developing asymmetry; furthermore, in the atypical 2-Cys Tpx, Glu150 exhibits conformation fluctuations suggesting involvement in a proton shuttle. The existence of a pathway of connected residues appears to propagate the asymmetry. The commonality of asymmetry and coupling pathways in both typical and atypical Prxs suggests a driving force toward dimer asymmetry as a common feature that plays a functional role in creating one active site with a lower cysteine pK(a).

The cellular response to DNA damage signaling by mismatch-repair (MMR) proteins is incompletely understood. It is generally accepted that MMR-dependent apoptosis pathway in response to DNA damage detection is independent of MMR's DNA repair function. In this study, we investigate correlated motions in response to the binding of mismatched and platinum cross-linked DNA fragments by MutSα, as derived from 50 ns molecular dynamics simulations. The protein dynamics in response to the mismatched and damaged DNA recognition suggests that MutSα signals their recognition through independent pathways providing evidence for the molecular origin of the MMR-dependent apoptosis. MSH2 subunit is indicated to play a key role in signaling both mismatched and damaged DNA recognition; localized and collective motions within the protein allow identifying sites on the MSH2 surface possible involved in recruiting proteins responsible for downstream events. Unlike in the mismatch complex, predicted key communication sites specific for the damage recognition are on the list of known cancer-causing mutations or deletions. This confirms MSH2's role in signaling DNA damage-induced apoptosis and suggests that defects in MMR alone is sufficient to trigger tumorigenesis, supporting the experimental evidence that MMR-damage response function could protect from the early occurrence of tumors. Identifying these particular communication sites may have implications for the treatment of cancers that are not defective for MMR, but are unable to function optimally for MMR-dependent responses following DNA damage such as the case of resistance to cisplatin.

We report, based on semi-empirical calculations, that Zn(2+) binds duplex DNA containing consecutive FdU-dA base pairs in the major groove with distorted trigonal bipyramidal geometry. In this previously uncharacterized binding motif, O4 and F5 on consecutive FdU are axial ligands while three water molecules complete the coordination sphere. NMR spectroscopy confirmed Zn(2+) complexation occurred with maintenance of base pairing while a slight hypsochromic shift in circular dichroism (CD) spectra indicated moderate structural distortion relative to B-form DNA. Zn(2+) complexation inhibited ethidium bromide (EtBr) intercalation and stabilized FdU-substituted duplex DNA (ΔT(m) > 15 °C). Mg(2+) neither inhibited EtBr complexation nor had as strong of a stabilizing effect. DNA sequences that did not contain consecutive FdU were not stabilized by Zn(2+). A lipofectamine preparation of the Zn(2+)-DNA complex displayed enhanced cytotoxicity toward prostate cancer cells relative to the individual components prepared as lipofectamine complexes indicating the potential utility of Zn(2+)-DNA complexes for cancer treatment.

Human mismatch repair proteins MSH2-MSH6 play an essential role in maintaining genetic stability and preventing disease. While protein functions have been extensively studied, the substantial amino-terminal region (NTR*) of MSH6 that is unique to eukaryotic proteins, has mostly evaded functional characterization. We demonstrate that a cluster of three nuclear localization signals (NLS) in the NTR direct nuclear import. Individual NLSs are capable of partially directing cytoplasmic protein into the nucleus; however only cooperative effects between all three NLSs efficiently transport MSH6 into the nucleus. In striking contrast to yeast and previous assumptions on required heterodimerization, human MSH6 does not determine localization of its heterodimeric partner, MSH2. A cancer-derived mutation localized between two of the three NLS significantly decreases nuclear localization of MSH6, suggesting altered protein localization can contribute to carcinogenesis. These results clarify the pending speculations on the functional role of the NTR in human MSH6 and identify a novel, cooperative nuclear localization signal.

Molecular dynamics simulations have become increasingly useful in studying biological systems of biomedical interest, and not just in the study of model or toy systems. In this article, the methods and principles of all-atom molecular dynamics will be elucidated with several examples provided of their utility to investigators interested on drug discovery., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

MSH proteins are capable of recognizing damage in DNA due to a common chemotheraputic, cisplatin, and consequently participate in the initiation of cell death pathways. While previous studies have used computational modeling and cell biology to demonstrate that there are specific structural responses to cisplatin damage that are critical to the initiation of apoptosis, this study demonstrates that there are also specific dynamical changes that are also associated with cisplatin binding. These changes further distinguish the undamaged MutS/DNA complex from the damaged MutS/DNA complex and suggest that there are dynamical aspects to the response of MSH proteins to the binding of DNA damaged by therapeutics; consideration of these responses may influence further drug design and development.

The beta-lactam-based antibiotics are among the most prescribed and effective antibacterial agents. Widespread use of these antibiotics, however, has created tremendous pressure for the emergence of resistance mechanisms in bacteria. The most common cause of antibiotic resistance is bacterial production of actamases that efficiently degrade antibiotics. The metallo-beta-lactamases are of particular clinical concern due to their transference between bacterial strains. We used molecular dynamics (MD) simulations to further study the conformational changes that occur due to binding of an inhibitor to the dicanzinc metallo-beta-lactamase from Bacteroides fragilis. Our studies confirm previous findings that the major flap is a major source of plasticity within the active site, therefore its dynamic response should be considered in drug development. However, our results also suggest the need for care in using MD simulations in evaluating loop mobility, both due to relaxation times and to the need to accurately model the zinc active site. Our study also reveals two new robust responses to ligand binding. First, there are specific localized changes in the zinc active site--a local loop flip--due to ligand intercalation that may be critical to the function of this enzyme. Second, inhibitor binding perturbs the dynamics throughout the protein, without otherwise perturbing the enzyme structure. These dynamic perturbations radiate outward from the active site and their existence suggests that long-range communication and dynamics may be important in the activity of this enzyme.

Avoidance of apoptosis is one of the hallmarks of cancer development and progression. Chemotherapeutic agents aim to initiate an apoptotic response, but often fail due to dysregulation. MSH proteins are capable of recognizing cisplatin damage in DNA and participate in the initiation of cell death. We have exploited this recognition and computationally simulated a MutS homolog (MSH) "death conformation". Screening and docking experiments based on this model determined that the MSH2-dependent cell-death pathway can be induced by a small molecule without DNA damage, reserpine. Reserpine was identified via virtual screening on structures obtained from molecular dynamics as a small molecule that selectively binds a protein "death" conformation. The virtual screening predicts that this small molecule binds in the absence of DNA. Cell biology confirmed that reserpine triggers the MSH2-dependent cell-death pathway. This result supports the hypothesis that the MSH2-dependent pathway is initiated by specific protein conformational changes triggered by binding to either DNA damage or small compound molecules. These findings have multiple implications for drug discovery and cell biology. Computational modeling may be used to identify and eventually design small molecules that selectively activate particular pathways through conformational control. Molecular dynamics simulations can be used to model the biologically relevant conformations and virtual screening can then be used to select for small molecules that bind specific conformations. The ability of a small molecule to induce the cell-death pathway suggests a broader role for MMR proteins in cellular events, such as cell-death pathways, than previously suspected.

Cysteine sulfenic acid (Cys-SOH), a reversible modification, is a catalytic intermediate at enzyme active sites, a sensor for oxidative stress, a regulator of some transcription factors, and a redox-signaling intermediate. This post-translational modification is not random: specific features near the cysteine control its reactivity. To identify features responsible for the propensity of cysteines to be modified to sulfenic acid, a list of 47 proteins (containing 49 known Cys-SOH sites) was compiled. Modifiable cysteines are found in proteins from most structural classes and many functional classes, but have no propensity for any one type of protein secondary structure. To identify features affecting cysteine reactivity, these sites were analyzed using both functional site profiling and electrostatic analysis. Overall, the solvent exposure of modifiable cysteines is not different from the average cysteine. The combined sequence, structure, and electrostatic approaches reveal mechanistic determinants not obvious from overall sequence comparison, including: (1) pKaS of some modifiable cysteines are affected by backbone features only; (2) charged residues are underrepresented in the structure near modifiable sites; (3) threonine and other polar residues can exert a large influence on the cysteine pKa; and (4) hydrogen bonding patterns are suggested to be important. This compilation of Cys-SOH modification sites and their features provides a quantitative assessment of previous observations and a basis for further analysis and prediction of these sites. Agreement with known experimental data indicates the utility of this combined approach for identifying mechanistic determinants at protein functional sites.

CheY is a response regulator protein involved in bacterial chemotaxis. Much is known about its active and inactive conformations, but little is known about the mechanisms underlying long-range interactions or correlated motions. To investigate these events, molecular dynamics simulations were performed on the unphosphorylated, inactive structure from Salmonella typhimurium and the CheY-BeF(3)(-) active mimic structure (with BeF(3)(-) removed) from Escherichia coli. Simulations utilized both sequences in each conformation to discriminate sequence- and structure-specific behavior. The previously identified conformational differences between the inactive and active conformations of the strand-4-helix-4 loop, which are present in these simulations, arise from the structural, and not the sequence, differences. The simulations identify previously unreported structure-specific flexibility features in this loop and sequence-specific flexibility features in other regions of the protein. Both structure- and sequence-specific long-range interactions are observed in the active and inactive ensembles. In the inactive ensemble, two distinct mechanisms based on Thr-87 or Ile-95 rotameric forms, are observed for the previously identified g+ and g- rotamer sampling by Tyr-106. These molecular dynamics simulations have thus identified both sequence- and structure-specific differences in flexibility, long-range interactions, and rotameric form of key residues. Potential biological consequences of differential flexibility and long-range correlated motion are discussed.

We determined the molecular mechanism of cell death response by MutS homologs in distinction to the repair event. Key protein-DNA contacts differ in the interaction of MutS homologs with cisplatinated versus mismatched DNA. Mutational analyses of protein-DNA contacts, which were predicted by molecular dynamics (MD) simulations, were performed. Mutations in suggested interaction sites can affect repair and cell death response independently, and to different extents. A glutamate residue is identified as the key contact with cisplatin-DNA. Mutation of the residue increases cisplatin resistance due to increased non-specific DNA binding. In contrast, the conserved phenylalanine that is instrumental and indispensable for mismatch recognition during repair is not required for cisplatin cytotoxicity. These differences in protein-DNA interactions are translated into localized conformational changes that affect nucleotide requirements and inter-subunit interactions. Specifically, the ability for ATP binding/hydrolysis has little consequence for the MMR-dependent damage response. As a consequence, intersubunit contacts are altered that most likely affect the interaction with downstream proteins. We here describe the interaction of MutS homologs with DNA damage, as it differs from the interaction with a mismatch, and its structural translation into all other functional regions of the protein as a mechanism to initiate cell death response and concomitantly inhibit repair.

We present a new hybrid explicit/implicit solvent method for dynamics simulations of macromolecular systems. The method models explicitly the hydration of the solute by either a layer or sphere of water molecules, and the generalized Born (GB) theory is used to treat the bulk continuum solvent outside the explicit simulation volume. To reduce the computational cost, we implemented a multigrid method for evaluating the pairwise electrostatic and GB terms. It is shown that for typical ion and protein simulations our method achieves similar equilibrium and dynamical observables as the conventional particle mesh Ewald (PME) method. Simulation timings are reported, which indicate that the hybrid method is much faster than PME, primarily due to a significant reduction in the number of explicit water molecules required to model hydration effects., ((c) 2004 Wiley Periodicals, Inc.)

In this work, we explore the question of whether pK(a) calculations based on a microscopic description of the protein and a macroscopic description of the solvent can be implemented to examine conformationally dependent proton shifts in proteins. To this end, we introduce a new method for performing constant-pH molecular dynamics (PHMD) simulations utilizing the generalized Born implicit solvent model. This approach employs an extended Hamiltonian in which continuous titration coordinates propagate simultaneously with the atomic motions of the system. The values adopted by these coordinates are modulated by potentials of mean force of isolated titratable model groups and the pH to control the proton occupation at particular sites in the polypeptide. Our results for four different proteins yield an absolute average error of approximately 1.6 pK units, and point to the role that thermally driven relaxation of the protein environment in the vicinity of titrating groups plays in modulating the local pK(a), thereby influencing the observed pK1/2 values. While the accuracy of our method is not yet equivalent to methods that obtain pK1/2 values through the ad hoc scaling of electrostatics, the present approach and constant pH methods in general provide a useful framework for studying pH-dependent phenomena. Further work to improve our model to approach quantitative agreement with experiment is outlined., (Copyright 2004 Wiley-Liss, Inc.)

The unusual temperature-dependent excited-state dynamics of a stilbene-antibody complex reported by Simeonov et al. are explored using theoretical methods. The anomalous temperature-dependent fluorescence emission and lifetime are shown to be the result of interplay among temperature-modulated protein flexibility, the excited-state potential surface for the stilbene central twist, and changes in the stilbene charge distribution upon excitation. Stilbene is found to possess a similar geometry and orientation within the antibody binding pocket at all temperatures in the ground state and at low temperatures (approximately 240 K) in the excited state. At higher temperatures (approximately 260 K), the excited-state conformation twists around the central double bond and adopts an alternate orientation within the binding pocket. These changes result from protein side chain and loop motions that are frozen out at lower temperatures and account for the observed red shift of the fluorescence emission spectrum (a calculated shift of 3.8 kcal mol-1 compares favorably with the approximately 5 kcal.mol-1 observed experimentally). Approximately 3.0 kcal mol-1 of this stabilization is global in nature and is not attributable to specific local interactions. Local interactions between stilbene and Tyr-B39 contribute approximately 0.8 kcal mol-1 to the fluorescence shift. The primary structural change in simulations of the high temperature excited state involves a decrease in the stilbene-tyrosine distance and a relative change in orientation of the aromatic rings. We identify several nearby charged residues that contribute to the fluorescence emission shift and provide targets for mutagenesis to probe the temperature-dependent dynamics of protein-chromophore interactions.

In a recent article (Lee, M. S.; Salsbury, F. R. Jr.; Brooks, C. L., III. J Chem Phys 2002, 116, 10606), we demonstrated that generalized Born (GB) theory provides a good approximation to Poisson electrostatic solvation energy calculations if one uses the same definitions of molecular volume for each. In this work, we present a new and improved analytic method for reproducing the Lee-Richards molecular volume, which is the most common volume definition for Poisson calculations. Overall, 1% errors are achieved for absolute solvation energies of a large set of proteins and relative solvation energies of protein conformations. We also introduce an accurate SASA approximation that uses the same machinery employed by our GB method and requires a small addition of computational cost. The combined methodology is shown to yield an efficient and accurate implicit solvent representation for simulations of biopolymers., (Copyright 2003 Wiley Periodicals, Inc.)

The structure and dynamics of an inhibitor-bound complex of the metallo-beta-lactamase from Bacteroides fragilis are studied by using molecular dynamics. A search of the conformational space was performed to obtain three distinct models of the complex, which were then subjected to solvated molecular dynamics. A solvated molecular dynamics study of the apo protein was performed to serve as a baseline for comparison with the bound simulations. We find loop conformation changes due to binding as well as a decrease in flexibility of the protein as a whole and especially in the major loop of the beta-lactamase. We report the structural and dynamical features of the inhibitor-bound and apo models, as well as experimentally measurable quantities, which should be capable of distinguishing the two binding modes we have determined., (Copyright 2001 Wiley-Liss, Inc.)