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

Impact of Solvent and Protonation State on Rotational Barriers in [s]-Triazines

Amine-substituted [s]-triazines display hindered rotation around the triazine-N bond. While this barrier, ΔG‡, has been measured to be between 15.1 and 17.7 kcal/mol for neutral triazines, the impacts that solvent and protonation state have not been addressed. Using a dimethylamine substituent as a reporter, ΔG‡ was measured to be 17.5-19.3 kcal/mol upon protonation across a range of solvents (D2O, DMSO-d6, MeCN-d3, MeOD-d4, tetrahydrofuran-d8, trifluoroethanol-d3). Furthermore, ΔG‡ increases as the solvent dielectric decreases (p < 0.01). This trend is consistent with the role that solvent plays in stabilizing the increased charge density on the triazine ring resulting from a loss of conjugation with the dimethylamine substituent. Across these solvents, ΔG‡ for the neutral molecule is smaller by ∼2-3 kcal/mol, ranging from 15.3-16.1 kcal/mol. In pyridine, ΔG‡ does not correlate with the solvent dielectric for the "protonated" model. The lower barrier is attributed to competitive protonation: the pKa of the protonated triazine (∼6) is similar to that of protonated pyridine-d5 (5.8). As additional acid is added, ΔG‡ increases. Adding additional acid to the protonated model in D2O or DMSO-d6 does not significantly affect ΔG‡.

Toward physics-based precision medicine: Exploiting protein dynamics to design new therapeutics and interpret variants

The goal of precision medicine is to utilize our knowledge of the molecular causes of disease to better diagnose and treat patients. However, there is a substantial mismatch between the small number of food and drug administration (FDA)-approved drugs and annotated coding variants compared to the needs of precision medicine. This review introduces the concept of physics-based precision medicine, a scalable framework that promises to improve our understanding of sequence-function relationships and accelerate drug discovery. We show that accounting for the ensemble of structures a protein adopts in solution with computer simulations overcomes many of the limitations imposed by assuming a single protein structure. We highlight studies of protein dynamics and recent methods for the analysis of structural ensembles. These studies demonstrate that differences in conformational distributions predict functional differences within protein families and between variants. Thanks to new computational tools that are providing unprecedented access to protein structural ensembles, this insight may enable accurate predictions of variant pathogenicity for entire libraries of variants. We further show that explicitly accounting for protein ensembles, with methods like alchemical free energy calculations or docking to Markov state models, can uncover novel lead compounds. To conclude, we demonstrate that cryptic pockets, or cavities absent in experimental structures, provide an avenue to target proteins that are currently considered undruggable. Taken together, our review provides a roadmap for the field of protein science to accelerate precision medicine.

Proton Affinity and Conformational Integrity of a 24-Atom Triazine Macrocycle across Physiologically Relevant pH

For 24-atom triazine macrocycles, protonation of the heterocycle leads to a rigid, folded structure presenting a network of hydrogen bonds. These molecules derive from dynamic covalent chemistry wherein triazine monomers bearing a protected hydrazine group and acetal tethered by the amino acid dimerize quantitatively in an acidic solution. Here, lysine is used, and the product is a tetracation. The primary amines of the lysine side chains do not interfere with quantitative yields of the desired bis(hydrazone) at concentrations of 5-125 mg/mL. Mathematical modeling of data derived from titration experiments of the macrocycle reveals that the pKa values of the protonated triazines are 5.6 and 6.7. Changes in chemical shifts of resonances in the 1H NMR spectra corroborate these values and further support assignment of the protonation sites. The pKa values of the lysine side chains are consistent with expectation. Upon deprotonation, the macrocycle enjoys greater conformational freedom as evident from the broadening of resonances in the 1H and 13C NMR spectra indicative of dynamic motion on the NMR time scale and the appearance of additional conformations at room temperature. While well-tempered metadynamics suggests only a modest difference in accessible conformational footprints of the protonated and deprotonated macrocycles, the shift in conformation(s) supports the stabilizing role that the protons adopt in the hydrogen-bonded network.

Addressing Challenges of Macrocyclic Conformational Sampling in Polar and Apolar Solvents: Lessons for Chameleonicity

We evaluated a workflow to reliably sample the conformational space of a set of 47 peptidic macrocycles. Starting from SMILES strings, we use accelerated molecular dynamics simulations to overcome high energy barriers, in particular, the cis-trans isomerization of peptide bonds. We find that our approach performs very well in polar solvents like water and dimethyl sulfoxide. Interestingly, the protonation state of a secondary amine in the ring only slightly influences the conformational ensembles of our test systems. For several of the macrocycles, determining the conformational distribution in chloroform turns out to be considerably more challenging. Especially, the choice of partial charges crucially influences the ensembles in chloroform. We address these challenges by modifying initial structures and the choice of partial charges. Our results suggest that special care has to be taken to understand the configurational distribution in apolar solvents, which is a key step toward a reliable prediction of membrane permeation of macrocycles and their chameleonic properties.

Supramolecular Polymerization of a Pyrene-Substituted Diamide and Its Ensemble of Kinetically Trapped Configurations

The prevalence of kinetically accessible states in supramolecular polymerization pathways has been exploited to control the growth of the polymer and thereby to obtain niche morphologies. Yet, these pathways themselves are not easily amenable for experimental delineation but could potentially be understood through molecular dynamics (MD) simulations. Herein, we report an extensive investigation of the self-assembly of pyrene-substituted diamide (PDA) monomers in solution, conducted using atomistic MD simulations and advanced sampling methods. We characterize such kinetic and thermodynamic states as well as the transition pathways and free energy barriers between them. PDA forms a dimeric segment with the N- to C-termini vectors of the diamide moieties arranged either in parallel or anti-parallel fashion. This characteristic, combined with the molecule's torsional flexibility and pyrene-solvent interactions, presents an ensemble of molecular configurations contributing to the kinetic state in the polymerization pathway. While this ensemble primarily comprises short oligomers containing a mix of anti-parallel and parallel dimeric segments, the thermodynamic state of the assembly is a right-handed polymer featuring parallel ones only. Our work thus offers an approach by which the landscape of any specific supramolecular polymerization can be deconstructed.

Adaptation of Empirical Methods to Predict the LogD of Triazine Macrocycles

Octanol/water partition coefficients guide drug design, but algorithms do not always accurately predict these values. For cationic triazine macrocycles that adopt a conserved folded shape in solution, common algorithms fall short. Here, the logD values for 12 macrocycles differing in amino acid choice were predicted and then measured experimentally. On average, AlogP, XlogP, and ChemAxon predictions deviate by 0.9, 2.8, and 3.9 log units, with XlogP overestimating lipophilicity and AlogP and ChemAxon underestimating lipophilicity. Importantly, however, a linear relationship (R2 > 0.98) exists between the values predicted by AlogP and the experimentally determined logD values, thus enabling more accurate predictions.

Controlling Swing Rates in Macrocyclic Molecular Mortise Hinges

Hinge motion is observed in macrocyclic, mortise-type molecular hinges using variable temperature NMR spectroscopy. The data is consistent with dynamic hinging from a folded-to-extended-to-folded enantiomeric state. Crystallographic and solution structures of the folded states are reported. Chemical shift predictions derived from crystallographic data corroborate fully revolute hinge motion. The rate of hinging is affected by steric congestion at the hinge axis. A macrocycle containing glycine, 1, hinges faster than one comprising aminoisobutyric acid, 2. The free energies of activation, ΔG≠ , for 1 and 2 were determined to be 13.3±0.3 kcal/mol and 16.3±0.3 kcal/mol, respectively. This barrier is largely independent of solvent across those surveyed (CD3 OD, CD3 CN, DMSO-d6 , pyridine-d5 , D2 O). Experiment and computation predict energy barriers that are consistent with disruption of an intramolecular network of hydrogen bonds. DFT calculations reveal a pathway for hinge motion.

A Model for the Rapid Assessment of Solution Structures for 24-Atom Macrocycles: The Impact of β-Branched Amino Acids on Conformation

Experiment and computation are used to develop a model to rapidly predict solution structures of macrocycles sharing the same Murcko framework. These 24-atom triazine macrocycles result from the quantitative dimerization of identical monomers presenting a hydrazine group and an acetal tethered to an amino acid linker. Monomers comprising glycine and the β-branched amino acids threonine, valine, and isoleucine yield macrocycles G-G, T-T, V-V, and I-I, respectively. Elements common to all members of the framework include the efficiency of macrocyclization (quantitative), the solution- and solid-state structures (folded), the site of protonation (opposite the auxiliary dimethylamine group), the geometry of the hydrazone (E), the C2 symmetry of the subunits (conserved), and the rotamer state adopted. In aggregate, the data reveal metrics predictive of the three-dimensional solution structure that derive from the fingerprint region of the 1D 1H spectrum and a network of rOes from a single resonance. The metrics also afford delineation of more nuanced structural features that allow subpopulations to be identified among the members of the framework. Well-tempered metadynamics provides free energy surfaces and population distributions of these macrocycles. The areas of the free energy surface decrease with increasing steric bulk (G-G > V-V ∼ T-T > I-I). In addition, the surfaces are increasingly isoenergetic with decreasing steric bulk (G-G > V-V ∼ T-T > I-I).

Computational and Experimental Evidence for Templated Macrocyclization: The Role of a Hydrogen Bond Network in the Quantitative Dimerization of 24-Atom Macrocycles

In the absence of preorganization, macrocyclization reactions are often plagued by oligomeric and polymeric side products. Here, a network of hydrogen bonds was identified as the basis for quantitative yields of macrocycles derived from the dimerization of monomers. Oligomers and polymers were not observed. Macrocyclization, the result of the formation of two hydrazones, was hypothesized to proceed in two steps. After condensation to yield the monohydrazone, a network of hydrogen bonds formed to preorganize the terminal acetal and hydrazine groups for cyclization. Experimental evidence for preorganization derived from macrocycles and acyclic models. Solution NMR spectroscopy and single-crystal X-ray diffraction revealed that the macrocycles isolated from the cyclization reaction were protonated twice. These protons contributed to an intramolecular network of hydrogen bonds that engaged distant carbonyl groups to realize a long-range order. DFT calculations showed that this network of hydrogen bonds contributed 8.7 kcal/mol to stability. Acyclic models recapitulated this network in solution. Condensation of an acetal and a triazinyl hydrazine, which adopted a number of conformational isomers, yielded a hydrazone that adopted a favored rotamer conformation in solution. The critical hydrogen-bonded proton was also evident. DFT calculations of acyclic models showed that the rotamers were isoenergetic when deprotonated. Upon protonation, however, energies diverged with one low-energy rotamer adopting the conformation observed in the macrocycle. This conformation anchored the network of hydrogen bonds of the intermediate. Computation revealed that the hydrogen-bonded network in the acyclic intermediate contributed up to 14 kcal/mol of stability and preorganized the acetal and hydrazine for cyclization.