Induced fit versus conformational selection: From rate constants to fluxes… and back to rate constants

Human Plasma Butyrylcholinesterase Hydrolyzes Atropine: Kinetic and Molecular Modeling Studies

The participation of butyrylcholinesterase (BChE) in the degradation of atropine has been recurrently addressed for more than 70 years. However, no conclusive answer has been provided for the human enzyme so far. In the present work, a steady-state kinetic analysis performed by spectrophotometry showed that highly purified human plasma BChE tetramer slowly hydrolyzes atropine at pH 7.0 and 25 °C. The affinity of atropine for the enzyme is weak, and the observed kinetic rates versus the atropine concentration was of the first order: the maximum atropine concentration in essays was much less than Km. Thus, the bimolecular rate constant was found to be kcat/Km = 7.7 × 104 M-1 min-1. Rough estimates of catalytic parameters provided slow kcat < 40 min-1 and high Km = 0.3-3.3 mM. Then, using a specific organophosphoryl agent, echothiophate, the time-dependent irreversible inhibition profiles of BChE for hydrolysis of atropine and the standard substrate butyrylthiocholine (BTC) were investigated. This established that both substrates are hydrolyzed at the same site, i.e., S198, as for all substrates of this enzyme. Lastly, molecular docking provided evidence that both atropine isomers bind to the active center of BChE. However, free energy perturbations yielded by the Bennett Acceptance Ratio method suggest that the L-atropine isomer is the most reactive enantiomer. In conclusion, the results provided evidence that plasma BChE slowly hydrolyzes atropine but should have no significant role in its metabolism under current conditions of medical use and even under administration of the highest possible doses of this antimuscarinic drug.

The complexities of ligand/receptor interactions: Exploring the role of molecular vibrations and quantum tunnelling

Molecular vibrations and quantum tunneling may link ligand binding to the function of pharmacological receptors. The well-established lock-and-key model explains a ligand's binding and recognition by a receptor; however, a general mechanism by which receptors translate binding into activation, inactivation, or modulation remains elusive. The Vibration Theory of Olfaction was proposed in the 1930s to explain this subset of receptor-mediated phenomena by correlating odorant molecular vibrations to smell, but a mechanism was lacking. In the 1990s, inelastic electron tunneling was proposed as a plausible mechanism for translating molecular vibration to odorant physiology. More recently, studies of ligands' vibrational spectra and the use of deuterated ligand analogs have provided helpful information to study this admittedly controversial hypothesis in metabotropic receptors other than olfactory receptors. In the present work, based in part on published experiments from our laboratory using planarians as an experimental organism, I will present a rationale and possible experimental approach for extending this idea to ligand-gated ion channels.

Glycan heterogeneity as a cause of the persistent fraction in HIV-1 neutralization

Neutralizing antibodies (NAbs) to multiple epitopes on the HIV-1-envelope glycoprotein (Env) have been isolated from infected persons. The potency of NAbs is measured more often than the size of the persistent fraction of infectivity at maximum neutralization, which may also influence preventive efficacy of active or passive immunization and the therapeutic outcome of the latter. Many NAbs neutralize HIV-1 CZA97.012, a clone of a Clade-C isolate, to ~100%. But here NAb PGT151, directed to a fusion-peptide epitope, left a persistent fraction of 15%. NAb PGT145, ligating the Env-trimer apex, left no detectable persistent fraction. The divergence in persistent fractions was further analyzed by depletion of pseudoviral populations of the most PGT151- and PGT145-reactive virions. Thereby, neutralization by the non-depleting NAb increased, whereas neutralization by the depleting NAb decreased. Furthermore, depletion by PGT151 increased sensitivity to autologous neutralization by sera from rabbits immunized with soluble native-like CZA97.012 trimer: substantial persistent fractions were reduced. NAbs in these sera target epitopes comprising residue D411 at the V4-β19 transition in a defect of the glycan shield on CZA97.012 Env. NAb binding to affinity-fractionated soluble native-like CZA97.012 trimer differed commensurately with neutralization in analyses by ELISA and surface plasmon resonance. Glycan differences between PGT151- and PGT145-purified trimer fractions were then demonstrated by mass spectrometry, providing one explanation for the differential antigenicity. These differences were interpreted in relation to a new structure at 3.4-Å resolution of the soluble CZA97.012 trimer determined by cryo-electron microscopy. The trimer adopted a closed conformation, refuting apex opening as the cause of reduced PGT145 binding to the PGT151-purified form. The evidence suggests that differences in binding and neutralization after trimer purification or pseudovirus depletion with PGT145 or PGT151 are caused by variation in glycosylation, and that some glycan variants affect antigenicity through direct effects on antibody contacts, whereas others act allosterically.

Competition in drug binding and … the race to equilibrium

Binding kinetics has become a popular topic in pharmacology due to its potential contribution to the selectivity and duration of drug action. Yet, the overall kinetic aspects of complex binding mechanisms are still merely described in terms of elaborate algebraic equations. Interestingly, it has been recommended some 10 years ago to examine such mechanisms in terms of binding fluxes instead of the conventional rate constants. Alike the velocity of product formation in enzymology, those fluxes refer to the velocity by which one target species converts into another one. Novel binding flux-based approaches are utilized to get a better visual insight into the "competition" between two drugs/ligands for a single target as well as between induced fit- and conformational selection pathways for a single ligand within a thermodynamic cycle. The present data were obtained by differential equation-based simulations. Early on, the ligand-binding steps "race" to equilibrium (i.e., when their forward and reverse fluxes are equal) at their individual pace. The overall/global equilibrium is only reached later on. For the competition association assays, this parting might produce a transient "overshoot" of one of the bound target species. A similar overshoot may also show up within a thermodynamic cycle and, at first glance, suggest that the induced fit pathway dominates. Yet, present findings show that under certain circumstances, it could rather be the other way round. Novel binding flux-based approaches offer visually attractive insights into crucial aspects of "complex" binding mechanisms under non-equilibrium conditions.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Elucidating the mechanisms underlying protein conformational switching using NMR spectroscopy

How proteins switch between various ligand-free and ligand-bound structures has been a key biophysical question ever since the postulation of the Monod-Wyman-Changeux and Koshland-Nemethy-Filmer models over six decades ago. The ability of NMR spectroscopy to provide structural and kinetic information on biomolecular conformational exchange places it in a unique position as an analytical tool to interrogate the mechanisms of biological processes such as protein folding and biomolecular complex formation. In addition, recent methodological developments in the areas of saturation transfer and relaxation dispersion have expanded the scope of NMR for probing the mechanics of transitions in systems where one or more states constituting the exchange process are sparsely populated and 'invisible' in NMR spectra. In this review, we highlight some of the strategies available from NMR spectroscopy for examining the nature of multi-site conformational exchange, using five case studies that have employed NMR, either in isolation, or in conjunction with other biophysical tools.