Rate Of Interconversion Research Articles

Evaluation of transition rates from nonequilibrium instantons

Equilibrium rate theories play a crucial role in understanding rare, reactive events. However, they are inapplicable to a range of irreversible processes in systems driven far from thermodynamic equilibrium like active and biological matter. Here we develop an efficient numerical method to compute the rate constant of rare nonequilibrium events in the weak-noise limit based on an instanton approximation to the stochastic path integral and illustrate its wide range of application. We demonstrate excellent agreement of the instanton rates with numerically exact results for a particle under a nonconservative force. We also study phase transitions in an active field theory. We elucidate how activity alters the stability of the two phases and their rates of interconversion in a manner that can be well described by modifying classical nucleation theory. Published by the American Physical Society 2024

Redox-Triggered Reversible Switching between Dynamic and Quasi-static α-Helical Peptides.

We report the reversible transformation between a singly stapled dynamic α-helical peptide and a doubly stapled quasi-static one through redox-triggered dithiol/disulfide conversions of a stapling moiety. This process allows the rate of interconversion between the right-handed (P) and left-handed (M) α-helices to be altered by a factor of approximately 103 before and after the transformation. An as-obtained doubly stapled α-helical peptide, which is composed of an achiral peptide having an l-valine carboxylic acid residue at the C-terminus, a disulfide-based reversible staple, and a biphenyl-based fixed staple, adopts an (M)-rich form as a kinetically trapped state. The (M)-rich helix was subsequently transformed into the thermodynamically stable (P)-rich form in 1,1,2,2-tetrachloroethane with the half-life time (t1/2) of approximately 44 days at 25 °C. Reduction of the doubly stapled peptide with tri-n-butylphosphine in tetrahydrofuran/water (10/1, v/v) produced the corresponding singly stapled dynamic α-helical peptide bearing two thiol groups at the side chains, which underwent solvent-induced reversible helicity inversion. The resulting dithiol of the singly stapled peptide could be reoxidized to form the original doubly stapled form using 4,4'-dithiodipyridine. Furthermore, the P/M interconversion of a doubly stapled peptide with two flexible hydrocarbon-based staples is considerably more rapid than that with more rigid staples.

19F NMR Reveals the Dynamics of Substrate Binding and Lid Closure for Iodotyrosine Deiodinase as a Complement to Steady-State Kinetics and Crystallography.

Active site lids are common features of enzymes and typically undergo conformational changes upon substrate binding to promote catalysis. Iodotyrosine deiodinase is no exception and contains a lid segment in all of its homologues from human to bacteria. The solution-state dynamics of the lid have now been characterized using 19F NMR spectroscopy with a CF3-labeled enzyme and CF3O-labeled ligands. From two-dimensional 19F-19F NMR exchange spectroscopy, interconversion rates between the free and bound states of a CF3O-substituted tyrosine (45 ± 10 s-1) and the protein label (40 ± 3 s-1) are very similar and suggest a correlation between ligand binding and conformational reorganization of the lid. Both occur at rates that are ∼100-fold faster than turnover, and therefore these steps do not limit catalysis. A simple CF3O-labeled phenol also binds to the active site and induces a conformational change in the lid segment that was not previously detectable by crystallography. Exchange rates of the ligand (130 ± 20 s-1) and protein (98 ± 8 s-1) in this example are faster than those above but remain self-consistent to affirm a correlation between ordering of the lid and binding of the ligand. Both ligands also protect the protein from limited proteolysis, as expected from their ability to stabilize a compact lid structure. However, the minimal turnover of simple phenol substrates indicates that such stabilization may be necessary but is not sufficient for efficient catalysis.

Lenacapavir Exhibits Atropisomerism-Mechanistic Pharmacokinetics and Disposition Studies of Lenacapavir Reveal Intestinal Excretion as a Major Clearance Pathway.

Lenacapavir (LEN), a long-acting injectable, is the first approved human immunodeficiency virus type 1 capsid inhibitor and one of a few Food and Drug Administration-approved drugs that exhibit atropisomerism. LEN exists as a mixture of two class 2 atropisomers that interconvert at a fast rate (half-life < 2 hours) with a ratio that is stable over time and unaffected by enzymes or binding to proteins in plasma. LEN exhibits low systemic clearance (CL) in nonclinical species and humans; however, in all species, the observed CL was higher than the in vitro predicted CL. The volume of distribution was moderate in nonclinical species and consistent with the tissue distribution observed by whole-body autoradiography in rats. LEN does not distribute to brain, consistent with being a P-glycoprotein (P-gp) substrate. Mechanistic drug disposition studies with [14C]LEN in intravenously dosed bile duct-cannulated rats and dogs showed a substantial amount of unchanged LEN (31%-60% of dose) excreted in feces, indicating that intestinal excretion (IE) was a major clearance pathway for LEN in both species. Coadministration of oral elacridar, a P-gp inhibitor, in rats decreased CL and IE of LEN. Renal excretion was < 1% of dose in both species. In plasma, almost all radioactivity was unchanged LEN. Low levels of metabolites in excreta included LEN conjugates with glutathione, pentose, and glucuronic acid, which were consistent with metabolites formed in vitro in Hμrel hepatocyte cocultures and those observed in human. Our studies highlight the importance of IE for efflux substrates that are highly metabolically stable compounds with slow elimination rates. SIGNIFICANCE STATEMENT: LEN is a long-acting injectable that exists as conformationally stable atropisomers. Due to an atropisomeric interconversion rate that significantly exceeds the in vivo elimination rate, the atropisomer ratio of LEN remains constant in circulation. The disposition of LEN highlights that intestinal excretion has a substantial part in the elimination of compounds that are metabolically highly stable and efflux transporter substrates.

Proton Transfer Kinetics in Histidine Side Chains Determined by pH-Dependent Multi-Nuclear NMR Relaxation.

Histidine is a key amino-acid residue in proteins with unique properties engendered by its imidazole side chain that can exist in three different states: two different neutral tautomeric forms and a protonated, positively charged one with a pKa value close to physiological pH. Commonly, two or all three states coexist and interchange rapidly, enabling histidine to act as both donor and acceptor of hydrogen bonds, coordinate metal ions, and engage in acid/base catalysis. Understanding the exchange dynamics among the three states is critical for assessing histidine's mechanistic role in catalysis, where the rate of proton exchange and interconversion among tautomers might be rate limiting for turnover. Here, we determine the exchange kinetics of histidine residues with pKa values representative of the accessible range from 5 to 9 by measuring pH-dependent 15N, 13C, and 1H transverse relaxation rate constants for 5 nuclei in each imidazole. Proton exchange between the imidazole and the solvent is mediated by hydronium ions at acidic and neutral pH, whereas hydroxide mediated exchange becomes the dominant mechanism at basic pH. Proton transfer is very fast and reaches the diffusion limit for pKa values near neutral pH. We identify a direct pathway between the two tautomeric forms, likely mediated by a bridging water molecule or, in the case of high pH, hydroxide ion. For histidines with pKa 7, we determine all rate constants (lifetimes) involving protonation over the entire pH range. Our approach should enable critical insights into enzymatic acid/base catalyzed reactions involving histidines in proteins.

A Transient, Excited Species of the Autoinhibited α-State of the Bacterial Transcription Factor RfaH May Represent an Early Intermediate on the Fold-Switching Pathway.

RfaH is a two-domain transcription factor in which the C-terminal domain switches fold from an α-helical hairpin to a β-roll upon binding the ops-paused RNA polymerase. To ascertain the presence of a sparsely populated excited state that may prime the autoinhibited resting state of RfaH for binding ops-paused RNA polymerase, we carried out a series of NMR-based exchange experiments to probe for conformational exchange on the millisecond time scale. Quantitative analysis of these data reveals exchange between major ground (∼95%) and sparsely populated excited (∼5%) states with an exchange lifetime of ∼3 ms involving residues at the interface between the N-terminal and C-terminal domains formed by the β3/β4 hairpin and helix α3 of the N-terminal domain and helices α4 and α5 of the C-terminal domain. The largest 15N backbone chemical shift differences are associated with the β3/β4 hairpin, leading us to suggest that the excited state may involve a rigid body lateral displacement/rotation away from the C-terminal domain to adopt a position similar to that seen in the active RNA polymerase-bound state. Such a rigid body reorientation would result in a reduction in the interface between the N- and C-terminal domains with the possible introduction of a cavity or cavities. This hypothesis is supported by the observation that the population of the excited species and the exchange rate of interconversion between ground and excited states are reduced at a high (2.5 kbar) pressure. Mechanistic implications for fold switching of the C-terminal domain in the context of RNA polymerase binding are discussed.

Quantum Nonlocality: Multicopy Resource Interconvertibility and Their Asymptotic Inequivalence.

Quantum nonlocality, pioneered in Bell's seminal work and subsequently verified through a series of experiments, has drawn substantial attention due to its practical applications in various protocols. Evaluating and comparing the extent of nonlocality within distinct quantum correlations holds significant practical relevance. Within the resource theoretic framework this can be achieved by assessing the interconversion rate among different nonlocal correlations under free local operations and shared randomness. In this study we, however, present instances of quantum nonlocal correlations that are incomparable in the strongest sense. Specifically, when starting with an arbitrary many copies of one nonlocal correlation, it becomes impossible to obtain even a single copy of the other correlation, and this incomparability holds in both directions. Such incomparable quantum correlations can be obtained even in the simplest Bell scenario, which involves two parties, each having two dichotomic measurements setups. Notably, there exist an uncountable number of such incomparable correlations. Our result challenges the notion of a "unique gold coin," often referred to as the "maximally resourceful state," within the framework of the resource theory of quantum nonlocality. To this end, we provide examples of isotropic quantum correlations that cannot be distilled up to the Tsirelson point, and thus partially answer a long-standing open question in nonlocality distillation.

Quadratic energy–momentum squared gravity: Constraints from big bang nucleosynthesis

In this work, we extend the standard cosmological model within the quadratic energy–momentum squared gravity (qEMSG) framework, introducing a nonminimal interaction between the usual material field (Tμν) and its accompanying partner field (qEMSF, TμνqEMSF), defined by f(T2)=αT2 with T2=TμνTμν. Adopting an analytical approach within the qEMSG framework, we present a comprehensive exploration of Big Bang Nucleosynthesis (BBN) dynamics. Our analysis selects the radiation-dominated universe solution compatible with the standard cosmological model limit as α→0 and reveals that qEMSF interaction model can modify the radiation energy density’s evolution, potentially altering neutron–proton interconversion rates and consequently affecting 4He abundance in various ways. By explicitly defining modifications to the predicted primordial 4He mass fraction, Yp, we establish the most stringent cosmological constraints on the parameter α based on recent measurements of Yp: (−8.81≤α≤8.14)×10−27eV−4 (68% CL) from Aver et al.’s primordial 4He abundance measurements, aligning with α=0. Additionally, (3.48≤α≤4.43)×10−27eV−4 (68% CL) from Fields et al.’s estimates, utilizing the Planck-CMB estimated baryon density within the standard cosmological model framework, diverges from α=0, thereby lending support to the qEMSF interaction model. The study also highlights the bidirectional nature of energy–momentum/entropy transfer in qEMSF interaction model, depending on the sign of α. The implications of qEMSF in the presence of additional relativistic relics are also explored, showcasing the model’s potential to accommodate deviations from standard cosmology and the Standard Model of particle physics.

Watching a signaling protein function: What has been learned over four decades of time-resolved studies of photoactive yellow protein.

Photoactive yellow protein (PYP) is a signaling protein whose internal p-coumaric acid chromophore undergoes reversible, light-induced trans-to-cis isomerization, which triggers a sequence of structural changes that ultimately lead to a signaling state. Since its discovery nearly 40 years ago, PYP has attracted much interest and has become one of the most extensively studied proteins found in nature. The method of time-resolved crystallography, pioneered by Keith Moffat, has successfully characterized intermediates in the PYP photocycle at near atomic resolution over 12 decades of time down to the sub-picosecond time scale, allowing one to stitch together a movie and literally watch a protein as it functions. But how close to reality is this movie? To address this question, results from numerous complementary time-resolved techniques including x-ray crystallography, x-ray scattering, and spectroscopy are discussed. Emerging from spectroscopic studies is a general consensus that three time constants are required to model the excited state relaxation, with a highly strained ground-state cis intermediate formed in less than 2.4 ps. Persistent strain drives the sequence of structural transitions that ultimately produce the signaling state. Crystal packing forces produce a restoring force that slows somewhat the rates of interconversion between the intermediates. Moreover, the solvent composition surrounding PYP can influence the number and structures of intermediates as well as the rates at which they interconvert. When chloride is present, the PYP photocycle in a crystal closely tracks that in solution, which suggests the epic movie of the PYP photocycle is indeed based in reality.

Beyond slow two-state protein conformational exchange using CEST: applications to three-state protein interconversion on the millisecond timescale.

Although NMR spectroscopy is routinely used to study the conformational dynamics of biomolecules, robust analyses of the data are challenged in cases where exchange is more complex than two-state, such as when a 'visible' major conformer exchanges with two 'invisible' minor states on the millisecond timescale. It is becoming increasingly clear that chemical exchange saturation transfer (CEST) NMR experiments that were initially developed to study systems undergoing slow interconversion are also sensitive to intermediate-fast timescale biomolecular conformational exchange. Here we investigate the utility of the amide 15N CEST experiment to characterise protein three-state exchange occurring on the millisecond timescale by studying the interconversion between the folded (F) state of the FF domain from human HYPA/FBP11 (WT FF) and two of its folding intermediates I1 and I2. Although 15N CPMG experiments are consistent with the F state interconverting with a single minor state on the millisecond timescale, 15N CEST data clearly establish an exchange process between F and a pair of minor states. A unique three-state exchange model cannot be obtained by analysis of 15N CEST data recorded at a single temperature. However, including the relative sign of the difference in the chemical shifts of the two minor states based on a simple two-state analysis of CEST data recorded at multiple temperatures, results in a robust three-state model in which the F, I1 and I2 states interconvert with each other on the millisecond timescale ( ~ 550s-1, ~ 1200s-1, ~ 5000s-1), with I1 and I2 sparsely populated at ~ 0.15% and ~ 0.35%, respectively, at 15°C. A computationally demanding grid-search of exchange parameter space is not required to extract the best-fit exchange parameters from the CEST data. The utility of the CEST experiment, thus, extends well beyond studies of conformers in slow exchange on the NMR chemical shift timescale, to include systems with interconversion rates on the order of thousands/second.

Use of dislodgeable foliar residue studies for the evaluation of isomerization potential of active substances for non-dietary risk assessment.

When assessing the non-dietary risks for operators, workers, bystanders and residents to active substances in plant protection products (PPPs) that contain stereoisomers, the EFSA guidance on stereoisomers recommends the application of uncertainty factors when the initial ratio of stereoisomers undergoes a significant stereoisomeric excess change of more than 10%. This precautionary approach may be over-conservative in cases where the change in isomers is due to differences in their degradation rates rather than conversion of one isomer to a potentially more toxic isomer. Therefore, the impact of isomeric composition in non-dietary risk assessments of PPPs was evaluated, with particular emphasis on dislodgeable foliar residues (DFR) data and potential enantioselective degradation processes. Hypothetical outcomes are discussed as well as an evaluation of data for 5 compounds from a total of 35 DFR studies conducted under GLP, aimed at deriving DFR values for higher tier assessments. The findings indicate that possible isomerization of the active substance on leaves can be detected in DFR studies, which is essential for assessing risks to workers and residents. The results of this evaluation highlight the potential for over-conservatism of the current EFSA method for calculating uncertainty factors in non-dietary risk assessments. The EFSA method's reliance on worst-case assumptions, coupled with the lack of a formal statistical basis, can lead to overestimation of exposure, as illustrated in our theoretical considerations and supported by empirical data from DFR studies. To address this, we propose an alternative approach to evaluating risk considering both degradation and interconversion rates.

Relevance of the calculation of the diffusion coefficient in a capillary electrophoresis experiment.

This paper analyzes the role of the diffusion coefficient in the movement of analytes that can reversibly react with a selector given a product in the presence of drift. The problem mimics the movement of enantiomers in a capillary electrophoresis experiment. As is well known, the signal in the capillary must be sharp enough to make a good determination of the effective mobility of the analytes being analyzed. The essence of the technique is based on fast interconversion rates. Therefore, the effective diffusion coefficient must be negligible during the experiment. In the present work, an exact expression for both the apparent mobility and the diffusion coefficient is obtained. This is done by writing the rate equations governing the process and solving them using the generating function technique. The effective mobility coincides with the Wren and Rowe equation, whereas the diffusion coefficient allows us to determine the values of the parameters to be taken into account so that this quantity is minimal or close to zero. On the other hand, the numerical solution of the kinetic equations and Monte Carlo simulations allow us to follow the signal in the capillary and to determine its space-time evolution.

Stapling strategy for slowing helicity interconversion of α-helical peptides and isolating chiral auxiliary-free one-handed forms

In nature, α-helical peptides adopt right-handed conformations that are dictated by L-amino acids. Isolating one-handed α-helical peptides composed of only achiral components remains a significant challenge. Here, this goal is achieved by optical resolution of the corresponding racemic (quasi-)static α-helical peptide with double stapling, which effectively freezes the interconversion between the right-handed (P)- and left-handed (M)-α-helices. An as-obtained doubly stapled analogue having an unprotected L-valine residue at the C-terminus transforms from a kinetically trapped (M)-α-helix to a thermodynamically stable (P)-α-helix upon heating. In contrast, the corresponding singly stapled α-helical peptide undergoes an acid/base-triggered and solvent-induced reversible inversion of its preferred helicity within minutes. The interconversion rates of the singly and doubly stapled α-helical peptide foldamers are approximately 106 and 1012 times slower, respectively, than that of a non-stapled dynamic helical peptide. Therefore, the enantiopure doubly-stapled (quasi-)static α-helical peptide would retain its optical activity for several years at 25 °C.

Kinetic Resolution of the Atomic 3D Structures Formed by Ground and Excited Conformational States in an RNA Dynamic Ensemble.

Knowing the 3D structures formed by the various conformations populating the RNA free-energy landscape, their relative abundance, and kinetic interconversion rates is required to obtain a quantitative and predictive understanding of how RNAs fold and function at the atomic level. While methods integrating ensemble-averaged experimental data with computational modeling are helping define the most abundant conformations in RNA ensembles, elucidating their kinetic rates of interconversion and determining the 3D structures of sparsely populated short-lived RNA excited conformational states (ESs) remains challenging. Here, we developed an approach integrating Rosetta-FARFAR RNA structure prediction with NMR residual dipolar couplings and relaxation dispersion that simultaneously determines the 3D structures formed by the ground-state (GS) and ES subensembles, their relative abundance, and kinetic rates of interconversion. The approach is demonstrated on HIV-1 TAR, whose six-nucleotide apical loop was previously shown to form a sparsely populated (∼13%) short-lived (lifetime ∼ 45 μs) ES. In the GS, the apical loop forms a broad distribution of open conformations interconverting on the pico-to-nanosecond time scale. Most residues are unpaired and preorganized to bind the Tat-superelongation protein complex. The apical loop zips up in the ES, forming a narrow distribution of closed conformations, which sequester critical residues required for protein recognition. Our work introduces an approach for determining the 3D ensemble models formed by sparsely populated RNA conformational states, provides a rare atomic view of an RNA ES, and kinetically resolves the atomic 3D structures of RNA conformational substates, interchanging on time scales spanning 6 orders of magnitude, from picoseconds to microseconds.

Optimizing molecular potential models by imposing kinetic constraints with path reweighting.

Empirical force fields employed in molecular dynamics simulations of complex systems are often optimized to reproduce experimentally determined structural and thermodynamic properties. In contrast, experimental knowledge about the interconversion rates between metastable states in such systems is hardly ever incorporated in a force field due to a lack of an efficient approach. Here, we introduce such a framework based on the relationship between dynamical observables, such as rate constants, and the underlying molecular model parameters using the statistical mechanics of trajectories. Given a prior ensemble of molecular dynamics trajectories produced with imperfect force field parameters, the approach allows for the optimal adaption of these parameters such that the imposed constraint of equally predicted and experimental rate constant is obeyed. To do so, the method combines the continuum path ensemble maximum caliber approach with path reweighting methods for stochastic dynamics. When multiple solutions are found, the method selects automatically the combination that corresponds to the smallest perturbation of the entire path ensemble, as required by the maximum entropy principle. To show the validity of the approach, we illustrate the method on simple test systems undergoing rare event dynamics. Next to simple 2D potentials, we explore particle models representing molecular isomerization reactions and protein-ligand unbinding. Besides optimal interaction parameters, the methodology gives physical insights into what parts of the model are most sensitive to the kinetics. We discuss the generality and broad implications of the methodology.

Machine learning of atomic dynamics and statistical surface identities in gold nanoparticles

It is known that metal nanoparticles (NPs) may be dynamic and atoms may move within them even at fairly low temperatures. Characterizing such complex dynamics is key for understanding NPs’ properties in realistic regimes, but detailed information on, e.g., the stability, survival, and interconversion rates of the atomic environments (AEs) populating them are non-trivial to attain. In this study, we decode the intricate atomic dynamics of metal NPs by using a machine learning approach analyzing high-dimensional data obtained from molecular dynamics simulations. Using different-shape gold NPs as a representative example, an AEs’ dictionary allows us to label step-by-step the individual atoms in the NPs, identifying the native and non-native AEs and populating them along the MD simulations at various temperatures. By tracking the emergence, annihilation, lifetime, and dynamic interconversion of the AEs, our approach permits estimating a “statistical equivalent identity” for metal NPs, providing a comprehensive picture of the intrinsic atomic dynamics that shape their properties.

Supramolecular assembly guided by photolytic redox cycling

AbstractIn living systems, the formation of structures relies on balancing kinetic and thermodynamic influences powered by reversible covalent bond chemistry. Although synthetic efforts have replicated these processes to some extent, elucidating their combination is necessary to identify mechanisms that confer nature’s structural precision and flexibility within a complex environment. Here we design a photolytic reaction cascade where competing redox pathways control the transience, interconversion and production rates of thiol/disulfide supramolecular monomers in situ. In contrast to direct assembly by dissolution, cascade generation of the same monomers formed hierarchical assemblies with different structural order. Redox-induced cycling between thiol–disulfide formation led to the emergence of new secondary structures and chirality within the final assemblies. These multiple structural states found within the same molecular system demonstrate the concept of assembly plasticity engaged frequently in biology. We demonstrate the importance of reaction complexity in controlling supramolecular propagation and in expanding the library of nanoarchitectures that can be created.

A pyrazolone-based dinuclear Cu(II) Schiff-base complex: DFT studies on the rate-determining steps of the tautomerism in the ligand and molecular docking modelling with COVID-19 main protease (6LU7)

A new dinuclear Schiff-base complex, [Cu2L'2(μ-OAc)2], was synthesized and characterized by various spectroscopic methods and single-crystal X-ray diffraction (SCXRD). The symmetrical H2L ligand was synthesized from the reaction between 4-acetyl-3-methyl-1-phenyl-2-pyrazoline-5-one (or acylpyrazolone) with 1,3-propanediamine, but upon reaction with a metal salt, it gave the unsymmetrical HL'. The transition states and rates of interconversion between different tautomeric forms of HL' were studied by DFT computations. The imine-one (Z) in the gas phase and imine-ol (Z) in the solution were the most stable forms, most probably due to better hydrogen bonding. We also applied molecular docking modeling on the various tautomeric forms of the ligand as well as the complex with the main protease (6LU7) of the SARS-CoV-2 virus. Based on the obtained results, the complex had better interactions with the receptor. The order of the interactions was [Cu2L2(μ-OAc)2] > amine-one (Z) > imine-one (E) > imine-ol (Z) > amine-one (E) > imine-one (Z) > imine-ol (E).

Atropisomerism in Drug Discovery: A Medicinal Chemistry Perspective Inspired by Atropisomeric Class I PI3K Inhibitors.

Atropisomerism is a type of axial chirality resulting from hindered rotation about a σ bond that gives rise to nonsuperimposable stereoisomers (termed "atropisomers"). The inversion of chirality of an atropisomeric axis is a time- and temperature-dependent dynamic process occurring by simple bond rotation. For this reason, the rotational energy barrier (ΔErot) and the interconversion rate between an atropisomeric pair of biologically active molecules are important parameters to consider in drug discovery.Many compounds with atropisomeric axes advance into development every year. The vast majority of them have low rotational energy barriers (ΔErot lower than 20 kcal/mol), meaning they are rapidly equilibrating conformers and considered achiral (class 1 atropisomers). Compounds in class 2 (ΔErot = 20 to 30 kcal/mol) can be challenging to develop given that the stereochemical integrity of the atropisomeric axes can be compromised over time. It has been recommended that small molecule drug candidates containing one or more atropisomeric axes with rotational energy barriers greater than 30 kcal/mol (class 3 atropisomers) should be developed as single atropisomers.In medicinal chemistry, a σ bond with restricted rotation is engineered into a bioactive molecule primarily to limit its number of accessible conformations, thereby minimizing entropic and/or enthalpic energy penalties associated with biological target binding. In addition to enhanced pharmacology, potential positive outcomes of introducing atropisomerism include improved physicochemical properties and superior pharmacokinetics/ADME profiles. The application of atropisomerism in medicinal chemistry has become increasingly enabled due to recent advances in synthesis, purification, and analysis, as described in this special issue and recent review articles.Herein, we discuss two case studies from our own work in which restricting rotation about axes of atropisomerism led to significant improvements in pharmacological, physicochemical, and ADME properties for different series of PI3K inhibitors. In the first instance, a restricted axis of rotation was designed to mitigate an acid-mediated hydrolytic degradation pathway observed in a series of PI3Kδ inhibitors. The conformational constraint disrupts conjugation between a quinazolinone and a pyridine, leading to improved chemical stability under acidic conditions. In the second case study, introduction of a restricted axis of rotation between two heteroaromatic systems in a series of PI3Kβ inhibitors generated pairs of atropisomeric compounds with significantly different biological activities. Advanced profiling also demonstrated clear substrate stereospecificity in regard to metabolism by aldehyde oxidase. Gratifyingly, the eutomer (more active atropisomer) shows significantly less susceptibility for oxidative metabolism relative to the distomer (less active atropisomer). The improvements in potency, selectivity, chemical stability, and metabolic stability discussed in this manuscript are all directly related to the concept of atropisomerism.

Well-Tempered Metadynamics Simulations Predict the Structural and Dynamic Properties of a Chiral 24-Atom Macrocycle in Solution.

Inspired by therapeutic potential, the molecular engineering of macrocycles is garnering increased interest. Exercising control with design, however, is challenging due to the dynamic behavior that these molecules must demonstrate in order to be bioactive. Herein, the value of metadynamics simulations is demonstrated: the free-energy surfaces calculated reveal folded and flattened accessible conformations of a 24-atom macrocycle separated by barriers of ∼6 kT under experimentally relevant conditions. Simulations reveal that the dominant conformer is folded-an observation consistent with a solid-state structure determined by X-ray crystallography and a network of rOes established by 1H NMR. Simulations suggest that the macrocycle exists as a rapidly interconverting pair of enantiomeric, folded structures. Experimentally, 1H NMR shows a single species at room temperature. However, at lower temperature, the interconversion rate between these enantiomers becomes markedly slower, resulting in the decoalescence of enantiotopic methylene protons into diastereotopic, distinguishable resonances due to the persistence of conformational chirality. The emergence of conformational chirality provides critical experimental support for the simulations, revealing the dynamic nature of the scaffold-a trait deemed critical for oral bioactivity.