Discovery of Novel Small-Molecule Calcium Sensitizers for Cardiac Troponin C: A Combined Virtual and Experimental Screening Approach

Biosensing Platforms for Cardiac Biomarker Detection

Myocardial infarction (MI) is a cardiovascular disease that occurs when there is an elevated demand for myocardial oxygen as a result of the rupture or erosion of atherosclerotic plaques. Globally, the mortality rates associated with MI are steadily on the rise. Traditional diagnostic biomarkers employed in clinical settings for MI diagnosis have various drawbacks, prompting researchers to investigate fast, precise, and highly sensitive biosensor platforms and technologies. Biosensors are analytical devices that combine biological elements with physicochemical transducers to detect and quantify specific compounds or analytes. These devices play a crucial role in various fields including healthcare, environmental monitoring, food safety, and biotechnology. Biosensors developed for the detection of cardiac biomarkers are typically electrochemical, mass, and optical biosensors. Nanomaterials have emerged as revolutionary components in the field of biosensing, offering unique properties that significantly enhance the sensitivity and specificity of the detection systems. This review provides a comprehensive overview of the advancements and applications of nanomaterial-based biosensing systems. Beginning with an exploration of the fundamental principles governing nanomaterials, we delve into their diverse properties, including but not limited to electrical, optical, magnetic, and thermal characteristics. The integration of these nanomaterials as transducers in biosensors has paved the way for unprecedented developments in analytical techniques. Moreover, the principles and types of biosensors and their applications in cardiovascular disease diagnosis are explained in detail. The current biosensors for cardiac biomarker detection are also discussed, with an elaboration of the pros and cons of existing platforms and concluding with future perspectives.

Deep Learning-Based Chemical Similarity for Accelerated Organic Light-Emitting Diode Materials Discovery

Thermally activated delayed fluorescence (TADF) material has attracted great attention as a promising metal-free organic light-emitting diode material with a high theoretical efficiency. To accelerate the discovery of novel TADF materials, computer-aided material design strategies have been developed. However, they have clear limitations due to the accessibility of only a few computationally tractable properties. Here, we propose TADF-likeness, a quantitative score to evaluate the TADF potential of molecules based on a data-driven concept of chemical similarity to existing TADF molecules. We used a deep autoencoder to characterize the common features of existing TADF molecules with common chemical descriptors. The score was highly correlated with the four essential electronic properties of TADF molecules and had a high success rate in large-scale virtual screening of millions of molecules to identify promising candidates at almost no cost, validating its feasibility for accelerating TADF discovery. The concept of TADF-likeness can be extended to other fields of materials discovery.

Discovery of Nanomolar Inhibitors for Human Dihydroorotate Dehydrogenase Using Structure-Based Drug Discovery Methods

We used a structure-based drug discovery approach to identify novel inhibitors of human dihydroorotate dehydrogenase (DHODH), which is a therapeutic target for treating cancer and autoimmune and inflammatory diseases. In the case of acute myeloid leukemia, no previously discovered DHODH inhibitors have yet succeeded in this clinical application. Thus, there remains a strong need for new inhibitors that could be used as alternatives to the current standard-of-care. Our goal was to identify novel inhibitors of DHODH. We implemented prefiltering steps to omit PAINS and Lipinski violators at the earliest stages of this project. This enriched compounds in the data set that had a higher potential of favorable oral druggability. Guided by Glide SP docking scores, we found 20 structurally unique compounds from the ChemBridge EXPRESS-pick library that inhibited DHODH with IC50, DHODH values between 91 nM and 2.7 μM. Ten of these compounds reduced MOLM-13 cell viability with IC50, MOLM-13 values between 2.3 and 50.6 μM. Compound 16 (IC50, DHODH = 91 nM) inhibited DHODH more potently than the known DHODH inhibitor, teriflunomide (IC50, DHODH = 130 nM), during biochemical characterizations and presented a promising scaffold for future hit-to-lead optimization efforts. Compound 17 (IC50, MOLM-13 = 2.3 μM) was most successful at reducing survival in MOLM-13 cell lines compared with our other hits. The discovered compounds represent excellent starting points for the development and optimization of novel DHODH inhibitors.

Mouse Models of Cardiomyopathies Caused by Mutations in Troponin C

Cardiac muscle contraction is regulated via Ca2+ exchange with the hetero-trimeric troponin complex located on the thin filament. Binding of Ca2+ to cardiac troponin C, a Ca2+ sensing subunit within the troponin complex, results in a series of conformational re-arrangements among the thin filament components, leading to an increase in the formation of actomyosin cross-bridges and muscle contraction. Ultimately, a decline in intracellular Ca2+ leads to the dissociation of Ca2+ from troponin C, inhibiting cross-bridge cycling and initiating muscle relaxation. Therefore, troponin C plays a crucial role in the regulation of cardiac muscle contraction and relaxation. Naturally occurring and engineered mutations in troponin C can lead to altered interactions among components of the thin filament and to aberrant Ca2+ binding and exchange with the thin filament. Mutations in troponin C have been associated with various forms of cardiac disease, including hypertrophic, restrictive, dilated, and left ventricular noncompaction cardiomyopathies. Despite progress made to date, more information from human studies, biophysical characterizations, and animal models is required for a clearer understanding of disease drivers that lead to cardiomyopathies. The unique use of engineered cardiac troponin C with the L48Q mutation that had been thoroughly characterized and genetically introduced into mouse myocardium clearly demonstrates that Ca2+ sensitization in and of itself should not necessarily be considered a disease driver. This opens the door for small molecule and protein engineering strategies to help boost impaired systolic function. On the other hand, the engineered troponin C mutants (I61Q and D73N), genetically introduced into mouse myocardium, demonstrate that Ca2+ desensitization under basal conditions may be a driving factor for dilated cardiomyopathy. In addition to enhancing our knowledge of molecular mechanisms that trigger hypertrophy, dilation, morbidity, and mortality, these cardiomyopathy mouse models could be used to test novel treatment strategies for cardiovascular diseases. In this review, we will discuss (1) the various ways mutations in cardiac troponin C might lead to disease; (2) relevant data on mutations in cardiac troponin C linked to human disease, and (3) all currently existing mouse models containing cardiac troponin C mutations (disease-associated and engineered).

Targeting Troponin C with Small Molecules Containing Diphenyl Moieties: Calcium Sensitivity Effects on Striated Muscles and Structure-Activity Relationship

Despite large investments from academia and industry, heart failure, which results from a disruption of the contractile apparatus, remains a leading cause of death. Cardiac muscle contraction is a calcium-dependent mechanism, which is regulated by the troponin protein complex (cTn) and specifically by the N-terminal domain of its calcium-binding subunit (cNTnC). There is an increasing need for the development of small molecules that increase calcium sensitivity without altering the systolic calcium concentration, thereby strengthening the cardiac function. Here, we examined the effect of our previously identified calcium-sensitizing small molecule, ChemBridge compound 7930079, in the context of several homologous muscle systems. The effect of this molecule on force generation in isolated cardiac trabeculae and slow skeletal muscle fibers was measured. Furthermore, we explored the use of Gaussian accelerated molecular dynamics in sampling highly predictive receptor conformations based on NMR-derived starting structures. Additionally, we took a rational computational approach for lead optimization based on lipophilic diphenyl moieties. This integrated structural-biochemical-physiological approach led to the identification of three novel low-affinity binders, which had similar binding affinities to the known positive inotrope trifluoperazine. The most potent identified calcium sensitizer was compound 16 with an apparent affinity of 117 ± 17 μM.

Umbrella Sampling Simulations of Cardiac Thin Filament Reveal Thermodynamic Consequences of Troponin I Inhibitory Peptide Mutations

The cardiac thin filament comprises F-actin, tropomyosin, and troponin (cTn). cTn is composed of three subunits: troponin C (cTnC), troponin I (cTnI), and troponin T (cTnT). To computationally study the effect of the thin filament on cTn activation events, we employed targeted molecular dynamics followed by umbrella sampling using a model of the thin filament to measure the thermodynamics of cTn transition events. Our simulations revealed that the thin filament causes an increase in the free energy required to open the cTnC hydrophobic patch and causes a more favorable interaction between this region and the cTnI switch peptide. Mutations to the cTn complex can lead to cardiomyopathy, a collection of diseases that present clinically with symptoms of hypertrophy or dilation of the cardiac muscle, leading to impairment of the heart's ability to function normally and ultimately myocardial infarction or heart failure. Upon introduction of cardiomyopathic mutations to R145 of cTnI, we observed a general decrease in the free energy of opening the cTnC hydrophobic patch, which is on par with previous experimental results. These mutations also exhibited a decrease in electrostatic interactions between cTnI-R145 and actin-E334. After introduction of a small molecule to the wild-type cTnI-actin interface to intentionally disrupt intersubunit contacts, we successfully observed similar thermodynamic consequences and disruptions to the same protein-protein contacts as observed with the cardiomyopathic mutations. Computational studies utilizing the cTn complex in isolation would have been unable to observe these effects, highlighting the importance of using a more physiologically relevant thin-filament model to investigate the global consequences of cardiomyopathic mutations to the cTn complex.

Targeting Troponin C with Small Molecules Containing Diphenyl Moieties: Calcium Sensitivity Effects on Striated Muscle and Structure Activity Relationship

Despite large investments from academia and industry, heart failure, which results from a disruption of the contractile apparatus, remains a leading cause of death. Cardiac muscle contraction is a calcium-dependent mechanism, which is regulated by the troponin protein complex (cTn) and specifically by the N-terminal domain of its calcium binding subunit (cNTnC). There is an increasing need for the development of small molecules that increase calcium sensitivity without altering systolic calcium concentration, thereby strengthening cardiac function. Here, we examined the effect of our previously identified calcium sensitizing small molecule, ChemBridge compound 7930079, in the context of several homologous muscle systems. The effect of this molecule on force generation in isolated cardiac trabeculae and slow skeletal muscle fibers was measured. Furthermore, we explored the use of Gaussian accelerated molecular dynamics in sampling highly predictive receptor conformations based on NMR derived starting structures. Additionally, we took a rational computational approach for lead optimization based on lipophilic diphenyl moieties. This led to the identification of three novel low affinity binders, which had similar binding affinities to known positive inotrope trifluoperazine. The most potent identified calcium sensitizer was compound 16 with an apparent affinity of 117 ± 17 μM .

Computational Exploration and Characterization of Potential Calcium Sensitizing Mutations in Cardiac Troponin C

Calcium-dependent heart muscle contraction is regulated by the cardiac troponin protein complex (cTn) and specifically by the N-terminal domain of its calcium binding subunit (cNTnC). cNTnC contains one calcium binding site (site II), and altered calcium binding in this site has been studied for decades. It has been previously shown that cNTnC mutants, which increase calcium sensitization may have therapeutic benefits, such as restoring cardiac muscle contractility and functionality post-myocardial infarction events. Here, we computationally characterized eight mutations for their potential effects on calcium binding affinity in site II of cNTnC. We utilized two distinct methods to estimate calcium binding: adaptive steered molecular dynamics (ASMD) and thermodynamic integration (TI). We observed a sensitizing trend for all mutations based on the employed ASMD methodology. The TI results showed excellent agreement with experimentally known calcium binding affinities in wild-type cNTnC. Based on the TI results, five mutants were predicted to increase calcium sensitivity in site II. This study presents an interesting comparison of the two computational methods, which have both been shown to be valuable tools in characterizing the impacts of calcium sensitivity in mutant cNTnC systems.

The interplay between cardiac dyads and mitochondria regulated the calcium handling in cardiomyocytes

Calcium mishandling and mitochondrial dysfunction have been increasingly recognized as significant factors involved in the progression procedure of cardiomyopathy. Ca2+ mishandling could cause calcium-triggered arrhythmias, which could enhance force development and ATP consumption. Mitochondrial disorganization and dysfunction in cardiomyopathy could disturb the balance of energy catabolic and anabolic procedure. Close spatial localization and arrangement of structural among T-tubule, sarcoplasmic reticulum, mitochondria are important for Ca2+ handling. So that, we illustrate the regulating network between calcium handling and mitochondrial homeostasis, as well as its intracellular mechanisms in this review, which would be worthy to develop novel therapeutic strategy and restore the function of injured cardiomyocytes.

Actives-Based Receptor Selection Strongly Increases the Success Rate in Structure-Based Drug Design and Leads to Identification of 22 Potent Cancer Inhibitors

Computer-aided drug design, an important component of the early stages of the drug discovery pipeline, routinely identifies large numbers of false positive hits that are subsequently confirmed to be experimentally inactive compounds. We have developed a methodology to improve true positive prediction rates in structure-based drug design and have successfully applied the protocol to twenty target systems and identified the top three performing conformers for each of the targets. Receptor performance was evaluated based on the area under the curve of the receiver operating characteristic curve for two independent sets of known actives. For a subset of five diverse cancer-related disease targets, we validated our approach through experimental testing of the top 50 compounds from a blind screening of a small molecule library containing hundreds of thousands of compounds. Our methods of receptor and compound selection resulted in the identification of 22 novel inhibitors in the low μM-nM range, with the most potent being an EGFR inhibitor with an IC₅₀ value of 7.96 nM. Additionally, for a subset of five independent target systems, we demonstrated the utility of Gaussian accelerated molecular dynamics to thoroughly explore a target system's potential energy surface and generate highly predictive receptor conformations.

Umbrella Sampling Simulations Measure Switch Peptide Binding and Hydrophobic Patch Opening Free Energies in Cardiac Troponin

The cardiac troponin (cTn) complex is an important regulatory protein in heart contraction. Upon binding of Ca2+, cTn undergoes a conformational shift that allows the troponin I switch peptide (cTnISP) to be released from the actin filament and bind to the troponin C hydrophobic patch (cTnCHP). Mutations and modifications to this complex can change its sensitivity to Ca2+ and alter the energetics of the transition from the Ca2+-unbound, cTnISP-unbound form to the Ca2+-bound, cTnISP-bound form. We utilized targeted molecular dynamics (TMD) to obtain a trajectory of this transition pathway, followed by umbrella sampling to estimate the free energy associated with the cTnISP-cTnCHP binding and the cTnCHP opening events for wild-type (WT) cTn. We were able to reproduce experimental values for the cTnISP-cTnCHP binding event and obtain cTnCHP opening free energies in agreement with previous computational measurements of smaller cTnC systems. This excellent agreement for WT cTn demonstrated the strength of computational methods in studying the dynamics and energetics of the cTn complex. We then introduced mutations to the cTn complex that cause cardiomyopathy or alter its Ca2+ sensitivity and observed a general decrease in the free energy of opening the cTnCHP. For these same mutations, we observed no general trend in the effect on the cTnISP-cTnCHP binding event. Our method sets the stage for future computational studies on this system that predict the consequences of yet uncharacterized mutations on cTn dynamics and energetics.

Caffeoyl malic acid is a potential dual inhibitor targeting TNFα/IL-4 evaluated by a combination strategy of network analysis-deep learning-molecular simulation

Atopic dermatitis (AD) is a common inflammatory skin disease involving multiple signaling pathways. One of the effective treatment strategies of AD is to develop a new drug capable of regulating the key therapeutic targets. Here we report the combination use of network analysis, deep learning, and molecular simulation for the identification of key therapeutic targets for AD and screening of potential multi-target drugs. From the TCM@Taiwan database, we identify a small molecule, namely caffeoyl malic acid (CMA), to inhibit the key therapeutic targets (TNFα and IL-4) for AD. CMA is further identified as a TNFα inhibitor by a deep learning model based on convolutional neural network. Molecular simulations demonstrate that CMA can stably bind to TNFα and IL-4, thereby producing diverse effects on the structural fluctuation, structural flexibility, looseness, and motion strength of each protein. Furthermore, conformation alignments reveal that CMA makes the distance between chain A and C of TNFα become wider and the slit between the two α helices of IL-4 get narrow obviously. CMA leads to the change of protein conformation, which hinders the formation of the protein-receptor complex. Collectively, our findings suggest that CMA is a potential dual TNFα/IL-4 inhibitor for the treatment of AD.

In Silico and Experimental ADAM17 Kinetic Modeling as Basis for Future Screening System for Modulators

Understanding the mechanisms of modulators’ action on enzymes is crucial for optimizing and designing pharmaceutical substances. The acute inflammatory response, in particular, is regulated mainly by a disintegrin and metalloproteinase (ADAM) 17. ADAM17 processes several disease mediators such as TNFα and APP, releasing their soluble ectodomains (shedding). A malfunction of this process leads to a disturbed inflammatory response. Chemical protease inhibitors such as TAPI-1 were used in the past to inhibit ADAM17 proteolytic activity. However, due to ADAM17′s broad expression and activity profile, the development of active-site-directed ADAM17 inhibitor was discontinued. New ‘exosite’ (secondary substrate binding site) inhibitors with substrate selectivity raised the hope of a substrate-selective modulation as a promising approach for inflammatory disease therapy. This work aimed to develop a high-throughput screen for potential ADAM17 modulators as therapeutic drugs. By combining experimental and in silico methods (structural modeling and docking), we modeled the kinetics of ADAM17 inhibitor. The results explain ADAM17 inhibition mechanisms and give a methodology for studying selective inhibition towards the design of pharmaceutical substances with higher selectivity.

Applications of machine learning in computer-aided drug discovery

Machine learning (ML) has revolutionised the field of structure-based drug design (SBDD) in recent years. During the training stage, ML techniques typically analyse large amounts of experimentally determined data to create predictive models in order to inform the drug discovery process. Deep learning (DL) is a subfield of ML, that relies on multiple layers of a neural network to extract significantly more complex patterns from experimental data, and has recently become a popular choice in SBDD. This review provides a thorough summary of the recent DL trends in SBDD with a particular focus on de novo drug design, binding site prediction, and binding affinity prediction of small molecules.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Molecular dynamics provides new insights into the mechanism of calcium signal transduction and interdomain interactions in cardiac troponin

Understanding the regulation of cardiac muscle contraction at a molecular level is crucial for the development of therapeutics for heart conditions. Despite the availability of atomic structures of the protein components of cardiac muscle thin filaments, detailed insights into their dynamics and response to calcium are yet to be fully depicted. In this study, we used molecular dynamics simulations of the core domains of the cardiac muscle protein troponin to characterize the equilibrium dynamics of its calcium-bound and calcium-free forms, with a focus on elements of cardiac muscle contraction activation and deactivation, that is, calcium binding to the cardiac troponin Ca²⁺ -binding subunit (TnC) and the release of the switch region of the troponin inhibitory subunit (TnI) from TnC. The process of calcium binding to the TnC binding site is described as a three-step process commencing with calcium capture by the binding site residues, followed by cooperative residue interplay bringing the calcium ion to the binding site, and finally, calcium-water exchange. Furthermore, we uncovered a set of TnC-TnI interdomain interactions that are critical for TnC N-lobe hydrophobic pocket dynamics. Absence of these interactions allows the closure of the TnC N-lobe hydrophobic pocket while the TnI switch region remains expelled, whereas if the interactions are maintained, the hydrophobic pocket remains open. Modification of these interactions may fine-tune the ability of the TnC N-lobe hydrophobic pocket to close or remain open, modulate cardiac contractility and present potential therapy-relevant targets.

High Throughput Screen Identifies Small Molecule Effectors That Modulate Thin Filament Activation in Cardiac Muscle

Current therapeutic interventions for both heart disease and heart failure are largely insufficient and associated with undesired side effects. Biomedical research has emphasized the role of sarcomeric protein function for the normal performance and energy efficiency of the heart, suggesting that directly targeting the contractile myofilaments themselves using small molecule effectors has therapeutic potential and will likely result in greater drug efficacy and selectivity. In this study, we developed a robust and highly reproducible fluorescence polarization-based high throughput screening (HTS) assay that directly targets the calcium-dependent interaction between cardiac troponin C (cTnC) and the switch region of cardiac troponin I (cTnI_SP), with the aim of identifying small molecule effectors of the cardiac thin filament activation pathway. We screened a commercially available small molecule library and identified several hit compounds with both inhibitory and activating effects. We used a range of biophysical and biochemical methods to characterize hit compounds and identified fingolimod, a sphingosin-1-phosphate receptor modulator, as a new troponin-based small molecule effector. Fingolimod decreased the ATPase activity and calcium sensitivity of demembranated cardiac muscle fibers in a dose-dependent manner, suggesting that the compound acts as a calcium desensitizer. We investigated fingolimod's mechanism of action using a combination of computational studies, biophysical methods, and synthetic chemistry, showing that fingolimod bound to cTnC repels cTnI_SP via mainly electrostatic repulsion of its positively charged tail. These results suggest that fingolimod is a potential new lead compound/scaffold for the development of troponin-directed heart failure therapeutics.