Ligand Binding Thermodynamics in Drug Discovery: Still a Hot Tip?

Computational Analysis of Heat Capacity Effects in Protein-Ligand Binding

Heat capacity effects in protein-ligand binding as measured by calorimetric experiments have recently attracted considerable attention, particularly in the field of enzyme inhibitor design. A significant negative heat capacity change upon ligand binding implies a marked temperature dependence of the binding enthalpy, which is of high relevance for attempts to optimize protein-ligand interactions. In this work, we address the question of how well such heat capacity changes can be predicted by computer simulations. We examine a series of human thrombin inhibitors that all bind with ΔCp values of about -0.4 kcal/mol/K and calculate heat capacity changes from plain molecular dynamics simulations of the bound and free states of the enzyme and ligand. The results show that accurate ΔCp estimates within a few tenths of a kcal/mol/K of the experimental values can be obtained with this approach. This allows us to address the structural and energetic origin of the negative heat capacity changes for the thrombin inhibitors, and it is found that conformational equilibria of the free ligands in solution make a major contribution to the observed effect.

Water network chemistry to exploit the nature of catalytic water molecules in Mtb DNA gyrase: a computational study to understand the binding mechanism of fluoroquinolones

The dynamics of DNA gyrase and mutants of DNA gyrA such as G88A, A90V, S91P, D94A, D94G, D94N, D94Y; and double-point mutant (S91P-D94G), are meticulously investigated using computational approaches. Molecular dynamics (MD) and hydration thermodynamics have shed light on the fundamental, mechanistic basis of mutations on the conformational stability of Quinolone Binding Pocket (QBP) of DNA gyrase. Analysis of MD results revealed the displacement of a single crystal water molecule (HOH201) from the catalytic site of wild-type (WT) and mutants of DNA gyrA. This prompted our research group to probe the five crystal water molecules present in the QBP of the enzyme using water thermodynamics. Hydration thermodynamics analysis revealed the displacement of HOH201 due to unstable thermodynamic signatures. Further, the analysis highlighted significant changes in thermodynamic signatures and locations of five crystal water hydration sites upon mutation. Integrated MD simulations and water thermodynamics provided promising insights into the conformational changes and inaccessibility of the catalytic water molecule that can influence the design of DNA gyrase inhibitors.Communicated by Ramaswamy H. Sarma.

Target-ligand binding affinity from single point enthalpy calculation and elemental composition

Reliable target-ligand binding thermodynamics data are essential for successful drug design and molecular engineering projects. Besides experimental methods, a number of theoretical approaches have been introduced for the generation of binding thermodynamics data. However, available approaches often neglect electronic effects or explicit water molecules influencing target-ligand interactions. To handle electronic effects within a reasonable time frame, we introduce a fast calculator QMH-L using a single target-ligand complex structure pre-optimized at the molecular mechanics level. QMH-L is composed of the semi-empirical quantum mechanics calculation of binding enthalpy with predicted explicit water molecules at the complex interface, and a simple descriptor based on the elemental composition of the ligand. QMH-L estimates the target-ligand binding free energy with a root mean square error (RMSE) of 0.94 kcal mol-1. The calculations also provide binding enthalpy values and they were compared with experimental binding thermodynamics data collected from the most reliable isothermal titration calorimetry studies of systems including various protein targets and challenging, large peptide ligands with a molecular weight of up to 2-3 thousand. The single point enthalpy calculations of QMH-L require modest computational resources and are based on short runs with open source and/or free software like Gromacs, Mopac, MobyWat, and Fragmenter. QMH-L can be applied for fast, automated scoring of drug candidates during a virtual screen, enthalpic engineering of new ligands or thermodynamic explanation of complex interactions.

Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems

Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area.

Improving binding entropy by higher ligand symmetry? - A case study with human matriptase

Understanding different contributions to the binding entropy of ligands is of utmost interest to better predict affinity and the thermodynamic binding profiles of protein-ligand interactions and to develop new strategies for ligand optimization. To these means, the largely neglected effects of introducing higher ligand symmetry, thereby reducing the number of energetically distinguishable binding modes on binding entropy using the human matriptase as a model system, were investigated. A set of new trivalent phloroglucinol-based inhibitors that address the roughly symmetric binding site of the enzyme was designed, synthesized, and subjected to isothermal titration calorimetry. These highly symmetric ligands that can adopt multiple indistinguishable binding modes exhibited high entropy-driven affinity in line with affinity-change predictions.

Elucidation of a Thermodynamical Feature Attributed to Substrate Binding to the Prokaryotic H+/Oligopeptide Cotransporter YdgR with Calorimetric Analysis: The Substrate Binding Driven by the Change in Entropy Implies the Release of Bound Water Molecules from the Binding Pocket

Here, we have elucidated the substrate recognition mechanism by a prokaryotic H⁺/oligopeptide cotransporter, YdgR, using isothermal titration calorimetry. Under acidic conditions (pH 6.0), the binding of a dipeptide, Val-Ala, to YdgR elicited endothermic enthalpy, which compensated for the increase in entropy due to dipeptide binding. A series of dipeptides were used in the binding titration. The dipeptides represent Val-X and X-Val, where X is Ala, Ser, Val, Tyr, or Phe. Most dipeptides revealed endothermic enthalpy, which was completely compensated by the increase in entropy due to dipeptide binding. The change in enthalpy due to binding correlated well with the change in entropy, whereas the Gibbs free energy involved in the binding of the dipeptide to YdgR remained unchanged irrespective of dipeptide sequences, implying that the binding reaction was driven by entropy, that is, the release of bound water molecules in the binding pocket. It is also important to clarify that, based on the prediction of water molecules in the ligand-binding pocket of YdgR, the release of three bound water molecules in the putative substrate binding pocket occurred through binding to YdgR. In the comparison of Val-X and X-Val dipeptides, the N-terminal region of the binding pocket might contain more bound water molecules than the C-terminal region. In light of these findings, we suggest that bound water molecules might play an important role in substrate recognition and binding by YdgR.

Eukaryotic Expression of the Cytochrome c Oxidase Subunit I of Sitophilus zeamais and Its Interaction with Allyl Isothiocyanate

Sitophilus zeamais Motschulsky (Coleoptera: Curculionidae) is a destructive pest of stored grains around the world. Allyl isothiocyanate (AITC) was shown to have good bioactivity in the control of S. zeamais. In this study, the interaction of AITC on cytochrome c oxidase core subunits I (COX I) and their binding mechanism were determined using spectroscopic, isothermal titration calorimetry and molecular docking techniques. The results indicate the binding constant (Ka) of AITC and COX I was 6.742 × 103 L/mol. Analysis of spectroscopic revealed that the binding of COX I to reduced Cyt c induced conformational changes of reduced Cyt c, while AITC could competitively bind and inhibit the activity of the COX I protein. Moreover, molecular docking results suggested a sulfur atom in the AITC structure could form a hydrogen bond having a length of 3.3 Å with the Gly- 27 of COX I.

Isothermal titration calorimetry determination of thermodynamics of binding of cocaine and its metabolites to humanized h2E2 anti-cocaine mAb

We analyzed the thermodynamics of binding of cocaine and several cocaine metabolites to a humanized anti-cocaine mAb (h2E2), which is under development for the treatment of cocaine use disorders, using isothermal titration calorimetry. The calculated equilibrium dissociation (binding) constants were consistent with previous findings using other methods. All three ligands that display high affinity (nM) binding to the mAb (cocaine, cocaethylene, and benzoylecgonine) displayed similar enthalpically driven binding with substantial enthalpy-entropy compensation. The increased affinity of the cocaethylene metabolite compared to cocaine and benzoylecgonine is mostly attributable to a substantially less negative entropic binding component for cocaethylene, resulting in a more favorable binding energy, and thus, a higher affinity. The much lower affinity cocaine metabolites, norcocaine and ecgonine methyl ester, have much lower binding enthalpies than the high affinity ligands, and in contrast to the three high affinity ligands, have favorable (positive) entropic thermodynamic components of binding. Surprisingly, approximately 3.7 molecules of norcocaine are bound per mAb Fab site, as determined by isothermal titration calorimetry. This is in contrast to the three high affinity ligands, which bound with the expected stoichiometry of one drug molecule bound per one mAb Fab site. The results are discussed in relation to the previously published Fab:benzoylecgonine crystal structure for this h2E2 mAb, and compared to the isothermal titration calorimetry results published previously using an unrelated anti-cocaine mAb, mAb08.

•

ITC was used to measure the thermodynamics of ligand binding to anti-cocaine h2E2 mAb.

•

Binding of high and low affinity cocaine metabolites were compared to cocaine binding.

•

Low affinity metabolites have favorable/positive entropic binding components.

•

h2E2 mAb ITC results differ from those published for mAb08 anti-cocaine mAb.

•

Surprisingly, ≈3.7 molecules of low affinity norcocaine are bound per mAb Fab site.

Electrostatic Contributions to the Binding Free Energy of Nicotine to the Acetylcholine Binding Protein

Biomolecular binding relies on specific attractive interactions between two partner molecules, including electrostatics, dispersion, hydrophobicity, and solvation. Assessing the contributions of electrostatic interactions to binding is key to the understanding of ligand binding mechanisms and the design of improved biomolecular binders. For example, nicotine is a well-known agonist of nicotinic acetylcholine receptors (nAChRs), but the molecular mechanisms for the differential action of nicotine on brain and muscle nAChRs remain elusive. In this work, we have chosen the acetylcholine binding protein (AChBP) in complex with nicotine as a model system to interrogate the electrostatic contributions to nicotine binding. Our absolute binding free energy simulations confirm that nicotine binds AChBP predominantly in its protonated (charged) form. By comparing energetic contributions from decomposed interactions for either neutral or charged nicotine, our calculations shed light on the nature of the binding of nicotine to the AChBP. The preferred binding of charged nicotine over neutral nicotine originates from its stronger electrostatic interactions with AChBP, a cation-π interaction to a tryptophan residue and a hydrogen bond between nicotine and the backbone carbonyl of the tryptophan, whereas the major force driving the binding process appears to be van der Waals interactions. The various nonelectrostatic terms can also indirectly modulate the electrostatic interactions through fine-tuning the binding pose of the ligand in the binding site, providing an explanation of why the binding specificity of nicotine to the brain versus muscle nAChRs is driven by electrostatic interaction, given that the immediate binding site residues, including the key tryptophan residue, are identical in the two receptors.

Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization

Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying “magic bullets” amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand–target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural–energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.

Calculation of Heat Capacity Changes in Enzyme Catalysis and Ligand Binding

It has been suggested that heat capacity changes in enzyme catalysis may be the underlying reason for temperature optima that are not related to unfolding of the enzyme. If this were to be a common phenomenon, it would have major implications for our interpretation of enzyme kinetics. In most cases, the support for the possible existence of a nonzero (negative) activation heat capacity, however, only relies on fitting such a kinetic model to experimental data. It is therefore of fundamental interest to try to use computer simulations to address this issue. One way is simply to calculate the temperature dependence of the activation free energy and determine whether the relationship is linear or not. An alternative approach is to calculate the absolute heat capacities of the reactant and transition states from plain molecular dynamics simulations using either the temperature derivative or fluctuation formula for the enthalpy. Here, we examine these different approaches for a designer enzyme with a temperature optimum that is not caused by unfolding. Benchmark calculations for the heat capacity of liquid water are first carried out using different thermostats. It is shown that the derivative formula for the heat capacity is generally the most robust and insensitive to the thermostat used and its parameters. The enzyme calculations using this method give results in agreement with direct calculations of activation free energies and show no sign of a negative activation heat capacity. We also provide a simple scheme for the calculation of binding heat capacity changes, which is of clear interest in ligand design, and demonstrate it for substrate binding to the designer enzyme. Neither in that case do the simulations predict any negative heat capacity change.

Insight into Thermodynamic and Kinetic Profiles in Small-Molecule Optimization

Structure-activity relationships (SARs) and structure-property relationships (SPRs) have been considered the most important factors during the drug optimization process. For medicinal chemists, improvements in the potencies and druglike properties of small molecules are regarded as their major goals. Among them, the binding affinity and selectivity of small molecules on their targets are the most important indicators. In recent years, there has been growing interest in using thermodynamic and kinetic profiles to analyze ligand-receptor interactions, which could provide not only binding affinities but also detailed binding parameters for small-molecule optimization. In this perspective, we are trying to provide an insight into thermodynamic and kinetic profiles in small-molecule optimization. Through a highlight of strategies on the small-molecule optimization with specific cases, we aim to put forward the importance of structure-thermodynamic relationships (STRs) and structure-kinetic relationships (SKRs), which could provide more guidance to find safe and effective small-molecule drugs.

Importance of Binding Site Hydration and Flexibility Revealed When Optimizing a Macrocyclic Inhibitor of the Keap1-Nrf2 Protein-Protein Interaction

Upregulation of the transcription factor Nrf2 by inhibition of the interaction with its negative regulator Keap1 constitutes an opportunity for the treatment of disease caused by oxidative stress. We report a structurally unique series of nanomolar Keap1 inhibitors obtained from a natural product-derived macrocyclic lead. Initial exploration of the structure-activity relationship of the lead, followed by structure-guided optimization, resulted in a 100-fold improvement in inhibitory potency. The macrocyclic core of the nanomolar inhibitors positions three pharmacophore units for productive interactions with key residues of Keap1, including R415, R483, and Y572. Ligand optimization resulted in the displacement of a coordinated water molecule from the Keap1 binding site and a significantly altered thermodynamic profile. In addition, minor reorganizations of R415 and R483 were accompanied by major differences in affinity between ligands. This study therefore indicates the importance of accounting both for the hydration and flexibility of the Keap1 binding site when designing high-affinity ligands.

Agrobacterium fabrum atu0526-Encoding Protein Is the Only Chemoreceptor That Regulates Chemoattraction toward the Broad Antibacterial Agent Formic Acid

Soil-born plant pathogens, especially Agrobacterium, generally navigate their way to hosts through recognition of the root exudates by chemoreceptors. However, there is still a lack of appropriate identification of chemoreceptors and their ligands in Agrobacterium. Here, Atu0526, a sCache-type chemoreceptor from Agrobacterium fabrum C58, was confirmed as the receptor of a broad antibacterial agent, formic acid. The binding of formic acid to Atu0526 was screened using a thermo shift assay and verified using isothermal titration calorimetry. Inconsistent with the previously reported antimicrobial properties, formic acid was confirmed to be a chemoattractant to A. fabrum and could promote its growth. The chemotaxis of A. fabrum C58 toward formic acid was completely lost with the knock-out of atu0526, and regained with the complementation of the gene, indicating that Atu0526 is the only chemoreceptor for formic acid in A. fabrum C58. The affinity of formic acid to Atu0526^LBD significantly increased after the arginine at position 115 was replaced by alanine. However, in vivo experiments showed that the R115A mutation fully abolished the chemotaxis of A. fabrum toward formic acid. Molecular docking based on a predicted 3D structure of Atu0526 suggested that the arginine may provide "an anchorage" for formic acid to pull the minor loop, thereby forming a conformational change that generates the ligand-binding signal. Collectively, our findings will promote an understanding of sCache-type chemoreceptors and their signal transduction mechanism.

Protein-Small Molecule Interactions in Native Mass Spectrometry

Small molecule drug discovery has been propelled by the continual development of novel scientific methodologies to occasion therapeutic advances. Although established biophysical methods can be used to obtain information regarding the molecular mechanisms underlying drug action, these approaches are often inefficient, low throughput, and ineffective in the analysis of heterogeneous systems including dynamic oligomeric assemblies and proteins that have undergone extensive post-translational modification. Native mass spectrometry can be used to probe protein-small molecule interactions with unprecedented speed and sensitivity, providing unique insights into polydisperse biomolecular systems that are commonly encountered during the drug discovery process. In this review, we describe potential and proven applications of native MS in the study of interactions between small, drug-like molecules and proteins, including large multiprotein complexes and membrane proteins. Approaches to quantify the thermodynamic and kinetic properties of ligand binding are discussed, alongside a summary of gas-phase ion activation techniques that have been used to interrogate the structure of protein-small molecule complexes. We additionally highlight some of the key areas in modern drug design for which native mass spectrometry has elicited significant advances. Future developments and applications of native mass spectrometry in drug discovery workflows are identified, including potential pathways toward studying protein-small molecule interactions on a whole-proteome scale.

Uncertainty in protein-ligand binding constants: asymmetric confidence intervals versus standard errors

Equilibrium binding constants (K_b) between chemical compounds and target proteins or between interacting proteins provide a quantitative understanding of biological interaction mechanisms. Reported uncertainties of measured experimental parameters are critical for decision-making in many scientific areas, e.g., in lead compound discovery processes and in comparing computational predictions with experimental results. Uncertainties in measured K_b values are commonly represented by a symmetric normal distribution, often quoted in terms of the experimental value plus-minus the standard deviation. However, in general, the distributions of measured K_b (and equivalent K_d) values and the corresponding free energy change ΔG_b are all asymmetric to varying degree. Here, using a simulation approach, we illustrate the effect of asymmetric K_b distributions within the realm of isothermal titration calorimetry (ITC) experiments. Further we illustrate the known, but perhaps not widely appreciated, fact that when distributions of any of K_b, K_d and ΔG_b are transformed into each other, their degree of asymmetry is changed. Consequently, we recommend that a more accurate way of expressing the uncertainties of K_b, K_d, and ΔG_b values is to consistently report 95% confidence intervals, in line with other authors' suggestions. The ways to obtain such error ranges are discussed in detail and exemplified for a binding reaction obtained by ITC.

Inhibitor development of MTH1 via high-throughput screening with fragment based library and MTH1 substrate binding cavity

MutT Homolog 1 (MTH1) has been proven to hydrolyze oxidized nucleotide triphosphates during DNA repair. It can prevent the incorporation of wrong nucleotides during DNA replication and mitigate cell apoptosis. In a cancer cell, abundant reactive oxygen species can lead to substantial DNA damage and DNA mutations by base-pairing mismatch. MTH1 could eliminate oxidized dNTP and prevent cancer cells from entering cell death. Therefore, inhibition of MTH1 activity is considered to be an anti-cancer therapeutic target. In this study, high-throughput screening techniques were combined with a fragment-based library containing 2,313 compounds, which were used to screen for lead compounds with MTH1 inhibitor activity. Four compounds with MTH1 inhibitor ability were selected, and compound MI0639 was found to have the highest effective inhibition. To discover the selectivity and specificity of this action, several derivatives based on the MTH1 and MI0639 complex structure were synthesized. We compared 14 complex structures of MTH1 and the various compounds in combination with enzymatic inhibition and thermodynamic analysis. Nanomolar-range IC₅₀ inhibition abilities by enzyme kinetics and K_d values by thermodynamic analysis were obtained for two compounds, named MI1020 and MI1024. Based on structural information and compound optimization, we aim to provide a strategy for the development of MTH1 inhibitors with high selectivity and specificity.

Understanding the Interaction of Polyelectrolyte Architectures with Proteins and Biosystems

The counterions neutralizing the charges on polyelectrolytes such as DNA or heparin may dissociate in water and greatly influence the interaction of such polyelectrolytes with biomolecules, particularly proteins. In this Review we give an overview of studies on the interaction of proteins with polyelectrolytes and how this knowledge can be used for medical applications. Counterion release was identified as the main driving force for the binding of proteins to polyelectrolytes: Patches of positive charge become multivalent counterions of the polyelectrolyte and lead to the release of counterions from the polyelectrolyte and a concomitant increase in entropy. This is shown from investigations on the interaction of proteins with natural and synthetic polyelectrolytes. Special emphasis is paid to sulfated dendritic polyglycerols (dPGS). The Review demonstrates that we are moving to a better understanding of charge-charge interactions in systems of biological relevance. Research along these lines will aid and promote the design of synthetic polyelectrolytes for medical applications.

Role of cardiac drugs and flavonoids on the IRE1-JNK pathway as revealed by re-ranked molecular docking scores, MM/PBSA and umbrella sampling

One of the important causes of cardiac dysfunction is the triggering of apoptosis through the IRE1-JNK signaling pathway due to excessive ER stress (endoplasmic reticulum stress). Although there are various studies on beneficial or harmful side effects of cardiac drugs, knowledge about the molecular mechanism of their interactions on this pathway is very limited. In this study, we investigated interactions of statins, ace inhibitors, antiarrhythmic drugs and flavonoids in IRE1, ASK1(apoptosis signal-regulating kinase 1) and JNK1 at an atomic level in comparison with their well-known inhibitors. The rank of scores obtained from four different docking algorithms (Autodock 4, Autodock Vina, iGEMDOCK and GOLD) were combined so that they could be compared with each other and evaluated together. According to combined results, the most potent compound for each compound group was selected for molecular dynamics simulations, MM/PBSA (molecular mechanics/Poisson-Boltzmann surface area) and umbrella sampling calculations. We observed that the statin group drugs had the best affinity by interacting with ASK1 and JNK1 by having a similar effect with their inhibitors, and atorvastatin and pitavastatin came to the fore. Norizalpinine from the flavonoid group had a strong binding interaction with IRE1, and amiodarone from the antiarrhythmic drug group had high binding affinities with IRE1, ASK1 and JNK1. Our study has shown that atorvastatin, pitavastatin, norizalpinine and amiodarone may have a role in preventing cardiac dysfunctions caused by ER stress and may shed light on further in vitro and in vivo research.Communicated by Ramaswamy H. Sarma.

Hidden Conformational States and Strange Temperature Optima in Enzyme Catalysis

The existence of temperature optima in enzyme catalysis that occur before protein melting sets in can be described by different types of kinetic models. Such optima cause distinctly curved Arrhenius plots and have, for example, been observed in several cold-adapted enzymes from psychrophilic species. The two main explanations proposed for this behavior either invoke conformational equilibria with inactive substrate-bound states or postulate differences in heat capacity between the reactant and transition states. Herein, we analyze the implications of the different types of kinetic models in terms of apparent activation enthalpies, entropies, and heat capacities, using the catalytic reaction of a cold-adapted α-amylase as a prototypic example. We show that the behavior of these thermodynamic activation parameters is fundamentally different between equilibrium and heat capacity models, and in the α-amylase case, computer simulations have shown the former model to be correct. A few other enzyme-catalyzed reactions are also discussed in this context.

Distinct binding of cetirizine enantiomers to human serum albumin and the human histamine receptor H1

Cetirizine, a major metabolite of hydroxyzine, became a marketed second-generation H₁ antihistamine that is orally active and has a rapid onset of action, long duration of effects and a very good safety record at recommended doses. The approved drug is a racemic mixture of (S)-cetirizine and (R)-cetirizine, the latter being the levorotary enantiomer that also exists in the market as a third-generation, non-sedating and highly selective antihistamine. Both enantiomers bind tightly to the human histamine H₁ receptor (hH₁R) and behave as inverse agonists but the affinity and residence time of (R)-cetirizine are greater than those of (S)-cetirizine. In blood plasma, cetirizine exists in the zwitterionic form and more than 90% of the circulating drug is bound to human serum albumin (HSA), which acts as an inactive reservoir. Independent X-ray crystallographic work has solved the structure of the hH₁R:doxepin complex and has identified two drug-binding sites for cetirizine on equine serum albumin (ESA). Given this background, we decided to model a membrane-embedded hH₁R in complex with either (R)- or (S)-cetirizine and also the complexes of both ESA and HSA with these two enantiomeric drugs to analyze possible differences in binding modes between enantiomers and also among targets. The ensuing molecular dynamics simulations in explicit solvent and additional computational chemistry calculations provided structural and energetic information about all of these complexes that is normally beyond current experimental possibilities. Overall, we found very good agreement between our binding energy estimates and extant biochemical and pharmacological evidence. A much higher degree of solvent exposure in the cetirizine-binding site(s) of HSA and ESA relative to the more occluded orthosteric binding site in hH₁R is translated into larger positional fluctuations and considerably lower affinities for these two nonspecific targets. Whereas it is demonstrated that the two known pockets in ESA provide enough stability for cetirizine binding, only one such site does so in HSA due to a number of amino acid replacements. At the histamine-binding site in hH₁R, the distinct interactions established between the phenyl and chlorophenyl moieties of the two enantiomers with the amino acids lining up the pocket and between their free carboxylates and Lys179 in the second extracellular loop account for the improved pharmacological profile of (R)-cetirizine.

Modulation of the Self-Assembly of π-Amphiphiles in Water from Enthalpy- to Entropy-Driven by Enwrapping Substituents

Depending on the connectivity of solubilizing oligoethylene glycol (OEG) side chains to the π-cores of amphiphilic naphthalene and perylene bisimide dyes, self-assembly in water occurs either upon heating or cooling. Herein, we show that this effect originates from differences in the enwrapping capability of the π-cores by the OEG chains. Rylene bisimides bearing phenyl substituents with three OEG chains attached directly to the hydrophobic π-cores are strongly sequestered by the OEG chains. These molecules self-assemble at elevated temperatures in an entropy-driven process according to temperature- and concentration-dependent UV/Vis spectroscopy and calorimetric dilution studies. In contrast, for rylene bisimides in which phenyl substituents with three OEG chains are attached via a methylene spacer, leading to much weaker sequestration, self-assembly originates upon cooling in an enthalpy-driven process. Our explanation for this controversial behavior is that the aggregation in the latter case is dictated by the release of "high energy water" from the hydrophobic π-surfaces as well as dispersion interactions between the π-scaffolds which drive the self-assembly in an enthalpically driven process. In contrast, for the former case we suggest that in addition to the conventional explanation of a dehydration of hydrogen-bonded water molecules from OEG units it is in particular the increase in conformational entropy of back-folded OEG side chains upon aggregation that provides the pronounced gain in entropy that drives the aggregation process. Thus, our studies revealed that a subtle change in the attachment of solubilizing substituents can switch the thermodynamic signature for the self-assembly of amphiphilic dyes in water from enthalpy- to entropy-driven.

Mapping the S1 and S1' subsites of cysteine proteases with new dipeptidyl nitrile inhibitors as trypanocidal agents

The cysteine protease cruzipain is considered to be a validated target for therapeutic intervention in the treatment of Chagas disease. A series of 26 new compounds were designed, synthesized, and tested against the recombinant cruzain (Cz) to map its S1/S1´ subsites. The same series was evaluated on a panel of four human cysteine proteases (CatB, CatK, CatL, CatS) and Leishmania mexicana CPB, which is a potential target for the treatment of cutaneous leishmaniasis. The synthesized compounds are dipeptidyl nitriles designed based on the most promising combinations of different moieties in P1 (ten), P2 (six), and P3 (four different building blocks). Eight compounds exhibited a Ki smaller than 20.0 nM for Cz, whereas three compounds met these criteria for LmCPB. Three inhibitors had an EC50 value of ca. 4.0 μM, thus being equipotent to benznidazole according to the antitrypanosomal effects. Our mapping approach and the respective structure-activity relationships provide insights into the specific ligand-target interactions for therapeutically relevant cysteine proteases.

[Ligand-receptor binding kinetics in drug design]

Traditionally, the thermodynamic values of affinity are considered as the main criterion for the development of new drugs. Usually, these values for drugs are measured <i>in vitro</i> at steady concentrations of the receptor and ligand, which are differed from <i>in vivo</i> environment. Recent studies have shown that the kinetics of the process of drug binding to its receptor make significant contribution in the drug effectiveness. This has increased attention in characterizing and predicting the rate constants of association and dissociation of the receptor ligand at the stage of preclinical studies of drug candidates. A drug with a long residence time can determine ligand-receptor selectivity (kinetic selectivity), maintain pharmacological activity of the drug at its low concentration in vivo. The paper discusses the theoretical basis of protein-ligand binding, molecular determinants that control the kinetics of the drug-receptor binding. Understanding the molecular features underlying the kinetics of receptor-ligand binding will contribute to the rational design of drugs with desired properties.

Thermodynamic profiling for fragment-based lead discovery and optimization

Introduction: The enthalpic and entropic components of the ligand-protein binding free energy reflect the type and quality of the interactions and relate to the physicochemical properties of the ligands. These findings have significance in medicinal chemistry optimizations since they suggest that the thermodynamic profiling of the binding may help monitor and control the unfavorable size and hydrophobicity increase typically accompanying affinity improvements and leading to suboptimal pharmacokinetic properties.Areas covered: This review describes the ligand-protein binding event in terms of elementary steps, their associated interactions, and their enthalpic and entropic consequences. The relationships among the breaking and forming interactions, the binding thermodynamic profile, and the physicochemical properties of the ligands are also discussed.Expert opinion: Analysis of the size dependence of available affinity and favorable enthalpy highlights the limitation of the simultaneous optimization of these quantities. Indeed, moderate, rather than very high affinities can be conciliated with favorable physicochemical and pharmacokinetic profiles as it is supported by the affinity range of historical oral drugs. Although thermodynamic quantities are not suitable endpoints for medicinal chemistry optimizations owing to the complexity of the binding thermodynamics, thermodynamic profiling together with structural studies can be advantageously used to understand the details of the binding process and to optimize it.

Thermodynamic, kinetic, and structural parameterization of human carbonic anhydrase interactions toward enhanced inhibitor design

The aim of rational drug design is to develop small molecules using a quantitative approach to optimize affinity. This should enhance the development of chemical compounds that would specifically, selectively, reversibly, and with high affinity interact with a target protein. It is not yet possible to develop such compounds using computational (i.e., in silico) approach and instead the lead molecules are discovered in high-throughput screening searches of large compound libraries. The main reason why in silico methods are not capable to deliver is our poor understanding of the compound structure-thermodynamics and structure-kinetics correlations. There is a need for databases of intrinsic binding parameters (e.g., the change upon binding in standard Gibbs energy (ΔGint), enthalpy (ΔHint), entropy (ΔSint), volume (ΔVintr), heat capacity (ΔCp,int), association rate (ka,int), and dissociation rate (kd,int)) between a series of closely related proteins and a chemically diverse, but pharmacophoric group-guided library of compounds together with the co-crystal structures that could help explain the structure-energetics correlations and rationally design novel compounds. Assembly of these data will facilitate attempts to provide correlations and train data for modeling of compound binding. Here, we report large datasets of the intrinsic thermodynamic and kinetic data including over 400 primary sulfonamide compound binding to a family of 12 catalytically active human carbonic anhydrases (CA). Thermodynamic parameters have been determined by the fluorescent thermal shift assay, isothermal titration calorimetry, and by the stopped-flow assay of the inhibition of enzymatic activity. Kinetic measurements were performed using surface plasmon resonance. Intrinsic thermodynamic and kinetic parameters of binding were determined by dissecting the binding-linked protonation reactions of the protein and sulfonamide. The compound structure-thermodynamics and kinetics correlations reported here helped to discover compounds that exhibited picomolar affinities, hour-long residence times, and million-fold selectivities over non-target CA isoforms. Drug-lead compounds are suggested for anticancer target CA IX and CA XII, antiglaucoma CA IV, antiobesity CA VA and CA VB, and other isoforms. Together with 85 X-ray crystallographic structures of 60 compounds bound to six CA isoforms, the database should be of help to continue developing the principles of rational target-based drug design.

Protein-Polyelectrolyte Complexes and Micellar Assemblies

In this review, we highlight the recent progress in our understanding of the structure, properties and applications of protein-polyelectrolyte complexes in both bulk and micellar assemblies. Protein-polyelectrolyte complexes form the basis of the genetic code, enable facile protein purification, and have emerged as enterprising candidates for simulating protocellular environments and as efficient enzymatic bioreactors. Such complexes undergo self-assembly in bulk due to a combined influence of electrostatic interactions and entropy gains from counterion release. Diversifying the self-assembly by incorporation of block polyelectrolytes has further enabled fabrication of protein-polyelectrolyte complex micelles that are multifunctional carriers for therapeutic targeted delivery of proteins such as enzymes and antibodies. We discuss research efforts focused on the structure, properties and applications of protein-polyelectrolyte complexes in both bulk and micellar assemblies, along with the influences of amphoteric nature of proteins accompanying patchy distribution of charges leading to unique phenomena including multiple complexation windows and complexation on the wrong side of the isoelectric point.

Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions

The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( K_a or K_d) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.

Interaction of Proteins with Polyelectrolytes: Comparison of Theory to Experiment

We discuss recent investigations of the interaction of polyelectrolytes with proteins. In particular, we review our recent studies on the interaction of simple proteins such as human serum albumin (HSA) and lysozyme with linear polyelectrolytes, charged dendrimers, charged networks, and polyelectrolyte brushes. In all cases discussed here, we combined experimental work with molecular dynamics (MD) simulations and mean-field theories. In particular, isothermal titration calorimetry (ITC) has been employed to obtain the respective binding constants K_b and the Gibbs free energy of binding. MD simulations with explicit counterions but implicit water demonstrate that counterion release is the main driving force for the binding of proteins to strongly charged polyelectrolytes: patches of positive charges located on the surface of the protein become multivalent counterions of the polyelectrolyte, thereby releasing a number of counterions condensed on the polyelectrolyte. The binding Gibbs free energy due to counterion release is predicted to scale with the logarithm of the salt concentration in the system, which is verified by both simulations and experiment. In several cases, namely, for the interaction of proteins with linear polyelectrolytes and highly charged hydrophilic dendrimers, the binding constant could be calculated from simulations to very good approximation. This finding demonstrated that in these cases explicit hydration effects do not contribute to the Gibbs free energy of binding. The Gibbs free energy can also be used to predict the kinetics of protein uptake by microgels for a given system by applying dynamic density functional theory. The entire discussion demonstrates that the direct comparison of theory with experiments can lead to a full understanding of the interaction of proteins with charged polymers. Possible implications for applications, such as drug design, are discussed.

A protein interaction free energy model based on amino acid residue contributions: Assessment of point mutation stability of T4 lysozyme

Here we present a model to estimate the interaction free energy contribution of each amino acid residue of a given protein. Protein interaction energy is described in terms of per-residue interaction factors, μ. Multibody interactions are implicitly captured in μ through the combination of amino acid terms (γ) guided by local conformation indices (σ). The model enables construction of an interaction factor heat map for a protein in a given fold, allows prima facie assessment of the degree of residue-residue interaction, and facilitates a qualitative and quantitative evaluation of protein association properties. The model was used to compute thermal stability of T4 bacteriophage lysozyme mutants across seven sites. Qualitative assessment of mutational effects provides a straightforward rationale regarding whether a particular site primarily perturbs native or non-native states, or both. The presented model was found to be in good agreement with experimental mutational data (R 2 = 0.73) and suggests an approach by which to convert structure space into energy space.

Broad-scale analysis of thermodynamic signatures in medicinal chemistry: are enthalpy-favored binders the better development option?

Thermodynamic profiles of ligand binding, particularly enthalpically favored binding signatures, have been suggested as a criterion to support the decision-making process around which compounds to select for further optimization in drug development. The concept was enthusiastically taken up, but turned out to be too superficial, either because many aspects determining thermodynamic profiles are insufficiently appreciated or because it is difficult to compare such data on a global scale. The impact of water, changes in protonation states, along with buffer dependencies and incompatible measurement conditions that are far from standard conditions hamper such broad-scale comparisons. However, thermodynamic signatures can make us aware of the impact of these aspects and provide important hints for improving our understanding of the binding process and defining criteria for drug optimization.

Role of Hydrophobicity and Charge of Amyloid-Beta Oligomer Eliminating d-Peptides in the Interaction with Amyloid-Beta Monomers

Inhibition of the self-assembly process of amyloid-beta and even more the removal of already existing toxic amyloid-beta assemblies represent promising therapeutic strategies against Alzheimer's disease. To approach this aim, we selected a d-enantiomeric peptide by phage-display based on the interaction with amyloid-beta monomers. This lead compound was successfully optimized by peptide microarrays with respect to its affinity and specificity to the target resulting in d-peptides with both increased hydrophobicity and net charge. Here, we present a detailed biophysical characterization of the interactions between these optimized d-peptides and amyloid-beta monomers in comparison to the original lead compound in order to obtain a more thorough understanding of the physicochemical determinants of the interactions. Kinetics and apparent stoichiometry of complex formation were studied using surface plasmon resonance. Potential modes of binding to amyloid-beta were identified, and the influences of ionic strength on complex stability, as well as on the inhibitory effect on amyloid-beta aggregation were investigated. The results reveal a very different mode of interaction of the optimized d-peptides based on a combination of electrostatic and hydrophobic interactions as compared to the mostly electrostatically driven interaction of the lead compound. These conclusions were supported by the thermodynamic profiles of the interaction between optimized d-peptides and Aβ monomers, which indicate an increase in binding entropy with respect to the lead compound.

Simple Entropy Terms for End-Point Binding Free Energy Calculations

We introduce a number of computationally inexpensive modifications to the MM/PBSA and MM/GBSA estimators for binding free energies, which are based on average receptor-ligand interaction energies in simulations of a noncovalent complex, to improve the treatment of entropy: second- and higher-order terms in a cumulant expansion and a confining potential on ligand external degrees of freedom. We also consider a filter for snapshots where ligands have drifted from the initial binding pose. The variations were tested on six sets of systems for which binding modes and free energies have previously been experimentally determined. For some data sets, none of the tested estimators led to results significantly correlated with measured free energies. In data sets with nontrivial correlation, a ligand RMSD cutoff of 3 Å and a second-order truncation of the cumulant expansion was found to be comparable or better than the average interaction energy by several statistical metrics.

Application of ITC-Based Characterization of Thermodynamic and Kinetic Association of Ligands With Proteins in Drug Design

A comprehensive characterization of the thermodynamic and kinetic profiling of ligands binding to a given target protein is crucial for the hit selection as well as the hit-to-lead-to-drug evolution. Isothermal titration calorimetry (ITC), widely known as an invaluable tool to measure the thermodynamic data, has recently found its way to determine the binding kinetics too. The extensive application of ITC in measurement of both thermodynamic and kinetic data manifests unique roles of ITC in drug discovery and development. This mini-review concentrates on elaborating how to gain the thermodynamic and kinetic data using ITC, highlighting the importance of these data in lead discovery and optimization, and intends to provide an overview of the technical and conceptual advances that offer unprecedented access to protein–ligand recognition by ITC measurement.

How Ligand Protonation State Controls Water in Protein-Ligand Binding

The role of water in protein-ligand binding has been an intensely studied topic in recent years; however, how ligand protonation state change perturbs water has not been considered. Here we show that water dynamics and interactions can be controlled by the protonation state of ligand using continuous constant pH molecular dynamics simulations of two closely related model systems, β-secretase 1 and 2 (BACE1 and BACE2), in complex with a small-molecule inhibitor. Simulations revealed that, upon binding, the inhibitor pyrimidine ring remains deprotonated in BACE1 but becomes protonated in BACE2. Pyrimidine protonation results in water displacement, rigidification of the binding pocket, and shift in the ligand binding mode from water-mediated to direct hydrogen bonding. These findings not only support but also rationalize the most recent structure-selectivity data in BACE1 drug design. Binding-induced protonation state changes are likely common; our work offers a glimpse at how modeling protein-ligand binding while allowing ligand titration can further advance the understanding of water and structure-based drug design.

Thermodynamics of the Binding of Lysozyme to a Dendritic Polyelectrolyte: Electrostatics Versus Hydration

The interaction between dendritic polyglycerol sulfate (dPGS) of the second generation and lysozyme was studied by isothermal titration calorimetry (ITC) at different temperatures and salt concentrations. Analysis by ITC showed that 2-3 lysozyme molecules were bound to each dPGS. The resulting binding constant K _b and the Gibbs free energy ΔG ^o decreased markedly with increasing salt concentration but were nearly independent of temperature. The salt dependence of K _b led to the conclusion that ca. 3 counterions bound to dPGS were released upon complex formation. The gain in entropy ΔG _ci by this counterion-release scales logarithmically with salt concentration and is the main driving force for binding. The temperature dependence of ΔG ^o was analyzed by the nonlinear van't Hoff plot, taking into account a finite heat capacity change ΔC _p,vH. This evaluation led to the binding enthalpy ΔH _vH and the binding entropy ΔS _vH. Both quantities varied strongly with temperature and even changed sign, but they compensated each other throughout the entire range of temperature. Coarse-grained computer simulations with explicit salt and implicit water were used to obtain the binding free energies that agreed with ITC results. Thus, electrostatic factors were the driving forces for binding whereas all hydration contributions leading to the strongly varying ΔH _vH and ΔS _vH canceled out. The calorimetric enthalpy ΔH _ITC measured directly by ITC differed largely from ΔH _vH. ITC measurements done in two buffer systems with different ionization enthalpies revealed that binding was linked to buffer ionization and a partial protonation of the protein.

Biophysical, Biochemical, and Cell Based Approaches Used to Decipher the Role of Carbonic Anhydrases in Cancer and to Evaluate the Potency of Targeted Inhibitors

Carbonic anhydrases (CAs) are thought to be important for regulating pH in the tumor microenvironment. A few of the CA isoforms are upregulated in cancer cells, with only limited expression in normal cells. For these reasons, there is interest in developing inhibitors that target these tumor-associated CA isoforms, with increased efficacy but limited nonspecific cytotoxicity. Here we present some of the biophysical, biochemical, and cell based techniques and approaches that can be used to evaluate the potency of CA targeted inhibitors and decipher the role of CAs in tumorigenesis, cancer progression, and metastatic processes. These techniques include esterase activity assays, stop flow kinetics, and mass inlet mass spectroscopy (MIMS), all of which measure enzymatic activity of purified protein, in the presence or absence of inhibitors. Also discussed is the application of X-ray crystallography and Cryo-EM as well as other structure-based techniques and thermal shift assays to the studies of CA structure and function. Further, large-scale genomic and proteomic analytical methods, as well as cell based techniques like those that measure cell growth, apoptosis, clonogenicity, and cell migration and invasion, are discussed. We conclude by reviewing approaches that test the metastatic potential of CAs and how the aforementioned techniques have contributed to the field of CA cancer research.

Crystal structure correlations with the intrinsic thermodynamics of human carbonic anhydrase inhibitor binding

The structure-thermodynamics correlation analysis was performed for a series of fluorine- and chlorine-substituted benzenesulfonamide inhibitors binding to several human carbonic anhydrase (CA) isoforms. The total of 24 crystal structures of 16 inhibitors bound to isoforms CA I, CA II, CA XII, and CA XIII provided the structural information of selective recognition between a compound and CA isoform. The binding thermodynamics of all structures was determined by the analysis of binding-linked protonation events, yielding the intrinsic parameters, i.e., the enthalpy, entropy, and Gibbs energy of binding. Inhibitor binding was compared within structurally similar pairs that differ by para- or meta-substituents enabling to obtain the contributing energies of ligand fragments. The pairs were divided into two groups. First, similar binders—the pairs that keep the same orientation of the benzene ring exhibited classical hydrophobic effect, a less exothermic enthalpy and a more favorable entropy upon addition of the hydrophobic fragments. Second, dissimilar binders—the pairs of binders that demonstrated altered positions of the benzene rings exhibited the non-classical hydrophobic effect, a more favorable enthalpy and variable entropy contribution. A deeper understanding of the energies contributing to the protein-ligand recognition should lead toward the eventual goal of rational drug design where chemical structures of ligands could be designed based on the target protein structure.

Dissecting the Structure-Activity Relationship of Galectin-Ligand Interactions

Galectins are β-galactoside-binding proteins. As carbohydrate-binding proteins, they participate in intracellular trafficking, cell adhesion, and cell-cell signaling. Accumulating evidence indicates that they play a pivotal role in numerous physiological and pathological activities, such as the regulation on cancer progression, inflammation, immune response, and bacterial and viral infections. Galectins have drawn much attention as targets for therapeutic interventions. Several molecules have been developed as galectin inhibitors. In particular, TD139, a thiodigalactoside derivative, is currently examined in clinical trials for the treatment of idiopathic pulmonary fibrosis. Herein, we provide an in-depth review on the development of galectin inhibitors, aiming at the dissection of the structure-activity relationship to demonstrate how inhibitors interact with galectin(s). We especially integrate the structural information established by X-ray crystallography with several biophysical methods to offer, not only in-depth understanding at the molecular level, but also insights to tackle the existing challenges.

Paclitaxel inhibited lysozyme fibrillation by increasing colloidal stability through formation of "off-pathway" oligomers

Protein fibrillation is a challenging issue in medicine, causing many diseases, and an impediment to pharmaceutics and protein industry. Many chemicals, especially polyphenol compounds and aromatic small molecules, have been widely used as an effective strategy to combat protein fibril formation. Hence, understanding mechanisms of fibrillation inhibition and contributing forces in this process are significant. In this study, the inhibitory effect of paclitaxel on lysozyme fibrillation was investigated with respect to thermal and colloidal stability. Fibrillation was monitored with ThT fluorescence, circular dichroism, and AFM; paclitaxel-lysozyme interaction with isothermal titration calorimetry and docking; thermal and colloidal stability with differential scanning calorimetry and zeta-pulse, respectively. Paclitaxel inhibited lysozyme fibrillation, and interacted with lysozyme through hydrogen bonds and van der Waals' interactions. The viability of PC12 cells retrieved as a result of fibrillation inhibition by paclitaxel. Hydrophobic forces dominantly shielded the aggregation-prone region of lysozyme and suppressed the effective interactions between lysozyme monomers. Although paclitaxel did not affect lysozyme's thermal stability, it increased lysozyme's colloidal stability by either increasing the surface charge density or charge distribution on lysozyme. In conclusion, our results suggest a model for paclitaxel's inhibitory role through two complementary steps driving to "off-pathway" oligomer formation and attenuation of fibril formation.