De novo design of drug-binding proteins with predictable binding energy and specificity.

The de novo design of small molecule-binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase-1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule-binding proteins with tuned interaction energies is feasible entirely from computation.

Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics.

Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.

Electron transfer rate modulation with mid-IR in butadiyne-bridged donor-bridge-acceptor compounds.

Controlling electron transfer (ET) processes in donor-bridge-acceptor (DBA) compounds by mid-IR excitation can enhance our understanding of the ET dynamics and may find practical applications in molecular sensing and molecular-scale electronics. Alkyne moieties are attractive to serve as ET bridges, as they offer the possibility of fast ET and present convenient vibrational modes to perturb the ET dynamics. Yet, these bridges introduce complexity because of the strong torsion angle dependence of the ET rates and transition dipoles among electronic states and a shallow torsion barrier. In this study, we implemented ultrafast 3-pulse laser spectroscopy to investigate how the ET from the dimethyl aniline (D) electron donor to the N-isopropyl-1,8-napthalimide (NAP) electron acceptor can be altered by exciting the CC stretching mode (νCC) of the butadiyne bridge linking the donor and acceptor. The electron transfer was initiated by electronically exciting the acceptor moiety at 400 nm, followed by vibrational excitation of the alkyne, νCC, and detecting the changes in the absorption spectrum in the visible spectral region. The experiments were performed at different delay times t1 and t2, which are the delays between UV-mid-IR and mid-IR-Vis pulses, respectively. Two sets of torsion-angle conformers were identified, one featuring a very fast mean ET time of 0.63 ps (group A) and another featuring a slower mean ET time of 4.3 ps (group B), in the absence of the mid-IR excitation. TD-DFT calculations were performed to determine key torsion angle dependent molecular parameters, including the electronic and vibrational transition dipoles, transition frequencies, and electronic couplings. To describe the 3-pulse data, we developed a kinetic model that includes a locally excited, acceptor-based S2 state, a charge separated S1 state, and their vibrationally excited counterparts, with either excited νCC (denoted as S1Atr, S1Btr, S2Atr, and S2Btr, where tr stands for the excited triplet bond, νCC) or excited daughter modes of the νCC relaxation (S1Ah, S1Bh, S2Ah, and S2Bh, where h stands for vibrationally hot species). The kinetic model was solved analytically, and the species-associated spectra (SAS) were determined numerically using a matrix approach, treating first the experiments with longer t1 delays and then using the already determined SAS for modeling the experiments with shorter t1 delays. Strong vibronic coupling of νCC and of vibrationally hot states makes the analysis complicated. Nevertheless, the SAS were identified and the ET rates of the vibrationally excited species, S2Atr, S2Btr and S2Bh, were determined. The results show that the ET rate for the S2A species is ca. 1.2-fold slower when the νCC mode is excited. The ET rate for species S2B is slower by ca. 1.3-fold if the compound is vibrationally hot and is essentially unchanged when the νCC mode is excited. The SAS determined for the tr and h species resemble the SAS for their respective precursor species in the 2-pulse transient absorption experiments, which validates the procedure used and the results.

Intra-cluster Charge Migration upon Hydration of Protonated Formic Acid Revealed by Anharmonic Analysis of Cold Ion Vibrational Spectra.

The rates of many chemical reactions are accelerated when carried out in micron-sized droplets, but the molecular origin of the rate acceleration remains unclear. One example is the condensation reaction of 1,2-diaminobenzene with formic acid to yield benzimidazole. The observed rate enhancements have been rationalized by invoking enhanced acidity at the surface of methanol solvent droplets with low water content to enable protonation of formic acid to generate a cationic species (protonated formic acid or PFA) formed by attachment of a proton to the neutral acid. Because PFA is the key feature in this reaction mechanism, vibrational spectra of cryogenically cooled, microhydrated PFA·(H2O)n=1-6 were acquired to determine how the extent of charge localization depends on the degree of hydration. Analysis of these highly anharmonic spectra with path integral ab initio molecular dynamics simulations reveals the gradual displacement of the excess proton onto the water network in the microhydration regime at low temperatures with n = 3 as the tipping point for intra-cluster proton transfer.

De novo design of drug-binding proteins with predictable binding energy and specificity.

The de novo design of small-molecule-binding proteins has seen exciting recent progress; however, the ability to achieve exquisite affinity for binding small molecules while tuning specificity has not yet been demonstrated directly from computation. Here, we develop a computational procedure that results in the highest affinity binders to date with predetermined relative affinities, targeting a series of PARP1 inhibitors. Two of four designed proteins bound with affinities ranging from < 5 nM to low µM, in a predictable manner. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free-energy calculations performed directly on the designed models are in excellent agreement with the experimentally measured affinities, suggesting that the de novo design of small-molecule-binding proteins with tuned interaction energies is now feasible entirely from computation. We expect these methods to open many opportunities in biomedicine, including rapid sensor development, antidote design, and drug delivery vehicles.

Synthetic Control of Exciton Dynamics in Bioinspired Cofacial Porphyrin Dimers.

Understanding how the complex interplay among excitonic interactions, vibronic couplings, and reorganization energy determines coherence-enabled transport mechanisms is a grand challenge with both foundational implications and potential payoffs for energy science. We use a combined experimental and theoretical approach to show how a modest change in structure may be used to modify the exciton delocalization, tune electronic and vibrational coherences, and alter the mechanism of exciton transfer in covalently linked cofacial Zn-porphyrin dimers (meso-beta linked ABm-ß and meso-meso linked AAm-m). While both ABm-ß and AAm-m feature zinc porphyrins linked by a 1,2-phenylene bridge, differences in the interporphyrin connectivity set the lateral shift between macrocycles, reducing electronic coupling in ABm-ß and resulting in a localized exciton. Pump-probe experiments show that the exciton dynamics is faster by almost an order of magnitude in the strongly coupled AAm-m dimer, and two-dimensional electronic spectroscopy (2DES) identifies a vibronic coherence that is absent in ABm-ß. Theoretical studies indicate how the interchromophore interactions in these structures, and their system-bath couplings, influence excitonic delocalization and vibronic coherence-enabled rapid exciton transport dynamics. Real-time path integral calculations reproduce the exciton transfer kinetics observed experimentally and find that the linking-modulated exciton delocalization strongly enhances the contribution of vibronic coherences to the exciton transfer mechanism, and that this coherence accelerates the exciton transfer dynamics. These benchmark molecular design, 2DES, and theoretical studies provide a foundation for directed explorations of nonclassical effects on exciton dynamics in multiporphyrin assemblies.

Electron ratcheting in self-assembled soft matter.

Ratcheted multi-step hopping electron transfer systems can plausibly produce directional charge transport over very large distances without requiring a source-drain voltage bias. We examine molecular strategies to realize ratcheted charge transport based on multi-step charge hopping, and we illustrate two ratcheting mechanisms with examples based on DNA structures. The charge transport times and currents that may be generated in these assemblies are also estimated using kinetic simulations. The first ratcheting mechanism described for nanoscale systems requires local electric fields on the 109 V/m scale to realize nearly 100% population transport. The second ratcheting mechanism for even larger systems, based on electrochemical gating, is estimated to generate currents as large as 0.1 pA for DNA structures that are a few µm in length with a gate voltage of about 5 V, a magnitude comparable to currents measured in DNA wires at the nanoscale when a source-drain voltage bias of similar magnitude is applied, suggesting an approach to considerably extend the distance range over which DNA charge transport devices may operate.

Delocalization-Assisted Transport through Nucleic Acids in Molecular Junctions.

The flow of charge through molecules is central to the function of supramolecular machines, and charge transport in nucleic acids is implicated in molecular signaling and DNA repair. We examine the transport of electrons through nucleic acids to understand the interplay of resonant and nonresonant charge carrier transport mechanisms. This study reports STM break junction measurements of peptide nucleic acids (PNAs) with a G-block structure and contrasts the findings with previous results for DNA duplexes. The conductance of G-block PNA duplexes is much higher than that of the corresponding DNA duplexes of the same sequence; however, they do not display the strong even-odd dependence conductance oscillations found in G-block DNA. Theoretical analysis finds that the conductance oscillation magnitude in PNA is suppressed because of the increased level of electronic coupling interaction between G-blocks in PNA and the stronger PNA-electrode interaction compared to that in DNA duplexes. The strong interactions in the G-block PNA duplexes produce molecular conductances as high as 3% G0, where G0 is the quantum of conductance, for 5 nm duplexes.

Nonexponential Length Dependence of Molecular Conductance in Acene-Based Molecular Wires.

In the nonresonant regime, molecular conductance decays exponentially with distance, limiting the fabrication of efficient molecular semiconductors at the nanoscale. In this work, we calculate the conductance of a series of acene derivatives connected to gold electrodes using density functional theory (DFT) combined with the nonequilibrium Green's function (NEGF) formalism. We show that these systems have near length-independent conductance and can exhibit a conductance increase with molecular length depending on the connection to the electrodes. The analysis of the molecular orbital energies and transmission functions attribute this behavior to the dramatic decrease of the highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gap with length, which shifts the transmission peaks near the Fermi level. These results demonstrate that the anchoring configuration determines the conductance behavior of acene derivatives, which are optimal building blocks to fabricate single-molecule devices with tunable charge transport properties.

Quantum interference enhances rectification behavior of molecular devices.

We present a theoretical and computational work and demonstrate that cross-conjugated molecules with electron-donating groups are efficient rectifiers with high conductance. The rectification ratios obtained are up to one order of magnitude at an applied bias voltage of 0.3 V. The use of electron-withdrawing groups to form donor-bridge-acceptor triads gives rectification ratios of the order of 102. We found that the high rectification results from localizing the Highest Occupied Molecular Orbital (HOMO) at one end of the molecular device. When the HOMO is localized, quantum interference effects substantially enhance rectification. Our observations rely on transport calculations of linearly-conjugated and cross-conjugated molecules using Non-Equilibrium Green's Function Technique and Density Functional Theory (NEGF-DFT). Analysis of transmission functions suggests a dependency of the rectification ratio on the anti-resonance peak position near the Fermi level of the electrode, allowing the possibility to modulate molecular rectification through electrochemical gating.

Symmetry controlled photo-selection and charge separation in butadiyne-bridged donor-bridge-acceptor compounds.

Electron transfer (ET) in donor-bridge-acceptor (DBA) compounds depends strongly on the structural and electronic properties of the bridge. Among the bridges that support donor-acceptor conjugation, alkyne bridges have attractive and unique properties: they are compact, possess linear structure permitting access to high symmetry DBA molecules, and allow torsional motion of D and A, especially for longer bridges. We report conformation dependent electron transfer dynamics in a set of novel DBA compounds featuring butadiyne (C4) bridge, N-isopropyl-1,8-napthalimide (NAP) acceptors, and donors that span a range of reduction potentials (trimethyl silane (Si-C4-NAP), phenyl (Ph-C4-NAP), and dimethyl aniline (D-C4-NAP)). Transient mid-IR absorption spectra of the C[triple bond, length as m-dash]C bridge stretching modes, transient spectra in the visible range, and TD-DFT calculations were used to decipher the ET mechanisms. We found that the electronic excited state energies and, especially, the transition dipoles (S0 → Sn) depend strongly on the dihedral angle (θ) between D and A and the frontier orbital symmetry, offering an opportunity to photo-select particular excited states with specific ranges of dihedral angles by exciting at chosen wavelengths. For example, excitation of D-C4-NAP at 400 nm predominantly prepares an S1 excited state in the planar conformations (θ ∼ 0) but selects an S2 state with θ ∼ 90°, indicating the dominant role of the molecular symmetry in the photophysics. Moreover, the symmetry of the frontier orbitals of such DBA compounds not only defines the photo-selection outcome, but also determines the rate of the S2 → S1 charge separation reaction. Unprecedented variation of the S2-S1 electronic coupling with θ by over four orders of magnitude results in slow ET at θ ca. 0° and 90° but extremely fast ET at θ of 20-60°. The unique features of high-symmetry alkyne bridged DBA structures enable frequency dependent ET rate selection and make this family of compounds promising targets for the vibrational excitation control of ET kinetics.

Unsymmetrical Bis-Alkynyl Complexes Based on Co(III)(cyclam): Synthesis, Ultrafast Charge Separation, and Analysis.

Donor-bridge-acceptor (D-B-A) systems with a polarizable bridge can afford rapid photoinduced electron transfer dynamics that may be susceptible to rate modulation by infrared excitation. We describe the synthesis, characterization, and electronic structure of a class of readily assembled D-B-A structures linked by a cobalt cyclam bridge. The reaction between [Co(cyclam)Cl2]Cl and 4-ethynyl-N-isopropyl-1,8-naphthalimide (HC2NAPiPr) yields [Co(cyclam)(C2NAPiPr)Cl]Cl (1), which reacts with LiC2Y at -78 °C to afford [Co(cyclam)(C2NAPiPr)(C2D)]Cl with D as C6H4-4-NMe2 (2a), NAPiPr (2b), Ph (2c), and C6H4-4-N(4-MeOPh)2 (2d). Molecular structures of 1 and 2a were established using single-crystal X-ray diffraction, while the redox properties and fluorescence profiles of compounds 1 and 2 were examined using voltammetric and steady-state emission techniques, respectively. The electronic structures and photophysical properties of these compounds were studied using density functional theory and time-dependent density functional theory methods. The excited-state dynamics of compounds 1, 2a, and 2d were explored using femtosecond transient absorption spectroscopy with 400 nm excitation and detection in both the visible and mid-IR spectral regions. Formation of a long-lived excited state was complete within 20 ps of excitation in all three compounds. Ultrafast spectral changes observed in 2a and 2d within the first 20 ps indicated the formation of a charge separated state (CS state, D+-B-A-) with characteristic times of less than 0.1 and 0.25 ps, respectively. The CS state undergoes rapid charge recombination (8 ps in 2a and 4 ps in 2d). The CS dynamics is facilitated by the Co center, which mixes the bright NAP-centered electronic state with a pure CS state. The mixing strength depends on the donor energetics and conformation, which significantly influences the charge transfer dynamics in 2a and 2d.