Trapping moving targets with small molecules

Single-Molecule Imaging of Integral Membrane Protein Dynamics and Function

Integral membrane proteins (IMPs) play central roles in cellular physiology and represent the majority of known drug targets. Single-molecule fluorescence and fluorescence resonance energy transfer (FRET) methods have recently emerged as valuable tools for investigating structure-function relationships in IMPs. This review focuses on the practical foundations required for examining polytopic IMP function using single-molecule FRET (smFRET) and provides an overview of the technical and conceptual frameworks emerging from this area of investigation. In this context, we highlight the utility of smFRET methods to reveal transient conformational states critical to IMP function and the use of smFRET data to guide structural and drug mechanism-of-action investigations. We also identify frontiers where progress is likely to be paramount to advancing the field.

Insight into de-regulation of amino acid feedback inhibition: a focus on structure analysis method

Regulation of amino acid's biosynthetic pathway is of significant importance to maintain homeostasis and cell functions. Amino acids regulate their biosynthetic pathway by end-product feedback inhibition of enzymes catalyzing committed steps of a pathway. Discovery of new feedback resistant enzyme variants to enhance industrial production of amino acids is a key objective in industrial biotechnology. Deregulation of feedback inhibition has been achieved for various enzymes using in vitro and in silico mutagenesis techniques. As enzyme's function, its substrate binding capacity, catalysis activity, regulation and stability are dependent on its structural characteristics, here, we provide detailed structural analysis of all feedback sensitive enzyme targets in amino acid biosynthetic pathways. Current review summarizes information regarding structural characteristics of various enzyme targets and effect of mutations on their structures and functions especially in terms of deregulation of feedback inhibition. Furthermore, applicability of various experimental as well as computational mutagenesis techniques to accomplish feedback resistance has also been discussed in detail to have an insight into various aspects of research work reported in this particular field of study.

Conformational Changes in Surface-Immobilized Proteins Measured Using Combined Atomic Force and Fluorescence Microscopy

Biological organisms rely on proteins to perform the majority of their functions. Most protein functions are based on their physical motions (conformational changes), which can be described as transitions between different conformational states in a multidimensional free-energy landscape. A comprehensive understanding of this free-energy landscape is therefore of paramount importance for understanding the biological functions of proteins. Protein dynamics includes both equilibrium and nonequilibrium motions, which typically exhibit a wide range of characteristic length and time scales. The relative probabilities of various conformational states in the energy landscape, the energy barriers between them, their dependence on external parameters such as force and temperature, and their connection to the protein function remain largely unknown for most proteins. In this paper, we present a multimolecule approach in which the proteins are immobilized at well-defined locations on Au substrates using an atomic force microscope (AFM)-based patterning method called nanografting. This method enables precise control over the protein location and orientation on the substrate, as well as the creation of biologically active protein ensembles that self-assemble into well-defined nanoscale regions (protein patches) on the gold substrate. We performed AFM-force compression and fluorescence experiments on these protein patches and measured the fundamental dynamical parameters such as protein stiffness, elastic modulus, and transition energies between distinct conformational states. Our results provide new insights into the processes that govern protein dynamics and its connection to protein function.

Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations

Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

Dual-Inhibition of Human N-Myristoyltransferase Subtypes Halts Common Cold Pathogenesis: Atomistic Perspectives from the Case of IMP-1088

The pharmacological inhibition of human N-myristoyltransferase (HsNMT) has emerged as an efficient strategy to completely prevent the replication process of rhinoviruses, a potential treatment for the common cold. This was corroborated by the recent discovery of compound IMP-1088, a novel inhibitor that demonstrated a dual-inhibitory activity against the two HsNMT subtypes 1 and 2 without inducing cytotoxicity. However, the molecular and structural basis for the dual-inhibitory potential of IMP-1088 has not been investigated. As such, we employ molecular modelling techniques to resolve the structural mechanisms that account for the dual-inhibitory prowess of IMP-1088. Sequence and nanosecond-based analyses identified Tyr296, Phe190, Tyr420, Leu453, Gln496, Val181, Leu474, Glu182, and Asn246 as residues common within the binding pockets of both HsNMT1 and HsNMT2 subtypes whose consistent interactions with IMP-1088 underpin the basis for its dual inhibitory potency. Nano-second-based assessment of interaction dynamics revealed that Tyr296 consistently elicited high-affinity π-π stacked interaction with IMP-1088, thus further highlighting its cruciality corroborating previous report. An exploration of resulting structural changes upon IMP-1088 binding further revealed a characteristic impeding of residue fluctuations, structural compactness, and a consequential burial of crucial hydrophobic residues, features required for HsNMT1/2 functionality. Findings present essential structural perspectives that augment previous experimental efforts and could also advance drug development for treating respiratory tract infections, especially those mediated by rhinoviruses.

Synthesis and spectral characterization of the first fluorescein-tagged iron(ii) clathrochelates, their supramolecular interactions with globular proteins, and cellular uptake

A fluorescein-tagged iron(ii) cage complex was obtained in a moderate total yield using a two-step synthetic procedure starting from its propargylamine-containing clathrochelate precursor. An 11-fold decrease in fluorescence quantum yield is observed in passing from the given fluorescein-based dye to its clathrochelate derivative. An excitation energy transfer from the terminal fluorescent group of the macrobicyclic molecule to its quasiaromatic highly π-conjugated clathrochelate framework can explain this effect. The kinetics of the hydrolysis of the acetyl groups of acetylated fluorescein azide and its clathrochelate derivative in the presence of one equivalent of BSA evidenced no strong supramolecular host-guest interactions between BSA and the tested compounds. Study of a chemical stability of the deacetylated iron(ii) clathrochelate suggested the formation of a supramolecular 1 : 1 BSA-clathrochelate assembly. Moreover, an addition of BSA or HSA to its solution caused the appearance of strong clathrochelate-based ICD outputs. The fluorescence emission anisotropy studies also evidenced the supramolecular binding of the fluorescein-tagged iron(ii) clathrochelate to the BSA macromolecule, leading to a high increase in this type of anisotropy. Subcellular uptake of the fluorescein-tagged molecules was visualized using fluorescence microscopy and showed its distribution to be mainly in the cytosol without entering the nucleus or accumulating in any other organelle. An X-rayed crystal of the above propargylamide macrobicyclic precursor with a reactive terminal C[triple bond, length as m-dash]C bond contains the clathrochelate molecules of two types, A and B. The encapsulated iron(ii) ion in these molecules is situated in the center of its FeN6-coordination polyhedron, the geometry of which is intermediate between a trigonal prism (TP) and a trigonal antiprism (TAP). The Fe-N distances vary from 1.8754(6) to 1.9286(4) Å and the heights h of their distorted TP-TAP polyhedra are very similar (2.30 and 2.31 Å); their values of φ are equal to 25.3 and 26.6°. In this crystal, the molecules of types A and B participate in different types of hydrogen bonding, giving H-bonded clathrochelate tetramers through their carboxylic and amide groups, respectively; these tetramers are connected to H-bonded chains.

Targeting protein self-association in drug design

Protein self-association is a universal phenomenon essential for stability and molecular recognition. Disrupting constitutive homomers constitutes an original and emerging strategy in drug design. Inhibition of homomeric proteins can be achieved through direct complex disruption, subunit intercalation, or by promoting inactive oligomeric states. Targeting self-interaction grants several advantages over active site inhibition because of the stimulation of protein degradation, the enhancement of selectivity, substoichiometric inhibition, and by-pass of compensatory mechanisms. This new landscape in protein inhibition is driven by the development of biophysical and biochemical tools suited for the study of homomeric proteins, such as differential scanning fluorimetry (DSF), native mass spectrometry (MS), Förster resonance energy transfer (FRET) spectroscopy, 2D nuclear magnetic resonance (NMR), and X-ray crystallography. In this review, we discuss the different aspects of this new paradigm in drug design.

Modeling of Hidden Structures Using Sparse Chemical Shift Data from NMR Relaxation Dispersion

NMR relaxation dispersion measurements report on conformational changes occurring on the μs-ms timescale. Chemical shift information derived from relaxation dispersion can be used to generate structural models of weakly populated alternative conformational states. Current methods to obtain such models rely on determining the signs of chemical shift changes between the conformational states, which are difficult to obtain in many situations. Here, we use a "sample and select" method to generate relevant structural models of alternative conformations of the C-terminal-associated region of Escherichia coli dihydrofolate reductase (DHFR), using only unsigned chemical shift changes for backbone amides and carbonyls (¹H, ¹⁵N, and ¹³C'). We find that CS-Rosetta sampling with unsigned chemical shift changes generates a diversity of structures that are sufficient to characterize a minor conformational state of the C-terminal region of DHFR. The excited state differs from the ground state by a change in secondary structure, consistent with previous predictions from chemical shift hypersurfaces and validated by the x-ray structure of a partially humanized mutant of E. coli DHFR (N23PP/G51PEKN). The results demonstrate that the combination of fragment modeling with sparse chemical shift data can determine the structure of an alternative conformation of DHFR sampled on the μs-ms timescale. Such methods will be useful for characterizing alternative states, which can potentially be used for in silico drug screening, as well as contributing to understanding the role of minor states in biology and molecular evolution.

Docking-based identification of small-molecule binding sites at protein-protein interfaces

Protein-protein interactions play an essential role in many biological processes, and their perturbation is a major cause of disease. The use of small molecules to modulate them is attracting increased attention, but protein interfaces generally do not have clear cavities for binding small compounds. A proposed strategy is to target interface hot-spot residues, but their identification through computational approaches usually require the complex structure, which is not often available. In this context, pyDock energy-based docking and scoring can predict hot-spots on the unbound proteins, thus not requiring the complex structure. Here, we have devised a new strategy to detect protein–protein inhibitor binding sites, based on the integration of molecular dynamics for the generation of transient cavities, and docking-based interface hot-spot prediction for the selection of the suitable cavities. This integrative approach has been validated on a test set formed by protein–protein complexes with known inhibitors for which complete structural data of unbound molecules and complexes is available. The results show that local conformational sampling with short molecular dynamics can generate transient cavities similar to the known inhibitor binding sites, and that docking simulations can identify the best cavities with similar predictive accuracy as when knowing the real interface. In a few cases, these predicted pockets are shown to be suitable for protein–ligand docking. The proposed strategy will be useful for many protein–protein complexes for which there is no available structure, as long as the the unbound proteins do not deviate dramatically from the bound conformations.

Protein Dynamics in F-like Bacterial Conjugation

Efficient in silico development of novel antibiotics requires high-resolution, dynamic models of drug targets. As conjugation is considered the prominent contributor to the spread of antibiotic resistance genes, targeted drug design to disrupt vital components of conjugative systems has been proposed to lessen the proliferation of bacterial antibiotic resistance. Advancements in structural imaging techniques of large macromolecular complexes has accelerated the discovery of novel protein-protein interactions in bacterial type IV secretion systems (T4SS). The known structural information regarding the F-like T4SS components and complexes has been summarized in the following review, revealing a complex network of protein-protein interactions involving domains with varying degrees of disorder. Structural predictions were performed to provide insight on the dynamicity of proteins within the F plasmid conjugative system that lack structural information.

Critical Review: digital resolution biomolecular sensing for diagnostics and life science research

One of the frontiers in the field of biosensors is the ability to quantify specific target molecules with enough precision to count individual units in a test sample, and to observe the characteristics of individual biomolecular interactions. Technologies that enable observation of molecules with "digital precision" have applications for in vitro diagnostics with ultra-sensitive limits of detection, characterization of biomolecular binding kinetics with a greater degree of precision, and gaining deeper insights into biological processes through quantification of molecules in complex specimens that would otherwise be unobservable. In this review, we seek to capture the current state-of-the-art in the field of digital resolution biosensing. We describe the capabilities of commercially available technology platforms, as well as capabilities that have been described in published literature. We highlight approaches that utilize enzymatic amplification, nanoparticle tags, chemical tags, as well as label-free biosensing methods.

Machine Learning Force Fields and Coarse-Grained Variables in Molecular Dynamics: Application to Materials and Biological Systems

Machine learning encompasses tools and algorithms that are now becoming popular in almost all scientific and technological fields. This is true for molecular dynamics as well, where machine learning offers promises of extracting valuable information from the enormous amounts of data generated by simulation of complex systems. We provide here a review of our current understanding of goals, benefits, and limitations of machine learning techniques for computational studies on atomistic systems, focusing on the construction of empirical force fields from ab initio databases and the determination of reaction coordinates for free energy computation and enhanced sampling.

The Regulatory Role of Key Metabolites in the Control of Cell Signaling

Robust biological systems are able to adapt to internal and environmental perturbations. This is ensured by a thick crosstalk between metabolism and signal transduction pathways, through which cell cycle progression, cell metabolism and growth are coordinated. Although several reports describe the control of cell signaling on metabolism (mainly through transcriptional regulation and post-translational modifications), much fewer information is available on the role of metabolism in the regulation of signal transduction. Protein-metabolite interactions (PMIs) result in the modification of the protein activity due to a conformational change associated with the binding of a small molecule. An increasing amount of evidences highlight the role of metabolites of the central metabolism in the control of the activity of key signaling proteins in different eukaryotic systems. Here we review the known PMIs between primary metabolites and proteins, through which metabolism affects signal transduction pathways controlled by the conserved kinases Snf1/AMPK, Ras/PKA and TORC1. Interestingly, PMIs influence also the mitochondrial retrograde response (RTG) and calcium signaling, clearly demonstrating that the range of this phenomenon is not limited to signaling pathways related to metabolism.

Dynamic design: manipulation of millisecond timescale motions on the energy landscape of cyclophilin A

Proteins need to interconvert between many conformations in order to function, many of which are formed transiently, and sparsely populated. Particularly when the lifetimes of these states approach the millisecond timescale, identifying the relevant structures and the mechanism by which they interconvert remains a tremendous challenge. Here we introduce a novel combination of accelerated MD (aMD) simulations and Markov state modelling (MSM) to explore these 'excited' conformational states. Applying this to the highly dynamic protein CypA, a protein involved in immune response and associated with HIV infection, we identify five principally populated conformational states and the atomistic mechanism by which they interconvert. A rational design strategy predicted that the mutant D66A should stabilise the minor conformations and substantially alter the dynamics, whereas the similar mutant H70A should leave the landscape broadly unchanged. These predictions are confirmed using CPMG and R1ρ solution state NMR measurements. By efficiently exploring functionally relevant, but sparsely populated conformations with millisecond lifetimes in silico, our aMD/MSM method has tremendous promise for the design of dynamic protein free energy landscapes for both protein engineering and drug discovery.

Protein-Ligand Complexes as Constrained Dynamical Systems

This study focuses on how the low-frequency end of the vibrational spectrum related to the functional motions changes as a protein binds to a small ligand(s). Our recently proposed residue-specific (RESPEC) elastic network model provides a natural laboratory for this aim due to its systematic mixed coarse-graining approach and parametrization. Current analysis on a large data set of protein-ligand complexes reveals a universal curve enclosing the frequency distributions, which bears the features of previous computational and experimental studies. We mostly observe positive frequency shifts in the collective modes of the protein upon ligand binding. This observation, conforming to the Rayleigh-Courant-Fisher theorem, points to a constraining effect imposed by ligands on protein dynamics, which may be accompanied by a negative vibrational entropy difference. Positive frequency shifts in the global modes can thus be linked to the harmonic well getting steeper, because of interactions with the ligand(s).

Exploring molecular mechanism of allosteric inhibitor to relieve drug resistance of multiple mutations in HIV-1 protease by enhanced conformational sampling

Recently, allosteric regulations of HIV-1 protease (PR) are suggested as a promising approach to relieve drug resistance of mutations toward inhibitors targeting the active site of PR. Replica-exchange molecular dynamics (REMD) simulations and normal mode analysis (NMA) are integrated to enhance conformational sampling of PR. Molecular mechanics generalized Born surface area (MM-GBSA) method was applied to calculate binding free energies of three inhibitors APV, DRV, and NIT to the wild-type (WT) and multidrug resistance (MDR) PRs. The results suggest that binding free energies of APV and DRV are decreased in the MDR PR relative to the WT PR, suggesting drug resistance of mutations on these two inhibitors. However, the binding ability of the allosteric inhibitor NIT is not impaired in the MDR PR. In addition, internal dynamics analysis based on REMD simulations proves that mutations hardly produce obvious effect on the conformation of the MDR PR in comparison to the WT PR. Scanning of hydrophobic contacts and hydrogen bond contacts of inhibitors with residues of PRs on the concatenated trajectories of REMD demonstrates that mutations change the symmetric interaction networks of APV and DRV with PR, but do not generate obvious influence on the asymmetric interaction network of NIT with PR. In summary, allosteric inhibitor NIT can adapt the MDR PR better than those inhibitors toward the active site of PR, thus allosteric inhibitors of PR may be a possible channel to overcome drug resistance of PR.

Promising new therapeutic targets for regulation of inflammation and immunity: RING-type E3 ubiquitin ligases

Ubiquitin-proteasome system (UPS) is a primary signaling pathway for regulation of protein turnover and removal of misfolded proteins in eukaryotic cells. Enzymes of the UPS pathway - E1 activating, E2 conjugating, E3 ligating - act together to covalently tag substrate proteins with a chain of ubiquitins, small regulatory proteins. The poly-ubiquitin chain then serves as a recognition motif for 26S proteasome to recognize and degrade the substrate. In recent years UPS has emerged as attractive enzymatic cascade for development of novel therapeutics against various human diseases. Building on the previous success of targeting this pathway in cancer - the broader scientific community is currently looking for ways to elucidate functions of E3 ligases, substrate-specific members of the UPS. RING-type E3 ubiquitin ligases, the largest class of E3s, represent prospective targets for small molecule modulation and their importance is reinforced by ever growing evidence of playing role in non-cancer diseases, primarily associated with inflammatory and immune disorders. In this review, we aim to briefly cover the current knowledge of biological functions of RING-type E3 ligases in inflammation and immunity.

Chemical and Biochemical Perspectives of Protein Lysine Methylation

Protein lysine methylation is a distinct posttranslational modification that causes minimal changes in the size and electrostatic status of lysine residues. Lysine methylation plays essential roles in regulating fates and functions of target proteins in an epigenetic manner. As a result, substrates and degrees (free versus mono/di/tri) of protein lysine methylation are orchestrated within cells by balanced activities of protein lysine methyltransferases (PKMTs) and demethylases (KDMs). Their dysregulation is often associated with neurological disorders, developmental abnormalities, or cancer. Methyllysine-containing proteins can be recognized by downstream effector proteins, which contain methyllysine reader domains, to relay their biological functions. While numerous efforts have been made to annotate biological roles of protein lysine methylation, limited work has been done to uncover mechanisms associated with this modification at a molecular or atomic level. Given distinct biophysical and biochemical properties of methyllysine, this review will focus on chemical and biochemical aspects in addition, recognition, and removal of this posttranslational mark. Chemical and biophysical methods to profile PKMT substrates will be discussed along with classification of PKMT inhibitors for accurate perturbation of methyltransferase activities. Semisynthesis of methyllysine-containing proteins will also be covered given the critical need for these reagents to unambiguously define functional roles of protein lysine methylation.

Small Molecule Modulators of RING-Type E3 Ligases: MDM and Cullin Families as Targets

Ubiquitin-proteasome system (UPS) is a primary signaling pathway for regulation of intracellular protein levels. E3 ubiquitin ligases, substrate-specific members of the UPS, represent highly attractive protein targets for drug discovery. The importance of E3 ligases as prospective targets for small molecule modulation is reinforced by ever growing evidence of their role in cancer and other diseases. To date the number of potent compounds targeting E3 ligases remains rather low and their rational design constitutes a challenging task. To successfully address this problem one must take into consideration the multi-subunit nature of many E3 ligases that implies multiple druggable pockets and protein-protein interfaces. In this review, we briefly cover the current state of drug discovery in the field of RING-type E3 ligases with focus on MDM and Cullin families as targets. We also provide an overview of small molecule chimeras that induce RING-type E3-mediated proteasomal degradation of substrate proteins of interest.

Dependence of prevalence of contiguous pathways in proteins on structural complexity

Allostery is a regulatory mechanism in proteins where an effector molecule binds distal from an active site to modulate its activity. Allosteric signaling may occur via a continuous path of residues linking the active and allosteric sites, which has been suggested by large conformational changes evident in crystal structures. An alternate possibility is that the signal occurs in the realm of ensemble dynamics via an energy landscape change. While the latter was first proposed on theoretical grounds, increasing evidence suggests that such a control mechanism is plausible. A major difficulty for testing the two methods is the ability to definitively determine that a residue is directly involved in allosteric signal transduction. Statistical Coupling Analysis (SCA) is a method that has been successful at predicting pathways, and experimental tests involving mutagenesis or domain substitution provide the best available evidence of signaling pathways. However, ascertaining energetic pathways which need not be contiguous is far more difficult. To date, simple estimates of the statistical significance of a pathway in a protein remain to be established. The focus of this work is to estimate such benchmarks for the statistical significance of contiguous pathways for the null model of selecting residues at random. We found that when 20% of residues in proteins are randomly selected, contiguous pathways at the 6 Å cutoff level were found with success rates of 51% in PDZ, 30% in p53, and 3% in MutS. The results suggest that the significance of pathways may have system specific factors involved. Furthermore, the possible existence of false positives for contiguous pathways implies that signaling could be occurring via alternate routes including those consistent with the energetic landscape model.

Recent excitements in protein NMR: Large proteins and biologically relevant dynamics

The advent of Transverse Relaxation Optimized SpectroscopY (TROSY) and perdeuteration allowed biomolecular NMR spectroscopists to overcome the size limitation barrier (approx. 20 kDa) in de novo structure determination of proteins. The utility of these techniques was immediately demonstrated on large proteins and protein complexes (e.g. GroELGroES, ClpP protease, Hsp90-p53, 20S proteasome, etc.). Further, recent methodological developments such as Residual Dipolar Couplings and Paramagnetic Relaxation Enhancement allowed accurate measurement of long-range structural restraints. Additionally, Carr-Purcell-Meiboom-Gill (CPMG), rotating frame relaxation experiments (R1(rho)) and saturation transfer experiments (CEST and DEST) created never-before accessibility to the (mu)s-ms timescale dynamic parameters that led to the deeper understanding of biological processes. Meanwhile, the excitement in the field continued with a series of developments in the fast data acquisition methods allowing rapid structural studies on less stable proteins. This review aims to discuss important developments in the field of biomolecular NMR spectroscopy in the recent past, i.e., in the post TROSY era. These developments not only gave access to the structural studies of large protein assemblies, but also revolutionized tools in the arsenal of today's biomolecular NMR and point to a bright future of biomolecular NMR spectroscopy.

Oncogenic K-Ras Binds to an Anionic Membrane in Two Distinct Orientations: A Molecular Dynamics Analysis

K-Ras is a membrane-associated GTPase that cycles between active and inactive conformational states to regulate a variety of cell signaling pathways. Somatic mutations in K-Ras are linked to 15-20% of all human tumors. K-Ras attaches to the inner leaflet of the plasma membrane via a farnesylated polybasic domain; however, the structural details of the complex remain poorly understood. Based on extensive (7.5 μs total) atomistic molecular dynamics simulations here we show that oncogenic mutant K-Ras interacts with a negatively charged lipid bilayer membrane in multiple orientations. Of these, two highly populated orientations account for ∼54% of the conformers whose catalytic domain directly interacts with the bilayer. In one of these orientation states, membrane binding involves helices 3 and 4 of the catalytic domain in addition to the farnesyl and polybasic motifs. In the other orientation, β-strands 1-3 and helix 2 on the opposite face of the catalytic domain contribute to membrane binding. Flexibility of the linker region was found to be important for the reorientation. The biological significance of these observations was evaluated by initial experiments in cells overexpressing mutant K-Ras as well as by an analysis of Ras-effector complex structures. The results suggest that only one of the two major orientation states is capable of effector binding. We propose that the different modes of membrane binding may be exploited in structure-based drug design efforts for cancer therapy.

Effect of allosteric molecules on structure and drug affinity of HIV-1 protease by molecular dynamics simulations

Recent experiments show that small molecules can bind onto the allosteric sites of HIV-1 protease (PR), which provides a starting point for developing allosteric inhibitors. However, the knowledge of the effect of such binding on the structural dynamics and binding free energy of the active site inhibitor and PR is still lacking. Here, we report 200ns long molecular dynamics simulation results to gain insight into the influences of two allosteric molecules (1H-indole-6-carboxylic acid, 1F1 and 2-methylcyclohexano, 4D9). The simulations demonstrate that both allosteric molecules change the PR conformation and stabilize the structures of PR and the inhibitor; the residues of the flaps are sensitive to the allosteric molecules and the flexibility of the residues is pronouncedly suppressed; the additions of the small molecules to the allosteric sites strengthen the binding affinities of 3TL-PR by about 12-15kal/mol in the binding free energy, which mainly arises from electrostatic term. Interestingly, it is found that the action mechanisms of 1F1 and 4D9 are different, the former behaviors like a doorman that keeps the inhibitor from escape and makes the flaps (door) partially open; the latter is like a wedge that expands the allosteric space and meanwhile closes the flaps. Our data provide a theoretical support for designing the allosteric inhibitor.

Ensembler: Enabling High-Throughput Molecular Simulations at the Superfamily Scale

The rapidly expanding body of available genomic and protein structural data provides a rich resource for understanding protein dynamics with biomolecular simulation. While computational infrastructure has grown rapidly, simulations on an omics scale are not yet widespread, primarily because software infrastructure to enable simulations at this scale has not kept pace. It should now be possible to study protein dynamics across entire (super)families, exploiting both available structural biology data and conformational similarities across homologous proteins. Here, we present a new tool for enabling high-throughput simulation in the genomics era. Ensembler takes any set of sequences—from a single sequence to an entire superfamily—and shepherds them through various stages of modeling and refinement to produce simulation-ready structures. This includes comparative modeling to all relevant PDB structures (which may span multiple conformational states of interest), reconstruction of missing loops, addition of missing atoms, culling of nearly identical structures, assignment of appropriate protonation states, solvation in explicit solvent, and refinement and filtering with molecular simulation to ensure stable simulation. The output of this pipeline is an ensemble of structures ready for subsequent molecular simulations using computer clusters, supercomputers, or distributed computing projects like Folding@home. Ensembler thus automates much of the time-consuming process of preparing protein models suitable for simulation, while allowing scalability up to entire superfamilies. A particular advantage of this approach can be found in the construction of kinetic models of conformational dynamics—such as Markov state models (MSMs)—which benefit from a diverse array of initial configurations that span the accessible conformational states to aid sampling. We demonstrate the power of this approach by constructing models for all catalytic domains in the human tyrosine kinase family, using all available kinase catalytic domain structures from any organism as structural templates. Ensembler is free and open source software licensed under the GNU General Public License (GPL) v2. It is compatible with Linux and OS X. The latest release can be installed via the conda package manager, and the latest source can be downloaded from https://github.com/choderalab/ensembler.

Proteins are the workhorses of the human body, and are involved in essentially every biological process. Many diseases are caused by proteins malfunctioning. To understand how a protein functions, it is necessary to know its physical properties. The field of structural biology provides many techniques for determining the three-dimensional structure of a protein. The dynamics of a protein, i.e. the way it moves, are of equal importance, but are more difficult to uncover with traditional experimental techniques. Computer simulations are an effective alternative method for understanding protein dynamics, but require experimental structural information as a starting point. While recent advances in genomics and experimental techniques have provided a wealth of such structural data, the appropriate software for using this data effectively has been lacking. To tackle this problem, we have developed a software package called Ensembler, which allows a user to automatically select appropriate experimentally derived structures for a given protein or family of proteins, and to use them to prepare a series of simulations. The resultant simulation data can then used to investigate the dynamics of the protein(s) in question, and their involvement in disease.

Mapping the Dynamics Landscape of Conformational Transitions in Enzyme: The Adenylate Kinase Case

Conformational transition describes the essential dynamics and mechanism of enzymes in pursuing their various functions. The fundamental and practical challenge to researchers is to quantitatively describe the roles of large-scale dynamic transitions for regulating the catalytic processes. In this study, we tackled this challenge by exploring the pathways and free energy landscape of conformational changes in adenylate kinase (AdK), a key ubiquitous enzyme for cellular energy homeostasis. Using explicit long-timescale (up to microseconds) molecular dynamics and bias-exchange metadynamics simulations, we determined at the atomistic level the intermediate conformational states and mapped the transition pathways of AdK in the presence and absence of ligands. There is clearly chronological operation of the functional domains of AdK. Specifically in the ligand-free AdK, there is no significant energy barrier in the free energy landscape separating the open and closed states. Instead there are multiple intermediate conformational states, which facilitate the rapid transitions of AdK. In the ligand-bound AdK, the closed conformation is energetically most favored with a large energy barrier to open it up, and the conformational population prefers to shift to the closed form coupled with transitions. The results suggest a perspective for a hybrid of conformational selection and induced fit operations of ligand binding to AdK. These observations, depicted in the most comprehensive and quantitative way to date, to our knowledge, emphasize the underlying intrinsic dynamics of AdK and reveal the sophisticated conformational transitions of AdK in fulfilling its enzymatic functions. The developed methodology can also apply to other proteins and biomolecular systems.

Antibody humanization by molecular dynamics simulations-in-silico guided selection of critical backmutations

Monoclonal antibodies represent the fastest growing class of biotherapeutic proteins. However, as they are often initially derived from rodent organisms, there is a severe risk of immunogenic reactions, hampering their applicability. The humanization of these antibodies remains a challenging task in the context of rational drug design. "Superhumanization" describes the direct transfer of the complementarity determining regions to a human germline framework, but this humanization approach often results in loss of binding affinity. In this study, we present a new approach for predicting promising backmutation sites using molecular dynamics simulations of the model antibody Ab2/3H6. The simulation method was developed in close conjunction with novel specificity experiments. Binding properties of mAb variants were evaluated directly from crude supernatants and confirmed using established binding affinity assays for purified antibodies. Our approach provides access to the dynamical features of the actual binding sites of an antibody, based solely on the antibody sequence. Thus we do not need structural data on the antibody-antigen complex and circumvent cumbersome methods to assess binding affinities. © 2016 The Authors Journal of Molecular Recognition Published by John Wiley & Sons Ltd.

Intracellular Generation of a Diterpene-Peptide Conjugate that Inhibits 14-3-3-Mediated Interactions

Synthetic agents that disrupt intracellular protein-protein interactions (PPIs) are highly desirable for elucidating signaling networks and developing new therapeutics. However, designing cell-penetrating large molecules equipped with the many functional groups necessary for binding to large interfaces remains challenging. Here, we describe a rational strategy for the intracellular oxime ligation-mediated generation of an amphipathic bivalent inhibitor composed of a peptide and diterpene natural product, fusicoccin, which binds 14-3-3 protein with submicromolar affinity. Our results demonstrate that co-treatment of cells with small module molecules, the aldehyde-containing fusicoccin 1 and the aminooxy-containing peptide 2, generates the corresponding conjugate 3 in cells, resulting in significant cytotoxicity. In contrast, chemically synthesized 3 is not cytotoxic, likely due to its inability to penetrate cells. Compound 3, but not 1 or 2, disrupts endogenous 14-3-3/cRaf interactions, suggesting that cell death is caused by inhibition of 14-3-3 activity. These results suggest that intracellular generation of large-sized molecules may serve as a new approach for modulating PPIs.

Inhibitors of protein methyltransferases as chemical tools

Protein methyltransferases (PMT)s play essential roles in many biological processes through methylation of histones and diverse nonhistone substrates. Dysregulation of these enzymes has been implicated in many diseases including cancers. While PMT-associated biology can be probed via genetic perturbation, this approach targets full-length PMTs rather than their methyltransferase activities and often lacks temporal, spatial and dose controls (timing, location and amount of dosed compounds). By contrast, small-molecule inhibitors of PMTs can be designed to specifically target the methyltransferase domains in a temporal, spatial and dose-dependent manner. This utility has motivated the development of hundreds of PMT inhibitors, but meanwhile can make it challenging to select the most suitable PMT inhibitors to interrogate PMT-associated biology. This perspective aims to provide timely guidance to evaluate these PMT inhibitors in their relevant biological contexts.

Targeting Cullin-RING E3 ubiquitin ligases for drug discovery: structure, assembly and small-molecule modulation

In the last decade, the ubiquitin–proteasome system has emerged as a valid target for the development of novel therapeutics. E3 ubiquitin ligases are particularly attractive targets because they confer substrate specificity on the ubiquitin system. CRLs [Cullin–RING (really interesting new gene) E3 ubiquitin ligases] draw particular attention, being the largest family of E3s. The CRLs assemble into functional multisubunit complexes using a repertoire of substrate receptors, adaptors, Cullin scaffolds and RING-box proteins. Drug discovery targeting CRLs is growing in importance due to mounting evidence pointing to significant roles of these enzymes in diverse biological processes and human diseases, including cancer, where CRLs and their substrates often function as tumour suppressors or oncogenes. In the present review, we provide an account of the assembly and structure of CRL complexes, and outline the current state of the field in terms of available knowledge of small-molecule inhibitors and modulators of CRL activity. A comprehensive overview of the reported crystal structures of CRL subunits, components and full-size complexes, alone or with bound small molecules and substrate peptides, is included. This information is providing increasing opportunities to aid the rational structure-based design of chemical probes and potential small-molecule therapeutics targeting CRLs.

Binding hotspots on K-ras: consensus ligand binding sites and other reactive regions from probe-based molecular dynamics analysis

We have used probe-based molecular dynamics (pMD) simulations to search for interaction hotspots on the surface of the therapeutically highly relevant oncogenic K-Ras G12D. Combining the probe-based query with an ensemble-based pocket identification scheme and an analysis of existing Ras-ligand complexes, we show that (i) pMD is a robust and cost-effective strategy for binding site identification, (ii) all four of the previously reported ligand binding sites are suitable for structure-based ligand design, and (iii) in some cases probe binding and expanded sampling of configurational space enable pocket expansion and increase the likelihood of site identification. Furthermore, by comparing the distribution of hotspots in nonpocket-like regions with known protein- and membrane-interacting interfaces, we propose that pMD has the potential to predict surface patches responsible for protein-biomolecule interactions. These observations have important implications for future drug design efforts and will facilitate the search for potential interfaces responsible for the proposed transient oligomerization or interaction of Ras with other biomolecules in the cellular milieu.

Identification of potential small molecule allosteric modulator sites on IL-1R1 ectodomain using accelerated conformational sampling method

The interleukin-1 receptor (IL-1R) is the founding member of the interleukin 1 receptor family which activates innate immune response by its binding to cytokines. Reports showed dysregulation of cytokine production leads to aberrant immune cells activation which contributes to auto-inflammatory disorders and diseases. Current therapeutic strategies focus on utilizing antibodies or chimeric cytokine biologics. The large protein-protein interaction interface between cytokine receptor and cytokine poses a challenge in identifying binding sites for small molecule inhibitor development. Based on the significant conformational change of IL-1R type 1 (IL-1R1) ectodomain upon binding to different ligands observed in crystal structures, we hypothesized that transient small molecule binding sites may exist when IL-1R1 undergoes conformational transition and thus suitable for inhibitor development. Here, we employed accelerated molecular dynamics (MD) simulation to efficiently sample conformational space of IL-1R1 ectodomain. Representative IL-1R1 ectodomain conformations determined from the hierarchy cluster analysis were analyzed by the SiteMap program which leads to identify small molecule binding sites at the protein-protein interaction interface and allosteric modulator locations. The cosolvent mapping analysis using phenol as the probe molecule further confirms the allosteric modulator site as a binding hotspot. Eight highest ranked fragment molecules identified from in silico screening at the modulator site were evaluated by MD simulations. Four of them restricted the IL-1R1 dynamical motion to inactive conformational space. The strategy from this study, subject to in vitro experimental validation, can be useful to identify small molecule compounds targeting the allosteric modulator sites of IL-1R and prevent IL-1R from binding to cytokine by trapping IL-1R in inactive conformations.

Perspectives on electrostatics and conformational motions in enzyme catalysis

Enzymes are essential for all living organisms, and their effectiveness as chemical catalysts has driven more than a half century of research seeking to understand the enormous rate enhancements they provide. Nevertheless, a complete understanding of the factors that govern the rate enhancements and selectivities of enzymes remains elusive, due to the extraordinary complexity and cooperativity that are the hallmarks of these biomolecules. We have used a combination of site-directed mutagenesis, pre-steady-state kinetics, X-ray crystallography, nuclear magnetic resonance (NMR), vibrational and fluorescence spectroscopies, resonance energy transfer, and computer simulations to study the implications of conformational motions and electrostatic interactions on enzyme catalysis in the enzyme dihydrofolate reductase (DHFR).

We have demonstrated that modest equilibrium conformational changes are functionally related to the hydride transfer reaction. Results obtained for mutant DHFRs illustrated that reductions in hydride transfer rates are correlated with altered conformational motions, and analysis of the evolutionary history of DHFR indicated that mutations appear to have occurred to preserve both the hydride transfer rate and the associated conformational changes. More recent results suggested that differences in local electrostatic environments contribute to finely tuning the substrate pK_a in the initial protonation step. Using a combination of primary and solvent kinetic isotope effects, we demonstrated that the reaction mechanism is consistent across a broad pH range, and computer simulations suggested that deprotonation of the active site Tyr100 may play a crucial role in substrate protonation at high pH.

Site-specific incorporation of vibrational thiocyanate probes into the ecDHFR active site provided an experimental tool for interrogating these microenvironments and for investigating changes in electrostatics along the DHFR catalytic cycle. Complementary molecular dynamics simulations in conjunction with mixed quantum mechanical/molecular mechanical calculations accurately reproduced the vibrational frequency shifts in these probes and provided atomic-level insight into the residues influencing these changes. Our findings indicate that conformational and electrostatic changes are intimately related and functionally essential. This approach can be readily extended to the study of other enzyme systems to identify more general trends in the relationship between conformational fluctuations and electrostatic interactions. These results are relevant to researchers seeking to design novel enzymes as well as those seeking to develop therapeutic agents that function as enzyme inhibitors.

Shedding light on protein folding, structural and functional dynamics by single molecule studies

The advent of advanced single molecule measurements unveiled a great wealth of dynamic information revolutionizing our understanding of protein dynamics and behavior in ways unattainable by conventional bulk assays. Equipped with the ability to record distribution of behaviors rather than the mean property of a population, single molecule measurements offer observation and quantification of the abundance, lifetime and function of multiple protein states. They also permit the direct observation of the transient and rarely populated intermediates in the energy landscape that are typically averaged out in non-synchronized ensemble measurements. Single molecule studies have thus provided novel insights about how the dynamic sampling of the free energy landscape dictates all aspects of protein behavior; from its folding to function. Here we will survey some of the state of the art contributions in deciphering mechanisms that underlie protein folding, structural and functional dynamics by single molecule fluorescence microscopy techniques. We will discuss a few selected examples highlighting the power of the emerging techniques and finally discuss the future improvements and directions.

A label-free untethered approach to single-molecule protein binding kinetics

Single molecule approaches provide rich real-time dynamics of molecular interactions that are not accessible to ensemble measurements. Previous single molecule studies have relied on labeling and tethering, which alters the natural state of the protein. Here we use the double-nanohole (DNH) optical tweezer approach to measure protein binding kinetics at the single molecule level in a label-free, free-solution (untethered) way. The binding kinetics of human serum albumin (HSA) to tolbutamide and to phenytoin are in quantitative agreement with previous measurements, and our single-molecule approach reveals a biexponential behavior characteristic of a multistep process. The DNH optical tweezer is an inexpensive platform for studying the real-time binding kinetics of protein-small molecule interactions in a label-free, free-solution environment, which will be of interest to future studies including drug discovery.

Broad-spectrum allosteric inhibition of herpesvirus proteases

Herpesviruses rely on a homodimeric protease for viral capsid maturation. A small molecule, DD2, previously shown to disrupt dimerization of Kaposi’s sarcoma-associated herpesvirus protease (KSHV Pr) by trapping an inactive monomeric conformation and two analogues generated through carboxylate bioisosteric replacement (compounds 2 and 3) were shown to inhibit the associated proteases of all three human herpesvirus (HHV) subfamilies (α, β, and γ). Inhibition data reveal that compound 2 has potency comparable to or better than that of DD2 against the tested proteases. Nuclear magnetic resonance spectroscopy and a new application of the kinetic analysis developed by Zhang and Poorman [Zhang, Z. Y., Poorman, R. A., et al. (1991) J. Biol. Chem. 266, 15591–15594] show DD2, compound 2, and compound 3 inhibit HHV proteases by dimer disruption. All three compounds bind the dimer interface of other HHV proteases in a manner analogous to binding of DD2 to KSHV protease. The determination and analysis of cocrystal structures of both analogues with the KSHV Pr monomer verify and elaborate on the mode of binding for this chemical scaffold, explaining a newly observed critical structure–activity relationship. These results reveal a prototypical chemical scaffold for broad-spectrum allosteric inhibition of human herpesvirus proteases and an approach for the identification of small molecules that allosterically regulate protein activity by targeting protein–protein interactions.

From a structural average to the conformational ensemble of a DNA bulge

Direct experimental measurements of conformational ensembles are critical for understanding macromolecular function, but traditional biophysical methods do not directly report the solution ensemble of a macromolecule. Small-angle X-ray scattering interferometry has the potential to overcome this limitation by providing the instantaneous distance distribution between pairs of gold-nanocrystal probes conjugated to a macromolecule in solution. Our X-ray interferometry experiments reveal an increasing bend angle of DNA duplexes with bulges of one, three, and five adenosine residues, consistent with previous FRET measurements, and further reveal an increasingly broad conformational ensemble with increasing bulge length. The distance distributions for the AAA bulge duplex (3A-DNA) with six different Au-Au pairs provide strong evidence against a simple elastic model in which fluctuations occur about a single conformational state. Instead, the measured distance distributions suggest a 3A-DNA ensemble with multiple conformational states predominantly across a region of conformational space with bend angles between 24 and 85 degrees and characteristic bend directions and helical twists and displacements. Additional X-ray interferometry experiments revealed perturbations to the ensemble from changes in ionic conditions and the bulge sequence, effects that can be understood in terms of electrostatic and stacking contributions to the ensemble and that demonstrate the sensitivity of X-ray interferometry. Combining X-ray interferometry ensemble data with molecular dynamics simulations gave atomic-level models of representative conformational states and of the molecular interactions that may shape the ensemble, and fluorescence measurements with 2-aminopurine-substituted 3A-DNA provided initial tests of these atomistic models. More generally, X-ray interferometry will provide powerful benchmarks for testing and developing computational methods.

Asymmetric perturbations of signalling oligomers

This review focuses on rapid and reversible noncovalent interactions for symmetric oligomers of signalling proteins. Symmetry mismatch, transient symmetry breaking and asymmetric perturbations via chemical (ligand binding) and physical (electric or mechanic) effects can initiate the signalling events. Advanced biophysical methods can reveal not only structural symmetries of stable membrane-bound signalling proteins but also asymmetric functional transition states. Relevant techniques amenable to distinguish between symmetric and asymmetric architectures are discussed including those with the capability of capturing low-populated transient conformational states. Typical examples of signalling proteins are overviewed for symmetry breaking in dimers (GPCRs, growth factor receptors, transcription factors); trimers (acid-sensing ion channels); tetramers (voltage-gated cation channels, ionotropic glutamate receptor, CNG and CHN channels); pentameric ligand-gated and mechanosensitive channels; higher order oligomers (gap junction channel, chaperonins, proteasome, virus capsid); as well as primary and secondary transporters. In conclusion, asymmetric perturbations seem to play important functional roles in a broad range of communicating networks.

Structural interrogation of phosphoproteome identified by mass spectrometry reveals allowed and disallowed regions of phosphoconformation

Background

High-throughput mass spectrometric (HT-MS) study is the method of choice for monitoring global changes in proteome. Data derived from these studies are meant for further validation and experimentation to discover novel biological insights. Here we evaluate use of relative solvent accessible surface area (rSASA) and DEPTH as indices to assess experimentally determined phosphorylation events deposited in PhosphoSitePlus.

Results

Based on accessibility, we map these identifications on allowed (accessible) or disallowed (inaccessible) regions of phosphoconformation. Surprisingly a striking number of HT-MS/MS derived events (1461/5947 sites or 24.6%) are present in the disallowed region of conformation. By considering protein dynamics, autophosphorylation events and/or the sequence specificity of kinases, 13.8% of these phosphosites can be moved to the allowed region of conformation. We also demonstrate that rSASA values can be used to increase the confidence of identification of phosphorylation sites within an ambiguous MS dataset.

Conclusion

While MS is a stand-alone technique for the identification of vast majority of phosphorylation events, identifications within disallowed region of conformation will benefit from techniques that independently probe for phosphorylation and protein dynamics. Our studies also imply that trapping alternate protein conformations may be a viable alternative to the design of inhibitors against mutation prone drug resistance kinases.

A general protocol for the generation of Nanobodies for structural biology

There is growing interest in using antibodies as auxiliary tools to crystallize proteins. Here we describe a general protocol for the generation of Nanobodies to be used as crystallization chaperones for the structural investigation of diverse conformational states of flexible (membrane) proteins and complexes thereof. Our technology has a competitive advantage over other recombinant crystallization chaperones in that we fully exploit the natural humoral response against native antigens. Accordingly, we provide detailed protocols for the immunization with native proteins and for the selection by phage display of in vivo-matured Nanobodies that bind conformational epitopes of functional proteins. Three representative examples illustrate that the outlined procedures are robust, making it possible to solve by Nanobody-assisted X-ray crystallography in a time span of 6-12 months.

INSIGHTS IN ENZYME FUNCTIONAL DYNAMICS AND ACTIVITY REGULATION BY SINGLE MOLECULE STUDIES

The advent of advanced single molecule measurements heralded the arrival of a wealth of dynamic information revolutionizing our understanding of protein dynamics and behavior in ways not deducible by conventional bulk assays. They offered the direct observation and quantification of the abundance and life time of multiple states and transient intermediates in the energy landscape that are typically averaged out in non-synchronized ensemble measurements, thus providing unprecedented insights into complex biological processes. Here we survey the current state of the art in single-molecule fluorescence microscopy methodology for studying the mechanism of enzymatic activity and the insights on protein functional dynamics. We will initially discuss the strategies employed to date, their limitations and possible ways to overcome them, and finally how single enzyme kinetics can advance our understanding on mechanisms underlying function and regulation of proteins.

[Formula: see text]Special Issue Comment: This review focuses on functional dynamics of individual enzymes and is related to the review on ion channels by Lu,44 the reviews on mathematical treatment of Flomenbom45 and Sach et al.,46 and review on FRET by Ruedas-Rama et al.41

Type 1 and Type 2 scenarios in hydrogen exchange mass spectrometry studies on protein–ligand complexes

Ligand binding to a protein can elicit a wide range of responses when studied by HDX mass spectrometry.