ABSTRACT
The molecular mechanisms governing the human voltage-gated proton channel hHv1 remain elusive. Here, we used membrane-enabled hybrid-solvent continuous constant pH molecular dynamics (CpHMD) simulations with pH replica exchange to further evaluate the structural models of hHv1 in the closed (hyperpolarized) and open (depolarized) states recently obtained with MD simulations and explore potential pH-sensing residues. The CpHMD titration at a set of symmetric pH conditions revealed three residues that can gain or lose protons upon channel depolarization. Among them, residue H168 at the intracellular end of the S3 helix switches from the deprotonated to the protonated state and its protonation is correlated with the increased tilting of the S3 helix during the transition from the closed to the open state. Thus, the simulation data suggest H168 as an interior pH sensor, in support of a recent finding based on electrophysiological experiments of Hv1 mutants. We propose that protonation of H168 acts as a key that unlocks the closed channel configuration by increasing the flexibility of the S2-S3 linker, which increases the tilt angle of S3 and enhances the mobility of the S4 helix, thus promoting channel opening. Our work represents an important step toward deciphering the pH-dependent gating mechanism of hHv1.
ABSTRACT
Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid-base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.
Subject(s)
Software , Molecular Dynamics Simulation , Genetic Variation , Algorithms , Thermodynamics , Proteins/chemistry , Crystallization , Nucleic Acids/chemistryABSTRACT
Chitosan-based materials have broad applications, from biotechnology to pharmaceutics. Recent experiments showed that the degree and pattern of acetylation along the chitosan chain modulate its biological and physicochemical properties; however, the molecular mechanism remains unknown. Here, we report, to the best of our knowledge, the first de novo all-atom molecular dynamics (MD) simulations to investigate chitosan's self-assembly process at different degrees and patterns of acetylation. Simulations revealed that 10 mer chitosan chains with 50% acetylation in either block or alternating patterns associate to form ordered nanofibrils comprised of mainly antiparallel chains in agreement with the fiber diffraction data of deacetylated chitosan. Surprisingly, regardless of the acetylation pattern, the same intermolecular hydrogen bonds mediate fibril sheet formation while water-mediated interactions stabilize sheet-sheet stacking. Moreover, acetylated units are involved in forming strong intermolecular hydrogen bonds (NH-O6 and O6H-O7), which offers an explanation for the experimental observation that increased acetylation lowers chitosan's solubility. Taken together, the present study provides atomic-level understanding the role of acetylation plays in modulating chitosan's physiochemical properties, contributing to the rational design of chitosan-based materials with the ability to tune by its degree and pattern of acetylation. Additionally, we disseminate the improved molecular mechanics parameters that can be applied in MD studies to further understand chitosan-based materials.
Subject(s)
Chitosan , Chitosan/chemistry , Acetylation , Molecular Dynamics SimulationABSTRACT
Twenty years ago, this journal published a review entitled "Biofabrication with Chitosan" based on the observations that (i) chitosan could be electrodeposited using low voltage electrical inputs (typically less than 5 V) and (ii) the enzyme tyrosinase could be used to graft proteins (via accessible tyrosine residues) to chitosan. Here, we provide a progress report on the coupling of electronic inputs with advanced biological methods for the fabrication of biopolymer-based hydrogel films. In many cases, the initial observations of chitosan's electrodeposition have been extended and generalized: mechanisms have been established for the electrodeposition of various other biological polymers (proteins and polysaccharides), and electrodeposition has been shown to allow the precise control of the hydrogel's emergent microstructure. In addition, the use of biotechnological methods to confer function has been extended from tyrosinase conjugation to the use of protein engineering to create genetically fused assembly tags (short sequences of accessible amino acid residues) that facilitate the attachment of function-conferring proteins to electrodeposited films using alternative enzymes (e.g., transglutaminase), metal chelation, and electrochemically induced oxidative mechanisms. Over these 20 years, the contributions from numerous groups have also identified exciting opportunities. First, electrochemistry provides unique capabilities to impose chemical and electrical cues that can induce assembly while controlling the emergent microstructure. Second, it is clear that the detailed mechanisms of biopolymer self-assembly (i.e., chitosan gel formation) are far more complex than anticipated, and this provides a rich opportunity both for fundamental inquiry and for the creation of high performance and sustainable material systems. Third, the mild conditions used for electrodeposition allow cells to be co-deposited for the fabrication of living materials. Finally, the applications have been expanded from biosensing and lab-on-a-chip systems to bioelectronic and medical materials. We suggest that electro-biofabrication is poised to emerge as an enabling additive manufacturing method especially suited for life science applications and to bridge communication between our biological and technological worlds.
Subject(s)
Chitosan , Chitosan/chemistry , Monophenol Monooxygenase/chemistry , Hydrogels , Proteins , BiopolymersABSTRACT
The ERK pathway is one of the most important signaling cascades involved in tumorigenesis. So far, eight noncovalent inhibitors of RAF and MEK kinases in the ERK pathway have been approved by the FDA for the treatment of cancers; however, their efficacies are limited due to various resistance mechanisms. There is an urgent need to develop novel targeted covalent inhibitors. Here we report a systematic study of the covalent ligandabilities of the ERK pathway kinases (ARAF, BRAF, CRAF, KSR1, KSR2, MEK1, MEK2, ERK1, and ERK2) using constant pH molecular dynamics titration and pocket analysis. Our data revealed that the hinge GK (gate keeper)+3 cysteine in RAF family kinases (ARAF, BRAF, CRAF, KSR1, and KSR2) and the back loop cysteine in MEK1 and MEK2 are reactive and ligandable. Structure analysis suggests that the type II inhibitors belvarafenib and GW5074 may be used as scaffolds for designing pan-RAF or CRAF-selective covalent inhibitors directed at the GK+3 cysteine, while the type III inhibitor cobimetinib may be modified to label the back loop cysteine in MEK1/2. The reactivities and ligandabilities of the remote cysteine in MEK1/2 and the DFG-1 cysteine in MEK1/2 and ERK1/2 are also discussed. Our work provides a starting point for medicinal chemists to design novel covalent inhibitors of the ERK pathway kinases. The computational protocol is general and can be applied to the systematic evaluation of covalent ligandabilities of the human cysteinome.
Subject(s)
MAP Kinase Kinase Kinases , MAP Kinase Signaling System , Humans , MAP Kinase Signaling System/physiology , MAP Kinase Kinase Kinases/metabolism , Proto-Oncogene Proteins B-raf/chemistry , Proto-Oncogene Proteins B-raf/metabolism , Cysteine/metabolism , Signal Transduction , raf Kinases/metabolismABSTRACT
Predictive modeling and understanding of chemical warhead reactivities have the potential to accelerate targeted covalent drug discovery. Recently, the carbanion formation free energies as well as other ground-state electronic properties from density functional theory (DFT) calculations have been proposed as predictors of glutathione reactivities of Michael acceptors; however, no clear consensus exists. By profiling the thiol-Michael reactions of a diverse set of singly- and doubly-activated olefins, including several model warheads related to afatinib, here we reexamined the question of whether low-cost electronic properties can be used as predictors of reaction barriers. The electronic properties related to the carbanion intermediate were found to be strong predictors, e.g., the change in the Cß charge accompanying carbanion formation. The least expensive reactant-only properties, the electrophilicity index, and the Cß charge also show strong rank correlations, suggesting their utility as quantum descriptors. A second objective of the work is to clarify the effect of the ß-dimethylaminomethyl (DMAM) substitution, which is incorporated in the warheads of several FDA-approved covalent drugs. Our data suggest that the ß-DMAM substitution is cationic at neutral pH in solution and promotes acrylamide's intrinsic reactivity by enhancing the charge accumulation at Cα upon carbanion formation. In contrast, the inductive effect of the ß-trimethylaminomethyl substitution is diminished due to steric hindrance. Together, these results reconcile the current views of the intrinsic reactivities of acrylamides and contribute to large-scale predictive modeling and an understanding of the structure-activity relationships of Michael acceptors for rational TCI design.
Subject(s)
Drug Discovery , Sulfhydryl Compounds , Structure-Activity Relationship , Afatinib , Glutathione/chemistryABSTRACT
The nation's opioid overdose deaths reached an all-time high in 2021. The majority of deaths are due to synthetic opioids represented by fentanyl. Naloxone, which is a FDA-approved reversal agent, antagonizes opioids through competitive binding at the µ-opioid receptor (mOR). Thus, knowledge of the opioid's residence time is important for assessing the effectiveness of naloxone. Here, we estimated the residence times (τ) of 15 fentanyl and 4 morphine analogs using metadynamics and compared them with the most recent measurement of the opioid kinetic, dissociation, and naloxone inhibitory constants (Mann et al. Clin. Pharmacol. Therapeut. 2022, 120, 1020-1232). Importantly, the microscopic simulations offered a glimpse at the common binding mechanism and molecular determinants of dissociation kinetics for fentanyl analogs. The insights inspired us to develop a machine learning approach to analyze the kinetic impact of fentanyl's substituents based on the interactions with mOR residues. This proof-of-concept approach is general; for example, it may be used to tune ligand residence times in computer-aided drug discovery.
Subject(s)
Analgesics, Opioid , Naloxone , Analgesics, Opioid/pharmacology , Naloxone/pharmacology , Naloxone/metabolism , Fentanyl/metabolism , Fentanyl/pharmacology , Morphine/chemistry , Receptors, Opioid, mu/metabolism , Narcotic AntagonistsABSTRACT
Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of Paxlovid, a drug approved by the U.S. Food and Drug Administration for high-risk COVID-19 patients. Recently, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir's inhibitory activity. As the COVID-19 cases skyrocket in China and the selective pressure of antiviral therapy builds in the US, there is an urgent need to characterize and understand how the H172Y mutation confers drug resistance. Here, we investigated the H172Y Mpro's conformational dynamics, folding stability, catalytic efficiency, and inhibitory activity using all-atom constant pH and fixed-charge molecular dynamics simulations, alchemical and empirical free energy calculations, artificial neural networks, and biochemical experiments. Our data suggest that the mutation significantly weakens the S1 pocket interactions with the N-terminus and perturbs the conformation of the oxyanion loop, leading to a decrease in the thermal stability and catalytic efficiency. Importantly, the perturbed S1 pocket dynamics weaken the nirmatrelvir binding in the P1 position, which explains the decreased inhibitory activity of nirmatrelvir. Our work demonstrates the predictive power of the combined simulation and artificial intelligence approaches, and together with biochemical experiments, they can be used to actively surveil continually emerging mutations of SARS-CoV-2 Mpro and assist the optimization of antiviral drugs. The presented approach, in general, can be applied to characterize mutation effects on any protein drug targets.
Subject(s)
COVID-19 , Humans , SARS-CoV-2/genetics , SARS-CoV-2/metabolism , Artificial Intelligence , Protease Inhibitors/chemistry , Antiviral Agents/chemistry , Molecular Dynamics Simulation , Mutation , Drug Resistance , Molecular Docking SimulationABSTRACT
Escherichia coli NhaA is a prototypical sodium-proton antiporter responsible for maintaining cellular ion and volume homeostasis by exchanging two protons for one sodium ion; despite two decades of research, the transport mechanism of NhaA remains poorly understood. Recent crystal structure and computational studies suggested Lys300 as a second proton-binding site; however, functional measurements of several K300 mutants demonstrated electrogenic transport, thereby casting doubt on the role of Lys300. To address the controversy, we carried out state-of-the-art continuous constant pH molecular dynamics simulations of NhaA mutants K300A, K300R, K300Q/D163N, and K300Q/D163N/D133A. Simulations suggested that K300 mutants maintain the electrogenic transport by utilizing an alternative proton-binding residue Asp133. Surprisingly, while Asp133 is solely responsible for binding the second proton in K300R, Asp133 and Asp163 jointly bind the second proton in K300A, and Asp133 and Asp164 jointly bind two protons in K300Q/D163N. Intriguingly, the coupling between Asp133 and Asp163 or Asp164 is enabled through the proton-coupled hydrogen-bonding network at the flexible intersection of two disrupted helices. These data resolve the controversy and highlight the intricacy of the compensatory transport mechanism of NhaA mutants. Alternative proton-binding site and proton sharing between distant aspartates may represent important general mechanisms of proton-coupled transport in secondary active transporters.
Subject(s)
Escherichia coli Proteins , Protons , Sodium-Hydrogen Exchangers , Aspartic Acid/chemistry , Aspartic Acid/genetics , Aspartic Acid/metabolism , Binding Sites , Escherichia coli Proteins/chemistry , Escherichia coli Proteins/genetics , Escherichia coli Proteins/metabolism , Hydrogen Bonding , Lysine/chemistry , Lysine/genetics , Lysine/metabolism , Molecular Dynamics Simulation , Mutation , Sodium-Hydrogen Exchangers/chemistry , Sodium-Hydrogen Exchangers/genetics , Sodium-Hydrogen Exchangers/metabolism , Static ElectricityABSTRACT
Malaria remains a global health threatâover 400,000 deaths occurred in 2019. Plasmepsins are promising targets of antimalarial therapeutics; however, no inhibitors have reached the clinic. To fuel the progress, a detailed understanding of the pH- and ligand-dependent conformational dynamics of plasmepsins is needed. Here we present the continuous constant pH molecular dynamics study of the prototypical plasmepsin II and its complexed form with a substrate analogue. The simulations revealed that the catalytic dyads D34 and D214 are highly coupled in the apo protein and that the pepstatin binding enhances the difference in proton affinity, making D34 the general base and D214 the general acid. The simulations showed that the flap adopts an open state regardless of pH; however, upon pepstatin binding the flap can close or open depending on the protonation state of D214. These and other data are discussed and compared with the off-targets human cathepsin D and renin. This study lays the groundwork for a systematic investigation of pH- and ligand-modulated dynamics of the entire family of plasmepsins to help design more potent and selective inhibitors.
Subject(s)
Aspartic Acid Endopeptidases , Malaria , Aspartic Acid Endopeptidases/chemistry , Humans , Hydrogen-Ion Concentration , Ligands , Malaria/drug therapy , Protein Conformation , Protozoan ProteinsABSTRACT
Renin is a pepsin-like aspartyl protease and an important drug target for the treatment of hypertension; despite three decades' research, its pH-dependent structure-function relationship remains poorly understood. Here, we employed continuous constant pH molecular dynamics (CpHMD) simulations to decipher the acid/base roles of renin's catalytic dyad and the conformational dynamics of the flap, which is a common structural feature among aspartyl proteases. The calculated pKa's suggest that catalytic Asp38 and Asp226 serve as the general base and acid, respectively, in agreement with experiment and supporting the hypothesis that renin's neutral optimum pH is due to the substrate-induced pKa shifts of the aspartic dyad. The CpHMD data confirmed our previous hypothesis that hydrogen bond formation is the major determinant of the dyad pKa order. Additionally, our simulations showed that renin's flap remains open regardless of pH, although a Tyr-inhibited state is occasionally formed above pH 5. These findings are discussed in comparison to the related aspartyl proteases, including ß-secretases 1 and 2, cathepsin D, and plasmepsin II. Our work represents a first step toward a systematic understanding of the pH-dependent structure-dynamics-function relationships of pepsin-like aspartyl proteases that play important roles in biology and human disease states.
Subject(s)
Molecular Dynamics Simulation , Renin , Catalysis , Humans , Hydrogen-Ion ConcentrationABSTRACT
The SARS coronavirus 2 (SARS-CoV-2) main protease (Mpro) is an attractive broad-spectrum antiviral drug target. Despite the enormous progress in structure elucidation, the Mpro's structure-function relationship remains poorly understood. Recently, a peptidomimetic inhibitor has entered clinical trial; however, small-molecule orally available antiviral drugs have yet to be developed. Intrigued by a long-standing controversy regarding the existence of an inactive state, we explored the proton-coupled dynamics of the Mpros of SARS-CoV-2 and the closely related SARS-CoV using a newly developed continuous constant pH molecular dynamics (MD) method and microsecond fixed-charge all-atom MD simulations. Our data supports a general base mechanism for Mpro's proteolytic function. The simulations revealed that protonation of His172 alters a conserved interaction network that upholds the oxyanion loop, leading to a partial collapse of the conserved S1 pocket, consistent with the first and controversial crystal structure of SARS-CoV Mpro determined at pH 6. Interestingly, a natural flavonoid binds SARS-CoV-2 Mpro in the close proximity to a conserved cysteine (Cys44), which is hyper-reactive according to the CpHMD titration. This finding offers an exciting new opportunity for small-molecule targeted covalent inhibitor design. Our work represents a first step toward the mechanistic understanding of the proton-coupled structure-dynamics-function relationship of CoV Mpros; the proposed strategy of designing small-molecule covalent inhibitors may help accelerate the development of orally available broad-spectrum antiviral drugs to stop the current pandemic and prevent future outbreaks.
Subject(s)
Antiviral Agents/chemistry , Antiviral Agents/pharmacology , Coronavirus 3C Proteases/chemistry , Coronavirus 3C Proteases/drug effects , SARS-CoV-2/drug effects , SARS-CoV-2/enzymology , Binding Sites , Cysteine/chemistry , Drug Design , Humans , Hydrogen-Ion Concentration , Models, Molecular , Molecular Dynamics Simulation , Protease Inhibitors/pharmacology , Protein Conformation , Protons , Small Molecule Libraries , Structure-Activity RelationshipABSTRACT
Broad-spectrum antiviral drugs are urgently needed to stop the Coronavirus Disease 2019 pandemic and prevent future ones. The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is related to the SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV), which have caused the previous outbreaks. The papain-like protease (PLpro) is an attractive drug target due to its essential roles in the viral life cycle. As a cysteine protease, PLpro is rich in cysteines and histidines, and their protonation/deprotonation modulates catalysis and conformational plasticity. Here, we report the pKa calculations and assessment of the proton-coupled conformational dynamics of SARS-CoV-2 in comparison to SARS-CoV and MERS-CoV PLpros using the recently developed graphical processing unit (GPU)-accelerated implicit-solvent continuous constant pH molecular dynamics method with a new asynchronous replica-exchange scheme, which allows computation on a single GPU card. The calculated pKa's support the catalytic roles of the Cys-His-Asp triad. We also found that several residues can switch protonation states at physiological pH among which is C270/271 located on the flexible blocking loop 2 (BL2) of SARS-CoV-2/CoV PLpro. Simulations revealed that the BL2 can open and close depending on the protonation state of C271/270, consistent with the most recent crystal structure evidence. Interestingly, despite the lack of an analogous cysteine, BL2 in MERS-CoV PLpro is also very flexible, challenging a current hypothesis. These findings are supported by the all-atom fixed-charge simulations and provide a starting point for more detailed studies to assist the structure-based design of broad-spectrum inhibitors against CoV PLpros.
Subject(s)
Antiviral Agents/pharmacology , Betacoronavirus/enzymology , Drug Design , Middle East Respiratory Syndrome Coronavirus/enzymology , Molecular Dynamics Simulation , Papain/chemistry , Papain/metabolism , Protons , Amino Acid Sequence , Histidine , Hydrogen-Ion Concentration , Papain/antagonists & inhibitors , Protein Domains , SARS-CoV-2ABSTRACT
A heme-dependent conformational rearrangement of the C-terminal domain of heme binding protein (PhuS) is required for interaction with the iron-regulated heme oxygenase (HemO). Herein, we further investigate the underlying mechanism of this conformational rearrangement and its implications for heme transfer via site-directed mutagenesis, resonance Raman (RR), hydrogen-deuterium exchange MS (HDX-MS) methods, and molecular dynamics (MD). HDX-MS revealed that the apo-PhuS C-terminal α6/α7/α8-helices are largely unstructured, whereas the apo-PhuS H212R variant showed an increase in structure within these regions. The increased rate of heme association with apo-PhuS H212R compared with the WT and lack of a detectable five-coordinate high-spin (5cHS) heme intermediate are consistent with a more folded and less dynamic C-terminal domain. HDX-MS and MD of holo-PhuS indicate an overall reduction in molecular flexibility throughout the protein, with significant structural rearrangement and protection of the heme binding pocket. We observed slow cooperative unfolding/folding events within the C-terminal helices of holo-PhuS and the N-terminal α1/α2-helices that are dampened or eliminated in the holo-PhuS H212R variant. Chemical cross-linking and MALDI-TOF MS mapped these same regions to the PhuS:HemO protein-protein interface. We previously proposed that the protein-protein interaction induces conformational rearrangement, promoting a ligand switch from His-209 to His-212 and triggering heme release to HemO. The reduced conformational freedom of holo-PhuS H212R combined with the increase in entropy and decrease in heme transfer on interaction with HemO further support this model. This study provides significant insight into the role of protein dynamics in heme binding and release in bacterial heme transport proteins.
Subject(s)
Bacterial Proteins/chemistry , Bacterial Proteins/metabolism , Carrier Proteins/chemistry , Carrier Proteins/metabolism , Heme Oxygenase (Decyclizing)/metabolism , Hemeproteins/chemistry , Hemeproteins/metabolism , Pseudomonas aeruginosa/metabolism , Allosteric Regulation , Bacterial Proteins/genetics , Carrier Proteins/genetics , Heme Oxygenase (Decyclizing)/chemistry , Heme Oxygenase (Decyclizing)/genetics , Heme-Binding Proteins , Hemeproteins/genetics , Ligands , Protein Binding , Protein Structure, Secondary , Pseudomonas aeruginosa/chemistry , Pseudomonas aeruginosa/enzymology , Pseudomonas aeruginosa/geneticsABSTRACT
Protein kinases are important cellular signaling molecules involved in cancer and a multitude of other diseases. It is well-known that inactive kinases display a remarkable conformational plasticity; however, the molecular mechanisms remain poorly understood. Conformational heterogeneity presents an opportunity but also a challenge in kinase drug discovery. The ability to predictively model various conformational states could accelerate selective inhibitor design. Here we performed a proton-coupled molecular dynamics study to explore the conformational landscape of a c-Src kinase. Starting from a completely inactive structure, the simulations captured all major types of conformational states without the use of a target structure, mutation, or bias. The simulations allowed us to test the experimental hypotheses regarding the mechanism of DFG flip, its coupling to the αC-helix movement, and the formation of regulatory spine. Perhaps the most significant finding is how key titratable residues, such as DFG-Asp, αC-Glu, and HRD-Asp, change protonation states dependent on the DFG, αC, and activation loop conformations. Our data offer direct evidence to support a long-standing hypothesis that protonation of Asp favors the DFG-out state and explain why DFG flip is also possible in simulations with deprotonated Asp. The simulations also revealed intermediate states, among which a unique DFG-out/α-C state formed as DFG-Asp is moved into a back pocket forming a salt bridge with catalytic Lys, which can be tested in selective inhibitor design. Our finding of how proton coupling enables the remarkable conformational plasticity may shift the paradigm of computational studies of kinases which assume fixed protonation states. Understanding proton-coupled conformational dynamics may hold a key to further innovation in kinase drug discovery.
Subject(s)
Protein-Tyrosine Kinases/chemistry , Humans , Hydrogen-Ion Concentration , Models, Molecular , Protein Conformation , Protein-Tyrosine Kinases/metabolism , Static ElectricityABSTRACT
Targeted covalent inhibitor design is gaining increasing interest and acceptance. A typical covalent kinase inhibitor design targets a reactive cysteine; however, this strategy is limited by the low abundance of cysteine and acquired drug resistance from point mutations. Inspired by the recent development of lysine-targeted chemical probes, we asked if nucleophilic (reactive) catalytic lysines are common on the basis of the published crystal structures of the human kinome. Using a newly developed p Ka prediction tool based on continuous constant pH molecular dynamics, the catalytic lysines of eight unique kinases from various human kinase groups were retrospectively and prospectively predicted to be nucleophilic, when kinase is in the rare DFG-out/αC-out type of conformation. Importantly, other reactive lysines as well as cysteines at various locations were also identified. On the basis of the findings, we proposed a new strategy in which selective type II reversible kinase inhibitors are modified to design highly selective, lysine-targeted covalent inhibitors. Traditional covalent drugs were discovered serendipitously; the presented tool, which can assess the reactivities of any potentially targetable residues, may accelerate the rational discovery of new covalent inhibitors. Another significant finding of the work is that lysines and cysteines in kinases may adopt neutral and charged states at physiological pH, respectively. This finding may shift the current paradigm of computational studies of kinases, which assume fixed solution protonation states.
Subject(s)
Computational Biology , Cysteine/metabolism , Drug Design , Lysine/metabolism , Protein Kinase Inhibitors/pharmacology , Protein Kinases/chemistry , Protein Kinases/metabolism , Biocatalysis , Humans , Molecular Dynamics Simulation , Molecular Targeted Therapy , Protein ConformationABSTRACT
We present a GPU implementation of the continuous constant pH molecular dynamics (CpHMD) based on the most recent generalized Born implicit-solvent model in the pmemd engine of the Amber molecular dynamics package. To test the accuracy of the tool for rapid pKa predictions, a series of 2 ns single-pH simulations were performed for over 120 titratable residues in 10 benchmark proteins that were previously used to test the various continuous CpHMD methods. The calculated pKa's showed a root-mean-square deviation of 0.80 and correlation coefficient of 0.83 with respect to experiment. Also, 90% of the pKa's were converged with estimated errors below 0.1 pH units. Surprisingly, this level of accuracy is similar to our previous replica-exchange simulations with 2 ns per replica and an exchange attempt frequency of 2 ps-1 (Huang, Harris, and Shen J. Chem. Inf. Model. 2018 , 58 , 1372 - 1383 ). Interestingly, for the linked titration sites in two enzymes, although residue-specific protonation state sampling in the single-pH simulations was not converged within 2 ns, the protonation fraction of the linked residues appeared to be largely converged, and the experimental macroscopic pKa values were reproduced to within 1 pH unit. Comparison with replica-exchange simulations with different exchange attempt frequencies showed that the splitting between the two macroscopic pKa's is underestimated with frequent exchange attempts such as 2 ps-1, while single-pH simulations overestimate the splitting. The same trend is seen for the single-pH vs replica-exchange simulations of a hydrogen-bonded aspartyl dyad in a much larger protein. A 2 ns single-pH simulation of a 400-residue protein takes about 1 h on a single NVIDIA GeForce RTX 2080 graphics card, which is over 1000 times faster than a CpHMD run on a single CPU core of a high-performance computing cluster node. Thus, we envision that GPU-accelerated continuous CpHMD may be used in routine pKa predictions for a variety of applications, from assisting MD simulations with protonation state assignment to offering pH-dependent corrections of binding free energies and identifying reactive hot spots for covalent drug design.
Subject(s)
Molecular Dynamics Simulation , Proteins/chemistry , Animals , Bacterial Proteins/chemistry , Humans , Hydrogen Bonding , Hydrogen-Ion Concentration , Protein Conformation , Protons , Software , Solvents/chemistry , ThermodynamicsABSTRACT
Biology is well-known for its ability to communicate through (i) molecularly-specific signaling modalities and (ii) a globally-acting electrical modality associated with ion flow across biological membranes. Emerging research suggests that biology uses a third type of communication modality associated with a flow of electrons through reduction/oxidation (redox) reactions. This redox signaling modality appears to act globally and has features of both molecular and electrical modalities: since free electrons do not exist in aqueous solution, the electrons must flow through molecular intermediates that can be switched between two states - with electrons (reduced) or without electrons (oxidized). Importantly, this global redox modality is easily accessible through its electrical features using convenient electrochemical instrumentation. In this review, we explain this redox modality, describe our electrochemical measurements, and provide four examples demonstrating that redox enables communication between biology and electronics. The first two examples illustrate how redox probing can acquire biologically relevant information. The last two examples illustrate how redox inputs can transduce biologically-relevant transitions for patterning and the induction of a synbio transceiver for two-hop molecular communication. In summary, we believe redox provides a unique ability to bridge bio-device communication because simple electrochemical methods enable global access to biologically meaningful information. Further, we envision that redox may facilitate the application of information theory to the biological sciences.
ABSTRACT
The growing importance of hydrogels in translational medicine has stimulated the development of top-down fabrication methods, yet often these methods lack the capabilities to generate the complex matrix architectures observed in biology. Here we show that temporally varying electrical signals can cue a self-assembling polysaccharide to controllably form a hydrogel with complex internal patterns. Evidence from theory and experiment indicate that internal structure emerges through a subtle interplay between the electrical current that triggers self-assembly and the electrical potential (or electric field) that recruits and appears to orient the polysaccharide chains at the growing gel front. These studies demonstrate that short sequences (minutes) of low-power (â¼1 V) electrical inputs can provide the program to guide self-assembly that yields hydrogels with stable, complex, and spatially varying structure and properties.