Molecular Insights into the Impact of Mutations on the Binding Affinity of Targeted Covalent Inhibitors of BTK.

Targeted covalent inhibitors (TCIs) have witnessed a significant resurgence in recent years, particularly in the kinase drug discovery field for treating diverse clinical indications. The inhibition of Bruton's tyrosine kinase (BTK) for treating B-cell cancers is a classic example where TCIs such as ibrutinib have had breakthroughs in targeted therapy. However, selectivity remains challenging, and the emergence of resistance mutations is a critical concern for clinical efficacy. Computational methods that can accurately predict the impact of mutations on inhibitor binding affinity could prove helpful in informing targeted approaches─providing insights into drug resistance mechanisms. In addition, such systems could help guide the systematic evaluation and impact of mutations in disease models for optimal experimental design. Here, we have employed in silico physics-based methods to understand the effects of mutations on the binding affinity and conformational dynamics of select TCIs of BTK. The TCIs studied include ibrutinib, acalabrutinib, and zanubrutinib─all of which are FDA-approved drugs for treating multiple forms of leukemia and lymphoma. Our results offer useful molecular insights into the structural determinants, thermodynamics, and conformational energies that impact ligand binding for this biological target of clinical relevance.

A discovery down under: decoding the draft genome sequence of Pantoea stewartii from Australia's Critically Endangered western ground parrot/kyloring (Pezoporus flaviventris).

Pantoea stewartii, a plant pathogen, is primarily transmitted through contaminated seeds and insect vectors, with the corn flea beetle (Chaetocnema pulicaria) being the primary carrier. P. stewartii is a bacterium belonging to the order Enterobacterales and can lead to crop diseases that have a significant economic impact worldwide. Due to its high potential for spread, P. stewartii is classified as a quarantine organism in numerous countries. Despite its impact on agriculture, the limited genome sequences of P. stewartii hamper understanding of its pathogenicity and host specificity, and the development of effective control strategies. In this study, a P. stewartii strain (C10109_Jinnung) was discovered in the faecal matter of the Critically Endangered western ground parrot/kyloring (Pezoporus flaviventris) in Australia, which to our knowledge is the first reported P. stewartii genome from a bird source. Whole-genome sequencing and phylogenomic analysis of strain C10109_Jinnung, obtained from a captive psittacine, provides new insights into the genetic diversity and potential transmission route for the spread of P. stewartii beyond insects and plants, where P. stewartii is typically studied. Our findings provide new insights into the potential transmission route for spread of P. stewartii and expand the known transmission agents beyond insects and plants. Expanding the catalogue of P. stewartii genomes is fundamental to improving understanding of the pathogenicity, evolution and dissemination, and to develop effective control strategies to reduce the substantial economic losses associated with P. stewartii in various crops and the potential impact of endangered animal species.

Investigating A Multi-Domain Polyketide Synthase in Amphidinium carterae.

Dinoflagellates are unicellular organisms that are implicated in harmful algal blooms (HABs) caused by potent toxins that are produced through polyketide synthase (PKS) pathways. However, the exact mechanisms of toxin synthesis are unknown due to a lack of genomic segregation of fat, toxins, and other PKS-based pathways. To better understand the underlying mechanisms, the actions and expression of the PKS proteins were investigated using the toxic dinoflagellate Amphidinium carterae as a model. Cerulenin, a known ketosynthase inhibitor, was shown to reduce acetate incorporation into all fat classes with the toxins amphidinol and sulpho-amphidinol. The mass spectrometry analysis of cerulenin-reacted synthetic peptides derived from ketosynthase domains of A. carterae multimodular PKS transcripts demonstrated a strong covalent bond that could be localized using collision-induced dissociation. One multi-modular PKS sequence present in all dinoflagellates surveyed to date was found to lack an AT domain in toxin-producing species, indicating trans-acting domains, and was shown by Western blotting to be post-transcriptionally processed. These results demonstrate how toxin synthesis in dinoflagellates can be differentiated from fat synthesis despite common underlying pathway.

Leveraging Advanced In Silico Techniques in Early Drug Discovery: A Study of Potent Small-Molecule YAP-TEAD PPI Disruptors.

Disruption of the YAP-TEAD protein-protein interaction is an attractive therapeutic strategy in oncology to suppress tumor progression and cancer metastasis. YAP binds to TEAD at a large flat binding interface (∼3500 Å2) devoid of a well-defined druggable pocket, so it has been difficult to design low-molecular-weight compounds to abrogate this protein-protein interaction directly. Recently, work by Furet and coworkers (ChemMedChem 2022, DOI: 10.1002/cmdc.202200303) reported the discovery of the first class of small molecules able to efficiently disrupt the transcriptional activity of TEAD by binding to a specific interaction site of the YAP-TEAD binding interface. Using high-throughput in silico docking, they identified a virtual screening hit from a hot spot derived from their previously rationally designed peptidic inhibitor. Structure-based drug design efforts led to the optimization of the hit compound into a potent lead candidate. Given advances in rapid high-throughput screening and rational approaches to peptidic ligand discovery for challenging targets, we analyzed the pharmacophore features involved in transferring from the peptidic to small-molecule inhibitor that could enable small-molecule discovery for such targets. Here, we show retrospectively that pharmacophore analysis augmented by solvation analysis of molecular dynamics trajectories can guide the designs, while binding free energy calculations provide greater insight into the binding conformation and energetics accompanying the association event. The computed binding free energy estimates agree well with experimental findings and offer useful insight into structural determinants that influence ligand binding to the TEAD interaction surface, even for such a shallow binding site. Taken together, our results demonstrates the utility of advanced in silico methods in structure-based design efforts for difficult-to-drug targets such as the YAP-TEAD transcription factor complex.

Benchmarking In Silico Tools for Cysteine pKa Prediction.

Accurate estimation of the pKa's of cysteine residues in proteins could inform targeted approaches in hit discovery. The pKa of a targetable cysteine residue in a disease-related protein is an important physiochemical parameter in covalent drug discovery, as it influences the fraction of nucleophilic thiolate amenable to chemical protein modification. Traditional structure-based in silico tools are limited in their predictive accuracy of cysteine pKa's relative to other titratable residues. Additionally, there are limited comprehensive benchmark assessments for cysteine pKa predictive tools. This raises the need for extensive assessment and evaluation of methods for cysteine pKa prediction. Here, we report the performance of several computational pKa methods, including single-structure and ensemble-based approaches, on a diverse test set of experimental cysteine pKa's retrieved from the PKAD database. The dataset consisted of 16 wildtype and 10 mutant proteins with experimentally measured cysteine pKa values. Our results highlight that these methods are varied in their overall predictive accuracies. Among the test set of wildtype proteins evaluated, the best method (MOE) yielded a mean absolute error of 2.3 pK units, highlighting the need for improvement of existing pKa methods for accurate cysteine pKa estimation. Given the limited accuracy of these methods, further development is needed before these approaches can be routinely employed to drive design decisions in early drug discovery efforts.

Novel irreversible peptidic inhibitors of transglutaminase 2.

Transglutaminase 2 (TG2), also referred to as tissue transglutaminase, plays crucial roles in both protein crosslinking and cell signalling. It is capable of both catalysing transamidation and acting as a G-protein, these activities being conformation-dependent, mutually exclusive, and tightly regulated. The dysregulation of both activities has been implicated in numerous pathologies. TG2 is expressed ubiquitously in humans and is localized both intracellularly and extracellularly. Targeted TG2 therapies have been developed but have faced numerous hurdles including decreased efficacy in vivo. Our latest efforts in inhibitor optimization involve the modification of a previous lead compound's scaffold by insertion of various amino acid residues into the peptidomimetic backbone, and derivatization of the N-terminus with substituted phenylacetic acids, resulting in 28 novel irreversible inhibitors. These inhibitors were evaluated for their ability to inhibit TG2 in vitro and their pharmacokinetic properties, and the most promising candidate 35 (k inact/K I = 760 × 103 M-1 min-1) was tested in a cancer stem cell model. Although these inhibitors display exceptional potency versus TG2, with k inact/K I ratios nearly ten-fold higher than their parent compound, their pharmacokinetic properties and cellular activity limit their therapeutic potential. However, they do serve as a scaffold for the development of potent research tools.

A Comparison of Dinoflagellate Thiolation Domain Binding Proteins Using In Vitro and Molecular Methods.

Dinoflagellates play important roles in ecosystems as primary producers and consumers making natural products that can benefit or harm environmental and human health but are also potential therapeutics with unique chemistries. Annotations of dinoflagellate genes have been hampered by large genomes with many gene copies that reduce the reliability of transcriptomics, quantitative PCR, and targeted knockouts. This study aimed to functionally characterize dinoflagellate proteins by testing their interactions through in vitro assays. Specifically, nine Amphidinium carterae thiolation domains that scaffold natural product synthesis were substituted into an indigoidine synthesizing gene from the bacterium Streptomyces lavendulae and exposed to three A. carterae phosphopantetheinyl transferases that activate synthesis. Unsurprisingly, several of the dinoflagellate versions inhibited the ability to synthesize indigoidine despite being successfully phosphopantetheinated. However, all the transferases were able to phosphopantetheinate all the thiolation domains nearly equally, defying the canon that transferases participate in segregated processes via binding specificity. Moreover, two of the transferases were expressed during growth in alternating patterns while the final transferase was only observed as a breakdown product common to all three. The broad substrate recognition and compensatory expression shown here help explain why phosphopantetheinyl transferases are lost throughout dinoflagellate evolution without a loss in a biochemical process.

Estimating the binding energetics of reversible covalent inhibitors of the SARS-CoV-2 main protease: an in silico study.

The main protease (Mpro) of the SARS-CoV-2 virus is an attractive therapeutic target for developing antivirals to combat COVID-19. Mpro is essential for the replication cycle of the SARS-CoV-2 virus, so inhibiting Mpro blocks a vital piece of the cell replication machinery of the virus. A promising strategy to disrupt the viral replication cycle is to design inhibitors that bind to the active site cysteine (Cys145) of the Mpro. Cysteine targeted covalent inhibitors are gaining traction in drug discovery owing to the benefits of improved potency and extended drug-target engagement. An interesting aspect of these inhibitors is that they can be chemically tuned to form a covalent, but reversible bond, with their targets of interest. Several small-molecule cysteine-targeting covalent inhibitors of the Mpro have been discovered-some of which are currently undergoing evaluation in early phase human clinical trials. Understanding the binding energetics of these inhibitors could provide new insights to facilitate the design of potential drug candidates against COVID-19. Motivated by this, we employed rigorous absolute binding free energy calculations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to estimate the energetics of binding of some promising reversible covalent inhibitors of the Mpro. We find that the inclusion of enhanced sampling techniques such as replica-exchange algorithm in binding free energy calculations can improve the convergence of predicted non-covalent binding free energy estimates of inhibitors binding to the Mpro target. In addition, our results indicate that binding free energy calculations coupled with multiscale simulations can be a useful approach to employ in ranking covalent inhibitors to their targets. This approach may be valuable in prioritizing and refining covalent inhibitor compounds for lead discovery efforts against COVID-19 and other coronavirus infections.

Dinoflagellate Phosphopantetheinyl Transferase (PPTase) and Thiolation Domain Interactions Characterized Using a Modified Indigoidine Synthesizing Reporter.

Photosynthetic dinoflagellates synthesize many toxic but also potential therapeutic compounds therapeutics via polyketide/non-ribosomal peptide synthesis, a common means of producing natural products in bacteria and fungi. Although canonical genes are identifiable in dinoflagellate transcriptomes, the biosynthetic pathways are obfuscated by high copy numbers and fractured synteny. This study focuses on the carrier domains that scaffold natural product synthesis (thiolation domains) and the phosphopantetheinyl transferases (PPTases) that thiolate these carriers. We replaced the thiolation domain of the indigoidine producing BpsA gene from Streptomyces lavendulae with those of three multidomain dinoflagellate transcripts and coexpressed these constructs with each of three dinoflagellate PPTases looking for specific pairings that would identify distinct pathways. Surprisingly, all three PPTases were able to activate all the thiolation domains from one transcript, although with differing levels of indigoidine produced, demonstrating an unusual lack of specificity. Unfortunately, constructs with the remaining thiolation domains produced almost no indigoidine and the thiolation domain for lipid synthesis could not be expressed in E. coli. These results combined with inconsistent protein expression for different PPTase/thiolation domain pairings present technical hurdles for future work. Despite these challenges, expression of catalytically active dinoflagellate proteins in E. coli is a novel and useful tool going forward.

Modeling the Binding and Conformational Energetics of a Targeted Covalent Inhibitor to Bruton's Tyrosine Kinase.

Targeted covalent inhibitors (TCIs) bind to their targets in both covalent and noncovalent modes, providing exceptionally high affinity and selectivity. These inhibitors have been effectively employed as inhibitors of protein kinases, with Taunton and coworkers (Nat. Chem. Biol. 2015, 11, 525-531) reporting a notable example of a TCI with a cyanoacrylamide warhead that forms a covalent thioether linkage to an active-site cysteine (Cys481) of Bruton's tyrosine kinase (BTK). The specific mechanism of the binding and the relative importance of the covalent and noncovalent interactions is difficult to determine experimentally, and established simulation methods for calculating the absolute binding affinity of an inhibitor cannot describe the covalent bond-forming steps. Here, an integrated approach using alchemical free-energy perturbation and QM/MM molecular dynamics methods was employed to model the complete Gibbs energy profile for the covalent inhibition of BTK by a cyanoacrylamide TCI. These calculations provide a rigorous and complete absolute Gibbs energy profile of the covalent modification binding process. Following a classic thiol-Michael addition mechanism, the target cysteine is deprotonated to form a nucleophilic thiolate, which then undergoes a facile conjugate addition to the electrophilic functional group to form a bond with the noncovalently bound ligand. This model predicts that the formation of the covalent linkage is highly exergonic relative to the noncovalent binding alone. Nevertheless, noncovalent interactions between the ligand and individual amino acid residues in the binding pocket of the enzyme are also essential for ligand binding, particularly van der Waals dispersion forces, which have a larger contribution to the binding energy than the covalent component in absolute terms. This model also shows that the mechanism of covalent modification of a protein occurs through a complex series of steps and that entropy, conformational flexibility, noncovalent interactions, and the formation of covalent linkage are all significant factors in the ultimate binding affinity of a covalent drug to its target.

A Global Approach to Estimating the Abundance and Duplication of Polyketide Synthase Domains in Dinoflagellates.

Many dinoflagellate species make toxins in a myriad of different molecular configurations but the underlying chemistry in all cases is presumably via modular synthases, primarily polyketide synthases. In many organisms modular synthases occur as discrete synthetic genes or domains within a gene that act in coordination thus forming a module that produces a particular fragment of a natural product. The modules usually occur in tandem as gene clusters with a syntenic arrangement that is often predictive of the resultant structure. Dinoflagellate genomes however are notoriously complex with individual genes present in many tandem repeats and very few synthetic modules occurring as gene clusters, unlike what has been seen in bacteria and fungi. However, modular synthesis in all organisms requires a free thiol group that acts as a carrier for sequential synthesis called a thiolation domain. We scanned 47 dinoflagellate transcriptomes for 23 modular synthase domain models and compared their abundance among 10 orders of dinoflagellates as well as their co-occurrence with thiolation domains. The total count of domain types was quite large with over thirty-thousand identified, 29 000 of which were in the core dinoflagellates. Although there were no specific trends in domain abundance associated with types of toxins, there were readily observable lineage specific differences. The Gymnodiniales, makers of long polyketide toxins such as brevetoxin and karlotoxin had a high relative abundance of thiolation domains as well as multiple thiolation domains within a single transcript. Orders such as the Gonyaulacales, makers of small polyketides such as spirolides, had fewer thiolation domains but a relative increase in the number of acyl transferases. Unique to the core dinoflagellates, however, were thiolation domains occurring alongside tetratricopeptide repeats that facilitate protein-protein interactions, especially hexa and hepta-repeats, that may explain the scaffolding required for synthetic complexes capable of making large toxins. Clustering analysis for each type of domain was also used to discern possible origins of duplication for the multitude of single domain transcripts. Single domain transcripts frequently clustered with synonymous domains from multi-domain transcripts such as the BurA and ZmaK like genes as well as the multi-ketosynthase genes, sometimes with a large degree of apparent gene duplication, while fatty acid synthesis genes formed distinct clusters. Surprisingly the acyl-transferases and ketoreductases involved in fatty acid synthesis (FabD and FabG, respectively) were found in very large clusters indicating an unprecedented degree of gene duplication for these genes. These results demonstrate a complex evolutionary history of core dinoflagellate modular synthases with domain specific duplications throughout the lineage as well as clues to how large protein complexes can be assembled to synthesize the largest natural products known.

Covalent and non-covalent binding free energy calculations for peptidomimetic inhibitors of SARS-CoV-2 main protease.

COVID-19, the disease caused by the newly discovered coronavirus-SARS-CoV-2, has created a global health, social, and economic crisis. As of mid-January 2021, there are over 90 million confirmed cases and more than 2 million reported deaths due to COVID-19. Currently, there are very limited therapeutics for the treatment or prevention of COVID-19. For this reason, it is important to find drug targets that will lead to the development of safe and effective therapeutics against the disease. The main protease (Mpro) of the virus is an attractive target for the development of effective antiviral therapeutics because it is required for proteolytic cleavage of viral polyproteins. Furthermore, the Mpro has no human homologues, so drugs designed to bind to this target directly have less risk for off-target effects. Recently, several high-resolution crystallographic structures of the Mpro in complex with inhibitors have been determined-to guide drug development and to spur efforts in structure-based drug design. One of the primary objectives of modern structure-based drug design is the accurate prediction of receptor-ligand binding affinities for rational drug design and discovery. Here, we perform rigorous alchemical absolute binding free energy calculations and QM/MM calculations to give insight into the total binding energy of two recently crystallized inhibitors of SARS-CoV-2 Mpro, namely, N3 and α-ketoamide 13b. The total binding energy consists of both covalent and non-covalent binding components since both compounds are covalent inhibitors of the Mpro. Our results indicate that the covalent and non-covalent binding free energy contributions of both inhibitors to the Mpro target differ significantly. The N3 inhibitor has more favourable non-covalent interactions, particularly hydrogen bonding, in the binding site of the Mpro than the α-ketoamide inhibitor. Also, the Gibbs energy of reaction for the Mpro-N3 covalent adduct is greater than the Gibbs reaction energy for the Mpro-α-ketoamide covalent adduct. These differences in the covalent and non-covalent binding free energy contributions for both inhibitors could be a plausible explanation for their in vitro differences in antiviral activity. Our findings are consistent with the reversible and irreversible character of both inhibitors as reported by experiment and highlight the importance of both covalent and non-covalent binding free energy contributions to the absolute binding affinity of a covalent inhibitor towards its target. This information could prove useful in the rational design, discovery, and evaluation of potent SARS-CoV-2 Mpro inhibitors for targeted antiviral therapy.

Quantum Chemical Methods for Modeling Covalent Modification of Biological Thiols.

Targeted covalent inhibitor drugs require computational methods that go beyond simple molecular-mechanical force fields in order to model the chemical reactions that occur when they bind to their targets. Here, several semiempirical and density-functional theory (DFT) methods are assessed for their ability to describe the potential energy surface and reaction energies of the covalent modification of a thiol by an electrophile. Functionals such as PBE and B3LYP fail to predict a stable enolate intermediate. This is largely due to delocalization error, which spuriously stabilizes the prereaction complex, in which excess electron density is transferred from the thiolate to the electrophile. Functionals with a high-exact exchange component, range-separated DFT functionals, and variationally optimized exact exchange (i.e., the LC-B05minV functional) correct this issue to various degrees. The large gradient behavior of the exchange enhancement factor is also found to significantly affect the results, leading to the improved performance of PBE0. While ωB97X-D and M06-2X were reasonably accurate, no method provided quantitative accuracy for all three electrophiles, making this a very strenuous test of functional performance. Additionally, one drawback of M06-2X was that molecular dynamics (MD) simulations using this functional were only stable if a fine integration grid was used. The low-cost semiempirical methods, PM3, AM1, and PM7, provide a qualitatively correct description of the reaction mechanism, although the energetics is not quantitatively reliable. As a proof of concept, the potential of mean force for the addition of methylthiolate to methylvinyl ketone was calculated using quantum mechanical/molecular mechanical MD in an explicit polarizable aqueous solvent. © 2019 Wiley Periodicals, Inc.

Mental Health Conditions and Hyperthyroidism.

OBJECTIVES: To evaluate the proportion of pediatric patients with concurrent diagnoses of hyperthyroidism and mental health conditions (MHCs) by using the Military Health System database. We hypothesized that the prevalence of mental health disorders would be higher in patients with hyperthyroidism compared with in the nonhyperthyroid population. METHODS: The prevalence of hyperthyroidism and MHCs was calculated by using data extracted from the Military Health System Data Repository on military beneficiaries between 10 and 18 years old who were eligible to receive care for at least 1 month during fiscal years 2008 through 2016. Prevalence ratios were used to compare MHC diagnoses in those with versus without a diagnosis of hyperthyroidism. RESULTS: There were 1894 female patients and 585 male patients diagnosed with hyperthyroidism during the study period. Prevalence ratios for MHCs in those with versus without hyperthyroidism ranged from 1.7 (attention-deficit/hyperactivity disorder [ADHD]) to 4.9 (bipolar disorder). Strikingly, suicidality was nearly 5 times more likely in patients diagnosed with hyperthyroidism than in patients who were never diagnosed with hyperthyroidism. For each of the MHCs examined, with the exception of suicidality, the MHC diagnosis was more commonly made before the diagnosis of hyperthyroidism, with the highest proportion of patients being diagnosed with ADHD before receiving a diagnosis of hyperthyroidism (68.3%). CONCLUSIONS: There is a clear association between hyperthyroidism and each of the following MHCs: ADHD, adjustment disorder, anxiety, bipolar disorder, depression, and suicidality. This study highlights the need to consider this association when evaluating patients with overlapping symptoms and for effective mental health screening tools and resources for clinicians.

The chemodiversity of algal dissolved organic matter from lysed Microcystis aeruginosa cells and its ability to form disinfection by-products during chlorination.

Algal-derived dissolved organic matter (ADOM) originating from lysed Microcystis aeruginosa cells was investigated as precursor material to form disinfection by-products upon disinfection with free chlorine. Non-targeted ultrahigh resolution 12 T negative mode electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) revealed high molecular diversity in solid-phase extracted and ionizable components of Microcystis aeruginosa ADOM. The toxin microcystin LR was effectively degraded by free chlorine, which was expected. However, we found a high diversity of disinfection by-products associated with the addition of free chlorine to the water-soluble and solid-phase extractable fraction of ADOM and of double-bond moieties in abundant and known unsaturated fatty acids. Aromatic DOM precursors were absent from known metabolites of Microcystis aeruginosa and no evidence for aromatic disinfection by-products (DBPs) was found, despite N-containing compounds. A large diversification of N-containing molecular formulas was observed after chlorination, which seems indicative for the breakdown and oxidation of larger peptides. Additionally, a diverse group of N-compounds with presumed chloramine functional groups was observed. This study highlights the importance to evaluate ADOM and its ability to form different DBPs when compared to allochthonous or terrestrially-derived DOM.

How Reactive are Druggable Cysteines in Protein Kinases?

Targeted covalent inhibitors (TCIs) have been successfully developed as high-affinity and selective inhibitors of enzymes of the protein kinase family. These drugs typically act by undergoing an electrophilic addition with an active-site cysteine residue, so design of a TCI begins with the identification of a "druggable" cysteine. These electrophilic additions generally require deprotonation of the thiol to form a reactive anionic thiolate, so the acidity of the residue is a critical factor. Few experimental measurements of the p Ka's of druggable cysteines have been reported, so computational prediction could prove to be very important in selecting reactive cysteine targets. Here we report the computed p Ka's of druggable cysteines in selected protein kinases that are of clinical relevance for targeted therapies. The p Ka's of the cysteines were calculated using advanced computational methods based on all-atom replica-exchange thermodynamic integration molecular dynamics simulations in explicit solvent. We found that the acidities of druggable cysteines within protein kinases are diverse and elevated, indicating enormous differences in their reactivity. Constant-pH molecular dynamics simulations were also performed on selected protein kinases, and the results confirmed this varied range in the acidities of druggable cysteines. Many of these active-site cysteines have low exposure to solvent molecules, elevating their p Ka values. Electrostatic interactions with nearby anionic residues also elevate the p Ka's of cysteine residues in the active site. The results suggest that some cysteine residues within kinase binding sites will be slow to react with a TCI because of their low acidity. Several oncogenic kinase mutations were also modeled and found to have p Ka's similar to that of the wild-type kinase.

The hydration structure of methylthiolate from QM/MM molecular dynamics.

Thiols are widely present in biological systems, most notably as the side chain of cysteine amino acids in proteins. Thiols can be deprotonated to form a thiolate which affords a diverse range of enzymatic activity and modes for chemical modification of proteins. Parameters for modeling thiolates using molecular mechanical force fields have not yet been validated, in part due to the lack of structural data on thiolate solvation. Here, the CHARMM36 and Amber models for thiolates in aqueous solutions are assessed using free energy perturbation and hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. The hydration structure of methylthiolate was calculated from 1 ns of QM/MM MD (PBE0-D3/def2-TZVP//TIP3P), which shows that the water-S- distances are approximately 2 Å with a coordination number near 6. The CHARMM thiolate parameters predict a thiolate S radius close to the QM/MM value and predict a hydration Gibbs energy of -329.2 kJ/mol, close to the experimental value of -318 kJ/mol. The cysteine thiolate model in the Amber force field underestimates the thiolate radius by 0.2 Å and overestimates the thiolate hydration energy by 119 kJ/mol because it uses the same Lennard-Jones parameters for thiolates as for thiols. A recent Drude polarizable model for methylthiolate with optimized thiolate parameters also performs well. SAPT2+ [Symmetry Adapted Perturbation Theory (SAPT)] analysis indicates that exchange repulsion is larger for the methylthiolate, consistent with it having a more diffuse electron density distribution in comparison with the parent thiol. These data demonstrate that it is important to define distinct non-bonded parameters for the protonated/deprotonated states of amino acid side chains in molecular mechanical force fields.

Limited Use of Outpatient Stress Testing in Young Patients With Atypical Chest Pain.

Low prevalence of coronary artery disease within this population suggests that younger patients may not require stress testing for chest pain evaluations as long as pretest likelihood is low.