Mechanism of Coupled Folding-upon-Binding of an Intrinsically Disordered Protein

The transition state for coupled folding and binding of a disordered DNA binding domain resembles the unbound state

The basic zippers (bZIPs) are one of two large eukaryotic families of transcription factors whose DNA binding domains are disordered in isolation but fold into stable α-helices upon target DNA binding. Here, we systematically disrupt pre-existing helical propensity within the DNA binding region of the homodimeric bZIP domain of cAMP-response element binding protein (CREB) using Ala-Gly scanning and examine the impact on target binding kinetics. We find that the secondary structure of the transition state strongly resembles that of the unbound state. The residue closest to the dimerization domain is largely folded within both unbound and transition states; dimerization apparently propagates additional helical propensity into the basic region. The results are consistent with electrostatically-enhanced DNA binding, followed by rapid folding from the folded zipper outwards. Fly-casting theory suggests that protein disorder can accelerate binding. Interestingly however, we did not observe higher association rate constants for mutants with lower levels of residual structure in the unbound state.

Alchemical Transformations and Beyond: Recent Advances and Real-World Applications of Free Energy Calculations in Drug Discovery

Computational methods constitute efficient strategies for screening and optimizing potential drug molecules. A critical factor in this process is the binding affinity between candidate molecules and targets, quantified as binding free energy. Among various estimation methods, alchemical transformation methods stand out for their theoretical rigor. Despite challenges in force field accuracy and sampling efficiency, advancements in algorithms, software, and hardware have increased the application of free energy perturbation (FEP) calculations in the pharmaceutical industry. Here, we review the practical applications of FEP in drug discovery projects since 2018, covering both ligand-centric and residue-centric transformations. We show that relative binding free energy calculations have steadily achieved chemical accuracy in real-world applications. In addition, we discuss alternative physics-based simulation methods and the incorporation of deep learning into free energy calculations.

Statistical accuracy of molecular dynamics-based methods for sampling conformational ensembles of disordered proteins

The characterization of the statistical ensemble of conformations of intrinsically disordered regions (IDRs) is a great challenge both from experimental and computational points of view. In this respect, a number of protocols have been developed using molecular dynamics (MD) simulations to sample the huge conformational space of the molecule. In this work, we consider one of the best methods available, replica exchange solute tempering (REST), as a reference to compare the results obtained using this method with the results obtained using other methods, in terms of experimentally measurable quantities. Along with the methods assessed, we propose here a novel protocol called probabilistic MD chain growth (PMD-CG), which combines the flexible-meccano and hierarchical chain growth methods with the statistical data obtained from tripeptide MD trajectories as the starting point. The system chosen for testing is a 20-residue region from the C-terminal domain of the p53 tumor suppressor protein (p53-CTD). Our results show that PMD-CG provides an ensemble of conformations extremely quickly, after suitable computation of the conformational pool for all peptide triplets of the IDR sequence. The measurable quantities computed on the ensemble of conformations agree well with those based on the REST conformational ensemble.

Extended Conformational Selection in the Antigen-Antibody Interaction of the PfAMA1 Protein

Plasmodium falciparum apical membrane antigen 1 (PfAMA1) is a surface protein found in two stages of the malaria life cycle. This is a protein involved in a reorientation movement of the parasite so that cell invasion occurs in the so-called "moving junction", relevant when the membranes of the parasite and the host are in contact. The structure of a conformational epitope of domain III of PfAMA1 in complex with the monoclonal antibody Fab F8.12.19 is experimentally known. Here, we used molecular dynamics with enhanced sampling by Hamiltonian replica exchange molecular dynamics (HREMD) to understand the effect of intermolecular interactions, conformational variability, and intrinsically disordered regions on the mechanism of antigen-antibody interaction. Clustering methods and the analysis of conformational variability were used in order to understand the influence of the presence of the partner protein in the complex. The free-state epitope accesses a broader conformational pool, including disordered conformations not seen in the bound state. The simulations suggest an extended conformational selection mechanism in which the antibody stabilizes a conformational set of the epitope existing in the free state. The stabilization of the active conformation occurs mainly through hydrogen bonds: Tyr(H33)-Asp493, His(L94)-Val510, Ser(L93)-Glu511, Tyr(H56)-Asp485, and Tyr(H35)-Asp493. The antibody has a structure with few flexible regions, and only the complementarity determining region (CDR) H3 shows greater plasticity in the presence of the epitope.

Predicting Conformational Ensembles of Intrinsically Disordered Proteins: From Molecular Dynamics to Machine Learning

Intrinsically disordered proteins and regions (IDP/IDRs) are ubiquitous across all domains of life. Characterized by a lack of a stable tertiary structure, IDP/IDRs populate a diverse set of transiently formed structural states that can promiscuously adapt upon binding with specific interaction partners and/or certain alterations in environmental conditions. This malleability is foundational for their role as tunable interaction hubs in core cellular processes such as signaling, transcription, and translation. Tracing the conformational ensemble of an IDP/IDR and its perturbation in response to regulatory cues is thus paramount for illuminating its function. However, the conformational heterogeneity of IDP/IDRs poses several challenges. Here, we review experimental and computational methods devised to disentangle the conformational landscape of IDP/IDRs, highlighting recent computational advances that permit proteome-wide scans of IDP/IDRs conformations. We briefly evaluate selected computational methods using the disordered N-terminal of the human copper transporter 1 as a test case and outline further challenges in IDP/IDRs ensemble prediction.

Deciphering the Mystery in p300 Taz2-p53 TAD2 Recognition

Intrinsically disordered proteins (IDPs) engage in various fundamental biological activities, and their behavior is of particular importance for a better understanding of the verbose but well-organized signal transduction in cells. IDPs exhibit uniquely paradoxical features with low affinity but simultaneously high specificity in recognizing their binding targets. The transcription factor p53 plays a crucial role in cancer suppression, carrying out some of its biological functions using its disordered regions, such as N-terminal transactivation domain 2 (TAD2). Exploration of the binding and unbinding processes between proteins is challenging, and the inherently disordered properties of these regions further complicate the issue. Computer simulations are a powerful tool to complement the experiments to fill gaps to explore the binding/unbinding processes between proteins. Here, we investigated the binding mechanism between p300 Taz2 and p53 TAD2 through extensive molecular dynamics (MD) simulations using the physics-based UNited RESidue (UNRES) force field with additional Go̅-like potentials. Distance restraints extracted from the NMR-resolved structures were imposed on intermolecular residue pairs to accelerate binding simulations, in which Taz2 was immobilized in a native-like conformation and disordered TAD2 was fully free. Starting from six structures with TAD2 placed at different positions around Taz2, we observed a metastable intermediate state in which the middle helical segment of TAD2 is anchored in the binding pocket, highlighting the significance of the TAD2 helix in directing protein recognition. Physics-based binding simulations show that successful binding is achieved after a series of stages, including (1) protein collisions to initiate the formation of encounter complexes, (2) partial attachment of TAD2, and finally (3) full attachment of TAD2 to the correct binding pocket of Taz2. Furthermore, machine-learning-based PathDetect-SOM was used to identify two binding pathways, the encounter complexes, and the intermediate states.

PepBinding: A Workflow for Predicting Peptide Binding Structures by Combining Peptide Docking and Peptide Gaussian Accelerated Molecular Dynamics Simulations

Predicting protein-peptide interactions is crucial for understanding peptide binding processes and designing peptide drugs. However, traditional computational modeling approaches face challenges in accurately predicting peptide-protein binding structures due to the slow dynamics and high flexibility of the peptides. Here, we introduce a new workflow termed "PepBinding" for predicting peptide binding structures, which combines peptide docking, all-atom enhanced sampling simulations using the Peptide Gaussian accelerated Molecular Dynamics (Pep-GaMD) method, and structural clustering. PepBinding has been demonstrated on seven distinct model peptides. In peptide docking using HPEPDOCK, the peptide backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures ranged from 3.8 to 16.0 Å, corresponding to the medium to inaccurate quality models according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. The Pep-GaMD simulations performed for only 200 ns significantly improved the docking models, resulting in five medium and two acceptable quality models. Therefore, PepBinding is an efficient workflow for predicting peptide binding structures and is publicly available at https://github.com/MiaoLab20/PepBinding.

Ligand Gaussian Accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional molecular dynamics (cMD), due to limited simulation time scales. Based on our previously developed ligand Gaussian accelerated molecular dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3″, in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding, and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as the model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 μs simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 were in agreement with the available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligands and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.

Transient Non-local Interactions Dominate the Dynamics of Measles Virus NTAIL

The RNA genome of measles virus is encapsidated by the nucleoprotein within a helical nucleocapsid that serves as template for both transcription and replication. The intrinsically disordered domain of the nucleoprotein (NTAIL), partly protruding outward from the nucleocapsid, is essential for binding the polymerase complex responsible for viral transcription and replication. As for many IDPs, binding of NTAIL occurs through a short molecular recognition element (MoRE) that folds upon binding, with the majority of NTAIL remaining disordered. Though NTAIL regions far from the MoRE influence the binding affinity, interactions between them and the MoRE have not been investigated in depth. Using an integrated approach, relying on photo-induced electron transfer (PET) experiments between tryptophan and cysteine pairs placed at different positions in the protein under varying salt and pH conditions, combined with simulations and analytical models, we identified transient interactions between two disordered regions distant in sequence, which dominate NTAIL dynamics, and regulate the conformational preferences of both the MoRE and the entire NTAIL domain. Co-evolutionary analysis corroborates our findings, and suggests an important functional role for the same intramolecular interactions. We propose mechanisms by which these non-local interactions may regulate binding to the phosphoprotein, polymerase recruitment, and ultimately viral transcription and replication. Our findings may be extended to other IDPs, where non-local intra-protein interactions affect the conformational preferences of intermolecular binding sites.

Impact of Nonnative Interactions on the Binding Kinetics of Intrinsically Disordered p53 with MDM2: Insights from All-Atom Simulation and Markov State Model Analysis

Intrinsically disordered proteins (IDPs) lack a well-defined tertiary structure but are essential players in various biological processes. Their ability to undergo a disorder-to-order transition upon binding to their partners, known as the folding-upon-binding process, is crucial for their function. One classical example is the intrinsically disordered transactivation domain (TAD) of the tumor suppressor protein p53, which quickly forms a structured α-helix after binding to its partner MDM2, with clinical significance for cancer treatment. However, the contribution of nonnative interactions between the IDP and its partner to the rapid binding kinetics, as well as their interplay with native interactions, is not well understood at the atomic level. Here, we used molecular dynamics simulation and Markov state model (MSM) analysis to study the folding-upon-binding mechanism between p53-TAD and MDM2. Our results suggest that the system progresses from the nascent encounter complex to the well-structured encounter complex and finally reaches the native complex, following an induced-fit mechanism. We found that nonnative hydrophobic and hydrogen bond interactions, combined with native interactions, effectively stabilize the nascent and well-structured encounter complexes. Among the nonnative interactions, Leu25p53-Leu54MDM2 and Leu25p53-Phe55MDM2 are particularly noteworthy, as their interaction strength is close to the optimum. Evidently, strengthening or weakening these interactions could both adversely affect the binding kinetics. Overall, our findings suggest that nonnative interactions are evolutionarily optimized to accelerate the binding kinetics of IDPs in conjunction with native interactions.

Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional Molecular Dynamics (cMD), due to limited simulation timescales. Based on our previously developed Ligand Gaussian accelerated Molecular Dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3", in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 microsecond simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 agreed with available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligand and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.

High-throughput assay exploiting disorder-to-order conformational switches: application to the proteasomal Rpn10:E6AP complex

Conformational switching is pervasively driven by protein interactions, particularly for intrinsically disordered binding partners. We developed a dually orthogonal fluorescence-based assay to monitor such events, exploiting environmentally sensitive fluorophores. This assay is applied to E3 ligase E6AP, as its AZUL domain induces a disorder-to-order switch in an intrinsically disordered region of the proteasome, the so-named Rpn10 AZUL-binding domain (RAZUL). By testing various fluorophores, we developed an assay appropriate for high-throughput screening of Rpn10:E6AP-disrupting ligands. We found distinct positions in RAZUL for fluorophore labeling with either acrylodan or Atto610, which had disparate spectral responses to E6AP binding. E6AP caused a hypsochromic shift with increased fluorescence of acrylodan-RAZUL while decreasing fluorescence intensity of Atto610-RAZUL. Combining RAZUL labeled with either acrylodan or Atto610 into a common sample achieved robust and orthogonal measurement of the E6AP-induced conformational switch. This approach is generally applicable to disorder-to-order (or vice versa) transitions mediated by molecular interactions.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm?

During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.

Dynamics and interactions of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are widespread in eukaryotes and participate in a variety of important cellular processes. Numerous studies using state-of-the-art experimental and theoretical methods have advanced our understanding of IDPs and revealed that disordered regions engage in a large repertoire of intra- and intermolecular interactions through their conformational dynamics, thereby regulating many intracellular functions in concert with folded domains. The mechanisms by which IDPs interact with their partners are diverse, depending on their conformational propensities, and include induced fit, conformational selection, and their mixtures. In addition, IDPs are implicated in many diseases, and progress has been made in designing inhibitors of IDP-mediated interactions. Here we review these recent advances with a focus on the dynamics and interactions of IDPs.

Biomolecular dynamics in the 21st century

The relevance of motions in biological macromolecules has been clear since the early structural analyses of proteins by X-ray crystallography. Computer simulations have been applied to provide a deeper understanding of the dynamics of biological macromolecules since 1976, and are now a standard tool in many labs working on the structure and function of biomolecules. In this mini-review we highlight some areas of current interest and active development for simulations, in particular all-atom molecular dynamics simulations.

Molecular simulations integrated with experiments for probing the interaction dynamics and binding mechanisms of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) exploit their plasticity to deploy a rich panoply of soft interactions and binding phenomena. Advances in tailoring molecular simulations for IDPs combined with experimental cross-validation offer an atomistic view of the mechanisms that control IDP binding, function, and dysfunction. The emerging theme is that unbound IDPs autonomously form transient local structures and self-interactions that determine their binding behavior. Recent results have shed light on whether and how IDPs fold, stay disordered or drive condensation upon binding; how they achieve binding specificity and select among competing partners. The disorder-binding paradigm is now being proactively used by researchers to target IDPs for rational drug design and engineer molecular responsive elements for biosensing applications.

Molecular Dynamics Simulations in Protein-Protein Docking

Concerted interactions between all the cell components form the basis of biological processes. Protein-protein interactions (PPIs) constitute a tremendous part of this interaction network. Deeper insight into PPIs can help us better understand numerous diseases and lead to the development of new diagnostic and therapeutic strategies. PPI interfaces, until recently, were considered undruggable. However, it is now believed that the interfaces contain "hot spots," which could be targeted by small molecules. Such a strategy would require high-quality structural data of PPIs, which are difficult to obtain experimentally. Therefore, in silico modeling can complement or be an alternative to in vitro approaches. There are several computational methods for analyzing the structural data of the binding partners and modeling of the protein-protein dimer/oligomer structure. The major problem with in silico structure prediction of protein assemblies is obtaining sufficient sampling of protein dynamics. One of the methods that can take protein flexibility and the effects of the environment into account is Molecular Dynamics (MD). While sampling of the whole protein-protein association process with plain MD would be computationally expensive, there are several strategies to harness the method to PPI studies while maintaining reasonable use of resources. This chapter reviews known applications of MD in the PPI investigation workflows.

Markov State Model of Solvent Features Reveals Water Dynamics in Protein-Peptide Binding

In this work, we investigate the role of solvent in the binding reaction of the p53 transactivation domain (TAD) peptide to its receptor MDM2. Previously, our group generated 831 μs of explicit-solvent aggregate molecular simulation trajectory data for the MDM2-p53 peptide binding reaction using large-scale distributed computing and subsequently built a Markov State Model (MSM) of the binding reaction (Zhou et al. 2017). Here, we perform a tICA analysis and construct an MSM with similar hyperparameters while using only solvent-based structural features. We find a remarkably similar landscape but accelerated implied timescales for the slowest motions. The solvent shells contributing most to the first tICA eigenvector are those centered on Lys24 and Thr18 of the p53 TAD peptide in the range of 3-6 Å. Important solvent shells were visualized to reveal solvation and desolvation transitions along the peptide-protein binding trajectories. Our results provide a solvent-centric view of the hydrophobic effect in action for a realistic peptide-protein binding scenario.

The C-terminal tail of CSNAP attenuates the CSN complex

Protein degradation is one of the essential mechanisms that enables reshaping of the proteome landscape in response to various stimuli. The largest E3 ubiquitin ligase family that targets proteins to degradation by catalyzing ubiquitination is the cullin-RING ligases (CRLs). Many of the proteins that are regulated by CRLs are central to tumorigenesis and tumor progression, and dysregulation of the CRL family is frequently associated with cancer. The CRL family comprises ∼300 complexes, all of which are regulated by the COP9 signalosome complex (CSN). Therefore, CSN is considered an attractive target for therapeutic intervention. Research efforts for targeted CSN inhibition have been directed towards inhibition of the complex enzymatic subunit, CSN5. Here, we have taken a fresh approach focusing on CSNAP, the smallest CSN subunit. Our results show that the C-terminal region of CSNAP is tightly packed within the CSN complex, in a groove formed by CSN3 and CSN8. We show that a 16 amino acid C-terminal peptide, derived from this CSN-interacting region, can displace the endogenous CSNAP subunit from the complex. This, in turn, leads to a CSNAP null phenotype that attenuates CSN activity and consequently CRLs function. Overall, our findings emphasize the potential of a CSNAP-based peptide for CSN inhibition as a new therapeutic avenue.

Structural Space of the Duffy Antigen/Receptor for Chemokines' Intrinsically Disordered Ectodomain 1 Explored by Temperature Replica-Exchange Molecular Dynamics Simulations

Plasmodium vivax malaria affects 14 million people each year. Its invasion requires interactions between the parasitic Duffy-binding protein (PvDBP) and the N-terminal extracellular domain (ECD1) of the host's Duffy antigen/receptor for chemokines (DARC). ECD1 is highly flexible and intrinsically disordered, therefore it can adopt different conformations. We computationally modeled the challenging ECD1 local structure. With T-REMD simulations, we sampled its dynamic behavior and collected its most representative conformations. Our results suggest that most of the DARC ECD1 domain remains in a disordered state during the simulated time. Globular local conformations are found in the analyzed local free-energy minima. These globular conformations share an α-helix spanning residues Ser18 to Ser29 and in many cases they comprise an antiparallel β-sheet, whose β-strands are formed around residues Leu10 and Ala49. The formation of a parallel β-sheet is almost negligible. So far, progress in understanding the mechanisms forming the basis of the P. vivax malaria infection of reticulocytes has been hampered by experimental difficulties, along with a lack of DARC structural information. Our collection of the most probable ECD1 structural conformations will help to advance modeling of the DARC structure and to explore DARC-ECD1 interactions with a range of physiological and pathological ligands.

Model-Dependent Solvation of the K-18 Domain of the Intrinsically Disordered Protein Tau

A known imbalance between intra-protein and protein-water interactions in many empirical force fields results in collapsed conformational ensembles of intrinsically disordered proteins in explicit solvent simulations that disagree with experiments. Multiple strategies have been introduced in the literature to modify protein-water interactions, which improve agreement between experiments and simulations. In this work, we combine simulations with standard and modified force fields with a spatially resolved analysis of solvation free energy contributions and compare the consequences of each strategy. We find that enhanced Lennard-Jones (LJ) interactions between protein atoms and water oxygens primarily improve the solvation of nonpolar functional groups of the protein. In contrast, modified electrostatics in the water model or strengthened LJ interactions between the protein and water hydrogens mainly affect the hydration of polar functional groups. Modified electrostatics further impact the average orientation of water molecules in the hydration shell. As a result, protein-water interactions with the first hydration layers are strengthened, while interactions with water molecules in higher hydration shells are weakened. Hence, distinct strategies to balance intra-protein and protein-water interactions in simulations have qualitatively different effects on protein solvation. These differences are not necessarily captured by comparisons to experiments that report on global parameters describing protein conformational ensembles, e.g., the radius of gyration, but will influence the tendency of a protein to form aggregates or phase-separated droplets.

Clustering Heterogeneous Conformational Ensembles of Intrinsically Disordered Proteins with t-Distributed Stochastic Neighbor Embedding

Intrinsically disordered proteins (IDPs) populate a range of conformations that are best described by a heterogeneous ensemble. Grouping an IDP ensemble into "structurally similar" clusters for visualization, interpretation, and analysis purposes is a much-desired but formidable task, as the conformational space of IDPs is inherently high-dimensional and reduction techniques often result in ambiguous classifications. Here, we employ the t-distributed stochastic neighbor embedding (t-SNE) technique to generate homogeneous clusters of IDP conformations from the full heterogeneous ensemble. We illustrate the utility of t-SNE by clustering conformations of two disordered proteins, Aβ42, and α-synuclein, in their APO states and when bound to small molecule ligands. Our results shed light on ordered substates within disordered ensembles and provide structural and mechanistic insights into binding modes that confer specificity and affinity in IDP ligand binding. t-SNE projections preserve the local neighborhood information, provide interpretable visualizations of the conformational heterogeneity within each ensemble, and enable the quantification of cluster populations and their relative shifts upon ligand binding. Our approach provides a new framework for detailed investigations of the thermodynamics and kinetics of IDP ligand binding and will aid rational drug design for IDPs.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long-time scale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We do not, however, find evidence for the existence of canonical conformational selection or induced fit binding pathways. We observe four kinetically separated native-like bound states that interconvert on time scales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

Binding-and-Folding Recognition of an Intrinsically Disordered Protein Using Online Learning Molecular Dynamics

Intrinsically disordered proteins participate in many biological processes by folding upon binding to other proteins. However, coupled folding and binding processes are not well understood from an atomistic point of view. One of the main questions is whether folding occurs prior to or after binding. Here we use a novel, unbiased, high-throughput adaptive sampling approach to reconstruct the binding and folding between the disordered transactivation domain of c-Myb and the KIX domain of the CREB-binding protein. The reconstructed long-term dynamical process highlights the binding of a short stretch of amino acids on c-Myb as a folded α-helix. Leucine residues, especially Leu298-Leu302, establish initial native contacts that prime the binding and folding of the rest of the peptide, with a mixture of conformational selection on the N-terminal region with an induced fit of the C-terminal.

Functional regulation of an intrinsically disordered protein via a conformationally excited state

A longstanding goal in the field of intrinsically disordered proteins (IDPs) is to characterize their structural heterogeneity and pinpoint the role of this heterogeneity in IDP function. Here, we use multinuclear chemical exchange saturation (CEST) nuclear magnetic resonance to determine the structure of a thermally accessible globally folded excited state in equilibrium with the intrinsically disordered native ensemble of a bacterial transcriptional regulator CytR. We further provide evidence from double resonance CEST experiments that the excited state, which structurally resembles the DNA-bound form of cytidine repressor (CytR), recognizes DNA by means of a "folding-before-binding" conformational selection pathway. The disorder-to-order regulatory switch in DNA recognition by natively disordered CytR therefore operates through a dynamical variant of the lock-and-key mechanism where the structurally complementary conformation is transiently accessed via thermal fluctuations.

Folding of Staphylococcal Nuclease Induced by Binding of Chemically Modified Substrate Analogues Sheds Light on Mechanisms of Coupled Folding/Binding Reactions

Several proteins have been shown to undergo a shift in the mechanism of ligand binding-induced folding from conformational selection (CS; folding precedes binding) to induced fit (IF; binding precedes folding) with increasing ligand concentration. In previous studies of the coupled folding/binding reaction of staphylococcal nuclease (SNase) in the presence of a substrate analogue, adenosine-3',5'-diphosphate (prAp), we found that the two phosphate groups make important energetic contributions toward stabilizing its complex with the native protein as well as transient conformational states encountered at high ligand concentrations favoring IF. However, the structural contributions of each phosphate group during the reaction remain unclear. To address this question, we relied on fluorescence, nuclear magnetic resonance (NMR), absorption, and isothermal titration calorimetry to study the effects of deletion of the phosphate groups of prAp on the kinetics of ligand-induced folding, using a strategy analogous to mutational ϕ-value analysis to interpret the results. Kinetic measurements over a wide range of ligand concentrations, together with structural characterization of a transient protein-ligand encounter complex using 2D NMR, indicated that, at high ligand concentrations favoring IF, (i) the 5'-phosphate group interacts weakly with denatured SNase during early stages of the reaction, resulting in loose docking of the two domains of SNase, and (ii) the 3'-phosphate group engages in some specific contacts with the polypeptide in the transition state prior to formation of the native SNase-prAp complex.

Uncovering the Contributions of Charge Regulation to the Stability of Single Alpha Helices

The single alpha helix (SAH) is a recurring motif in biology. The consensus sequence has a di-block architecture that includes repeats of four consecutive glutamate residues followed by four consecutive lysine residues. Measurements show that the overall helicity of sequences with consensus E4 K4 repeats is insensitive to a wide range of pH values. Here, we use the recently introduced q-canonical ensemble, which allows us to decouple measurements of charge state and conformation, to explain the observed insensitivity of SAH helicity to pH. We couple the outputs from separate measurements of charge and conformation with atomistic simulations to derive residue-specific quantifications of preferences for being in an alpha helix and for the ionizable residues to be charged vs. uncharged. We find a clear preference for accommodating uncharged Glu residues within internal positions of SAH-forming sequences. The stabilities of alpha helical conformations increase with the number of E4 K4 repeats and so do the numbers of accessible charge states that are compatible with forming conformations of high helical content. There is conformational buffering whereby charge state heterogeneity buffers against large-scale conformational changes thus making the overall helicity insensitive to large changes in pH. Further, the results clearly argue against a single, rod-like alpha helical conformation being the only or even dominant conformation in the ensembles of so-called SAH sequences.

Functional advantages of building nanosystems using multiple molecular components

Over half of all the natural nanomachines in living organisms are multimeric and likely exploit the self-assembly of their components to provide functional benefits. However, the advantages and disadvantages of building nanosystems using multiple molecular components remain relatively unexplored at the thermodynamic, kinetic and functional levels. In this study we used theory and a simple DNA-based model that forms the same nanostructures with different numbers of components to advance our knowledge in this area. Despite its lower assembly rate, we found that a system built with three components may undergo a more cooperative assembly transition from less preorganized components, which facilitates the emergence of functionalities. Using simple variations of its components, we also found that trimeric nanosystems display a much higher level of programmability than their dimeric counterparts because they can assemble with various levels of cooperativity, self-inhibition and time-dependent properties. We show here how two simple strategies (for example, cutting and adding components) can be employed to efficiently programme the regulatory function of a more complex, artificially selected, RNA-cleaving catalytic nanosystem.

Predicting Biomolecular Binding Kinetics: A Review

Biomolecular binding kinetics including the association (k_on) and dissociation (k_off) rates are critical parameters for therapeutic design of small-molecule drugs, peptides, and antibodies. Notably, the drug molecule residence time or dissociation rate has been shown to correlate with their efficacies better than binding affinities. A wide range of modeling approaches including quantitative structure-kinetic relationship models, Molecular Dynamics simulations, enhanced sampling, and Machine Learning has been developed to explore biomolecular binding and dissociation mechanisms and predict binding kinetic rates. Here, we review recent advances in computational modeling of biomolecular binding kinetics, with an outlook for future improvements.

Biophysical and Integrative Characterization of Protein Intrinsic Disorder as a Prime Target for Drug Discovery

Protein intrinsic disorder is increasingly recognized for its biological and disease-driven functions. However, it represents significant challenges for biophysical studies due to its high conformational flexibility. In addressing these challenges, we highlight the complementary and distinct capabilities of a range of experimental and computational methods and further describe integrative strategies available for combining these techniques. Integrative biophysics methods provide valuable insights into the sequence–structure–function relationship of disordered proteins, setting the stage for protein intrinsic disorder to become a promising target for drug discovery. Finally, we briefly summarize recent advances in the development of new small molecule inhibitors targeting the disordered N-terminal domains of three vital transcription factors.

Coupled folding-upon-binding of human tumor suppressor MIG6 to lung cancer EGFR kinase domain and molecular trimming/stapling of MIG6-derived β-hairpins to target the coupling event

Human epidermal growth factor receptor (EGFR) is involved in strong association with malignant proliferation, which has been shown to play a central role in the development and progression of non-small cell lung cancer and other solid tumors. The tumor-suppressor protein MIG6 is a negative regulator of EGFR kinase activity by binding at the activation interface of asymmetric dimer of EGFR kinase domain to disrupt EGFR dimerization and then inactivate the kinase. The protein adopts two discrete fragments 1 and 2 to directly interact with EGFR. It is revealed that the MIG6 fragment 2 is intrinsically disordered in free unbound state, but would fold into a well-structured β-hairpin when binding to EGFR, thus characterized by a so-called coupled folding-upon-binding process, which can be regarded as a compromise between favorable direct readout and unfavorable indirect readout. Here, a 23-mer F2P peptide was derived from MIG6 fragment 2, trimmed into a 17-mer tF2P peptide that contains the binding hotspot region of the fragment 2, and then constrained with an ordered hairpin conformation in free unbound state by disulfide stapling, finally resulting in a rationally stapled/trimmed stF2P peptide that largely minimizes the unfavorable indirect readout effect upon its binding to EGFR kinase domain, with affinity improved considerably upon the trimming and stapling/trimming. These rationally designed β-hairpin peptides may be further exploited as potent anti-lung cancer agents to target the activation event of EGFR dimerization.

Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

The binding affinity of human pediatric respiratory syncytial virus Phosphoprotein's C-terminal tail to nucleocapsid can be improved by a rationally designed halogen-bonded system

Human respiratory syncytial virus (hRSV) is a common contagious virus that causes infections of pediatric pneumonia and specifically impacts infants and small children. The hRSV phosphoprotein is a key component of the viral RNA polymerase, which can interact with nucleocapsid and other partners through its C-terminal tail (CTT) to promote the formation of viral transcriptase complex, where the Phe241 is a key anchor residue. Based on the crystal template-modeled complex structure of hRSV nucleocapsid with a peptidic segment derived from the phosphoprotein's CTT, we successfully introduced a rationally designed halogen-bonded system to the complex interface by substituting para (p)-position of the side-chain phenyl moiety of CTT Phe241 residue with a halogen atom X (X = F, Cl, Br or I). The halogen-bonded system consists of a halogen bond (X-bond) between nucleocapsid Ser131 residue and CTT Phe241 residue as well as a hydrogen bond (H-bond) between nucleocapsid Ser131 residue and nucleocapsid Glu128 residue; the X-bond and H-bond share a common hydroxyl group of nucleocapsid Ser131 residue. High-level theoretical calculations suggested that bromine Br is the best choice that can render strong potency for the X-bond and can confer high affinity to the nucleocapsid-CTT binding. Affinity analysis revealed that the p-brominated CTT ([p]bCTT) exhibited 6.3-fold affinity improvement relative to its nonhalogenated counterpart. In contrast, the Br-substitutions at ortho (o)- and meta (m)-positions, which resulted in two negative controls of o-brominated [o]bCTT and m-brominated [m]bCTT, respectively, were unable to form effective X-bond with nucleocapsid according to theoretical investigation and did not improve the binding affinity essentially relative to native CTT.

Structural modeling, energetic analysis and molecular design of a π-stacking system at the complex interface of pediatric respiratory syncytial virus nucleocapsid with the C-terminal peptide of phosphoprotein

Human respiratory syncytial virus (RSV) is a primary cause of lower respiratory tract infections and hospital visits during infancy and childhood. The RSV phosphoprotein (P) is a major polymerase cofactor that interacts with nucleoprotein (N) to promote the recognition of ribonucleoprotein complex (RNP) by viral RNA polymerase. The binding pocket of N protein is chemically diverse, in or around which a number of aromatic and charged amino acid residues are observed. Previously, a nonapeptide segment (P peptide, ²³³DNDLSLEDF²⁴¹) representing the C-terminal tail of P protein was identified to mediate the N-P interaction with a moderate affinity, in which the Phe241 at the end of P's C-terminus plays a critical role in the binding of P peptide to N protein. Here, we found that the side-chain aromatic phenyl moiety of P Phe241 residue can form short- and long-range cation-π interactions with N Arg132 and Arg150 residues, respectively, as well as T-shaped and parallel-displaced π-π stackings with N Tyr135 and His151 residues, respectively, which co-define a geometrically satisfactory π-stacking system at the complex interface of N protein with P peptide, thus largely stabilizing the complex architecture. The stacking effect was further optimized by systematically mutating the P Phe241 residue to other natural and non-natural aromatic amino acids with diverse chemical substitutions at the phenyl moiety to examine their structural and energetic effects on π-stacking system and on protein-peptide binding. The electron-donating mutations at the phenyl moiety of P Phe241 residue can effectively enhance the π-stacking system and then promote peptide binding, whereas the bulky and positively charged mutations would considerably impair the peptide potency by introducing steric hindrance and electrostatic repulsion. The [Tyr]P, [Thp]P and [Fph]P mutants were determined to have an increased affinity relative to wild-type P peptide, which could be used as self-inhibitory peptides to competitively disrupt the native interaction between N and P proteins.

The biophysics of disordered proteins from the point of view of single-molecule fluorescence spectroscopy

Intrinsically disordered proteins (IDPs) and regions (IDRs) have emerged as key players across many biological functions and diseases. Differently from structured proteins, disordered proteins lack stable structure and are particularly sensitive to changes in the surrounding environment. Investigation of disordered ensembles requires new approaches and concepts for quantifying conformations, dynamics, and interactions. Here, we provide a short description of the fundamental biophysical properties of disordered proteins as understood through the lens of single-molecule fluorescence observations. Single-molecule Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) provides an extensive and versatile toolbox for quantifying the characteristics of conformational distributions and the dynamics of disordered proteins across many different solution conditions, both in vitro and in living cells.

Spontaneous binding of single-stranded RNAs to RRM proteins visualized by unbiased atomistic simulations with a rescaled RNA force field

Recognition of single-stranded RNA (ssRNA) by RNA recognition motif (RRM) domains is an important class of protein-RNA interactions. Many such complexes were characterized using nuclear magnetic resonance (NMR) and/or X-ray crystallography techniques, revealing ensemble-averaged pictures of the bound states. However, it is becoming widely accepted that better understanding of protein-RNA interactions would be obtained from ensemble descriptions. Indeed, earlier molecular dynamics simulations of bound states indicated visible dynamics at the RNA-RRM interfaces. Here, we report the first atomistic simulation study of spontaneous binding of short RNA sequences to RRM domains of HuR and SRSF1 proteins. Using a millisecond-scale aggregate ensemble of unbiased simulations, we were able to observe a few dozen binding events. HuR RRM3 utilizes a pre-binding state to navigate the RNA sequence to its partially disordered bound state and then to dynamically scan its different binding registers. SRSF1 RRM2 binding is more straightforward but still multiple-pathway. The present study necessitated development of a goal-specific force field modification, scaling down the intramolecular van der Waals interactions of the RNA which also improves description of the RNA-RRM bound state. Our study opens up a new avenue for large-scale atomistic investigations of binding landscapes of protein-RNA complexes, and future perspectives of such research are discussed.

Intrinsically Disordered Proteins: An Overview

Many proteins and protein segments cannot attain a single stable three-dimensional structure under physiological conditions; instead, they adopt multiple interconverting conformational states. Such intrinsically disordered proteins or protein segments are highly abundant across proteomes, and are involved in various effector functions. This review focuses on different aspects of disordered proteins and disordered protein regions, which form the basis of the so-called "Disorder-function paradigm" of proteins. Additionally, various experimental approaches and computational tools used for characterizing disordered regions in proteins are discussed. Finally, the role of disordered proteins in diseases and their utility as potential drug targets are explored.

Anchor Residues Govern Binding and Folding of an Intrinsically Disordered Domain

Minimal protein mimics have yielded novel classes of protein-protein interaction inhibitors; however, this success has not been extended to targeting intrinsically disordered proteins, which represent a significant proportion of important therapeutic targets. We sought to determine the requirements for binding an intrinsically disordered region (IDR) by its native binding partner as a prelude to developing minimal protein mimics that regulate IDR interactions. Our analysis reinforces the hypothesis that IDRs reside on a fulcrum between unfolded and folded states and that a handful of key binding residues on partner protein surfaces dictate their folding. Our studies also suggest that minimal mimics of protein surfaces may not offer specific ligands for IDRs and that it would be more judicious to target the globular protein partners of IDRs.

A litmus test for classifying recognition mechanisms of transiently binding proteins

Partner recognition in protein binding is critical for all biological functions, and yet, delineating its mechanism is challenging, especially when recognition happens within microseconds. We present a theoretical and experimental framework based on straight-forward nuclear magnetic resonance relaxation dispersion measurements to investigate protein binding mechanisms on sub-millisecond timescales, which are beyond the reach of standard rapid-mixing experiments. This framework predicts that conformational selection prevails on ubiquitin’s paradigmatic interaction with an SH3 (Src-homology 3) domain. By contrast, the SH3 domain recognizes ubiquitin in a two-state binding process. Subsequent molecular dynamics simulations and Markov state modeling reveal that the ubiquitin conformation selected for binding exhibits a characteristically extended C-terminus. Our framework is robust and expandable for implementation in other binding scenarios with the potential to show that conformational selection might be the design principle of the hubs in protein interaction networks.

The authors provide a litmus test for the recognition mechanism of transiently binding proteins based on nuclear magnetic resonance and find a conformational selection binding mechanism through concentration-dependent kinetics of ubiquitin and SH3.

A Small Molecule Stabilizes the Disordered Native State of the Alzheimer's Aβ Peptide

The stabilization of native states of proteins is a powerful drug discovery strategy. It is still unclear, however, whether this approach can be applied to intrinsically disordered proteins. Here, we report a small molecule that stabilizes the native state of the Aβ42 peptide, an intrinsically disordered protein fragment associated with Alzheimer's disease. We show that this stabilization takes place by a disordered binding mechanism, in which both the small molecule and the Aβ42 peptide remain disordered. This disordered binding mechanism involves enthalpically favorable local π-stacking interactions coupled with entropically advantageous global effects. These results indicate that small molecules can stabilize disordered proteins in their native states through transient non-specific interactions that provide enthalpic gain while simultaneously increasing the conformational entropy of the proteins.

Elucidating the mechanisms underlying protein conformational switching using NMR spectroscopy

How proteins switch between various ligand-free and ligand-bound structures has been a key biophysical question ever since the postulation of the Monod-Wyman-Changeux and Koshland-Nemethy-Filmer models over six decades ago. The ability of NMR spectroscopy to provide structural and kinetic information on biomolecular conformational exchange places it in a unique position as an analytical tool to interrogate the mechanisms of biological processes such as protein folding and biomolecular complex formation. In addition, recent methodological developments in the areas of saturation transfer and relaxation dispersion have expanded the scope of NMR for probing the mechanics of transitions in systems where one or more states constituting the exchange process are sparsely populated and 'invisible' in NMR spectra. In this review, we highlight some of the strategies available from NMR spectroscopy for examining the nature of multi-site conformational exchange, using five case studies that have employed NMR, either in isolation, or in conjunction with other biophysical tools.

Minimal Increments of Hydrophobic Collapse within the N-Terminus of the Neuropeptide Galanin

The neuropeptide galanin has a 35-year history as an intriguing target in drug design owing to its implication as a potential anticonvulsant and neuronal trophic factor among many other therapeutically interesting functions including analgesia and mood alteration. In this study, we report the structural characterization of three synthetic fragments of the galanin N-terminus in buffered aqueous solution: hGal(2-12)KK, hGal(1-12)KK, and hGal(1-17)KK. High-field two-dimensional ¹H-¹H nuclear magnetic resonance (NMR) data were acquired for these fragments and used to derive distance restraints. We further utilized modified hydrogen bonding and dihedral restraints to reflect chemical shift patterns in the data, which revealed the signature of a weakly folded helix. Together, these sets of restraints were used to generate NMR structures of all three fragments, which depict a core of hydrophobic residues that cluster together regardless of the presence of a helical structure, and correspond to residues in the N-terminus of galanin that have been previously shown to be critical for binding its receptors. The helical structure only appears following the inclusion of Gly(1) in the sequence, and at longer sequence lengths, unlike many other peptides, the helix does not propagate. Rather, a few turns of poorly ordered helix appear to be a secondary consequence of clusters of hydrophobic sidechains that are conserved across all of the peptides in this study; the helices themselves appear ordered as a consequence of this clustering, and these clusters compare directly to those observed recently to make contacts between galanin and two of its receptor subtypes. Collapsed hydrophobic residues therefore organize and compose the functional core of human galanin and raise interesting questions about the nature of the conformational order in ligands that bind cell surface receptors.

Energetics of a protein disorder-order transition in small molecule recognition

Many proteins recognise other proteins via mechanisms that involve the folding of intrinsically disordered regions upon complex formation. Here we investigate how the selectivity of a drug-like small molecule arises from its modulation of a protein disorder-to-order transition. Binding of the compound AM-7209 has been reported to confer order upon an intrinsically disordered ‘lid’ region of the oncoprotein MDM2. Calorimetric measurements revealed that truncation of the lid region of MDM2 increases the apparent dissociation constant of AM-7209 250-fold. By contrast, lid truncation has little effect on the binding of the ligand Nutlin-3a. Insights into these differential binding energetics were obtained via a complete thermodynamic analysis that featured adaptive absolute alchemical free energy of binding calculations with enhanced-sampling molecular dynamics simulations. The simulations reveal that in apo MDM2 the ordered lid state is energetically disfavoured. AM-7209, but not Nutlin-3a, shows a significant energetic preference for ordered lid conformations, thus shifting the balance towards ordering of the lid in the AM-7209/MDM2 complex. The methodology reported herein should facilitate broader targeting of intrinsically disordered regions in medicinal chemistry.

Molecular simulations and biophysical measurements elucidate why the ligand AM-7209 orders a disordered region of the protein MDM2 on binding. This work expands strategies available to medicinal chemists for targeting disordered proteins.

Residues of the minimal sequence of G0S2 collectively contribute to ATGL inhibition while C-and N-terminal extensions promote binding to ATGL

The protein encoded by the G0/G1 switch gene 2 (G0S2) is a potent inhibitor of adipose triglyceride lipase (ATGL) and thus an important regulator of intracellular lipolysis. Since dysfunction of lipolysis is associated with metabolic diseases including diabetes and obesity, inhibition of ATGL is considered a therapeutic strategy. G0S2 interacts with ATGL's patatin-domain to mediate non-competitive inhibition, however atomic details of the inhibition mechanism are incompletely understood. Sequences of G0S2 from higher organisms show a highly conserved N-terminal part, including a hydrophobic region covering amino acids 27 to 42. We show that predicted G0S2 orthologs from platypus, chicken and Japanese rice-fish are able to inhibit human and mouse ATGL, emphasizing the contribution of conserved amino acid to ATGL inhibition. Our site directed mutagenesis and truncation studies give insights in the protein-protein interaction on a per-residue level. We determine that the minimal sequence required for ATGL inhibition ranges from amino acids 20 to 44. Residues Y27, V28, G30, A34 G37, V39 or L42 within this sequence play a substantial role in ATGL inhibition. Furthermore, we show that unspecific interactions of the N-terminal part (amino acids 20-27) of the minimal sequence facilitate the interaction to ATGL. Our studies also demonstrate that full-length G0S2 shows higher tolerance to specific single amino acid exchanges in the hydrophobic region due to the stronger contributions of unspecific interactions. However, exchanges of more than one amino-acid in the hydrophobic region also result in the loss of function as ATGL inhibitor even in the full-length protein.

Molecular Basis of Small-Molecule Binding to α-Synuclein

Intrinsically disordered proteins (IDPs) are implicated in many human diseases. They have generally not been amenable to conventional structure-based drug design, however, because their intrinsic conformational variability has precluded an atomic-level understanding of their binding to small molecules. Here we present long-time-scale, atomic-level molecular dynamics (MD) simulations of monomeric α-synuclein (an IDP whose aggregation is associated with Parkinson’s disease) binding the small-molecule drug fasudil in which the observed protein–ligand interactions were found to be in good agreement with previously reported NMR chemical shift data. In our simulations, fasudil, when bound, favored certain charge–charge and π-stacking interactions near the C terminus of α-synuclein but tended not to form these interactions simultaneously, rather breaking one of these interactions and forming another nearby (a mechanism we term dynamic shuttling). Further simulations with small molecules chosen to modify these interactions yielded binding affinities and key structural features of binding consistent with subsequent NMR experiments, suggesting the potential for MD-based strategies to facilitate the rational design of small molecules that bind with disordered proteins.

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives

Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.

Transient exposure of a buried phosphorylation site in an autoinhibited protein

Autoinhibition is a mechanism used to regulate protein function, often by making functional sites inaccessible through the interaction with a cis-acting inhibitory domain. Such autoinhibitory domains often display a substantial degree of structural disorder when unbound, and only become structured in the inhibited state. These conformational dynamics make it difficult to study the structural origin of regulation, including effects of regulatory post-translational modifications. Here, we study the autoinhibition of the Dbl Homology domain in the protein Vav1 by the so-called acidic inhibitory domain. We use molecular simulations to study the process by which a mostly unstructured inhibitory domain folds upon binding and how transient exposure of a key buried tyrosine residue makes it accessible for phosphorylation. We show that the inhibitory domain, which forms a helix in the bound and inhibited stated, samples helical structures already before binding and that binding occurs via a molten-globule-like intermediate state. Together, our results shed light on key interactions that enable the inhibitory domain to sample a finely tuned equilibrium between an inhibited and a kinase-accessible state.