The rise of the machines in chemistry.

The use of artificial intelligence and, more specifically, deep learning methods in chemistry is becoming increasingly common. Applications in informatics fields, such as cheminformatics and proteomics, structural biology, and spectroscopy, including NMR, are on the rise. Recent developments in model architectures, such as graph convolutional neural networks and transformers, have been enabled by advancements in computational hardware and software. However, model architectures with more predictive power often require larger amounts of training data, which can be challenging to acquire, but this requirement can be mitigated through techniques like pretraining and fine-tuning. In spite of these successes, challenges remain, such as normalization and scaling of data, availability of experimentally acquired data, and model explainability.

Detection of chemical exchange in methyl groups of macromolecules.

The zero- and double-quantum methyl TROSY Hahn-echo and the methyl 1H-1H dipole-dipole cross-correlation nuclear magnetic resonance experiments enable estimation of multiple quantum chemical exchange broadening in methyl groups in proteins. The two relaxation rate constants are established to be linearly dependent using molecular dynamics simulations and empirical analysis of experimental data. This relationship allows chemical exchange broadening to be recognized as an increase in the Hahn-echo relaxation rate constant. The approach is illustrated by analyzing relaxation data collected at three temperatures for E. coli ribonuclease HI and by analyzing relaxation data collected for different cofactor and substrate complexes of E. coli AlkB.

Dynamics of GCN4 facilitate DNA interaction: a model-free analysis of an intrinsically disordered region.

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) are known to play important roles in regulatory and signaling pathways. A critical aspect of these functions is the ability of IDP/IDRs to form highly specific complexes with target molecules. However, elucidation of the contributions of conformational dynamics to function has been limited by challenges associated with structural heterogeneity of IDP/IDRs. Using NMR spin relaxation parameters ((15)N R1, (15)N R2, and {(1)H}-(15)N heteronuclear NOE) collected at four static magnetic fields ranging from 14.1 to 21.1 T, we have analyzed the backbone dynamics of the basic leucine-zipper (bZip) domain of the Saccharomyces cerevisiae transcription factor GCN4, whose DNA binding domain is intrinsically disordered in the absence of DNA substrate. We demonstrate that the extended model-free analysis can be applied to proteins with IDRs such as apo GCN4 and that these results significantly extend previous NMR studies of GCN4 dynamics performed using a single static magnetic field of 11.74 T [Bracken, et al., J. Mol. Biol., 1999, 285, 2133-2146] and correlate well with molecular dynamics simulations [Robustelli, et al., J. Chem. Theory Comput., 2013, 9, 5190-5200]. In contrast to the earlier work, data at multiple static fields allows the time scales of internal dynamics of GCN4 to be reliably quantified. Large amplitude dynamic fluctuations in the DNA-binding region have correlation times (τs ≈ 1.4-2.5 ns) consistent with a two-step mechanism in which partially ordered bZip conformations of GCN4 form initial encounter complexes with DNA and then rapidly rearrange to the high affinity state with fully formed basic region recognition helices.

Protein dynamics control the progression and efficiency of the catalytic reaction cycle of the Escherichia coli DNA-repair enzyme AlkB.

A central goal of enzymology is to understand the physicochemical mechanisms that enable proteins to catalyze complex chemical reactions with high efficiency. Recent methodological advances enable the contribution of protein dynamics to enzyme efficiency to be explored more deeply. Here, we utilize enzymological and biophysical studies, including NMR measurements of conformational dynamics, to develop a quantitative mechanistic scheme for the DNA repair enzyme AlkB. Like other iron/2-oxoglutarate-dependent dioxygenases, AlkB employs a two-step mechanism in which oxidation of 2-oxoglutarate generates a highly reactive enzyme-bound oxyferryl intermediate that, in the case of AlkB, slowly hydroxylates an alkylated nucleobase. Our results demonstrate that a microsecond-to-millisecond time scale conformational transition facilitates the proper sequential order of substrate binding to AlkB. Mutations altering the dynamics of this transition allow generation of the oxyferryl intermediate but promote its premature quenching by solvent, which uncouples 2-oxoglutarate turnover from nucleobase oxidation. Therefore, efficient catalysis by AlkB depends upon the dynamics of a specific conformational transition, establishing another paradigm for the control of enzyme function by protein dynamics.

Multiplet-filtered and gradient-selected zero-quantum TROSY experiments for 13C1H3 methyl groups in proteins.

Multiplet-filtered and gradient-selected heteronuclear zero-quantum coherence (gsHZQC) TROSY experiments are described for measuring (1)H-(13)C correlations for (13)CH(3) methyl groups in proteins. These experiments provide improved suppression of undesirable, broad outer components of the heteronuclear zero-quantum multiplet in medium-sized proteins, or in flexible sites of larger proteins, compared to previously described HZQC sequences (Tugarinov et al. in J Am Chem Soc 126:4921-4925, 2004; Ollerenshaw et al. in J Biomol NMR 33:25-41, 2005). Hahn-echo versions of the gsHZQC experiment also are described for measuring zero- and double-quantum transverse relaxation rate constants for identification of chemical exchange broadening. Application of the proposed pulse sequences to Escherichia coli ribonuclease HI, with a molecular mass of 18 kD, indicates that improved multiplet suppression is obtained without substantial loss of sensitivity.

Crystallization and characterization of the thallium form of the Oxytricha nova G-quadruplex.

The crystal structure of the Tl+ form of the G-quadruplex formed from the Oxytricha nova telomere sequence, d(G4T4G4), has been solved to 1.55 A. This G-quadruplex contains five Tl+ ions, three of which are interspersed between adjacent G-quartet planes and one in each of the two thymine loops. The structure displays a high degree of similarity to the K+ crystal structure [Haider et al. (2002), J. Mol. Biol., 320, 189-200], including the number and location of the monovalent cation binding sites. The highly isomorphic nature of the two structures, which contain such a large number of monovalent binding sites (relative to nucleic acid content), verifies the ability of Tl+ to mimic K+ in nucleic acids. Information from this report confirms and extends the assignment of 205Tl resonances from a previous report [Gill et al. (2005), J. Am. Chem. Soc., 127, 16 723-16 732] where 205Tl NMR was used to study monovalent cation binding to this G-quadruplex. The assignment of these resonances provides evidence for the occurrence of conformational dynamics in the thymine loop region that is in slow exchange on the 205Tl timescale.

Survey of solution dynamics in Src kinase reveals allosteric cross talk between the ligand binding and regulatory sites.

The catalytic domain of protein tyrosine kinases can interconvert between active and inactive conformations in response to regulatory inputs. We recently demonstrated that Src kinase features an allosteric network that couples substrate-binding sites. However, the extent of conformational and dynamic changes that are propagated throughout the kinase domain remains poorly understood. Here, we monitor by NMR the effect of conformationally selective inhibitors on kinase backbone dynamics. We find that inhibitor binding and activation loop autophosphorylation induces dynamic changes across the entire kinase. We identify a highly conserved amino acid, Gly449, that is necessary for Src activation. Finally, we show for the first time how the SH3-SH2 domains perturb the dynamics of the kinase domain in the context of the full length protein. We provide experimental support for long-range communication in Src kinase that leads to the relative stabilization of active or inactive conformations and modulation of substrate affinity.

A fluorescence study of the structure and accessibility of plasmid DNA condensed with cationic gene delivery vehicles.

The cationic lipids 1,2-dioleoyl-3-trimethylammonium-propane and dimethyldioctadecylammonium bromide, with or without the helper lipids 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine or cholesterol, and the cationic polymer polyethyleneimine, were compared for their ability to displace fluorescent dyes from DNA. Differences in displacement of the intercalating dyes ethidium bromide and ethidium homodimer correlate with their relative affinities with DNA, with the extent of ethidium homodimer displacement significantly less. Differences in ethidium homodimer and ethidium bromide displacement as a function of the ratio of polycation to DNA and the charge density of the polycation suggest a greater sensitivity of the former to topological changes in condensed DNA. Marked differences in the ability of these cationic delivery systems to displace the minor groove binding dyes 4',6-diamidino-2-phenylindole and Hoechst 33258 upon interaction with DNA are also apparent, with the majority of Hoechst 33258 remaining bound to DNA. Changes in the spectral properties of Hoechst 33258 were further used to characterize polycation-induced changes in solvent accessibility of the DNA minor groove. Taken together, these studies demonstrate differences in the interaction of various cationic lipids and polyethyleneimine in terms of regional displacement of dyes, polycation-induced structural changes in DNA, as well as polycation-mediated changes in solvent accessibility of the minor groove. The relevance of these studies to current models of the structure and assembly of polycation/DNA complexes are discussed.

Local isotropic diffusion approximation for coupled internal and overall molecular motions in NMR spin relaxation.

The present work demonstrates that NMR spin relaxation rate constants for molecules interconverting between states with different diffusion tensors can be modeled theoretically by combining orientational correlation functions for exchanging spherical molecules with locally isotropic approximations for the diffusion anisotropic tensors. The resulting expressions are validated by comparison with correlation functions obtained by Monte Carlo simulations and are accurate for moderate degrees of diffusion anisotropy typically encountered in investigations of globular proteins. The results are complementary to an elegant, but more complex, formalism that is accurate for all degrees of diffusion anisotropy [Ryabov, Y.; Clore, G. M.; Schwieters, C. D. J. Chem. Phys. 2012, 136, 034108].

Using empirical phase diagrams to understand the role of intramolecular dynamics in immunoglobulin G stability.

Understanding the relationship between protein dynamics and stability is of paramount importance to the fields of biology and pharmaceutics. Clarifying this relationship is complicated by the large amount of experimental data that must be generated and analyzed if motions that exist over the wide range of timescales are to be included. To address this issue, we propose an approach that utilizes a multidimensional vector-based empirical phase diagram (EPD) to analyze a set of dynamic results acquired across a temperature-pH perturbation plane. This approach is applied to a humanized immunoglobulin G1 (IgG1), a protein of major biological and pharmaceutical importance whose dynamic nature is linked to its multiple biological roles. Static and dynamic measurements are used to characterize the IgG and to construct both static and dynamic EPDs. Between pH 5 and 8, a single, pH-dependent transition is observed that corresponds to thermal unfolding of the IgG. Under more acidic conditions, evidence exists for the formation of a more compact, aggregation resistant state of the immunoglobulin, known as A-form. The dynamics-based EPD presents a considerably more detailed pattern of apparent phase transitions over the temperature-pH plane. The utility and potential applications of this approach are discussed.

205Tl NMR methods for the characterization of monovalent cation binding to nucleic acids.

Monovalent cations play an important role in many biological functions. The guanine rich sequence, d(G4T4G4), requires monovalent cations for formation of the G-quadruplex, d(G4T4G4)2. This requirement can be satisfied by thallium (Tl+), a potassium (K+) surrogate. To verify that the structure of d(G4T4G4)2 in the presence of Tl+ is similar to the K+-form of the G-quadruplex, the solution structure of the Tl+-form of d(G4T4G4)2 was determined. The 10 lowest energy structures have an all atom RMSD of 0.76 +/- 0.16 A. Comparison of this structure to the identical G-quadruplex formed in the presence of K+ validates the isomorphous nature of Tl+ and K+. Using a 1H-205Tl spin-echo difference experiment we show that, in the Tl+-form of d(G4T4G4)2, small scalar couplings (<1 Hz) exist between 205Tl and protons in the G-quadruplex. These data comprise the first 1H-205Tl scalar couplings observed in a biological system and have the potential to provide important constraints for structure determination. These experiments can be applied to any system in which the substituted Tl+ cations are in slow exchange with the bulk ions in solution.

Conservation of mus-ms enzyme motions in the apo- and substrate-mimicked state.

Solution NMR spin-relaxation experiments were used to compare mus-ms dynamics in RNase A in the apo form and as complexed to the substrate-mimic, pTppAp. The crystal structure of the RNase A/pTppAp complex was determined and demonstrates that this ligand binds at the active site and utilizes established substrate binding sites in its interaction with RNase A. Relaxation-compensated CPMG experiments identify flexible residues in and around the active site in both the apo and pTppAp-bound enzyme. Quantitative analysis of the NMR spin-relaxation dispersion curves show that the time scale of motion in RNase A is unchanged when pTppAp binds and is similar to the time scale for the rate-determining step of the catalytic reaction. Temperature-dependent measurements provide an activation barrier for motion of 5.2 +/- 1.0 kcal/mol and 4.5 +/- 1.2 kcal/mol for the apo and pTppAp forms of RNase A, respectively. These data indicate very similar motion exists in the free and bound enzyme. Additionally, chemical shift data suggests that the magnitude of motion is also similar for these two forms and that it is likely that apo enzyme interconverts to a structure that resembles a ligand-bound form. Likewise, it appears that the bound conformation samples the apo enzyme form even when ligand is present. Taken together the data imply that RNase A is in a preexisting dynamic equilibrium between two conformations that represent the open and closed enzyme forms. These data suggest that ligand binding stabilizes the bound conformer but does not induce it.

Crystal structure of a group I intron splicing intermediate.

A recently reported crystal structure of an intact bacterial group I self-splicing intron in complex with both its exons provided the first molecular view into the mechanism of RNA splicing. This intron structure, which was trapped in the state prior to the exon ligation reaction, also reveals the architecture of a complex RNA fold. The majority of the intron is contained within three internally stacked, but sequence discontinuous, helical domains. Here the tertiary hydrogen bonding and stacking interactions between the domains, and the single-stranded joiner segments that bridge between them, are fully described. Features of the structure include: (1) A pseudoknot belt that circumscribes the molecule at its longitudinal midpoint; (2) two tetraloop-tetraloop receptor motifs at the peripheral edges of the structure; (3) an extensive minor groove triplex between the paired and joiner segments, P6-J6/6a and P3-J3/4, which provides the major interaction interface between the intron's two primary domains (P4-P6 and P3-P9.0); (4) a six-nucleotide J8/7 single stranded element that adopts a mu-shaped structure and twists through the active site, making critical contacts to all three helical domains; and (5) an extensive base stacking architecture that realizes 90% of all possible stacking interactions. The intron structure was validated by hydroxyl radical footprinting, where strong correlation was observed between experimental and predicted solvent accessibility. Models of the pre-first and pre-second steps of intron splicing are proposed with full-sized tRNA exons. They suggest that the tRNA undergoes substantial angular motion relative to the intron between the two steps of splicing.

The structural organization of cationic lipid-DNA complexes.

The interaction of cationic liposomes with supercoiled plasmid DNA results in a major rearrangement of each component to form compact multilamellar structures comprised of alternating layers of two-dimensional arrays of DNA sandwiched between lipid bilayers. Fluorescence resonance energy transfer was used to estimate the distance of closest approach of DNA to the lipid bilayers in these complexes. The effect of several compositional variables on this distance, including the ratio of cationic lipid to DNA, and the charge density, intrinsic curvature, and fluidity of the lipid bilayer were examined. Additionally, the effect of ionic strength was studied. For complexes prepared at or above a 3:1 charge ratio (+/-), the observed distance of closest approach was found to be in agreement with the intercalation of DNA between lipid bilayers. As the charge ratio was decreased, a monotonic increase in the distance was observed with a maximum observed at 0.5:1. Correlations between differences in the proximity of DNA to the lipid bilayer and the hydrodynamic size of the complexes were also found. A model based on these observations and previous reports suggests the formation of discrete populations of complexes below a charge ratio of 0.5:1 and above 3:1. The structure of the negatively charged complexes is consistent with DNA extending from the surface of the particles, whereas those possessing excess positive charge were multilamellar aggregates with the DNA effectively condensed between lipid bilayers. Complexes between these two states consist of weighted fractions of these two species.