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

Free fatty acid receptors beyond fatty acids: A computational journey to explore peptides as possible binders of GPR120

Free fatty acids receptors, with members among G protein-coupled receptors (GPCRs), are crucial for biological signaling, including the perception of the so called "fatty taste". In recent years, GPR120, a protein belonging to the GPCR family, drew attention as an interesting pharmacological target to cope with obesity, satiety and diabetes. Apart from long chain fatty acids, which are GPR120 natural agonists, other synthetic molecules were identified as agonists expanding the chemical space of GPR120's ligands. In this scenario, we unveiled peptides as possible GPR120 binders toward a better understanding of this multifaceted and relevant target. This study analyzed a virtual library collecting 531 441 low-polar hexapeptides, providing mechanistic insights on the GPR120 activation and further extending the possible chemical space of GPR120 agonists. The computational pipeline started with a narrow filtering of hexapeptides based on their chemical similarity with known GPR120 agonists. The best hits were tested through docking studies, molecular dynamics and umbrella sampling simulations, which pointed to G[I,L]FGGG as a promising GPR120 agonist sequence. The presence of both peptides in food-related proteins was thoroughly assessed, revealing they may occur in mushrooms, food-grade bacteria and rice. Simulations on the counterparts with D-amino acids were also performed. Umbrella sampling simulations described that GdIFGGG may have a better interaction compared to its all-L counterpart (-13 kCal/mol ΔG and -6 kCal/mol ΔG, respectively). Overall, we obtained a predictive model to better understand the underpinning mechanism of GPR120-hexapeptides interaction, hierarchizing novel potential agonist peptides for further analysis and describing promising food sources worth of further dedicated investigations.

Revolutionizing Peptide-Based Drug Discovery: Advances in the Post-AlphaFold Era

Peptide-based drugs offer high specificity, potency, and selectivity. However, their inherent flexibility and differences in conformational preferences between their free and bound states create unique challenges that have hindered progress in effective drug discovery pipelines. The emergence of AlphaFold (AF) and Artificial Intelligence (AI) presents new opportunities for enhancing peptide-based drug discovery. We explore recent advancements that facilitate a successful peptide drug discovery pipeline, considering peptides' attractive therapeutic properties and strategies to enhance their stability and bioavailability. AF enables efficient and accurate prediction of peptide-protein structures, addressing a critical requirement in computational drug discovery pipelines. In the post-AF era, we are witnessing rapid progress with the potential to revolutionize peptide-based drug discovery such as the ability to rank peptide binders or classify them as binders/non-binders and the ability to design novel peptide sequences. However, AI-based methods are struggling due to the lack of well-curated datasets, for example to accommodate modified amino acids or unconventional cyclization. Thus, physics-based methods, such as docking or molecular dynamics simulations, continue to hold a complementary role in peptide drug discovery pipelines. Moreover, MD-based tools offer valuable insights into binding mechanisms, as well as the thermodynamic and kinetic properties of complexes. As we navigate this evolving landscape, a synergistic integration of AI and physics-based methods holds the promise of reshaping the landscape of peptide-based drug discovery.

An in silico approach to determine inter-subunit affinities in human septin complexes

The septins are a conserved family of filament-forming guanine nucleotide binding proteins, often named the fourth component of the cytoskeleton. Correctly assembled septin structures are required for essential intracellular processes such as cytokinesis, vesicular transport, polarity establishment, and cellular adhesion. Structurally, septins belong to the P-Loop NTPases but they do not mediate signals to effectors through GTP binding and hydrolysis. GTP binding and hydrolysis are believed to contribute to septin complex integrity, but biochemical approaches addressing this topic are hampered by the stability of septin complexes after recombinant expression and the lack of nucleotide-depleted complexes. To overcome this limitation, we used a molecular dynamics-based approach to determine inter-subunit binding free energies in available human septin dimer structures and in their apo forms, which we generated in silico. The nucleotide in the GTPase active subunits SEPT2 and SEPT7, but not in SEPT6, was identified as a stabilizing element in the G interface. Removal of GDP from SEPT2 and SEPT7 results in flipping of a conserved Arg residue and disruption of an extensive hydrogen bond network in the septin unique element, concomitant with a decreased inter-subunit affinity. Based on these findings we propose a singular "lock-hydrolysis" mechanism stabilizing human septin filaments.

Uncovering the Kinetic Characteristics and Degradation Preference of PROTAC Systems with Advanced Theoretical Analyses

Proteolysis-targeting chimeras (PROTACs), which can selectively induce the degradation of target proteins, represent an attractive technology in drug discovery. A large number of PROTACs have been reported, but due to the complicated structural and kinetic characteristics of the target-PROTAC-E3 ligase ternary interaction process, the rational design of PROTACs is still quite challenging. Here, we characterized and analyzed the kinetic mechanism of MZ1, a PROTAC that targets the bromodomain (BD) of the bromodomain and extra terminal (BET) protein (Brd2, Brd3, or Brd4) and von Hippel-Lindau E3 ligase (VHL), from the kinetic and thermodynamic perspectives of view by using enhanced sampling simulations and free energy calculations. The simulations yielded satisfactory predictions on the relative residence time and standard binding free energy (r_p > 0.9) for MZ1 in different Brd^BD-MZ1-VHL ternary complexes. Interestingly, the simulation of the PROTAC ternary complex disintegration illustrates that MZ1 tends to remain on the surface of VHL with the BD proteins dissociating alone without a specific dissociation direction, indicating that the PROTAC prefers more to bind with E3 ligase at the first step in the formation of the target-PROTAC-E3 ligase ternary complex. Further exploration of the degradation difference of MZ1 in different Brd systems shows that the PROTAC with higher degradation efficiency tends to leave more lysine exposed on the target protein, which is guaranteed by the stability (binding affinity) and durability (residence time) of the target-PROTAC-E3 ligase ternary complex. It is quite possible that the underlying binding characteristics of the Brd^BD-MZ1-VHL systems revealed by this study may be shared by different PROTAC systems as a general rule, which may accelerate rational PROTAC design with higher degradation efficiency.

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.