Molecular Dynamics Simulations Based on Newly Developed Force Field Parameters for Cu2+ Spin Labels Provide Insights into Double-Histidine-Based Double Electron-Electron Resonance

The chromatin remodeler SMARCA5 binds to d-block metal supports: Characterization of affinities by IMAC chromatography and QM analysis

The ISWI family protein SMARCA5 contains the ATP-binding pocket that coordinates the catalytic Mg2+ ion and water molecules for ATP hydrolysis. In this study, we demonstrate that SMARCA5 can also possess an alternative metal-binding ability. First, we isolated SMARCA5 on the cobalt column (IMAC) to near homogeneity. Examination of the interactions of SMARCA5 with metal-chelating supports showed that, apart from Co2+, it binds to Cu2+, Zn2+ and Ni2+. The efficiency of the binding to the last-listed metal was influenced by the chelating ligand, resulting in a strong preference for Ni-NTA over the Ni-CM-Asp equivalent. To gain insight in the preferential affinity for the Ni-NTA ligand, QM calculations were performed on model systems and metal-ligand complexes with a limited protein fragment of SMARCA5 containing the double-histidine (dHis) motif. The calculations correlated the observed affinity with the relative stability of the d-block metals to tetradentate ligand coordination over tridentate, as well as their overall octahedral coordination capacity. Likewise, binding free energies derived from model imidazole complexes mirrored the observed Ni-NTA/Ni-CM-Asp preferential affinity. Finally, similar calculations on complexes with a SMARCA5 peptide fragment derived from the AlphaFold structural prediction, captured almost accurately the expected relative stability of the TM complexes, and produced a large energetic separation (~10 kcal∙mol-1) between Ni-NTA and Ni-CM-Asp in favour of the former.

Modeling of Cu(II)-based protein spin labels using rotamer libraries

The bifunctional spin label double-histidine copper-(II) capped with nitrilotriacetate [dHis-Cu(II)-NTA], used in conjunction with electron paramagnetic resonance (EPR) methods can provide high-resolution distance data for investigating protein structure and backbone conformational diversity. Quantitative utilization of this data is limited due to a lack of rapid and accurate dHis-Cu(II)-NTA modeling methods that can be used to translate experimental data into modeling restraints. Here, we develop two dHis-Cu(II)-NTA rotamer libraries using a set of recently published molecular dynamics simulations and a semi-empirical meta-dynamics-based conformational ensemble sampling tool for use with the recently developed chiLife bifunctional spin label modeling method. The accuracy of both the libraries and the modeling method are tested by comparing model predictions to experimentally determined distance distributions. We show that this method is accurate with absolute deviation between the predicted and experimental modes between 0.0-1.2 Å with an average of 0.6 Å over the test data used. In doing so, we also validate the generality of the chiLife bifunctional label modeling method. Taken together, the increased structural resolution and modeling accuracy of dHis-Cu(II)-NTA over other spin labels promise improvements in the accuracy and resolution of protein models by EPR.

Integrating Electron Paramagnetic Resonance Spectroscopy and Computational Modeling to Measure Protein Structure and Dynamics

Electron paramagnetic resonance (EPR) has become a powerful probe of conformational heterogeneity and dynamics of biomolecules. In this Review, we discuss different computational modeling techniques that enrich the interpretation of EPR measurements of dynamics or distance restraints. A variety of spin labels are surveyed to provide a background for the discussion of modeling tools. Molecular dynamics (MD) simulations of models containing spin labels provide dynamical properties of biomolecules and their labels. These simulations can be used to predict EPR spectra, sample stable conformations and sample rotameric preferences of label sidechains. For molecular motions longer than milliseconds, enhanced sampling strategies and de novo prediction software incorporating or validated by EPR measurements are able to efficiently refine or predict protein conformations, respectively. To sample large-amplitude conformational transition, a coarse-grained or an atomistic weighted ensemble (WE) strategy can be guided with EPR insights. Looking forward, we anticipate an integrative strategy for efficient sampling of alternate conformations by de novo predictions, followed by validations by systematic EPR measurements and MD simulations. Continuous pathways between alternate states can be further sampled by WE-MD including all intermediate states.

Differentiating between Label and Protein Conformers in Pulsed Dipolar EPR Spectroscopy with the dHis-Cu2+ (NTA) Motif

Pulsed dipolar EPR spectroscopy (PDS) in combination with site-directed spin labeling is a powerful tool in structural biology. However, the commonly used spin labels are conjugated to biomolecules via rather long and flexible linkers, which hampers the translation of distance distributions into biomolecular conformations. In contrast, the spin label copper(II)-nitrilotriacetic acid [Cu2+ (NTA)] bound to two histidines (dHis) is rigid and yields narrow distance distributions, which can be more easily translated into biomolecular conformations. Here, we use this label on the 71 kDa Yersinia outer protein O (YopO) to decipher whether a previously experimentally observed bimodal distance distribution is due to two conformations of the biomolecule or of the flexible spin labels. Two different PDS experiments, that is, pulsed electron-electron double resonance (PELDOR aka DEER) and relaxation-induced dipolar modulation enhancement (RIDME), yield unimodal distance distribution with the dHis-Cu2+ (NTA) motif; this result suggests that the α-helical backbone of YopO adopts a single conformation in frozen solution. In addition, we show that the Cu2+ (NTA) label preferentially binds to the target double histidine (dHis) sites even in the presence of 22 competing native histidine residues. Our results therefore suggest that the generation of a His-null background is not required for this spin labeling methodology. Together these results highlight the value of the dHis-Cu2+ (NTA) motif in PDS experiments.

Direct observation of negative cooperativity in a detoxification enzyme at the atomic level by Electron Paramagnetic Resonance spectroscopy and simulation

The catalytic activity of human glutathione S-transferase A1-1 (hGSTA1-1), a homodimeric detoxification enzyme, is dependent on the conformational dynamics of a key C-terminal helix α9 in each monomer. However, the structural details of how the two monomers interact upon binding of substrates is not well understood and the structure of the ligand-free state of the hGSTA1-1 homodimer has not been resolved. Here, we used a combination of electron paramagnetic resonance (EPR) distance measurements and weighted ensemble (WE) simulations to characterize the conformational ensemble of the ligand-free state at the atomic level. EPR measurements reveal a broad distance distribution between a pair of Cu(II) labels in the ligand-free state that gradually shifts and narrows as a function of increasing ligand concentration. These shifts suggest changes in the relative positioning of the two α9 helices upon ligand binding. WE simulations generated unbiased pathways for the seconds-timescale transition between alternate states of the enzyme, leading to the generation of atomically detailed structures of the ligand-free state. Notably, the simulations provide direct observations of negative cooperativity between the monomers of hGSTA1-1, which involve the mutually exclusive docking of α9 in each monomer as a lid over the active site. We identify key interactions between residues that lead to this negative cooperativity. Negative cooperativity may be essential for interaction of hGSTA1-1 with a wide variety of toxic substrates and their subsequent neutralization. More broadly, this work demonstrates the power of integrating EPR distances with WE rare-events sampling strategy to gain mechanistic information on protein function at the atomic level.

A novel approach to modeling side chain ensembles of the bifunctional spin label RX

We introduce a novel approach to modeling side chain ensembles of bifunctional spin labels. This approach utilizes rotamer libraries to generate side chain conformational ensembles. Because the bifunctional label is constrained by two attachment sites, the label is split into two monofunctional rotamers which are first attached to their respective sites, then rejoined by a local optimization in dihedral space. We validate this method against a set of previously published experimental data using the bifunctional spin label, RX. This method is relatively fast and can readily be used for both experimental analysis and protein modeling, providing significant advantages over modeling bifunctional labels with molecular dynamics simulations. Use of bifunctional labels for site directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy dramatically reduces label mobility, which can significantly improve resolution of small changes in protein backbone structure and dynamics. Coupling the use of bifunctional labels with side chain modeling methods allows for improved quantitative application of experimental SDSL EPR data to protein modeling.

"Store-bought is fine": Sensitivity considerations using shaped pulses for DEER measurements on Cu(II) labels

The narrow excitation bandwidth of monochromic pulses is a sensitivity limitation for pulsed dipolar spectroscopy on Cu(II)-based measurements. In response, frequency-swept pulses with large excitation bandwidths have been adopted to probe a greater range of the EPR spectrum. However, much of the work utilizing frequency-swept pulses in Cu(II) distance measurements has been carried out on home-built spectrometers and equipment. Herein, we carry out systematic Cu(II) based distance measurements to demonstrate the capability of chirp pulses on commercial instrumentation. More importantly we delineate sensitivity considerations under acquisition schemes that are necessary for robust distance measurements using Cu(II) labels for proteins. We show that a 200 MHz sweeping bandwidth chirp pulse can improve the sensitivity of long-range distance measurements by factors of three to four. The sensitivity of short-range distances only increases slightly due to special considerations for the chirp pulse duration relative to the period length of the modulated dipolar signal. Enhancements in sensitivity also dramatically reduce measurement collection times enabling rapid collection of orientationally averaged Cu(II) distance measurements in under two hours.

Efficient sampling of molecular orientations for Cu(II)-based DEER on protein labels

Combining rigid Cu(II) labels and pulsed-EPR techniques enables distance constraint measurements that are incisive probes of protein structure and dynamics. However, the labels can lead to a dipolar signal that is biased by the relative orientation of the two spins, which is typically unknown a priori in a bilabeled protein. This effect, dubbed orientational selectivity, becomes a bottleneck in measuring distances. This phenomenon also applies to other pulsed-EPR techniques that probe electron-nucleus interactions. In this work, we dissect orientational selectivity by generating an in silico sample of Cu(II)-labeled proteins to evaluate pulse excitation in the context of double electron-electron resonance (DEER) at Q-band frequencies. This approach enables the observation of the contribution of each protein orientation to the dipolar signal, which provides direct insights into optimizing acquisition schemes to mitigate orientational effects. Furthermore, we incorporate the excitation profile of realistic pulses to identify the excited spins. With this method, we show that rectangular pulses, despite their imperfect inversion capability, can sample similar spin orientations as other sophisticated pulses with the same bandwidth. Additionally, we reveal that the efficiency of exciting spin-pairs in DEER depends on the frequency offset of two pulses used in the experiment and the relative orientation of the two spins. Therefore, we systematically examine the frequency offset of the two pulses used in this double resonance experiment to determine the optimal frequency offset for optimal distance measurements. This procedure leads to a protocol where two measurements are sufficient to acquire orientational-independent DEER at Q-band. Notably, this procedure is feasible with any commercial pulsed-EPR spectrometer. Furthermore, we experimentally validate the computational results using DEER experiments on two different proteins. Finally, we show that increasing the amplitude of the rectangular pulse can increase the efficiency of DEER experiments by almost threefold. Overall, this work provides an attractive new approach for analyzing pulsed-EPR spectroscopy to obtain microscopic nuances that cannot be easily discerned from analytical or numerical calculations.

Reconciling membrane protein simulations with experimental DEER spectroscopy data

Spectroscopy experiments are crucial to study membrane proteins for which traditional structure determination methods still prove challenging. Double electron-electron resonance (DEER) spectroscopy experiments provide protein residue-pair distance distributions that are indicative of their conformational heterogeneity. Atomistic molecular dynamics (MD) simulations are another tool that have been proven to be vital to study the structural dynamics of membrane proteins such as to identify inward-open, occluded, and outward-open conformations of transporter membrane proteins, among other partially open or closed states of the protein. Yet, studies have reported that there is no direct consensus between the distributional data from DEER experiments and MD simulations, which has challenged validation of structures obtained from long-timescale simulations and using simulations to design experiments. Current coping strategies for comparisons rely on heuristics, such as mapping the nearest matching peaks between two ensembles or biased simulations. Here we examine the differences in residue-pair distance distributions arising due to the choice of membranes around the protein and covalent modification of a pair of residues to nitroxide spin labels in DEER experiments. Through comparing MD simulations of two proteins, PepT_So and LeuT-both of which have been characterized using DEER experiments previously-we show that the proteins' dynamics are similar despite the choice of the detergent micelle as a membrane mimetic in DEER experiments. On the other hand, covalently modified residues show slight local differences in their dynamics and a huge divergence when the oxygen atom pair distances between spin labeled residues are measured rather than protein backbone distances. Given the computational expense associated with pairwise MTSSL labeled MD simulations, we examine the use of biased simulations to explore the conformational dynamics of the spin labels only to reveal that such simulations alter the underlying protein dynamics. Our study identifies the main cause for the mismatch between DEER experiments and MD simulations and will accelerate the development of potential mitigation strategies to improve the match.

A new 13C trityl-based spin label enables the use of DEER for distance measurements

Triarylmethyl (TAM)-based labels, while still underutilized, are a powerful class of labels for pulsed-Electron Spin Resonance (ESR) distance measurements. They feature slow relaxation rates for long-lasting signals, high stability for cellular experiments, and narrow spectral features for efficient excitation of the spins. However, the typical narrow line shape limits the available distance measurements to only single-frequency experiments, such as Double Quantum Coherence (DQC) and Relaxation Induced Dipolar Modulation Enhancement (RIDME), which can be complicated to perform or hard to process. Therefore, widespread usage of TAM labels can be enhanced by the use of Double Electron-Electron Resonance (DEER) distance measurements. In this work, we developed a new spin label, 13C1-mOX063-d24, with a 13C isotope as the radical center. Due to the resolved hyperfine splitting, the spectrum is sufficiently broadened to permit DEER-based experiments at Q-band spectrometers. Additionally, this new label can be incorporated orthogonally with Cu(II)-based protein label. The orthogonal labeling scheme enables DEER distance measurement at X-band frequencies. Overall, the new trityl label allows for DEER-based distance measurements that complement existing TAM-label DQC and RIDME experiments.

The use of EPR spectroscopy to study transcription mechanisms

Electron paramagnetic resonance (EPR) spectroscopy has become a promising structural biology tool to resolve complex and dynamic biological mechanisms in-vitro and in-cell. Here, we focus on the advantages of continuous wave (CW) and pulsed EPR distance measurements to resolve transcription processes and protein-DNA interaction. The wide range of spin-labeling approaches that can be used to follow structural changes in both protein and DNA render EPR a powerful method to study protein-DNA interactions and structure-function relationships in other macromolecular complexes. EPR-derived data goes well beyond static structural information and thus serves as the method of choice if dynamic insight is needed. Herein, we describe the conceptual details of the theory and the methodology and illustrate the use of EPR to study the protein-DNA interaction of the copper-sensitive transcription factor, CueR.

Beyond structure: Deciphering site-specific dynamics in proteins from double histidine-based EPR measurements

Site-specific dynamics in proteins are at the heart of protein function. While electron paramagnetic resonance (EPR) has potential to measure dynamics in large protein complexes, the reliance on flexible nitroxide labels is limitating especially for the accurate measurement of site-specific β-sheet dynamics. Here, we employed EPR spectroscopy to measure site-specific dynamics across the surface of a protein, GB1. Through the use of the double Histidine (dHis) motif, which enables labeling with a Cu(II) - nitrilotriacetic acid (NTA) complex, dynamics information was obtained for both α-helical and β-sheet sites. Spectral simulations of the resulting CW-EPR report unique site-specific fluctuations across the surface of GB1. Additionally, we performed molecular dynamics (MD) simulations to complement the EPR data. The dynamics observed from MD agree with the EPR results. Furthermore, we observe small changes in gǁ values for different sites, which may be due to small differences in coordination geometry and/or local electrostatics of the site. Taken together, this work expands the utility of Cu(II)NTA-based EPR measurements to probe information beyond distance constraints.

An optimal acquisition scheme for Q-band EPR distance measurements using Cu2+-based protein labels

Recent advances in site-directed Cu²⁺ labeling of proteins and nucleic acids have added an attractive new methodology to measure the structure-function relationship in biomolecules. Despite the promise, accessing the higher sensitivity of Q-band Double Electron Electron Resonance (DEER) has been challenging for Cu²⁺ labels designed for proteins. Q-band DEER experiments on this label typically require many measurements at different magnetic fields, since the pulses can excite only a few orientations at a given magnetic field. Herein, we analyze such orientational effects through simulations and show that three DEER measurements, at strategically selected magnetic fields, are generally sufficient to acquire an orientational-averaged DEER time trace for this spin label at Q-band. The modeling results are experimentally verified on Cu²⁺ labeled human glutathione S-transferase (hGSTA1-1). The DEER distance distribution measured at the Q-band shows good agreement with the distance distribution sampled by molecular dynamics (MD) simulations and X-band experiments. The concordance of MD sampled distances and experimentally measured distances adds growing evidence that MD simulations can accurately predict distances for the Cu²⁺ labels, which remains a key bottleneck for the commonly used nitroxide label. In all, this minimal collection scheme reduces data collection time by as much as six-fold and is generally applicable to many octahedrally coordinated Cu²⁺ systems. Furthermore, the concepts presented here may be applied to other metals and pulsed EPR experiments.

dHis-troying Barriers: Deuteration Provides a Pathway to Increase Sensitivity and Accessible Distances for Cu2+ Labels

Recently, site-directed Cu²⁺ labeling has emerged as an incisive biophysical tool to directly report on distance constraints that pertain to the structure, conformational transitions, and dynamics of proteins and nucleic acids. However, short phase memory times inherent to the Cu²⁺ labels limit measurable distances to 4-5 nm. In this work we systematically examine different methods to dampen electron-nuclear and electron-electron coupled interactions to decrease rapid relaxation. We show that using Cu²⁺ spin concentrations up to ca. 800 μM has an invariant effect on relaxation and that increasing the cryoprotectant concentration reduces contributions of solvent protons to relaxation. On the other hand, the deuteration of protein and solvent dramatically increases the duration of the dipolar modulated signal by over 6-fold to 32 μs. Based on this increase in signal longevity, distances up to 9 nm and beyond can potentially be measured with Cu²⁺ labels.

A general model to optimise CuII labelling efficiency of double-histidine motifs for pulse dipolar EPR applications

Electron paramagnetic resonance (EPR) distance measurements are making increasingly important contributions to studies of biomolecules underpinning health and disease by providing highly accurate and precise geometric constraints. Combining double-histidine (dH) motifs with CuII spin labels shows promise for further increasing the precision of distance measurements, and for investigating subtle conformational changes. However, non-covalent coordination-based spin labelling is vulnerable to low binding affinity. Dissociation constants of dH motifs for CuII-nitrilotriacetic acid were previously investigated via relaxation induced dipolar modulation enhancement (RIDME), and demonstrated the feasibility of exploiting the dH motif for EPR applications at sub-μM protein concentrations. Herein, the feasibility of using modulation depth quantitation in CuII-CuII RIDME to simultaneously estimate a pair of non-identical independent KD values in such a tetra-histidine model protein is addressed. Furthermore, we develop a general speciation model to optimise CuII labelling efficiency, depending upon pairs of identical or disparate KD values and total CuII label concentration. We find the dissociation constant estimates are in excellent agreement with previously determined values, and empirical modulation depths support the proposed model.

Probing the Structure of Toxic Amyloid-β Oligomers with Electron Spin Resonance and Molecular Modeling

Structural models of the toxic species involved in the development of Alzheimer's disease are of utmost importance to understand the molecular mechanism and to describe early biomarkers of the disease. Among toxic species, soluble oligomers of amyloid-β (Aβ) peptides are particularly important, because they are responsible for spreading cell damages over brain regions, thus rapidly impairing brain functions. In this work we obtain structural information on a carefully prepared Aβ(1-42) sample, representing a toxic state for cell cultures, by combining electron spin resonance spectroscopy and computational models. We exploited the binding of Cu²⁺ to Aβ(1-42) and used copper as a probe for estimating Cu-Cu distances in the oligomers by applying double electron-electron resonance (DEER) pulse sequence. The DEER trace of this sample displays a unique feature that fits well with structural models of oligomers formed by Cu-cross-linked peptide dimers. Because Cu is bound to the Aβ(1-42) N-terminus, for the first time structural constraints that are missing in reported studies are provided at physiological conditions for the Aβ N-termini. These constraints suggest the Aβ(1-42) dimer as the building block of soluble oligomers, thus changing the scenario for any kinetic model of Aβ(1-42) aggregation.

Pulse Dipolar EPR Reveals Double-Histidine Motif CuII-NTA Spin-Labeling Robustness against Competitor Ions

Pulse-dipolar EPR is an appealing strategy for structural characterization of complex systems in solution that complements other biophysical techniques. Significantly, the emergence of genetically encoded self-assembling spin labels exploiting exogenously introduced double-histidine motifs in conjunction with Cu^II-chelates offers high precision distance determination in systems nonpermissive to thiol-directed spin labeling. However, the noncovalency of this interaction exposes potential vulnerabilities to competition from adventitious divalent metal ions, and pH sensitivity. Herein, a combination of room-temperature isothermal titration calorimetry (ITC) and cryogenic relaxation-induced dipolar modulation enhancement (RIDME) measurements are applied to the model protein Streptococcus sp. group G. protein G, B1 domain (GB1). Results demonstrate double-histidine motif spin labeling using Cu^II-nitrilotriacetic acid (Cu^II-NTA) is robust against the competitor ligand Zn^II-NTA at >1000-fold molar excess, and high nM binding affinity is surprisingly retained under acidic and basic conditions even though room temperature affinity shows a stronger pH dependence. This indicates the strategy is well-suited for diverse biological applications, with the requirement of other metal ion cofactors or slightly acidic pH not necessarily being prohibitive.

Pulse Dipolar EPR Reveals Double-Histidine Motif CuII-NTA Spin-Labeling Robustness against Competitor Ions

Pulse-dipolar EPR is an appealing strategy for structural characterization of complex systems in solution that complements other biophysical techniques. Significantly, the emergence of genetically encoded self-assembling spin labels exploiting exogenously introduced double-histidine motifs in conjunction with Cu^II-chelates offers high precision distance determination in systems nonpermissive to thiol-directed spin labeling. However, the noncovalency of this interaction exposes potential vulnerabilities to competition from adventitious divalent metal ions, and pH sensitivity. Herein, a combination of room-temperature isothermal titration calorimetry (ITC) and cryogenic relaxation-induced dipolar modulation enhancement (RIDME) measurements are applied to the model protein Streptococcus sp. group G. protein G, B1 domain (GB1). Results demonstrate double-histidine motif spin labeling using Cu^II-nitrilotriacetic acid (Cu^II-NTA) is robust against the competitor ligand Zn^II-NTA at >1000-fold molar excess, and high nM binding affinity is surprisingly retained under acidic and basic conditions even though room temperature affinity shows a stronger pH dependence. This indicates the strategy is well-suited for diverse biological applications, with the requirement of other metal ion cofactors or slightly acidic pH not necessarily being prohibitive.

Going the dHis-tance: Site-Directed Cu2+ Labeling of Proteins and Nucleic Acids

In this Account, we showcase site-directed Cu²⁺ labeling in proteins and DNA, which has opened new avenues for the measurement of the structure and dynamics of biomolecules using electron paramagnetic resonance (EPR) spectroscopy. In proteins, the spin label is assembled in situ from natural amino acid residues and a metal complex and requires no post-expression synthetic modification or purification procedures. The labeling scheme exploits a double histidine (dHis) motif, which utilizes endogenous or site-specifically mutated histidine residues to coordinate a Cu²⁺ complex. Pulsed EPR measurements on such Cu²⁺-labeled proteins potentially yield distance distributions that are up to 5 times narrower than the common protein spin label-the approach, thus, overcomes the inherent limitation of the current technology, which relies on a spin label with a highly flexible side chain. This labeling scheme provides a straightforward method that elucidates biophysical information that is costly, complicated, or simply inaccessible by traditional EPR labels. Examples include the direct measurement of protein backbone dynamics at β-sheet sites, which are largely inaccessible through traditional spin labels, and rigid Cu²⁺-Cu²⁺ distance measurements that enable higher precision in the analysis of protein conformations, conformational changes, interactions with other biomolecules, and the relative orientations of two labeled protein subunits. Likewise, a Cu²⁺ label has been developed for use in DNA, which is small, is nucleotide independent, and is positioned within the DNA helix. The placement of the Cu²⁺ label directly reports on the biologically relevant backbone distance. Additionally, for both of these labeling techniques, we have developed models for interpretation of the EPR distance information, primarily utilizing molecular dynamics (MD) simulations. Initial results using force fields developed for both protein and DNA labels have agreed with experimental results, which has been a major bottleneck for traditional spin labels. Looking ahead, we anticipate new combinations of MD and EPR to further our understanding of protein and DNA conformational changes, as well as working synergistically to investigate protein-DNA interactions.

Orientation and dynamics of Cu2+ based DNA labels from force field parameterized MD elucidates the relationship between EPR distance constraints and DNA backbone distances

Pulsed electron paramagnetic resonance (EPR) based distance measurements using the recently developed Cu2+-DPA label present a promising strategy for measuring DNA backbone distance constraints. Herein we develop force field parameters for Cu2+-DPA in order to understand the features of this label at an atomic level. We perform molecular dynamics (MD) simulations using the force field parameters of Cu2+-DPA on four different DNA duplexes. The distance between the Cu2+ centers, extracted from the 2 μs MD trajectories, agrees well with the experimental distance for all the duplexes. Further analyses of the trajectory provide insight into the orientation of the Cu2+-DPA inside the duplex that leads to such agreement with experiments. The MD results also illustrate the ability of the Cu2+-DPA to report on the DNA backbone distance constraints. Furthermore, measurement of fluctuations of individual residues showed that the flexibility of Cu2+-DPA in a DNA depends on the position of the label in the duplex, and a 2 μs MD simulation is not sufficient to fully capture the experimental distribution in some cases. Finally, the MD trajectories were utilized to understand the key aspects of the double electron electron resonance (DEER) results. The lack of orientational selectivity effects of the Cu2+-DPA at Q-band frequency is rationalized in terms of fluctuations in the Cu2+ coordination environment and rotameric fluctuations of the label linker. Overall, a combination of EPR and MD simulations based on the Cu2+-DPA labelling strategy can contribute towards understanding changes in DNA backbone conformations during protein-DNA interactions.

A Switch between Two Intrinsically Disordered Conformational Ensembles Modulates the Active Site of a Basic-Helix-Loop-Helix Transcription Factor

We report a conformational switch between two distinct intrinsically disordered subensembles within the active site of a transcription factor. This switch highlights an evolutionary benefit conferred by the high plasticity of intrinsically disordered domains, namely, their potential to dynamically sample a heterogeneous conformational space housing multiple states with tailored properties. We focus on proto-oncogenic basic-helix-loop-helix (bHLH)-type transcription factors, as these play key roles in cell regulation and function. Despite intense research efforts, the understanding of structure-function relations of these transcription factors remains incomplete as they feature intrinsically disordered DNA-interaction domains that are difficult to characterize, theoretically as well as experimentally. Here we characterize the structural dynamics of the intrinsically disordered region DNA-binding site of the vital MYC-associated transcription factor X (MAX). Integrating nuclear magnetic resonance (NMR) measurements, molecular dynamics (MD) simulations, and electron paramagnetic resonance (EPR) measurements, we show that, in the absence of DNA, the binding site of the free MAX₂ homodimer samples two intrinsically disordered conformational subensembles. These feature distinct structural properties: one subensemble consists of a set of highly flexible and spatially extended conformers, while the second features a set of "hinged" conformations. In this latter ensemble, the disordered N-terminal tails of MAX₂ fold back along the dimer, forming transient long-range contacts with the HLH-region and thereby exposing the DNA binding site to the solvent. The features of these divergent substates suggest two mechanisms by which protein conformational dynamics in MAX₂ might modulate DNA-complex formation: by enhanced initial recruitment of free DNA ligands, as a result of the wider conformational space sampled by the extended ensemble, and by direct exposure of the binding site and the corresponding strong electrostatic attractions presented while in the hinged conformations.

Development of Cu2+-Based Distance Methods and Force Field Parameters for the Determination of PNA Conformations and Dynamics by EPR and MD Simulations

Peptide nucleic acids (PNAs) are a promising group of synthetic analogues of DNA and RNA that offer several distinct advantages over the naturally occurring nucleic acids for applications in biosensing, drug delivery, and nanoelectronics. Because of its structural differences from DNA/RNA, methods to analyze and assess the structure, conformations, and dynamics are needed. In this work, we develop synergistic techniques for the study of the PNA conformation. We use CuQ₂, a Cu²⁺ complex with 8-hydroxyquinoline (HQ), as an alternative base pair and as a spin label in electron paramagnetic resonance (EPR) distance methods. We use molecular dynamics (MD) simulations with newly developed force field parameters for the spin labels to interpret the distance constraints determined by EPR. We complement these methods by UV-vis and circular dichroism measurements and assess the efficacy of the Cu²⁺ label on a PNA duplex whose backbone is based on aminoethylglycine and a duplex with a hydroxymethyl backbone modification. We show that the Cu²⁺ label functions efficiently within the standard PNA and the hydroxymethyl-modified PNA and that the MD parameters may be used to accurately reproduce our EPR findings. Through the combination of EPR and MD, we gain new insights into the PNA structure and conformations as well as into the mechanism of orientational selectivity in Cu²⁺ EPR at X-band. These results present for the first time a rigid Cu²⁺ spin label used for EPR distance measurements in PNA and the accompanying MD force fields for the spin label. Our studies also reveal that the spin labels have a low impact on the structure of the PNA duplexes. The combined MD and EPR approach represents an important new tool for the characterization of the PNA duplex structure and provides valuable information to aid in the rational application of PNA at large.

Robust Atomistic Modeling of Materials, Organometallic, and Biochemical Systems

Modern chemistry seems to be unlimited in molecular size and elemental composition. Metal‐organic frameworks or biological macromolecules involve complex architectures and a large variety of elements. Yet, a general and broadly applicable theoretical method to describe the structures and interactions of molecules beyond the 1000‐atom size regime semi‐quantitatively is not self‐evident. For this purpose, a generic force field named GFN‐FF is presented, which is completely newly developed to enable fast structure optimizations and molecular‐dynamics simulations for basically any chemical structure consisting of elements up to radon. The freely available computer program requires only starting coordinates and elemental composition as input from which, fully automatically, all potential‐energy terms are constructed. GFN‐FF outperforms other force fields in terms of generality and accuracy, approaching the performance of much more elaborate quantum‐mechanical methods in many cases.

Modeling of spin–spin distance distributions for nitroxide labeled biomacromolecules

Combining CREST and MD simulations based on GFN-FF for the automated computation of distance distributions for nitroxide labeled (metallo-) proteins.

Accessing distributions of exchange and dipolar couplings in stiff molecular rulers with Cu(ii) centres

Novel approaches to quantitatively analyse distributed exchange couplings are described and tested on experimental data sets for stiff synthetic molecules.