Folding and binding cascades: dynamic landscapes and population shifts

Effect of the Crystal Environment on Side-Chain Conformational Dynamics in Cyanovirin-N Investigated through Crystal and Solution Molecular Dynamics Simulations

Side chains in protein crystal structures are essential for understanding biochemical processes such as catalysis and molecular recognition. However, crystal packing could influence side-chain conformation and dynamics, thus complicating functional interpretations of available experimental structures. Here we investigate the effect of crystal packing on side-chain conformational dynamics with crystal and solution molecular dynamics simulations using Cyanovirin-N as a model system. Side-chain ensembles for solvent-exposed residues obtained from simulation largely reflect the conformations observed in the X-ray structure. This agreement is most striking for crystal-contacting residues during crystal simulation. Given the high level of correspondence between our simulations and the X-ray data, we compare side-chain ensembles in solution and crystal simulations. We observe large decreases in conformational entropy in the crystal for several long, polar and contacting residues on the protein surface. Such cases agree well with the average loss in conformational entropy per residue upon protein folding and are accompanied by a change in side-chain conformation. This finding supports the application of surface engineering to facilitate crystallization. Our simulation-based approach demonstrated here with Cyanovirin-N establishes a framework for quantitatively comparing side-chain ensembles in solution and in the crystal across a larger set of proteins to elucidate the effect of the crystal environment on protein conformations.

Effects of flexibility and electrostatic interactions on the coupled binding–folding mechanisms of Chz.core and H2A.z–H2B

The intrinsically disordered protein (IDP) Chz.core, which is the interaction core of Chz1, shows binding preference to histone variant H2A.z. The coupled folding–binding mechanism of the complex can be quantified by the free energy landscapes.

Molecular dynamics simulations elucidate conformational selection and induced fit mechanisms in the binding of PD-1 and PD-L1

PD-L1 binds PD-1 through a complex mechanism including both the conformational selection and induced fit pathways.

Information-Driven, Ensemble Flexible Peptide Docking Using HADDOCK

Modeling protein-peptide interactions remains a significant challenge for docking programs due to the inherent highly flexible nature of peptides, which often adopt different conformations whether in their free or bound forms. We present here a protocol consisting of a hybrid approach, combining the most frequently found peptide conformations in complexes with representative conformations taken from molecular dynamics simulations of the free peptide. This approach intends to broaden the range of conformations sampled during docking. The resulting ensemble of conformations is used as a starting point for information-driven flexible docking with HADDOCK. We demonstrate the performance of this protocol on six cases of increasing difficulty, taken from a protein-peptide benchmark set. In each case, we use knowledge of the binding site on the receptor to drive the docking process. In the majority of cases where MD conformations are added to the starting ensemble for docking, we observe an improvement in the quality of the resulting models.

Long-Range Communication between Different Functional Sites in the Picornaviral 3C Protein

The 3C protein is a master regulator of the picornaviral infection cycle, responsible for both cleaving viral and host proteins, and interacting with genomic RNA replication elements. Here we use nuclear magnetic resonance spectroscopy and molecular dynamics simulations to show that 3C is conformationally dynamic across multiple timescales. Binding of peptide and RNA lead to structural dynamics changes at both the protease active site and the RNA-binding site, consistent with these sites being dynamically coupled. Indeed, binding of RNA influences protease activity, and likewise, interactions at the active site affect RNA binding. We propose that RNA and peptide binding re-shapes the conformational energy landscape of 3C to regulate subsequent functions, including formation of complexes with other viral proteins. The observed channeling of the 3C energy landscape may be important for regulation of the viral infection cycle.

Predicting protein conformational changes for unbound and homology docking: learning from intrinsic and induced flexibility

Predicting protein conformational changes from unbound structures or even homology models to bound structures remains a critical challenge for protein docking. Here we present a study directly addressing the challenge by reducing the dimensionality and narrowing the range of the corresponding conformational space. The study builds on cNMA-our new framework of partner- and contact-specific normal mode analysis that exploits encounter complexes and considers both intrinsic and induced flexibility. First, we established over a CAPRI (Critical Assessment of PRedicted Interactions) target set that the direction of conformational changes from unbound structures and homology models can be reproduced to a great extent by a small set of cNMA modes. In particular, homology-to-bound interface root-mean-square deviation (iRMSD) can be reduced by 40% on average with the slowest 30 modes. Second, we developed novel and interpretable features from cNMA and used various machine learning approaches to predict the extent of conformational changes. The models learned from a set of unbound-to-bound conformational changes could predict the actual extent of iRMSD with errors around 0.6 Å for unbound proteins in a held-out benchmark subset, around 0.8 Å for unbound proteins in the CAPRI set, and around 1 Å even for homology models in the CAPRI set. Our results shed new insights into origins of conformational differences between homology models and bound structures and provide new support for the low-dimensionality of conformational adjustment during protein associations. The results also provide new tools for ensemble generation and conformational sampling in unbound and homology docking. Proteins 2017; 85:544-556. © 2016 Wiley Periodicals, Inc.

Conformational selection in amyloid-based immunotherapy: Survey of crystal structures of antibody-amyloid complexes

BACKGROUND

The dominant feature in neurodegenerative diseases is protein aggregations that lead to neuronal loss. Immunotherapies using antibodies or antibody fragments to target the aggregations are a highly perused approach. The molecular mechanisms underlying the amyloid-based immunotherapy are complex. Deciphering the properties of amyloidogenic proteins responsible for these diseases is essential to obtain insights into antibody recognition of the amyloid antigens.

SCOPE OF REVIEW

We systematically explore all available crystal structures of antibody-amyloid complexes related to neurodegenerative diseases, including antibodies that recognize the Aβ peptide, tau protein, prion protein, alpha-synuclein, huntingtin protein (mHTT), and polyglutamine.

MAJOR CONCLUSIONS

We found that antibodies mostly use the conformational selection mechanism to recognize the highly flexible amyloid antigens. In particular, solanezumab bound to Aβ12-28 tripeptide motif conformation (F19F20A21), which is shared with the Aβ42 fibril. This motif, which is trapped by the antibody, may provide the missing link in amyloid formation. Water molecules often bridge between the antibody and amyloid, contributing to the recognition.

GENERAL SIGNIFICANCE

This paper provides the structural basis for antibody recognition of amyloidogenic proteins. The analysis and discussion of known structures are expected to help in the design and optimization of antibodies in neurodegenerative diseases. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.

Investigation of Intradomain Motions of a Y-Family DNA Polymerase during Substrate Binding and Catalysis

DNA polymerases catalyze DNA synthesis through a stepwise kinetic mechanism that begins with binding to DNA, followed by selection, binding, and incorporation of a nucleotide into an elongating primer. It is hypothesized that subtle active site adjustments in a polymerase to align reactive moieties limit the rate of correct nucleotide incorporation. DNA damage can impede this process for many DNA polymerases, causing replication fork stalling, genetic mutations, and potentially cell death. However, specialized Y-family DNA polymerases are structurally evolved to efficiently bypass DNA damage in vivo, albeit at the expense of replication fidelity. Dpo4, a model Y-family polymerase from Sulfolobus solfataricus, has been well-studied kinetically, structurally, and computationally, which yielded a mechanistic understanding of how the Y-family DNA polymerases achieve their unique catalytic properties. We previously employed a real-time Förster resonance energy transfer (FRET) technique to characterize the global conformational motions of Dpo4 during DNA binding as well as nucleotide binding and incorporation by monitoring changes in distance between sites on the polymerase and DNA, and even between domains of Dpo4. Here, we extend the utility of our FRET methodology to observe conformational transitions within individual domains of Dpo4 during DNA binding and nucleotide incorporation. The results of this novel, intradomain FRET approach unify findings from many studies to fully clarify the complex DNA binding mechanism of Dpo4. Furthermore, intradomain motions in the Finger domain during nucleotide binding and incorporation, for the first time, report on the rate-limiting step of a single-nucleotide addition catalyzed by Dpo4.

Structure-guided residence time optimization of a dabigatran reversal agent

Novel oral anticoagulants are effective and safe alternatives to vitamin-K antagonists for anticoagulation therapy. However, anticoagulation therapy in general is associated with an elevated risk of bleeding. Idarucizumab is a reversal agent for the direct thrombin inhibitor, dabigatran etexilate (Pradaxa®) and is currently in Phase 3 studies. Here, we report data on the antibody fragment aDabi-Fab2, a putative backup molecule for idarucizumab. Although aDabi-Fab2 completely reversed effects of dabigatran in a rat model in vivo, we observed significantly reduced duration of action compared to idarucizumab. Rational protein engineering, based on the X-ray structure of aDabi-Fab2, led to the identification of mutant Y103W. The mutant had optimized shape complementarity to dabigatran while maintaining an energetically favored hydrogen bond. It displayed increased affinity for dabigatran, mainly driven by a slower off-rate. Interestingly, the increased residence time translated into longer duration of action in vivo. It was thus possible to further enhance the efficacy of aDabi-Fab2 based on rational design, giving it the potential to serve as a back-up candidate for idarucizumab.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Effect of osmolytes on the binding of EGR1 transcription factor to DNA

Osmolytes play a key role in maintaining protein stability and mediating macromolecular interactions within the intracellular environment of the cell. Herein, we show that osmolytes such as glycerol, sucrose, and polyethylene glycol 400 (PEG400) mitigate the binding of early growth response (protein) 1 (EGR1) transcription factor to DNA in a differential manner. Thus, while physiological concentrations of glycerol only moderately reduce the binding affinity, addition of sucrose and PEG400 is concomitant with a loss in the binding affinity by an order of magnitude. This salient observation suggests that EGR1 is most likely subject to conformational equilibrium and that the osmolytes exert their effect via favorable interactions with the unliganded conformation. Consistent with this notion, our analysis reveals that while EGR1 displays rather high structural stability in complex with DNA, the unliganded conformation becomes significantly destabilized in solution. In particular, while liganded EGR1 adopts a well-defined arc-like architecture, the unliganded protein samples a comparatively large conformational space between two distinct states that periodically interconvert between an elongated rod-like shape and an arc-like conformation on a submicrosecond time scale. Consequently, the ability of osmolytes to favorably interact with the unliganded conformation so as to stabilize it could account for the negative effect of osmotic stress on EGR1-DNA interaction observed here. Taken together, our study sheds new light on the role of osmolytes in modulating a key protein-DNA interaction.

Addressing the Role of Conformational Diversity in Protein Structure Prediction

Computational modeling of tertiary structures has become of standard use to study proteins that lack experimental characterization. Unfortunately, 3D structure prediction methods and model quality assessment programs often overlook that an ensemble of conformers in equilibrium populates the native state of proteins. In this work we collected sets of publicly available protein models and the corresponding target structures experimentally solved and studied how they describe the conformational diversity of the protein. For each protein, we assessed the quality of the models against known conformers by several standard measures and identified those models ranked best. We found that model rankings are defined by both the selected target conformer and the similarity measure used. 70% of the proteins in our datasets show that different models are structurally closest to different conformers of the same protein target. We observed that model building protocols such as template-based or ab initio approaches describe in similar ways the conformational diversity of the protein, although for template-based methods this description may depend on the sequence similarity between target and template sequences. Taken together, our results support the idea that protein structure modeling could help to identify members of the native ensemble, highlight the importance of considering conformational diversity in protein 3D quality evaluations and endorse the study of the variability of the native structure for a meaningful biological analysis.

Mechanism of allosteric propagation across a β-sheet structure investigated by molecular dynamics simulations

The bacterial adhesin FimH consists of an allosterically regulated mannose-binding lectin domain and a covalently linked inhibitory pilin domain. Under normal conditions, the two domains are bound to each other, and FimH interacts weakly with mannose. However, under tensile force, the domains separate and the lectin domain undergoes conformational changes that strengthen its bond with mannose. Comparison of the crystallographic structures of the low and the high affinity state of the lectin domain reveals conformational changes mainly in the regulatory inter-domain region, the mannose binding site and a large β sheet that connects the two distally located regions. Here, molecular dynamics simulations investigated how conformational changes are propagated within and between different regions of the lectin domain. It was found that the inter-domain region moves towards the high affinity conformation as it becomes more compact and buries exposed hydrophobic surface after separation of the pilin domain. The mannose binding site was more rigid in the high affinity state, which prevented water penetration into the pocket. The large central β sheet demonstrated a soft spring-like twisting. Its twisting motion was moderately correlated to fluctuations in both the regulatory and the binding region, whereas a weak correlation was seen in a direct comparison of these two distal sites. The results suggest a so called "population shift" model whereby binding of the lectin domain to either the pilin domain or mannose locks the β sheet in a rather twisted or flat conformation, stabilizing the low or the high affinity state, respectively. Proteins 2016; 84:990-1008. © 2016 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Knowledge-Based Methods To Train and Optimize Virtual Screening Ensembles

Ensemble docking can be a successful virtual screening technique that addresses the innate conformational heterogeneity of macromolecular drug targets. Yet, lacking a method to identify a subset of conformational states that effectively segregates active and inactive small molecules, ensemble docking may result in the recommendation of a large number of false positives. Here, three knowledge-based methods that construct structural ensembles for virtual screening are presented. Each method selects ensembles by optimizing an objective function calculated using the receiver operating characteristic (ROC) curve: either the area under the ROC curve (AUC) or a ROC enrichment factor (EF). As the number of receptor conformations, N, becomes large, the methods differ in their asymptotic scaling. Given a set of small molecules with known activities and a collection of target conformations, the most resource intense method is guaranteed to find the optimal ensemble but scales as O(2^N). A recursive approximation to the optimal solution scales as O(N²), and a more severe approximation leads to a faster method that scales linearly, O(N). The techniques are generally applicable to any system, and we demonstrate their effectiveness on the androgen nuclear hormone receptor (AR), cyclin-dependent kinase 2 (CDK2), and the peroxisome proliferator-activated receptor δ (PPAR-δ) drug targets. Conformations that consisted of a crystal structure and molecular dynamics simulation cluster centroids were used to form AR and CDK2 ensembles. Multiple available crystal structures were used to form PPAR-δ ensembles. For each target, we show that the three methods perform similarly to one another on both the training and test sets.

The interaction with gold suppresses fiber-like conformations of the amyloid β (16-22) peptide

Inorganic surfaces and nanoparticles can accelerate or inhibit the fibrillation process of proteins and peptides, including the biomedically relevant amyloid β peptide. However, the microscopic mechanisms that determine such an effect are still poorly understood. By means of large-scale, state-of-the-art enhanced sampling molecular dynamics simulations, here we identify an interaction mechanism between the segments 16-22 of the amyloid β peptide, known to be fibrillogenic by itself, and the Au(111) surface in water that leads to the suppression of fiber-like conformations from the peptide conformational ensemble. Moreover, thanks to advanced simulation analysis techniques, we characterize the conformational selection vs. induced fit nature of the gold effect. Our results disclose an inhibition mechanism that is rooted in the details of the microscopic peptide-surface interaction rather than in general phenomena such as peptide sequestration from the solution.

Rigid Residue Scan Simulations Systematically Reveal Residue Entropic Roles in Protein Allostery

Intra-protein information is transmitted over distances via allosteric processes. This ubiquitous protein process allows for protein function changes due to ligand binding events. Understanding protein allostery is essential to understanding protein functions. In this study, allostery in the second PDZ domain (PDZ2) in the human PTP1E protein is examined as model system to advance a recently developed rigid residue scan method combining with configurational entropy calculation and principal component analysis. The contributions from individual residues to whole-protein dynamics and allostery were systematically assessed via rigid body simulations of both unbound and ligand-bound states of the protein. The entropic contributions of individual residues to whole-protein dynamics were evaluated based on covariance-based correlation analysis of all simulations. The changes of overall protein entropy when individual residues being held rigid support that the rigidity/flexibility equilibrium in protein structure is governed by the La Châtelier's principle of chemical equilibrium. Key residues of PDZ2 allostery were identified with good agreement with NMR studies of the same protein bound to the same peptide. On the other hand, the change of entropic contribution from each residue upon perturbation revealed intrinsic differences among all the residues. The quasi-harmonic and principal component analyses of simulations without rigid residue perturbation showed a coherent allosteric mode from unbound and bound states, respectively. The projection of simulations with rigid residue perturbation onto coherent allosteric modes demonstrated the intrinsic shifting of ensemble distributions supporting the population-shift theory of protein allostery. Overall, the study presented here provides a robust and systematic approach to estimate the contribution of individual residue internal motion to overall protein dynamics and allostery.

Probing the Impact of Local Structural Dynamics of Conformational Epitopes on Antibody Recognition

Antibody-antigen interactions are governed by recognition of specific residues and structural complementarity between the antigen epitope and antibody paratope. While X-ray crystallography has provided detailed insights into static conformations of antibody-antigen complexes, factors such as conformational flexibility and dynamics, which are not readily apparent in the structures, can also have an impact on the binding event. Here we investigate the contribution of dynamics in the HIV-1 gp120 glycoprotein to antibody recognition of conserved conformational epitopes, including the CD4- and coreceptor-binding sites, and an inner domain site that is targeted by ADCC-active antibodies. Hydrogen/deuterium-exchange mass spectrometry (HDX-MS) was used to measure local structural dynamics across a panel of variable loop truncation mutants of HIV-1 gp120, including full-length gp120, ΔV3, ΔV1/V2, and extended core, which includes ΔV1/V2 and V3 loop truncations. CD4-bound full-length gp120 was also examined as a reference state. HDX-MS revealed a clear trend toward an increased level of order of the conserved subunit core resulting from loop truncation. Combined with biolayer interferometry and enzyme-linked immunosorbent assay measurements of antibody-antigen binding, we demonstrate that an increased level of ordering of the subunit core was associated with better recognition by an array of antibodies targeting complex conformational epitopes. These results provide detailed insight into the influence of structural dynamics on antibody-antigen interactions and suggest the importance of characterizing the structural stability of vaccine candidates to improve antibody recognition of complex epitopes.

Emerging Computational Methods for the Rational Discovery of Allosteric Drugs

Allosteric drug development holds promise for delivering medicines that are more selective and less toxic than those that target orthosteric sites. To date, the discovery of allosteric binding sites and lead compounds has been mostly serendipitous, achieved through high-throughput screening. Over the past decade, structural data has become more readily available for larger protein systems and more membrane protein classes (e.g., GPCRs and ion channels), which are common allosteric drug targets. In parallel, improved simulation methods now provide better atomistic understanding of the protein dynamics and cooperative motions that are critical to allosteric mechanisms. As a result of these advances, the field of predictive allosteric drug development is now on the cusp of a new era of rational structure-based computational methods. Here, we review algorithms that predict allosteric sites based on sequence data and molecular dynamics simulations, describe tools that assess the druggability of these pockets, and discuss how Markov state models and topology analyses provide insight into the relationship between protein dynamics and allosteric drug binding. In each section, we first provide an overview of the various method classes before describing relevant algorithms and software packages.

Long range recognition and selection in IDPs: the interactions of the C-terminus of p53

The C-terminal domain of p53 is an extensively studied IDP, interacting with different partners through multiple distinct conformations. To explore the interplay between preformed structural elements and intrinsic fluctuations in its folding and binding we combine extensive atomistic equilibrium and non-equilibrium simulations. We find that the free peptide segment rapidly interconverts between ordered and disordered states with significant populations of the conformations that are seen in the complexed states. The underlying global folding-binding landscape points to a synergistic mechanism in which recognition is dictated via long range electrostatic recognition which results in the formation of reactive structures as far away as 10 Å, and binding proceeds with the steering of selected conformations followed by induced folding at the target surface or within a close range.

Thermodynamic Basis of Selectivity in the Interactions of Tissue Inhibitors of Metalloproteinases N-domains with Matrix Metalloproteinases-1, -3, and -14

The four tissue inhibitors of metalloproteinases (TIMPs) are potent inhibitors of the many matrixins (MMPs), except that TIMP1 weakly inhibits some MMPs, including MMP14. The broad-spectrum inhibition of MMPs by TIMPs and their N-domains (NTIMPs) is consistent with the previous isothermal titration calorimetric finding that their interactions are entropy-driven but differ in contributions from solvent and conformational entropy (ΔSsolv, ΔSconf), estimated using heat capacity changes (ΔCp). Selective engineered NTIMPs have potential applications for treating MMP-related diseases, including cancer and cardiomyopathy. Here we report isothermal titration calorimetric studies of the effects of selectivity-modifying mutations in NTIMP1 and NTIMP2 on the thermodynamics of their interactions with MMP1, MMP3, and MMP14. The weak inhibition of MMP14 by NTIMP1 reflects a large conformational entropy penalty for binding. The T98L mutation, peripheral to the NTIMP1 reactive site, enhances binding by increasing ΔSsolv but also reduces ΔSconf However, the same mutation increases NTIMP1 binding to MMP3 in an interaction that has an unusual positive ΔCp This indicates a decrease in solvent entropy compensated by increased conformational entropy, possibly reflecting interactions involving alternative conformers. The NTIMP2 mutant, S2D/S4A is a selective MMP1 inhibitor through electrostatic effects of a unique MMP-1 arginine. Asp-2 increases reactive site polarity, reducing ΔCp, but increases conformational entropy to maintain strong binding to MMP1. There is a strong negative correlation between ΔSsolv and ΔSconf for all characterized interactions, but the data for each MMP have characteristic ranges, reflecting intrinsic differences in the structures and dynamics of their free and inhibitor-bound forms.

Salt bridge as a gatekeeper against partial unfolding

Salt bridges are frequently observed in protein structures. Because the energetic contribution of salt bridges is strongly dependent on the environmental context, salt bridges are believed to contribute to the structural specificity rather than the stability. To test the role of salt bridges in enhancing structural specificity, we investigated the contribution of a salt bridge to the energetics of native-state partial unfolding in a cysteine-free version of Escherichia coli ribonuclease H (RNase H*). Thermolysin cleaves a protruding loop of RNase H(*) through transient partial unfolding under native conditions. Lys86 and Asp108 in RNase H(*) form a partially buried salt bridge that tethers the protruding loop. Investigation of the global stability of K86Q/D108N RNase H(*) showed that the salt bridge does not significantly contribute to the global stability. However, K86Q/D108N RNase H(*) is greatly more susceptible to proteolysis by thermolysin than wild-type RNase H(*) is. The free energy for partial unfolding determined by native-state proteolysis indicates that the salt bridge significantly increases the energy for partial unfolding by destabilizing the partially unfolded form. Double mutant cycles with single and double mutations of the salt bridge suggest that the partially unfolded form is destabilized due to a significant decrease in the interaction energy between Lys86 and Asp108 upon partial unfolding. This study demonstrates that, even in the case that a salt bridge does not contribute to the global stability, the salt bridge may function as a gatekeeper against partial unfolding that disturbs the optimal geometry of the salt bridge.

Ras Conformational Ensembles, Allostery, and Signaling

Ras proteins are classical members of small GTPases that function as molecular switches by alternating between inactive GDP-bound and active GTP-bound states. Ras activation is regulated by guanine nucleotide exchange factors that catalyze the exchange of GDP by GTP, and inactivation is terminated by GTPase-activating proteins that accelerate the intrinsic GTP hydrolysis rate by orders of magnitude. In this review, we focus on data that have accumulated over the past few years pertaining to the conformational ensembles and the allosteric regulation of Ras proteins and their interpretation from our conformational landscape standpoint. The Ras ensemble embodies all states, including the ligand-bound conformations, the activated (or inactivated) allosteric modulated states, post-translationally modified states, mutational states, transition states, and nonfunctional states serving as a reservoir for emerging functions. The ensemble is shifted by distinct mutational events, cofactors, post-translational modifications, and different membrane compositions. A better understanding of Ras biology can contribute to therapeutic strategies.

Insights into Protein-Ligand Interactions: Mechanisms, Models, and Methods

Molecular recognition, which is the process of biological macromolecules interacting with each other or various small molecules with a high specificity and affinity to form a specific complex, constitutes the basis of all processes in living organisms. Proteins, an important class of biological macromolecules, realize their functions through binding to themselves or other molecules. A detailed understanding of the protein–ligand interactions is therefore central to understanding biology at the molecular level. Moreover, knowledge of the mechanisms responsible for the protein-ligand recognition and binding will also facilitate the discovery, design, and development of drugs. In the present review, first, the physicochemical mechanisms underlying protein–ligand binding, including the binding kinetics, thermodynamic concepts and relationships, and binding driving forces, are introduced and rationalized. Next, three currently existing protein-ligand binding models—the “lock-and-key”, “induced fit”, and “conformational selection”—are described and their underlying thermodynamic mechanisms are discussed. Finally, the methods available for investigating protein–ligand binding affinity, including experimental and theoretical/computational approaches, are introduced, and their advantages, disadvantages, and challenges are discussed.

Interplay between the hydrophobic effect and dipole interactions in peptide aggregation at interfaces

Amphipathic octapeptide aggregation at hydrophobic–hydrophilic interfaces is largely driven by backbone dipole interactions in peptide aggregation at interfaces.

Using THz Spectroscopy, Evolutionary Network Analysis Methods, and MD Simulation to Map the Evolution of Allosteric Communication Pathways in c-Type Lysozymes

It is now widely accepted that protein function is intimately tied with the navigation of energy landscapes. In this framework, a protein sequence is not described by a distinct structure but rather by an ensemble of conformations. And it is through this ensemble that evolution is able to modify a protein's function by altering its landscape. Hence, the evolution of protein functions involves selective pressures that adjust the sampling of the conformational states. In this work, we focus on elucidating the evolutionary pathway that shaped the function of individual proteins that make-up the mammalian c-type lysozyme subfamily. Using both experimental and computational methods, we map out specific intermolecular interactions that direct the sampling of conformational states and accordingly, also underlie shifts in the landscape that are directly connected with the formation of novel protein functions. By contrasting three representative proteins in the family we identify molecular mechanisms that are associated with the selectivity of enhanced antimicrobial properties and consequently, divergent protein function. Namely, we link the extent of localized fluctuations involving the loop separating helices A and B with shifts in the equilibrium of the ensemble of conformational states that mediate interdomain coupling and concurrently moderate substrate binding affinity. This work reveals unique insights into the molecular level mechanisms that promote the progression of interactions that connect the immune response to infection with the nutritional properties of lactation, while also providing a deeper understanding about how evolving energy landscapes may define present-day protein function.

Structural studies on Mycobacterium tuberculosis RecA: molecular plasticity and interspecies variability

Structures of crystals of Mycobacterium tuberculosis RecA, grown and analysed under different conditions, provide insights into hitherto underappreciated details of molecular structure and plasticity. In particular, they yield information on the invariant and variable features of the geometry of the P-loop, whose binding to ATP is central for all the biochemical activities of RecA. The strengths of interaction of the ligands with the P-loop reveal significant differences. This in turn affects the magnitude of the motion of the 'switch' residue, Gln195 in M. tuberculosis RecA, which triggers the transmission of ATP-mediated allosteric information to the DNA binding region. M. tuberculosis RecA is substantially rigid compared with its counterparts from M. smegmatis and E. coli, which exhibit concerted internal molecular mobility. The interspecies variability in the plasticity of the two mycobacterial proteins is particularly surprising as they have similar sequence and 3D structure. Details of the interactions of ligands with the protein, characterized in the structures reported here, could be useful for design of inhibitors against M. tuberculosis RecA.

cNMA: a framework of encounter complex-based normal mode analysis to model conformational changes in protein interactions

MOTIVATION

It remains both a fundamental and practical challenge to understand and anticipate motions and conformational changes of proteins during their associations. Conventional normal mode analysis (NMA) based on anisotropic network model (ANM) addresses the challenge by generating normal modes reflecting intrinsic flexibility of proteins, which follows a conformational selection model for protein-protein interactions. But earlier studies have also found cases where conformational selection alone could not adequately explain conformational changes and other models have been proposed. Moreover, there is a pressing demand of constructing a much reduced but still relevant subset of protein conformational space to improve computational efficiency and accuracy in protein docking, especially for the difficult cases with significant conformational changes.

METHOD AND RESULTS

With both conformational selection and induced fit models considered, we extend ANM to include concurrent but differentiated intra- and inter-molecular interactions and develop an encounter complex-based NMA (cNMA) framework. Theoretical analysis and empirical results over a large data set of significant conformational changes indicate that cNMA is capable of generating conformational vectors considerably better at approximating conformational changes with contributions from both intrinsic flexibility and inter-molecular interactions than conventional NMA only considering intrinsic flexibility does. The empirical results also indicate that a straightforward application of conventional NMA to an encounter complex often does not improve upon NMA for an individual protein under study and intra- and inter-molecular interactions need to be differentiated properly. Moreover, in addition to induced motions of a protein under study, the induced motions of its binding partner and the coupling between the two sets of protein motions present in a near-native encounter complex lead to the improved performance. A study to isolate and assess the sole contribution of intermolecular interactions toward improvements against conventional NMA further validates the additional benefit from induced-fit effects. Taken together, these results provide new insights into molecular mechanisms underlying protein interactions and new tools for dimensionality reduction for flexible protein docking.

AVAILABILITY AND IMPLEMENTATION

Source codes are available upon request.

Molecular dynamics simulations elucidate the mode of protein recognition by Skp1 and the F-box domain in the SCF complex

Polyubiquitination of the target protein by a ubiquitin transferring machinery is key to various cellular processes. E3 ligase Skp1-Cul1-F-box (SCF) is one such complex which plays crucial role in substrate recognition and transfer of the ubiquitin molecule. Previous computational studies have focused on S-phase kinase-associated protein 2 (Skp2), cullin, and RING-finger proteins of this complex, but the roles of the adapter protein Skp1 and F-box domain of Skp2 have not been determined. Using sub-microsecond molecular dynamics simulations of full-length Skp1, unbound Skp2, Skp2-Cks1 (Cks1: Cyclin-dependent kinases regulatory subunit 1), Skp1-Skp2, and Skp1-Skp2-Cks1 complexes, we have elucidated the function of Skp1 and the F-box domain of Skp2. We found that the L16 loop of Skp1, which was deleted in previous X-ray crystallography studies, can offer additional stability to the ternary complex via its interactions with the C-terminal tail of Skp2. Moreover, Skp1 helices H6, H7, and H8 display vivid conformational flexibility when not bound to Skp2, suggesting that these helices can recognize and lock the F-box proteins. Furthermore, we observed that the F-box domain could rotate (5°-129°), and that the binding partner determined the degree of conformational flexibility. Finally, Skp1 and Skp2 were found to execute a domain motion in Skp1-Skp2 and Skp1-Skp2-Cks1 complexes that could decrease the distance between ubiquitination site of the substrate and the ubiquitin molecule by 3 nm. Thus, we propose that both the F-box domain of Skp2 and Skp1-Skp2 domain motions displaying preferential conformational control can together facilitate polyubiquitination of a wide variety of substrates.

Towards a Molecular Understanding of the Link between Imatinib Resistance and Kinase Conformational Dynamics

Due to its inhibition of the Abl kinase domain in the BCR-ABL fusion protein, imatinib is strikingly effective in the initial stage of chronic myeloid leukemia with more than 90% of the patients showing complete remission. However, as in the case of most targeted anti-cancer therapies, the emergence of drug resistance is a serious concern. Several drug-resistant mutations affecting the catalytic domain of Abl and other tyrosine kinases are now known. But, despite their importance and the adverse effect that they have on the prognosis of the cancer patients harboring them, the molecular mechanism of these mutations is still debated. Here by using long molecular dynamics simulations and large-scale free energy calculations complemented by in vitro mutagenesis and microcalorimetry experiments, we model the effect of several widespread drug-resistant mutations of Abl. By comparing the conformational free energy landscape of the mutants with those of the wild-type tyrosine kinases we clarify their mode of action. It involves significant and complex changes in the inactive-to-active dynamics and entropy/enthalpy balance of two functional elements: the activation-loop and the conserved DFG motif. What is more the T315I gatekeeper mutant has a significant impact on the binding mechanism itself and on the binding kinetics.

Conformational selection and induced fit for RNA polymerase and RNA/DNA hybrid backtracked recognition

RNA polymerase catalyzes transcription with a high fidelity. If DNA/RNA mismatch or DNA damage occurs downstream, a backtracked RNA polymerase can proofread this situation. However, the backtracked mechanism is still poorly understood. Here we have performed multiple explicit-solvent molecular dynamics (MD) simulations on bound and apo DNA/RNA hybrid to study backtracked recognition. MD simulations at room temperature suggest that specific electrostatic interactions play key roles in the backtracked recognition between the polymerase and DNA/RNA hybrid. Kinetics analysis at high temperature shows that bound and apo DNA/RNA hybrid unfold via a two-state process. Both kinetics and free energy landscape analyses indicate that bound DNA/RNA hybrid folds in the order of DNA/RNA contracting, the tertiary folding and polymerase binding. The predicted Φ-values suggest that C7, G9, dC12, dC15, and dT16 are key bases for the backtracked recognition of DNA/RNA hybrid. The average RMSD values between the bound structures and the corresponding apo ones and Kolmogorov-Smirnov (KS) P-test analyses indicate that the recognition between DNA/RNA hybrid and polymerase might follow an induced fit mechanism for DNA/RNA hybrid and conformation selection for polymerase. Furthermore, this method could be used to relative studies of specific recognition between nucleic acid and protein.

R248Q mutation--Beyond p53-DNA binding

R248 in the DNA binding domain (DBD) of p53 interacts directly with the minor groove of DNA. Earlier nuclear magnetic resonance (NMR) studies indicated that the R248Q mutation resulted in conformation changes in parts of DBD far from the mutation site. However, how information propagates from the mutation site to the rest of the DBD is still not well understood. We performed a series of all-atom molecular dynamics (MD) simulations to dissect sterics and charge effects of R248 on p53-DBD conformation: (i) wild-type p53 DBD; (ii) p53 DBD with an electrically neutral arginine side-chain; (iii) p53 DBD with R248A; (iv) p53 DBD with R248W; and (v) p53 DBD with R248Q. Our results agree well with experimental observations of global conformational changes induced by the R248Q mutation. Our simulations suggest that both charge- and sterics are important in the dynamics of the loop (L3) where the mutation resides. We show that helix 2 (H2) dynamics is altered as a result of a change in the hydrogen bonding partner of D281. In turn, neighboring L1 dynamics is altered: in mutants, L1 predominantly adopts the recessed conformation and is unable to interact with the major groove of DNA. We focused our attention the R248Q mutant that is commonly found in a wide range of cancer and observed changes at the zinc-binding pocket that might account for the dominant negative effects of R248Q. Furthermore, in our simulations, the S6/S7 turn was more frequently solvent exposed in R248Q, suggesting that there is a greater tendency of R248Q to partially unfold and possibly lead to an increased aggregation propensity. Finally, based on the observations made in our simulations, we propose strategies for the rescue of R248Q mutants.

Dynamics Govern Specificity of a Protein-Protein Interface: Substrate Recognition by Thrombin

Biomolecular recognition is crucial in cellular signal transduction. Signaling is mediated through molecular interactions at protein-protein interfaces. Still, specificity and promiscuity of protein-protein interfaces cannot be explained using simplistic static binding models. Our study rationalizes specificity of the prototypic protein-protein interface between thrombin and its peptide substrates relying solely on binding site dynamics derived from molecular dynamics simulations. We find conformational selection and thus dynamic contributions to be a key player in biomolecular recognition. Arising entropic contributions complement chemical intuition primarily reflecting enthalpic interaction patterns. The paradigm “dynamics govern specificity” might provide direct guidance for the identification of specific anchor points in biomolecular recognition processes and structure-based drug design.

Directed evolution of the Escherichia coli cAMP receptor protein at the cAMP pocket

The Escherichia coli cAMP receptor protein (CRP) requires cAMP binding to undergo a conformational change for DNA binding and transcriptional regulation. Two CRP residues, Thr(127) and Ser(128), are known to play important roles in cAMP binding through hydrogen bonding and in the cAMP-induced conformational change, but the connection between the two is not completely clear. Here, we simultaneously randomized the codons for these two residues and selected CRP mutants displaying high CRP activity in a cAMP-producing E. coli. Many different CRP mutants satisfied the screening condition for high CRP activity, including those that cannot form any hydrogen bonds with the incoming cAMP at the two positions. In vitro DNA-binding analysis confirmed that these selected CRP mutants indeed display high CRP activity in response to cAMP. These results indicate that the hydrogen bonding ability of the Thr(127) and Ser(128) residues is not critical for the cAMP-induced CRP activation. However, the hydrogen bonding ability of Thr(127) and Ser(128) was found to be important in attaining high cAMP affinity. Computational analysis revealed that most natural cAMP-sensing CRP homologs have Thr/Ser, Thr/Thr, or Thr/Asn at positions 127 and 128. All of these pairs are excellent hydrogen bonding partners and they do not elevate CRP activity in the absence of cAMP. Taken together, our analyses suggest that CRP evolved to have hydrogen bonding residues at the cAMP pocket residues 127 and 128 for performing dual functions: preserving high cAMP affinity and keeping CRP inactive in the absence of cAMP.

Protein Flexibility in Drug Discovery: From Theory to Computation

Nowadays it is widely accepted that the mechanisms of biomolecular recognition are strongly coupled to the intrinsic dynamic of proteins. In past years, this evidence has prompted the development of theoretical models of recognition able to describe ligand binding assisted by protein conformational changes. On a different perspective, the need to take into account protein flexibility in structure-based drug discovery has stimulated the development of several and extremely diversified computational methods. Herein, on the basis of a parallel between the major recognition models and the simulation strategies used to account for protein flexibility in ligand binding, we sort out and describe the most innovative and promising implementations for structure-based drug discovery.

Standardization and quality control in quantifying non-enzymatic oxidative protein modifications in relation to ageing and disease: Why is it important and why is it hard?

Post-translational modifications (PTM) of proteins determine the activity, stability, specificity, transportability and lifespan of a protein. Some PTM are highly specific and regulated involving various enzymatic pathways, but there are other non-enzymatic PTM (nePTM), which occur stochastically, depend on the ternary structure of proteins and can be damaging. It is often observed that inactive and abnormal proteins accumulate in old cells and tissues. The nature, site and extent of nePTM give rise to a population of that specific protein with alterations in structure and function ranging from being fully active to totally inactive molecules. Determination of the type and the amount (abundance) of nePTM is essential for establishing connection between specific protein structure and specific biological role. This article summarizes analytical demands for reliable quantification of nePTM, including requirements for the assay performance, standardization and quality control, and points to the difficulties, uncertainties and un-resolved issues.

Identifying sequential substrate binding at the single-molecule level by enzyme mechanical stabilization

Enzyme-substrate binding is a dynamic process intimately coupled to protein structural changes, which in turn changes the unfolding energy landscape. By the use of single-molecule force spectroscopy (SMFS), we characterize the open-to-closed conformational transition experienced by the hyperthermophilic adenine diphosphate (ADP)-dependent glucokinase from Thermococcus litoralis triggered by the sequential binding of substrates. In the absence of substrates, the mechanical unfolding of TlGK shows an intermediate 1, which is stabilized in the presence of Mg·ADP(-), the first substrate to bind to the enzyme. However, in the presence of this substrate, an additional unfolding event is observed, intermediate 1*. Finally, in the presence of both substrates, the unfolding force of intermediates 1 and 1* increases as a consequence of the domain closure. These results show that SMFS can be used as a powerful experimental tool to investigate binding mechanisms of different enzymes with more than one ligand, expanding the repertoire of protocols traditionally used in enzymology.

Global ITC fitting methods in studies of protein allostery

Allostery is a nearly ubiquitous feature of biological systems in which ligand binding or covalent modification at one site alters the activities of distant sites in a macromolecule or macromolecular complex. The molecular mechanisms underlying this phenomenon have been studied for decades. Nevertheless there are many aspects that remain poorly understood. ITC yields detailed information on the thermodynamics of biomacromolecular interactions and their coupling to additional equilibria, therefore in principle it is a powerful tool for better understanding how allostery is achieved. A particularly powerful approach involves simultaneously fitting multiple ITC data sets together with those of complementary techniques, especially nuclear magnetic resonance and circular dichroism spectroscopies. In this review, we describe several group-fitting methods for discriminating between different binding models and for improving the accuracy of thermodynamic parameters extracted from variable-temperature ITC data. The techniques were applied to the antibiotic resistance-causing enzyme aminoglycoside-6'-acetyltransferase Ii, uncovering the existence of competition between opposing mechanisms and ligand-dependent switching of the underlying mechanism. These novel observations underline the potential of combining ITC and spectroscopic techniques to study allostery.

Relaxation dispersion NMR spectroscopy for the study of protein allostery

Allosteric transmission of information between distant sites in biological macromolecules often involves collective transitions between active and inactive conformations. Nuclear magnetic resonance (NMR) spectroscopy can yield detailed information on these dynamics. In particular, relaxation dispersion techniques provide structural, dynamic, and mechanistic information on conformational transitions occurring on the millisecond to microsecond timescales. In this review, we provide an overview of the theory and analysis of Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiments and briefly describe their application to the study of allosteric dynamics in the homeodomain from the PBX transcription factor (PBX-HD). CPMG NMR data show that local folding (helix/coil) transitions in one part of PBX-HD help to communicate information between two distant binding sites. Furthermore, the combination of CPMG and other spin relaxation data show that this region can also undergo local misfolding, reminiscent of conformational ensemble models of allostery.

Role of indirect readout mechanism in TATA box binding protein-DNA interaction

Gene expression generally initiates from recognition of TATA-box binding protein (TBP) to the minor groove of DNA of TATA box sequence where the DNA structure is significantly different from B-DNA. We have carried out molecular dynamics simulation studies of TBP-DNA system to understand how the DNA structure alters for efficient binding. We observed rigid nature of the protein while the DNA of TATA box sequence has an inherent flexibility in terms of bending and minor groove widening. The bending analysis of the free DNA and the TBP bound DNA systems indicate presence of some similar structures. Principal coordinate ordination analysis also indicates some structural features of the protein bound and free DNA are similar. Thus we suggest that the DNA of TATA box sequence regularly oscillates between several alternate structures and the one suitable for TBP binding is induced further by the protein for proper complex formation.