Theoretical and computational advances in biomolecular NMR spectroscopy

Solvent Dependent Nuclear Magnetic Resonance Molecular Parameters Based on a Polarization Consistent Screened Range Separated Hybrid Density Functional Theory Framework

Nuclear magnetic resonance (NMR) properties of solvated molecules are significantly affected by the solvent. We, therefore, employ a polarization consistent framework that efficiently addresses the solvent polarizing environment effects. Toward this goal a dielectric screened range separated hybrid (SRSH) functional is invoked with a polarizable continuum model (PCM) to properly represent the orbital gap in the condensed phase. We build on the success of range separated hybrid (RSH) functionals to address the erroneous tendency of traditional density functional theory (DFT) to collapse the orbital gap. Recently, the impact of RSH that properly opens up the orbital gap in gas-phase calculations on NMR properties has been assessed. Here, we report the use of SRSH-PCM that produces properly solute orbital gaps in calculating isotropic nuclear magnetic shielding and chemical shift parameters of molecular systems in the condensed phase. We show that in contrast to simpler DFT-PCM approaches, SRSH-PCM successfully follows expected dielectric constant trends.

Combination of pure shift NMR and chemical shift selective filters for analysis of Fischer-Tropsch waste-water

Fischer-Tropsch (F-T) process is an important synthesis route to acquire clean liquid fuels through modern coal chemical industry, which converts syngas (CO and H₂) into hydrocarbon, and also generates oxygenates discharged as the F-T waste-water. These oxygen-containing compounds in F-T waste-water have the similar molecular weight and some are even isomers of each other. Hence, it is necessary to develop rapid and efficient analysis tools to obtain identification and quantitative information of the F-T waste-water. The pure shift NMR techniques provided only chemical shift information in one-dimension ¹H NMR spectra, without homonuclear J_H-H coupling. In this work, we tested and compared three pure shift NMR techniques (including Zangger-Sterk, PSYCHE and TSE-PSYCHE methods) in the analysis of two F-T waste-water model mixtures, genuine waste-water and two alcohol isomer mixtures. The results show that J_H-H coupling multiplicities are collapsed into singlets corresponding to individual chemically distinct protons of the compound. For some severely overlapped signals in the pure shift NMR spectra, the chemical shift selective filters with TOCSY (CSSF-TOCSY) experiments were conducted to assist the signal assignment. Thus, pure shift NMR approaches can identify most signals of components, and CSSF-TOCSY can extract the signal of a specific compound. The combination of these two NMR techniques offers a powerful tool to analyze the F-T waste-water or other complex mixtures including isomer mixtures.

PyShifts: A PyMOL Plugin for Chemical Shift-Based Analysis of Biomolecular Ensembles

Here, we present PyShifts-a PyMOL plugin for chemical shift-based analysis of biomolecular ensembles. With PyShifts, users can compare and visualize differences between experimentally measured and computationally predicted chemical shifts. When analyzing multiple conformations of a biomolecule with PyShifts, users can also sort a set of conformations based on chemical shift differences and identify the conformers that exhibit the best agreement between measured and predicted chemical shifts. Although we have integrated PyShifts with the chemical shift predictors LARMOR^D and LARMOR^Cα, PyShifts can read in chemical shifts from any source, and so, users can employ PyShifts to analyze biomolecular structures using chemical shifts computed by any chemical shift predictor. We envision, therefore, that PyShifts (https://github.com/atfrank/PyShifts) will find utility as a general-purpose tool for exploring chemical shift-structure relationships in biomolecular ensembles.

Dynamic Descriptions of Highly Flexible Molecules from NMR Dipolar Couplings: Physical Basis and Limitations

Biomolecules that control physiological function by changing their conformation play key roles in biology and remain poorly characterized. NMR dipolar couplings (DCs) depend intrinsically on both molecular shape and structural fluctuations, thereby providing the enticing prospect of tracking these conformational changes at atomic detail. Although this dual dependence has until now severely complicated analysis of DCs from highly dynamic systems, general approaches have recently been proposed that simplify interpretation of experimental DCs, by entirely eliminating molecular alignment from the analysis. Using simple and intuitive simulation of target ensembles, we investigate the impact of such approaches on the resulting descriptions of the conformational energy landscape. We find that ensemble descriptions of highly flexible systems derived from DCs without explicit consideration of the alignment properties of the constituent conformations can be compromised and inaccurate, despite exhibiting high correlation with experimental measurement.

The PRE-Derived NMR Model of the 38.8-kDa Tri-Domain IsdH Protein from Staphylococcus aureus Suggests That It Adaptively Recognizes Human Hemoglobin

Staphylococcus aureus is a medically important bacterial pathogen that, during infections, acquires iron from human hemoglobin (Hb). It uses two closely related iron-regulated surface determinant (Isd) proteins to capture and extract the oxidized form of heme (hemin) from Hb, IsdH and IsdB. Both receptors rapidly extract hemin using a conserved tri-domain unit consisting of two NEAT (near iron transporter) domains connected by a helical linker domain. To gain insight into the mechanism of extraction, we used NMR to investigate the structure and dynamics of the 38.8-kDa tri-domain IsdH protein (IsdH(N2N3), A326-D660 with a Y642A mutation that prevents hemin binding). The structure was modeled using long-range paramagnetic relaxation enhancement (PRE) distance restraints, dihedral angle, small-angle X-ray scattering, residual dipolar coupling and inter-domain NOE nuclear Overhauser effect data. The receptor adopts an extended conformation wherein the linker and N3 domains pack against each other via a hydrophobic interface. In contrast, the N2 domain contacts the linker domain via a hydrophilic interface and, based on NMR relaxation data, undergoes inter-domain motions enabling it to reorient with respect to the body of the protein. Ensemble calculations were used to estimate the range of N2 domain positions compatible with the PRE data. A comparison of the Hb-free and Hb-bound forms reveals that Hb binding alters the positioning of the N2 domain. We propose that binding occurs through a combination of conformational selection and induced-fit mechanisms that may promote hemin release from Hb by altering the position of its F helix.

Modeling helical proteins using residual dipolar couplings, sparse long-range distance constraints and a simple residue-based force field

We present a fast and simple protocol to obtain moderate-resolution backbone structures of helical proteins. This approach utilizes a combination of sparse backbone NMR data (residual dipolar couplings and paramagnetic relaxation enhancements) or EPR data with a residue-based force field and Monte Carlo/simulated annealing protocol to explore the folding energy landscape of helical proteins. By using only backbone NMR data, which are relatively easy to collect and analyze, and strategically placed spin relaxation probes, we show that it is possible to obtain protein structures with correct helical topology and backbone RMS deviations well below 4 Å. This approach offers promising alternatives for the structural determination of proteins in which nuclear Overha-user effect data are difficult or impossible to assign and produces initial models that will speed up the high-resolution structure determination by NMR spectroscopy.

Interplay between conformational selection and induced fit in multidomain protein-ligand binding probed by paramagnetic relaxation enhancement

The binding of ligands and substrates to proteins has been extensively studied for many years and can be described, in its simplest form, by two limiting mechanisms: conformational selection and induced fit. Conformational selection involves the binding of ligand to a pre-existing sparsely-populated conformation of the free protein that is the same as that in the final protein-ligand complex. In the case of induced fit, the ligand binds to the major conformation of the free protein and only subsequent to binding undergoes a conformational change to the final protein-ligand complex. While these two mechanisms can be dissected and distinguished by transient kinetic measurements, direct direction, characterization and visualization of transient, sparsely-populated states of proteins are experimentally challenging. Unless trapped, sparsely-populated states are generally invisible to conventional structural and biophysical techniques, including crystallography and most NMR measurements. In this review we summarize some recent developments in the use of paramagnetic relaxation enhancement to directly study sparsely-populated states of proteins and illustrate the application of this approach to two proteins, maltose binding protein and calmodulin, both of which undergo large rigid body conformational rearrangements upon ligand binding from an open apo state to a closed ligand-bound holo state. We show that the apo state ensemble comprises a small population of partially-closed configurations that are similar but not identical to that of the holo state. These results highlight the complementarity and interplay of induced fit and conformational selection and suggest that the existence of partially-closed states in the absence of ligand facilitates the transition to the closed ligand-bound state.

Paramagnetic nuclear magnetic resonance relaxation and molecular mechanics studies of the chloroperoxidase-indole complex: insights into the mechanism of chloroperoxidase-catalyzed regioselective oxidation of indole

To unravel the mechanism of chloroperoxidase (CPO)-catalyzed regioselective oxidation of indole, we studied the structure of the CPO-indole complex using nuclear magnetic resonance (NMR) relaxation measurements and computational techniques. The dissociation constant (KD) of the CPO-indole complex was calculated to be approximately 21 mM. The distances (r) between protons of indole and the heme iron calculated via NMR relaxation measurements and molecular docking revealed that the pyrrole ring of indole is oriented toward the heme with its 2-H pointing directly at the heme iron. Both KD and r values are independent of pH in the range of 3.0-6.5. The stability and structure of the CPO-indole complex are also independent of the concentration of chloride or iodide ion. Molecular docking suggests the formation of a hydrogen bond between the NH group of indole and the carboxyl O of Glu 183 in the binding of indole to CPO. Simulated annealing of the CPO-indole complex using r values from NMR experiments as distance restraints reveals that the van der Waals interactions were much stronger than the Coulomb interactions in the binding of indole to CPO, indicating that the association of indole with CPO is primarily governed by hydrophobic rather than electrostatic interactions. This work provides the first experimental and theoretical evidence of the long-sought mechanism that leads to the "unexpected" regioselectivity of the CPO-catalyzed oxidation of indole. The structure of the CPO-indole complex will serve as a lighthouse in guiding the design of CPO mutants with tailor-made activities for biotechnological applications.

The structural basis for homotropic and heterotropic cooperativity of midazolam metabolism by human cytochrome P450 3A4

Human cytochrome P450 3A4 (CYP3A4) metabolizes a significant portion of clinically relevant drugs and often exhibits complex steady-state kinetics that can involve homotropic and heterotropic cooperativity between bound ligands. In previous studies, the hydroxylation of the sedative midazolam (MDZ) exhibited homotropic cooperativity via a decrease in the ratio of 1'-OH-MDZ to 4-OH-MDZ at higher drug concentrations. In this study, MDZ exhibited heterotropic cooperativity with the antiepileptic drug carbamazepine (CBZ) with characteristic decreases in the 1'-OH-MDZ to 4-OH-MDZ ratios. To unravel the structural basis of MDZ cooperativity, we probed MDZ and CBZ bound to CYP3A4 using longitudinal T(1) nuclear magnetic resonance (NMR) relaxation and molecular docking with AutoDock 4.2. The distances calculated from longitudinal T(1) NMR relaxation were used during simulated annealing to constrain the molecules to the substrate-free X-ray crystal structure of CYP3A4. These simulations revealed that either two MDZ molecules or an MDZ molecule and a CBZ molecule assume a stacked configuration within the CYP3A4 active site. In either case, the proton at position 4 of the MDZ molecule was closer to the heme than the protons of the 1'-CH(3) group. In contrast, molecular docking of a single molecule of MDZ revealed that the molecule was preferentially oriented with the 1'-CH(3) position closer to the heme than position 4. This study provides the first detailed molecular analysis of heterotropic and homotropic cooperativity of a human cytochrome P450 from an NMR-based model. Cooperativity of ligand binding through direct interaction between stacked molecules may represent a common motif for homotropic and heterotropic cooperativity.

Mammalian testis-determining factor SRY and the enigma of inherited human sex reversal: frustrated induced fit in a bent protein-DNA complex

Mammalian testis-determining factor SRY contains a high mobility group box, a conserved eukaryotic motif of DNA bending. Mutations in SRY cause XY gonadal dysgenesis and somatic sex reversal. Although such mutations usually arise de novo in spermatogenesis, some are inherited and so specify male development in one genetic background (the father) but not another (the daughter). Here, we describe the biophysical properties of a representative inherited mutation, V60L, within the minor wing of the L-shaped domain (box position 5). Although the stability and DNA binding properties of the mutant domain are similar to those of wild type, studies of SRY-induced DNA bending by subnanosecond time-resolved fluorescence resonance energy transfer (FRET) revealed enhanced conformational fluctuations leading to long range variation in bend angle. (1)H NMR studies of the variant protein-DNA complex demonstrated only local perturbations near the mutation site. Because the minor wing of SRY folds on DNA binding, the inherited mutation presumably hinders induced fit. Stopped-flow FRET studies indicated that such frustrated packing leads to accelerated dissociation of the bent complex. Studies of SRY-directed transcriptional regulation in an embryonic gonadal cell line demonstrated partial activation of downstream target Sox9. Our results have demonstrated a nonlocal coupling between DNA-directed protein folding and protein-directed DNA bending. Perturbation of this coupling is associated with a genetic switch poised at the threshold of activity.

NMR-derived models of amidopyrine and its metabolites in complexes with rabbit cytochrome P450 2B4 reveal a structural mechanism of sequential N-dealkylation

To understand the molecular basis of sequential N-dealkylation by cytochrome P450 2B enzymes, we studied the binding of amidopyrine (AP) as well as the metabolites of this reaction, desmethylamidopyrine (DMAP) and aminoantipyrine (AAP), using the X-ray crystal structure of rabbit P450 2B4 and two nuclear magnetic resonance (NMR) techniques: saturation transfer difference (STD) spectroscopy and longitudinal (T(1)) relaxation NMR. Results of STD NMR of AP and its metabolites bound to P450 2B4 were similar, suggesting that they occupy similar niches within the enzyme's active site. The model-dependent relaxation rates (R(M)) determined from T(1) relaxation NMR of AP and DMAP suggest that the N-linked methyl is closest to the heme. To determine the orientation(s) of AP and its metabolites within the P450 2B4 active site, we used distances calculated from the relaxation rates to constrain the metabolites to the X-ray crystal structure of P450 2B4. Simulated annealing of the complex revealed that the metabolites do indeed occupy similar hydrophobic pockets within the active site, while the N-linked methyls are free to rotate between two binding modes. From these bound structures, a model of N-demethylation in which the N-linked methyl functional groups rotate between catalytic and noncatalytic positions was developed. This study is the first to provide a structural model of a drug and its metabolites complexed to a cytochrome P450 based on NMR and to provide a structural mechanism for how a drug can undergo sequential oxidations without unbinding. The rotation of the amide functional group might represent a common structural mechanism for N-dealkylation reactions for other drugs such as the local anesthetic lidocaine.

Biomolecular Structure and Modeling: Historical Perspective

physics, chemistry, and biology have been connected by a web of causal explanation organized by induction-based theories that telescope into one another. … Thus, quantum theory underlies atomic physics, which is the foundation of reagent chemistry and its specialized offshoot biochemistry, which interlock with molecular biology — essentially, the chemistry of organic macromolecules — and hence, through successively higher levels of organization, cellular, organismic, and evolutionary biology. … Such is the unifying and highly productive understanding of the world that has evolved in the natural sciences.

Improving NMR protein structure quality by Rosetta refinement: a molecular replacement study

The structure of human protein HSPC034 has been determined by both solution nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography. Refinement of the NMR structure ensemble, using a Rosetta protocol in the absence of NMR restraints, resulted in significant improvements not only in structure quality, but also in molecular replacement (MR) performance with the raw X-ray diffraction data using MOLREP and Phaser. This method has recently been shown to be generally applicable with improved MR performance demonstrated for eight NMR structures refined using Rosetta (Qian et al., Nature 2007;450:259-264). Additionally, NMR structures of HSPC034 calculated by standard methods that include NMR restraints have improvements in the RMSD to the crystal structure and MR performance in the order DYANA, CYANA, XPLOR-NIH, and CNS with explicit water refinement (CNSw). Further Rosetta refinement of the CNSw structures, perhaps due to more thorough conformational sampling and/or a superior force field, was capable of finding alternative low energy protein conformations that were equally consistent with the NMR data according to the Recall, Precision, and F-measure (RPF) scores. On further examination, the additional MR-performance shortfall for NMR refined structures as compared with the X-ray structure were attributed, in part, to crystal-packing effects, real structural differences, and inferior hydrogen bonding in the NMR structures. A good correlation between a decrease in the number of buried unsatisfied hydrogen-bond donors and improved MR performance demonstrates the importance of hydrogen-bond terms in the force field for improving NMR structures. The superior hydrogen-bond network in Rosetta-refined structures demonstrates that correct identification of hydrogen bonds should be a critical goal of NMR structure refinement. Inclusion of nonbivalent hydrogen bonds identified from Rosetta structures as additional restraints in the structure calculation results in NMR structures with improved MR performance.

NMR structure and dynamics of human ephrin-B2 ectodomain: the functionally critical C-D and G-H loops are highly dynamic in solution

Eph receptors and ephrins constitute the largest family of receptor tyrosine kinases with 15 individual receptors and nine ligands. Its ectodomains represent attractive targets not only for understanding fundamental mechanisms underlying axon guidance, cell migration, segmentation, tumorigenesis, and bone remodeling, but also for drug screening/design to treat cancers, bone diseases and viral infection. So far no NMR study on the ephrin ectodomains is available and as such their properties in solution still remain unknown. In this study, we presented the first NMR structure and dynamics of the human ephrin-B2 ectodomain as well as its interaction with the receptor EphB2. Strikingly, the NMR study reveals a picture different from those previously obtained by X-ray crystallography. Although in solution it still adopts the same Greek key fold, with the central beta-barrel ( approximately 30% of the molecule) highly similar to that in crystal structures, the other regions are highly dynamic and accessible to the bulk solvent. In particular, the functionally critical C-D and G-H loops of the ephrin-B2 ectodomain are highly flexible as reflected by several NMR probes including hydrogen exchange and (15)N backbone relaxation data. Nevertheless, as revealed by ITC and NMR, the ephrin-B2 ectodomain binds to EphB2 with a K(d) of 22.3 nM to form a tight complex in which the tip of the C-D loop and the C-terminus still remain largely flexible. The present results may bear critical implications in understanding the molecular details as well as designing antagonists of therapeutic interest for Eph-ephrin interactions.

Structure and dynamics of Ca2+-binding domain 1 of the Na+/Ca2+ exchanger in the presence and in the absence of Ca2+

The Na(+)/Ca(2+) exchanger is the major exporter of Ca(2+) across the cell membrane of cardiomyocytes. The activity of the exchanger is regulated by a large intracellular loop that contains two Ca(2+)-binding domains, calcium-binding domain (CBD) 1 and CBD2. CBD1 binds Ca(2+) with much higher affinity than CBD2 and is considered to be the primary Ca(2+) sensor. The effect of Ca(2+) on the structure and dynamics of CBD1 has been characterized by NMR spectroscopy using chemical shifts, residual dipolar couplings, and spin relaxation. Residual dipolar couplings are used in a new way for residue selection in the determination of the anisotropic rotational diffusion tensor from spin relaxation data. The results provide a highly consistent description across these complementary data sets and show that Ca(2+) binding is accompanied by a selective conformational change among the binding site residues. Residues that exhibit a significant conformational change are also sites of altered dynamics. In particular, Ca(2+) binding restricts the mobility of the major acidic segment and affects the dynamics of several nearby binding loops. These observations indicate that Ca(2+) elicits a local transition to a well-ordered coordination geometry in the CBD1-binding site.

Protein-ligand NOE matching: a high-throughput method for binding pose evaluation that does not require protein NMR resonance assignments

Given the three-dimensional (3D) structure of a protein, the binding pose of a ligand can be determined using distance restraints derived from assigned intra-ligand and protein-ligand nuclear Overhauser effects (NOEs). A primary limitation of this approach is the need for resonance assignments of the ligand-bound protein. We have developed an approach that utilizes data from 3D 13C-edited, 13C/15N-filtered HSQC-NOESY spectra for evaluating ligand binding poses without requiring protein NMR resonance assignments. Only the 1H NMR assignments of the bound ligand are essential. Trial ligand binding poses are generated by any suitable method (e.g., computational docking). For each trial binding pose, the 3D 13C-edited, 13C/15N-filtered HSQC-NOESY spectrum is predicted, and the predicted and observed patterns of protein-ligand NOEs are matched and scored using a fast, deterministic bipartite graph matching algorithm. The best scoring (lowest "cost") poses are identified. Our method can incorporate any explicit restraints or protein assignment data that are available, and many extensions of the basic procedure are feasible. Only a single sample is required, and the method can be applied to both slowly and rapidly exchanging ligands. The method was applied to three test cases: one complex involving muscle fatty acid-binding protein (mFABP) and two complexes involving the leukocyte function-associated antigen 1 (LFA-1) I-domain. Without using experimental protein NMR assignments, the method identified the known binding poses with good accuracy. The addition of experimental protein NMR assignments improves the results. Our "NOE matching" approach is expected to be widely applicable; i.e., it does not appear to depend on a fortuitous distribution of binding pocket residues.

The neurobiologist's guide to structural biology: a primer on why macromolecular structure matters and how to evaluate structural data

Structural biology now plays a prominent role in addressing questions central to understanding how excitable cells function. Although interest in the insights gained from the definition and dissection of macromolecular anatomy is high, many neurobiologists remain unfamiliar with the methods employed. This primer aims to help neurobiologists understand approaches for probing macromolecular structure and where the limits and challenges remain. Using examples of macromolecules with neurobiological importance, the review covers X-ray crystallography, electron microscopy (EM), small-angle X-ray scattering (SAXS), and nuclear magnetic resonance (NMR) and biophysical methods with which these approaches are often paired: isothermal titration calorimetry (ITC), equilibrium analytical ultracentifugation, and molecular dynamics (MD).

The prediction of the wild-type telomerase RNA pseudoknot structure and the pivotal role of the bulge in its formation

In this study, the three-dimensional structure of the wild-type human telomerase RNA pseudoknot was predicted via molecular modeling. The wild-type pseudoknot structure is then compared to the recent NMR solution structure of the telomerase pseudoknot, which does not contain the U177 bulge. The removal of the bulge from the pseudoknot structure results in higher stability and significant reduction of activity of telomerase. We show that the effect of the bulge on the structure results in a significant transformation of the pseudoknot junction region where the starting base pairs are disrupted and unique triple base pairs are formed. We found that the formation of the junction region is greatly influenced by interactions of the U177 bulge with loop residues and rotation of residue A174. Moreover, this is the first study to our knowledge where a structure as complex as the pseudoknot has been solved by purely theoretical methods.

Molecular Fragment Replacement Approach to Protein Structure Determination by Chemical Shift and Dipolar Homology Database Mining

A novel approach is described for determining backbone structures of proteins that is based on finding fragments in the protein data bank (PDB). For each fragment in the target protein, usually chosen to be 7-10 residues in length, PDB fragments are selected that best fit to experimentally determined one-bond heteronuclear dipolar couplings and that show agreement between chemical shifts predicted for the PDB fragment and experimental values for the target fragment. These fragments are subsequently refined by simulated annealing to improve agreement with the experimental data. If the lowest-energy refined fragments form a unique structural cluster, this structure is accepted and side chains are added on the basis of a conformational database potential. The sequential backbone assembly process extends the chain by translating an accepted fragment onto it. For several small proteins, with extensive sets of dipolar couplings measured in two alignment media, a unique final structure is obtained that agrees well with structures previously solved by conventional methods. With less dipolar input data, large, oriented fragments of each protein are obtained, but their relative positioning requires either a small set of translationally restraining nuclear Overhauser enhancements (NOEs) or a protocol that optimizes burial of hydrophobic groups and pairing of beta-strands.

Separation and detection methods for covalent drug–protein adducts

Covalent binding of reactive metabolites of drugs to proteins has been a predominant hypothesis for the mechanism of toxicity caused by numerous drugs. The development of efficient and sensitive analytical methods for the separation, identification, quantification of drug-protein adducts have important clinical and toxicological implications. In the last few decades, continuous progress in analytical methodology has been achieved with substantial increase in the number of new, more specific and more sensitive methods for drug-protein adducts. The methods used for drug-protein adduct studies include those for separation and for subsequent detection and identification. Various chromatographic (e.g., affinity chromatography, ion-exchange chromatography, and high-performance liquid chromatography) and electrophoretic techniques [e.g., sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), two-dimensional SDS-PAGE, and capillary electrophoresis], used alone or in combination, offer an opportunity to purify proteins adducted by reactive drug metabolites. Conventionally, mass spectrometric (MS), nuclear magnetic resonance, and immunological and radioisotope methods are used to detect and identify protein targets for reactive drug metabolites. However, these methods are labor-intensive, and have provided very limited sequence information on the target proteins adducted, and thus the identities of the protein targets are usually unknown. Moreover, the antibody-based methods are limited by the availability, quality, and specificity of antibodies to protein adducts, which greatly hindered the identification of specific protein targets of drugs and their clinical applications. Recently, the use of powerful MS technologies (e.g., matrix-assisted laser desorption/ionization time-of-flight) together with analytical proteomics have enabled one to separate, identify unknown protein adducts, and establish the sequence context of specific adducts by offering the opportunity to search for adducts in proteomes containing a large number of proteins with protein adducts and unmodified proteins. The present review highlights the separation and detection technologies for drug-protein adducts, with an emphasis on methodology, advantages and limitations to these techniques. Furthermore, a brief discussion of the application of these techniques to individual drugs and their target proteins will be outlined.

NMR and structural genomics

The role of NMR in structural genomics is outlined, with particular emphasis on using protein domains as targets. Strategies for domain expression, characterization, and labeling are presented.

Molecular dynamics of biological macromolecules: a brief history and perspective

A description of the origin of my interest in and the development of molecular dynamics simulations of biomolecules is presented with a historical overview, including the role of my interactions with Shneior Lifson and his group in Israel. Some early applications of the methodology by members of my group are summarized, followed by a description of examples of recent applications and some discussion of possible future directions.

Improving the accuracy of NMR structures of RNA by means of conformational database potentials of mean force as assessed by complete dipolar coupling cross-validation

The description of the nonbonded contact terms used in simulated annealing refinement can have a major impact on nucleic acid structures generated from NMR data. Using complete dipolar coupling cross-validation, we demonstrate that substantial improvements in coordinate accuracy of NMR structures of RNA can be obtained by making use of two conformational database potentials of mean force: a nucleic acid torsion angle database potential consisting of various multidimensional torsion angle correlations; and an RNA specific base-base positioning potential that provides a simple geometric, statistically based, description of sequential and nonsequential base-base interactions. The former is based on 416 nucleic acid crystal structures solved at a resolution of </=2 A and an R-factor </=25%; the latter is based on 131 RNA crystal structures solved at a resolution of </=3 A and an R-factor of </=25%, and includes both the large and small subunits of the ribosome. The application of these two database potentials is illustrated for the structure refinement of an RNA aptamer/theophylline complex for which extensive NOE and residual dipolar coupling data have been measured in solution.