Backbone and side chain dynamics of mutant calmodulin-peptide complexes

Biomolecular NMR spectroscopy in the era of artificial intelligence

Biomolecular nuclear magnetic resonance (NMR) spectroscopy and artificial intelligence (AI) have a burgeoning synergy. Deep learning-based structural predictors have forever changed structural biology, yet these tools currently face limitations in accurately characterizing protein dynamics, allostery, and conformational heterogeneity. We begin by highlighting the unique abilities of biomolecular NMR spectroscopy to complement AI-based structural predictions toward addressing these knowledge gaps. We then highlight the direct integration of deep learning approaches into biomolecular NMR methods. AI-based tools can dramatically improve the acquisition and analysis of NMR spectra, enhancing the accuracy and reliability of NMR measurements, thus streamlining experimental processes. Additionally, deep learning enables the development of novel types of NMR experiments that were previously unattainable, expanding the scope and potential of biomolecular NMR spectroscopy. Ultimately, a combination of AI and NMR promises to further revolutionize structural biology on several levels, advance our understanding of complex biomolecular systems, and accelerate drug discovery efforts.

Extending the Detection Limit in Fragment Screening of Proteins Using Reverse Micelle Encapsulation

Detection of very weak (K_d > 10 mM) interactions of proteins with small molecules has been elusive. This is particularly important for fragment-based drug discovery, where it is suspected that the majority of potentially useful fragments will be invisible to current screening methodologies. We describe an NMR approach that permits detection of protein-fragment interactions in the very low affinity range and extends the current detection limit of ∼10 mM up to ∼200 mM and beyond. Reverse micelle encapsulation is leveraged to effectively reach very high fragment and protein concentrations, a principle that is validated by binding model fragments to E. coli dihydrofolate reductase. The method is illustrated by target-detected screening of a small polar fragment library against interleukin-1β, which lacks a known ligand-binding pocket. Evaluation of binding by titration and structural context allows for validation of observed hits using rigorous structural and statistical criteria. The 21 curated hit molecules represent a remarkable hit rate of nearly 10% of the library. Analysis shows that fragment binding involves residues comprising two-thirds of the protein's surface. Current fragment screening methods rely on detection of relatively tight binding to ligand binding pockets. The method presented here illustrates a potential to faithfully discover starting points for development of small molecules that bind to a desired region of the protein, even if the targeted region is defined by a relatively flat surface.

Allostery in Orai1 binding to calmodulin revealed from conformational thermodynamics

Here, we study microscopic mechanism of complex formation between Ca²⁺-bound calmodulin (holoCaM) and Orai1 that regulates Ca²⁺-dependent inactivation process in eukaryotic cells. We compute conformational thermodynamic changes in holoCaM with respect to complex of Orai1 bound to C-terminal domain of holoCaM using histograms of dihedral angles of the proteins over trajectories from molecular dynamics simulations. Our analysis shows that the N-terminal domain residues L4, T5, Q41, N42, T44 and E67 of holoCaM get destabilized and disordered due to Orai1 binding to C-terminal domain of calmodulin affect the N-terminal domain residues. Among these residues, polar T44, having maximum destabilization and disorder via backbone fluctuations, shows the largest change in solvent exposure. This suggests that N-terminal domain is allosterically regulated via T44 by the binding of Orai1 to the C-terminal domain.

Binding Mechanism of the N-Terminal SH3 Domain of CrkII and Proline-Rich Motifs in cAbl

The N-terminal Src homology 3 (nSH3) domain of a signaling adaptor protein, CT-10 regulator of kinase II (CrkII), recognizes proline-rich motifs (PRMs) of binding partners, such as cAbl kinase. The interaction between CrkII and cAbl kinase is involved in the regulation of cell spreading, microbial pathogenesis, and cancer metastasis. Here, we report the detailed biophysical characterizations of the interactions between the nSH3 domain of CrkII and PRMs in cAbl. We identified that the nSH3 domain of CrkII binds to three PRMs in cAbl with virtually identical affinities. Structural studies, by using x-ray crystallography and NMR spectroscopy, revealed that the binding modes of all three nSH3:PRM complexes are highly similar to each other. Van 't Hoff analysis revealed that nSH3:PRM interaction is associated with favorable enthalpy and unfavorable entropy change. The combination of experimentally determined thermodynamic parameters, structure-based calculations, and (15)N NMR relaxation analysis highlights the energetic contribution of conformational entropy change upon the complex formation, and water molecules structured in the binding interface of the nSH3:PRM complex. Understanding the molecular basis of nSH3:PRM interaction will provide, to our knowledge, new insights for the rational design of small molecules targeting the interaction between CrkII and cAbl.

Disorder transitions and conformational diversity cooperatively modulate biological function in proteins

Structural differences between conformers sustain protein biological function. Here, we studied in a large dataset of 745 intrinsically disordered proteins, how ordered-disordered transitions modulate structural differences between conformers as derived from crystallographic data. We found that almost 50% of the proteins studied show no transitions and have low conformational diversity while the rest show transitions and a higher conformational diversity. In this last subset, 60% of the proteins become more ordered after ligand binding, while 40% more disordered. As protein conformational diversity is inherently connected with protein function our analysis suggests differences in structure-function relationships related to order-disorder transitions.

The ensemble nature of allostery

Allostery is the process by which biological macromolecules (mostly proteins) transmit the effect of binding at one site to another, often distal, functional site, allowing for regulation of activity. Recent experimental observations demonstrating that allostery can be facilitated by dynamic and intrinsically disordered proteins have resulted in a new paradigm for understanding allosteric mechanisms, which focuses on the conformational ensemble and the statistical nature of the interactions responsible for the transmission of information. Analysis of allosteric ensembles reveals a rich spectrum of regulatory strategies, as well as a framework to unify the description of allosteric mechanisms from different systems.

Structure and dynamics of calmodulin (CaM) bound to nitric oxide synthase peptides: effects of a phosphomimetic CaM mutation

Nitric oxide synthase (NOS) plays a major role in a number of key physiological and pathological processes. Knowledge of how this is regulated is important. The small acidic calcium binding protein, calmodulin (CaM), is required to fully activate the enzyme. The exact mechanism of how CaM activates NOS is not fully understood. Studies have shown CaM to act like a switch that causes a conformational change in NOS to allow for the transfer of an electron between the reductase and oxygenase domains through a process that is thought to be highly dynamic. To investigate the dynamic properties of CaM-NOS interactions, we determined the solution structure of CaM bound to the inducible NOS (iNOS) and endothelial NOS (eNOS) CaM binding region peptides. In addition, we investigated the effect of CaM phosphorylation. Tyrosine 99 (Y99) of CaM is reported to be phosphorylated in vivo. We have produced a phosphomimetic Y99E CaM to investigate the structural and functional effects that the phosphorylation of this residue may have on nitric oxide production. All three mammalian NOS isoforms were included in the investigation. Our results show that a phosphomimetic Y99E CaM significantly reduces the maximal synthase activity of eNOS by 40% while having little effect on nNOS or iNOS activity. A comparative nuclear magnetic resonance study between phosphomimetic Y99E CaM and wild-type CaM bound to the eNOS CaM binding region peptide was performed. This investigation provides important insights into how the increased electronegativity of a phosphorylated CaM protein affects the binding, dynamics, and activation of the NOS enzymes.

Changes in Lysozyme Flexibility upon Mutation Are Frequent, Large and Long-Ranged

We investigate changes in human c-type lysozyme flexibility upon mutation via a Distance Constraint Model, which gives a statistical mechanical treatment of network rigidity. Specifically, two dynamical metrics are tracked. Changes in flexibility index quantify differences within backbone flexibility, whereas changes in the cooperativity correlation quantify differences within pairwise mechanical couplings. Regardless of metric, the same general conclusions are drawn. That is, small structural perturbations introduced by single point mutations have a frequent and pronounced affect on lysozyme flexibility that can extend over long distances. Specifically, an appreciable change occurs in backbone flexibility for 48% of the residues, and a change in cooperativity occurs in 42% of residue pairs. The average distance from mutation to a site with a change in flexibility is 17–20 Å. Interestingly, the frequency and scale of the changes within single point mutant structures are generally larger than those observed in the hen egg white lysozyme (HEWL) ortholog, which shares 61% sequence identity with human lysozyme. For example, point mutations often lead to substantial flexibility increases within the β-subdomain, which is consistent with experimental results indicating that it is the nucleation site for amyloid formation. However, β-subdomain flexibility within the human and HEWL orthologs is more similar despite the lowered sequence identity. These results suggest compensating mutations in HEWL reestablish desired properties.

The functional importance of protein dynamics is universally accepted, making the study of dynamical similarities and differences among proteins of the same function an intriguing problem. While some metrics are likely to be conserved across family, differences are also very common. In previous works we have used a Distance Constraint Model to quantify flexibility differences across sets of orthologous proteins, which reproduce this diversity. In the same manner, this work investigates changes occurring upon individual point mutations. Somewhat surprisingly, the small structural perturbations caused by mutation lead to changes throughout the protein. These changes can be quite large, actually surpassing the scale for differences between ortholog pairs. Moreover, changes in flexibility frequently occur at sites far from the mutation site. These results underscore the sensitivity of protein dynamics in connection with allostery, and help explain why differences across protein families are so common.

Exploring the role of structure and dynamics in the function of chymotrypsin inhibitor 2

Increasing awareness of the possible role of internal dynamics in protein function has led to the development of new methods for experimentally characterizing protein dynamics across multiple time scales, especially using NMR spectroscopy. A few analyses of the conformational dynamics of proteins ranging from nonallosteric single domains to multidomain allosteric enzymes are now available; however, demonstrating a connection between dynamics and function remains difficult on account of the comparative lack of studies examining both changes in dynamics and changes in function in response to the same perturbations. In previous work, we characterized changes in structure and dynamics on the ps–ns time scale resulting from hydrophobic core mutations in chymotrypsin inhibitor 2 and found that there are moderate, persistent global changes in dynamics in the absence of gross structural changes (Whitley et al., Biochemistry 2008;47:8566–8576). Here, we assay those and additional mutants for inhibitory ability toward the serine proteases elastase and chymotrypsin to determine the effects of mutation on function. Results indicate that core mutation has only a subtle effect on CI2 function. Using chemical shifts, we also studied the effect of complex formation on CI2 structure and found that perturbations are greatest at the complex interface but also propagate toward CI2's hydrophobic core. The structure–dynamics–function data set completed here suggests that dynamics plays a limited role in the function of this small model system, although we do observe a correlation between nanosecond-scale reactive loop motions and inhibitory ability for mutations at one key position in the hydrophobic core.

Probing slow protein dynamics by adiabatic R(1rho) and R(2rho) NMR experiments

Slow micros/ms dynamics involved in protein folding, binding, catalysis, and allostery are currently detected using NMR dispersion experiments such as CPMG (Carr-Purcell-Meiboom-Gill) or spin-lock R(1rho). In these methods, protein dynamics are obtained by analyzing relaxation dispersion curves obtained from either changing the time spacing between 180 degree pulses or by changing the effective spin-locking field strength. In this Communication, we introduce a new method to induce a dispersion of relaxation rates. Our approach relies on altering the shape of the adiabatic full passage pulse and is conceptually different from existing approaches. By changing the nature of the adiabatic radiofrequency irradiation, we are able to obtain rotating frame R(1rho) and R(2rho) dispersion curves that are sensitive to slow micros/ms protein dynamics (demonstrated with ubiquitin). The strengths of this method are to (a) extend the dynamic range of the relaxation dispersion analysis, (b) avoid the need for multiple magnetic field strengths to extract dynamic parameters, (c) measure accurate relaxation rates that are independent of frequency offset, and (d) reduce the stress to NMR hardware (e.g., cryoprobes).

Detection of native-state nonadditivity in double mutant cycles via hydrogen exchange

Proteins have evolved to exploit long-range structural and dynamic effects as a means of regulating function. Understanding communication between sites in proteins is therefore vital to our comprehension of such phenomena as allostery, catalysis, and ligand binding/ejection. Double mutant cycle analysis has long been used to determine the existence of communication between pairs of sites, proximal or distal, in proteins. Typically, nonadditivity (or "thermodynamic coupling") is measured from global transitions in concert with a single probe. Here, we have applied the atomic resolution of NMR in tandem with native-state hydrogen exchange (HX) to probe the structure/energy landscape for information transduction between a large number of distal sites in a protein. Considering the event of amide proton exchange as an energetically quantifiable structural perturbation, m n-dimensional cycles can be constructed from mutation of n-1 residues, where m is the number of residues for which HX data is available. Thus, efficient mapping of a large number of couplings is made possible. We have applied this technique to one additive and two nonadditive double mutant cycles in a model system, eglin c. We find heterogeneity of HX-monitored couplings for each cycle, yet averaging results in strong agreement with traditionally measured values. Furthermore, long-range couplings observed at locally exchanging residues indicate that the basis for communication can occur within the native state ensemble, a conclusion not apparent from traditional measurements. We propose that higher-order couplings can be obtained and show that such couplings provide a mechanistic basis for understanding lower-order couplings via "spheres of perturbation". The method is presented as an additional tool for identifying a large number of couplings with greater coverage of the protein of interest.

Protein functional landscapes, dynamics, allostery: a tortuous path towards a universal theoretical framework

Energy landscape theories have provided a common ground for understanding the protein folding problem, which once seemed to be overwhelmingly complicated. At the same time, the native state was found to be an ensemble of interconverting states with frustration playing a more important role compared to the folding problem. The landscape of the folded protein – the native landscape – is glassier than the folding landscape; hence, a general description analogous to the folding theories is difficult to achieve. On the other hand, the native basin phase volume is much smaller, allowing a protein to fully sample its native energy landscape on the biological timescales. Current computational resources may also be used to perform this sampling for smaller proteins, to build a ‘topographical map’ of the native landscape that can be used for subsequent analysis. Several major approaches to representing this topographical map are highlighted in this review, including the construction of kinetic networks, hierarchical trees and free energy surfaces with subsequent structural and kinetic analyses. In this review, we extensively discuss the important question of choosing proper collective coordinates characterizing functional motions. In many cases, the substates on the native energy landscape, which represent different functional states, can be used to obtain variables that are well suited for building free energy surfaces and analyzing the protein's functional dynamics. Normal mode analysis can provide such variables in cases where functional motions are dictated by the molecule's architecture. Principal component analysis is a more expensive way of inferring the essential variables from the protein's motions, one that requires a long molecular dynamics simulation. Finally, the two popular models for the allosteric switching mechanism, ‘preexisting equilibrium’ and ‘induced fit’, are interpreted within the energy landscape paradigm as extreme points of a continuum of transition mechanisms. Some experimental evidence illustrating each of these two models, as well as intermediate mechanisms, is presented and discussed.

The role of conformational entropy in molecular recognition by calmodulin

The physical basis for high affinity interactions involving proteins is complex and potentially involves a range of energetic contributions. Among these are changes in protein conformational entropy, which cannot yet be reliably computed from molecular structures. We have recently employed changes in conformational dynamics as a proxy for changes in conformational entropy of calmodulin upon association with domains from regulated proteins. The apparent change in conformational entropy was linearly related to the overall binding entropy. This view warrants a more quantitative foundation. Here we calibrate an “entropy meter” employing an experimental dynamical proxy based on NMR relaxation and show that changes in the conformational entropy of calmodulin are a significant component of the energetics of binding. Furthermore, the distribution of motion at the interface between the target domain and calmodulin are surprisingly non-complementary. These observations promote modification of our understanding of the energetics of protein-ligand interactions.

Nuclear magnetic resonance study of the role of M42 in the solution dynamics of Escherichia coli dihydrofolate reductase

It is widely recognized that key positions throughout a protein's structure contribute unequally to function. In light of recent studies that suggest protein dynamics are required for function, a number of these residues may serve to promote motions required for ligand binding and catalysis. In this nuclear magnetic resonance (NMR) study, the conformational dynamics of the dihydrofolate reductase (DHFR) mutant M42W, in the presence of methotrexate and NADPH, are characterized and compared to those of the wild-type enzyme. M42 is distal to the active site, yet the M42W substitution regulates catalysis and ligand affinity and is therefore analogous to an allosteric modulator of DHFR function. To gain understanding of how this mutation regulates activity, we employ a "pandynamic" strategy by measuring conformational fluctuations of backbone amide and side-chain methyl groups on multiple time scales. Changes in pico- to nanosecond dynamics indicate that the mutational effects are propagated throughout a network of interacting residues within DHFR, consistent with a role for M42 as a dynamic communication hub. On the micro- to millisecond time scale, mutation increases the rate of switching in the catalytic core. Mutation also introduces switching in the adenosine binding subdomain that occurs at a higher frequency than in the catalytic core and which correlates with the rate of product release for M42W-DHFR. Finally, a structurally inferred analysis of side-chain dynamics suggests that the M42W mutation dampens motional contributions from nonlocal sources. These data show that the M42W mutation alters the dynamics of DHFR and are consistent with theoretical analysis that suggests this mutation disrupts motion that promotes catalysis.

Global and local mobility of apocalmodulin monitored through fast-field cycling relaxometry

Calmodulin (CaM) is a ubiquitous eukaryotic protein with two conformationally independent domains that can bind up to two calcium ions each. In the calcium-bound state, CaM is able to regulate a vast number of cellular activities by binding to a multiplicity of target proteins in different modes. Its versatility has been ascribed to its anomalously high flexibility. The calcium-free form (apoCaM), which is the resting state of CaM in cells, is also able to functionally bind a number of protein targets, but its dynamics has received less attention. At variance with the calcium-bound form, the crystal structure of apoCaM shows a compact organization of the two domains, but NMR measurements could not detect any contact between them, thus indicating the presence of mobility in solution. The mobility of apoCaM is here investigated through protein proton relaxation rate measurements performed with a high-sensitivity fast-field cycling relaxometer. Such measurements provide direct access to the spectral density function and show that 1), the reorientation time is in agreement with a closed form of the protein; but 2), the collective order parameter is much smaller than for other well folded compact proteins, indicating that a remarkably large side-chain mobility must be present.

Detection of nanosecond time scale side-chain jumps in a protein dissolved in water/glycerol solvent

In solution, the correlation time of the overall protein tumbling, tau(R), plays a role of a natural dynamics cutoff-internal motions with correlation times on the order of tau ( R ) or longer cannot be reliably identified on the basis of spin relaxation data. It has been proposed some time ago that the 'observation window' of solution experiments can be expanded by changing the viscosity of solvent to raise the value of tau(R). To further explore this concept, we prepared a series of samples of alpha-spectrin SH3 domain in solvent with increasing concentration of glycerol. In addition to the conventional (15)N labeling, the protein was labeled in the Val, Leu methyl positions ((13)CHD(2) on a deuterated background). The collected relaxation data were used in asymmetric fashion: backbone (15)N relaxation rates were used to determine tau(R) across the series of samples, while methyl (13)C data were used to probe local dynamics (side-chain motions). In interpreting the results, it has been initially suggested that addition of glycerol leads only to increases in tau(R), whereas local motional parameters remain unchanged. Thus the data from multiple samples can be analyzed jointly, with tau(R) playing the role of experimentally controlled variable. Based on this concept, the extended model-free model was constructed with the intent to capture the effect of ns time-scale rotameric jumps in valine and leucine side chains. Using this model, we made a positive identification of nanosecond dynamics in Val-23 where ns motions were already observed earlier. In several other cases, however, only tentative identification was possible. The lack of definitive results was due to the approximate character of the model-contrary to what has been assumed, addition of glycerol led to a gradual 'stiffening' of the protein. This and other observations also shed light on the interaction of the protein with glycerol, which is one of the naturally occurring osmoprotectants. In particular, it has been found that the overall protein tumbling is controlled by the bulk solvent, and not by a thin solvation layer which contains a higher proportion of water.

Frameworks for understanding long-range intra-protein communication

The phenomenon of intra-protein communication is fundamental to such processes as allostery and signaling, yet comparatively little is understood about its physical origins despite notable progress in recent years. This review introduces contemporary but distinct frameworks for understanding intra-protein communication by presenting both the ideas behind them and a discussion of their successes and shortcomings. The first framework holds that intra-protein communication is accomplished by the sequential mechanical linkage of residues spanning a gap between distal sites. According to the second framework, proteins are best viewed as ensembles of distinct structural microstates, the dynamical and thermodynamic properties of which contribute to the experimentally observable macroscale properties. Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying intra-protein communication, and the insights into both frameworks it provides are presented through a discussion of numerous examples from the literature. Distinct from mechanical and thermodynamic considerations of intra-protein communication are recently applied graph and network theoretic analyses. These computational methods reduce complex three dimensional protein architectures to simple maps comprised of nodes (residues) connected by edges (inter-residue "interactions"). Analysis of these graphs yields a characterization of the protein's topology and network characteristics. These methods have shown proteins to be "small world" networks with moderately high local residue connectivities existing concurrently with a small but significant number of long range connectivities. However, experimental studies of the tantalizing idea that these putative long range interaction pathways facilitate one or several macroscopic protein characteristics are unfortunately lacking at present. This review concludes by comparing and contrasting the presented frameworks and methodologies for studying intra-protein communication and suggests a manner in which they can be brought to bear simultaneously to further enhance our understanding of this important fundamental phenomenon.

Dynamic dysfunction in dihydrofolate reductase results from antifolate drug binding: modulation of dynamics within a structural state

The arduous task of rationally designing small-molecule enzyme inhibitors is complicated by the inherent flexibility of the protein scaffold. To gain insight into the changes in dynamics associated with small-molecule-based inhibition, we have characterized, using NMR spectroscopy, Escherichia coli dihydrofolate reductase in complex with two drugs: methotrexate and trimethoprim. The complexes allowed the intrinsic dynamic effects of drug binding to be revealed within the context of the "closed" structural ensemble. Binding of both drugs results in an identical decoupling of global motion on the micro- to millisecond timescale. Consistent with a change in overall dynamic character, the drugs' perturbations to pico- to nanosecond backbone and side-chain methyl dynamics are also highly similar. These data show that the inhibitors simultaneously modulate slow concerted switching and fast motions at distal regions of dihydrofolate reductase, providing a dynamic link between the substrate binding site and distal loop residues known to affect catalysis.

Allostery and cooperativity revisited

Although phenomenlogical models that account for cooperativity in allosteric systems date back to the early and mid-60's (e.g., the KNF and MWC models), there is resurgent interest in the topic due to the recent experimental and computational studies that attempted to reveal, at an atomistic level, how allostery actually works. In this review, using systems for which atomistic simulations have been carried out in our groups as examples, we describe the current understanding of allostery, how the mechanisms go beyond the classical MWC/Pauling-KNF descriptions, and point out that the "new view" of allostery, emphasizing "population shifts," is, in fact, an "old view." The presentation offers not only an up-to-date description of allostery from a theoretical/computational perspective, but also helps to resolve several outstanding issues concerning allostery.

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.

Compensatory and long-range changes in picosecond-nanosecond main-chain dynamics upon complex formation: 15N relaxation analysis of the free and bound states of the ubiquitin-like domain of human plexin-B1 and the small GTPase Rac1

The formation of a complex between Rac1 and the cytoplasmic domain of plexin-B1 is one of the first documented cases of a direct interaction between a small guanosine 5'-triphosphatase (GTPase) and a transmembrane receptor. Structural studies have begun to elucidate the role of this interaction for the signal transduction mechanism of plexins. Mapping of the Rac1 GTPase surface that contacts the Rho GTPase binding domain of plexin-B1 by solution NMR spectroscopy confirms the plexin domain as a GTPase effector protein. Regions neighboring the GTPase switch I and II regions are also involved in the interaction and there is considerable interest to examine the changes in protein dynamics that take place upon complex formation. Here we present main-chain nitrogen-15 relaxation measurements for the unbound proteins as well as for the Rho GTPase binding domain and Rac1 proteins in their complexed state. Derived order parameters, S2, show that considerable motions are maintained in the bound state of plexin. In fact, some of the changes in S2 on binding appear compensatory, exhibiting decreased as well as increased dynamics. Fluctuations in Rac1, already a largely rigid protein on the picosecond-nanosecond timescale, are overall diminished, but isomerization dynamics in the switch I and II regions of the GTPase are retained in the complex and appear to be propagated to the bound plexin domain. Remarkably, fluctuations in the GTPase are attenuated at sites, including helices alpha6 (the Rho-specific insert helix), alpha7 and alpha8, that are spatially distant from the interaction region with plexin. This effect of binding on long-range dynamics appears to be communicated by hinge sites and by subtle conformational changes in the protein. Similar to recent studies on other systems, we suggest that dynamical protein features are affected by allosteric mechanisms. Altered protein fluctuations are likely to prime the Rho GTPase-plexin complex for interactions with additional binding partners.

A Computational Model for the Prediction of Aqueous Solubility That Includes Crystal Packing, Intrinsic Solubility, and Ionization Effects

The optimization of aqueous solubility is an important step along the route to bringing a new therapeutic to market. We describe the development of an empirical computational model to rank the pH-dependent aqueous solubility of drug candidates. The model consists of three core components to describe aqueous solubility. The first is a multivariate QSAR model for the prediction of the intrinsic solubility of the neutral solute. The second facet of the approach is the consideration of ionization using a predicted pKa and the Henderson-Hasselbalch equation. The third aspect of the model is a novel method for assessing the effects of crystal packing on solubility through a series of short molecular dynamics simulations of an actual or hypothetical small molecule crystal structure at escalating temperatures. The model also includes a Monte Carlo error function that considers the variability of each of the underlying components of the model to estimate the 90% confidence interval of estimation.

Estimation of hERG inhibition of drug candidates using multivariate property and pharmacophore SAR

We describe the development of a computational model for the prediction of the inhibition of K(+) flow through the hERG ion channel. Using a collection of 1075 discovery compounds with hERG inhibition measured in our standard patch-clamp electrophysiology assay, molecular features important for drug-induced inhibition were identified using a combination of statistical inference algorithms and manual hypothesis generation and testing. While many of the features used in the model reflect those referenced in the literature, several aspects of the model provide new insight into the role of physicochemical properties, electrostatics, and novel pharmacophores in hERG inhibition. Coefficients for these 10 features were then determined by least median squares regression, resulting in a model with an R(2) approximately 0.66 and RMS error (RMSe) of 0.47 log units for an external test set. Significant additional validation performed using a large collection of subsequent discovery data has been very encouraging with an R(2)=0.54 and an RMSe of 0.63 log units. The performance of the model across several different chemotypes is demonstrated and discussed.

Pump-probe molecular dynamics as a tool for studying protein motion and long range coupling

A new method for analyzing the dynamics of proteins is developed and tested. The method, pump-probe molecular dynamics, excites selected atoms or residues with a set of oscillating forces, and the transmission of the impulse to other parts of the protein is probed using Fourier transform of the atomic motions. From this analysis, a coupling profile can be determined which quantifies the degree of interaction between pump and probe residues. Various physical properties of the method such as reciprocity and speed of transmission are examined to establish the soundness of the method. The coupling strength can be used to address questions such as the degree of interaction between different residues at the level of dynamics, and identify propagation of influence of one part of the protein on another via "pathways" through the protein. The method is illustrated by analysis of coupling between different secondary structure elements in the allosteric protein calmodulin, and by analysis of pathways of residue-residue interaction in the PDZ domain protein previously elucidated by genomics and mutational studies.

A mutation in the S-switch region of the Runt domain alters the dynamics of an allosteric network responsible for CBFbeta regulation

The Runt domain is the DNA binding domain of the core binding factor (CBF) Runx subunits. The CBFs are transcription factors that play critical roles in hematopoiesis, bone, and neuron development in mammals. A common non-DNA binding CBFbeta subunit heterodimerizes with the Runt domain of the Runx proteins and allosterically regulates its affinity for DNA. Previous NMR dynamics studies suggested a model whereby CBFbeta allosterically regulates DNA binding by quenching conformational exchange in the Runt domain, particularly in the S-switch region and the betaE'-F loop. We sought to test this model, and to this end introduced all possible single amino acid substitutions into the S-switch region and the betaE'-F loop, and screened for mutations that enhanced DNA-binding. We demonstrate that one Runt domain mutant, R164N, binds both DNA and CBFbeta with higher affinity, but it is less sensitive to allosteric regulation by CBFbeta. Analysis of NMR relaxation data shows that the chemical exchange exhibited by the wild-type Runt domain is largely quenched by the R164N substitution. These data support a model in which the dynamic behavior of a network of residues connecting the CBFbeta and DNA binding sites on the Runt domain plays a critical role in the mechanism of allosteric regulation. This study provides an important functional link between dynamic behavior and protein allosteric function, consistent with results on other allosterically regulated proteins.

Evaluation of energetic and dynamic coupling networks in a PDZ domain protein

A number of computational and experimental studies have identified intramolecular communication "pathways" or "networks" important for transmitting allostery. Here, we have used mutagenesis and NMR relaxation methods to investigate the scope and nature of the communication networks found in the second post-synaptic density-95/discs large/zonula occludens-1 (PDZ) domain of the human protein tyrosine phosphatase 1E protein (hPTP1E) (PDZ2). It was found that most mutations do not have a significant energetic contribution to peptide ligand binding. Three mutants that showed significant changes in binding also displayed context-dependent dynamic effects. Both a mutation at a partially exposed site (H71Y) and a buried core position (I35V) had a limited response in side-chain (2)H-based dynamics when compared to wild-type PDZ2. In contrast, a change at a second core position (I20F) that had previously been shown to be part of an energetic and dynamic network, resulted in extensive changes in side-chain dynamics. This response is reminiscent to that seen previously upon peptide ligand binding. These results shed light on the nature of the PDZ2 dynamic network and suggest that position 20 in PDZ2 acts as a "hub" that is energetically and dynamically critical for transmitting changes in dynamics throughout the PDZ domain.