Slow Conformational Dynamics Research Articles

Massively parallel free energy calculations for in silico affinity maturation of designed miniproteins.

Computational protein design efforts continue to make remarkable advances, yet the discovery of high-affinity binders typically requires large-scale experimental screening of site-saturated mutant (SSM) libraries. Here, we explore how massively parallel free energy methods can be used for in silico affinity maturation of de novo designed binding proteins. Using an expanded ensemble (EE) approach, we perform exhaustive relative binding free energy calculations for SSM variants of three miniproteins designed to bind influenza A H1 hemagglutinin by Chevalier et al. (2017). We compare our predictions to experimental ΔΔ G values inferred from a Bayesian analysis of the high-throughput sequencing data, and to state-of-the-art predictions made using the Flex ddG Rosetta protocol. A systematic comparison reveals prediction accuracies around 2 kcal/mol, and identifies net charge changes, large numbers of alchemical atoms, and slow side chain conformational dynamics as key contributors to the uncertainty of the EE predictions. Flex ddG predictions are more accurate on average, but highly conservative. In contrast, EE predictions can better classify stabilizing and destabilizing mutations. We also explored the ability of SSM scans to rationalize known affinity-matured variants containing multiple mutations, which are non-additive due to epistatic effects. Simple electrostatic models fail to explain non-additivity, but observed mutations are found at positions with higher Shannon entropies. Overall, this work suggests that simulation-based free energy methods can provide predictive information for in silico affinity maturation of designed miniproteins, with many feasible improvements to the efficiency and accuracy within reach.

Evaluating Polarizable Biomembrane Simulations against Experiments.

Owing to the increase of available computational capabilities and the potential for providing a more accurate description, polarizable molecular dynamics force fields are gaining popularity in modeling biomolecular systems. It is, however, crucial to evaluate how much precision is truly gained with increasing cost and complexity of the simulation. Here, we leverage the NMRlipids open collaboration and Databank to assess the performance of available polarizable lipid models─the CHARMM-Drude and the AMOEBA-based parameters─against high-fidelity experimental data and compare them to the top-performing nonpolarizable models. While some improvement in the description of ion binding to membranes is observed in the most recent CHARMM-Drude parameters, and the conformational dynamics of AMOEBA-based parameters are excellent, the best nonpolarizable models tend to outperform their polarizable counterparts for each property we explored. The identified shortcomings range from inaccuracies in describing the conformational space of lipids to excessively slow conformational dynamics. Our results provide valuable insights for the further refinement of polarizable lipid force fields and for selecting the best simulation parameters for specific applications.

Acceleration of generalized replica exchange with solute tempering simulations of large biological systems on massively parallel supercomputer.

Generalized replica exchange with solute tempering (gREST) is one of the enhanced sampling algorithms for proteins or other systems with rugged energy landscapes. Unlike the replica-exchange molecular dynamics (REMD) method, solvent temperatures are the same in all replicas, while solute temperatures are different and are exchanged frequently between replicas for exploring various solute structures. Here, we apply the gREST scheme to large biological systems containing over one million atoms using a large number of processors in a supercomputer. First, communication time on a multi-dimensional torus network is reduced by matching each replica to MPI processors optimally. This is applicable not only to gREST but also to other multi-copy algorithms. Second, energy evaluations, which are necessary for the multistate bennet acceptance ratio (MBAR) method for free energy estimations, are performed on-the-fly during the gREST simulations. Using these two advanced schemes, we observed 57.72 ns/day performance in 128-replica gREST calculations with 1.5 million atoms system using 16,384 nodes in Fugaku. These schemes implemented in the latest version of GENESIS software could open new possibilities to answer unresolved questions on large biomolecular complex systems with slow conformational dynamics.

Slow conformational changes in the rigid and highly stable chymotrypsin inhibitor 2.

Slow conformational changes are often directly linked to protein function. It is however less clear how such processes may perturb the overall folding stability of a protein. We previously found that the stabilizing double mutant L49I/I57V in the small protein chymotrypsin inhibitor 2 from barley led to distributed increased nanosecond and faster dynamics. Here we asked what effects the L49I and I57V substitutions, either individually or together, have on the slow conformational dynamics of CI2. We used 15 N CPMG spin relaxation dispersion experiments to measure the kinetics, thermodynamics and structural changes associated with slow conformational change in CI2. These changes result in an excited state that is populated to 4.3% at 1 °C. As the temperature is increased the population of the excited state decreases. Structural changes in the excited state are associated with residues that interact with water molecules that have well defined positions and are found at these positions in all crystal structures of CI2. The substitutions in CI2 have only little effect on the structure of the excited state whereas the stability of the excited state to some extent follows the stability of the main state. The minor state is thus most populated for the most stable CI2 variant and least populated for the least stable variant. We hypothesize that the interactions between the substituted residues and the well-ordered water molecules links subtle structural changes around the substituted residues to the region in the protein that experience slow conformational changes. This article is protected by copyright. All rights reserved.

Contryphan sequence diversity: Messy N-terminus processing, effects on chromatographic behaviour and mass spectrometric fragmentation

Contryphans, peptides containing a single disulfide bond, are found abundantly in cone snail venom. The analysis of a large dataset of available contryphan sequences permits a classification based on the occurrence of proline residues at positions 2 and 5 within the macrocyclic 23-membered disulfide loop. Further sequence diversity is generated by variable proteolytic processing of the contryphan precursor proteins. In the majority of contryphans, presence of Pro at position 2 and a D-residue at position 3 leads to a slow conformational dynamics, manifesting as anomalous chromatographic profiles during LC analysis. LC-MS analysis of diverse contryphans suggests that elution profiles may be used as a rapid diagnostic for the presence of the Pro2-DXxx3 motif. Natural sequences from C.inscriptus and C.frigidus together with synthetic analogs permit the delineation of the features necessary for abnormal chromatographic behaviour. A diagnostic for the presence of Pro at position 5 is obtained by the observation of non-canonical fragment ions, generated by N-Cα bond cleavage at the dehydroalanine residue formed by disulfide cleavage. Anomalous LC profiles supports Pro at position 2, while non-canonical mass spectral fragments established Pro at position 5, providing a rapid method for contryphan analysis from LC-ESI-MS/MS profiles of crude Conus venom. SignificanceContryphans are peptides, widely distributed in cone snail venom, which display extensive sequence diversity. Heterogeneity of proteolytic processing of contryphan precursor proteins, together with post-translational modifications contributes to contryphan diversity. Contryphans, identified by a combination of mass spectrometry and transcriptomic analysis, are classified on the basis of sequence features, primarily the number of proline residues within the disulfide loop. Conformational diversity arises in contryphans by cis-trans isomerization of Cys-Pro bonds, resulting in characteristic chromatographic profiles, permitting identification even in crude venom mixtures. Rapid identification of contryphans in cone snail peptide libraries is also facilitated by diagnostic mass spectral fragments arising by non-canonical cleavage of the N-Cα bond at Cys(7).

Different Forms of Disorder in NMDA-Sensitive Glutamate Receptor Cytoplasmic Domains Are Associated with Differences in Condensate Formation.

The N-methyl-D-aspartate (NMDA)-sensitive glutamate receptor (NMDAR) helps assemble downstream signaling pathways through protein interactions within the postsynaptic density (PSD), which are mediated by its intracellular C-terminal domain (CTD). The most abundant NMDAR subunits in the brain are GluN2A and GluN2B, which are associated with a developmental switch in NMDAR composition. Previously, we used single molecule fluorescence resonance energy transfer (smFRET) to show that the GluN2B CTD contained an intrinsically disordered region with slow, hop-like conformational dynamics. The CTD from GluN2B also undergoes liquid-liquid phase separation (LLPS) with synaptic proteins. Here, we extend these observations to the GluN2A CTD. Sequence analysis showed that both subunits contain a form of intrinsic disorder classified as weak polyampholytes. However, only GluN2B contained matched patterning of arginine and aromatic residues, which are linked to LLPS. To examine the conformational distribution, we used discrete molecular dynamics (DMD), which revealed that GluN2A favors extended disordered states containing secondary structures while GluN2B favors disordered globular states. In contrast to GluN2B, smFRET measurements found that GluN2A lacked slow conformational dynamics. Thus, simulation and experiments found differences in the form of disorder. To understand how this affects protein interactions, we compared the ability of these two NMDAR isoforms to undergo LLPS. We found that GluN2B readily formed condensates with PSD-95 and SynGAP, while GluN2A failed to support LLPS and instead showed a propensity for colloidal aggregation. That GluN2A fails to support this same condensate formation suggests a developmental switch in LLPS propensity.

Slow conformational dynamics of the human A2A adenosine receptor are temporally ordered

A more complete depiction of protein energy landscapes includes the identification of different function-related conformational states and the determination of the pathways connecting them. We used total internal reflection fluorescence (TIRF) imaging to investigate the conformational dynamics of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), at the single-molecule level. Slow, reversible conformational exchange was observed among three different fluorescence emission states populated for agonist-bound A2AAR. Transitions among these states predominantly occurred in a specific order, and exchange between inactive and active-like conformations proceeded through an intermediate state. Models derived from molecular dynamics simulations with available A2AAR structures rationalized the relative fluorescence emission intensities for the highest and lowest emission states but not the transition state. This suggests that the functionally critical intermediate state required to achieve activation is not currently visualized among available A2AAR structures.

Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions.

Enzymes employ a wide range of protein motions to achieve efficient catalysis of chemical reactions. While the role of collective protein motions in substrate binding, product release, and regulation of enzymatic activity is generally understood, their roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations, enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the roles of catalytic residues in catalysis and the conformational change in the enzyme. This study reveals that the motions in the active site, which occur on a time scale of picoseconds to nanoseconds, link the catalytic reaction to the slow conformational dynamics of the enzyme by modulating the free energy landscapes of subdomain motions. In particular, substantial conformational rearrangement occurs in the active site following the catalytic reaction. This rearrangement not only affects the reaction barrier but also promotes a more open conformation of the enzyme after the reaction, which then results in an accelerated opening of the enzyme compared to that of the reactant state. The results illustrate a linkage between enzymatic catalysis and collective protein motions, whereby the disparate time scales between the two processes are bridged by a cascade of intermediate-scale motion of catalytic residues modulating the free energy landscapes of the catalytic and conformational change processes.

Ligand modulation of the conformational dynamics of the A2A adenosine receptor revealed by single-molecule fluorescence

G protein-coupled receptors (GPCRs) are the largest class of transmembrane proteins, making them an important target for therapeutics. Activation of these receptors is modulated by orthosteric ligands, which stabilize one or several states within a complex conformational ensemble. The intra- and inter-state dynamics, however, is not well documented. Here, we used single-molecule fluorescence to measure ligand-modulated conformational dynamics of the adenosine A2A receptor (A2AR) on nanosecond to millisecond timescales. Experiments were performed on detergent-purified A2R in either the ligand-free (apo) state, or when bound to an inverse, partial or full agonist ligand. Single-molecule Förster resonance energy transfer (smFRET) was performed on detergent-solubilized A2AR to resolve active and inactive states via the separation between transmembrane (TM) helices 4 and 6. The ligand-dependent changes of the smFRET distributions are consistent with conformational selection and with inter-state exchange lifetimes ≥ 3 ms. Local conformational dynamics around residue 2296.31 on TM6 was measured using fluorescence correlation spectroscopy (FCS), which captures dynamic quenching due to photoinduced electron transfer (PET) between a covalently-attached dye and proximal aromatic residues. Global analysis of PET-FCS data revealed fast (150–350 ns), intermediate (50–60 μs) and slow (200–300 μs) conformational dynamics in A2AR, with lifetimes and amplitudes modulated by ligands and a G-protein mimetic (mini-Gs). Most notably, the agonist binding and the coupling to mini-Gs accelerates and increases the relative contribution of the sub-microsecond phase. Molecular dynamics simulations identified three tyrosine residues (Y112, Y2887.53, and Y2907.55) as being responsible for the dynamic quenching observed by PET-FCS and revealed associated helical motions around residue 2296.31 on TM6. This study provides a quantitative description of conformational dynamics in A2AR and supports the idea that ligands bias not only GPCR conformations but also the dynamics within and between distinct conformational states of the receptor.

Conformational heterogeneity of Savinase from NMR, HDX-MS and X-ray diffraction analysis.

BackgroundSeveral examples have emerged of enzymes where slow conformational changes are of key importance for function and where low populated conformations in the resting enzyme resemble the conformations of intermediate states in the catalytic process. Previous work on the subtilisin protease, Savinase, from Bacillus lentus by NMR spectroscopy suggested that this enzyme undergoes slow conformational dynamics around the substrate binding site. However, the functional importance of such dynamics is unknown.MethodsHere we have probed the conformational heterogeneity in Savinase by following the temperature dependent chemical shift changes. In addition, we have measured changes in the local stability of the enzyme when the inhibitor phenylmethylsulfonyl fluoride is bound using hydrogen-deuterium exchange mass spectrometry (HDX-MS). Finally, we have used X-ray crystallography to compare electron densities collected at cryogenic and ambient temperatures and searched for possible low populated alternative conformations in the crystals.ResultsThe NMR temperature titration shows that Savinase is most flexible around the active site, but no distinct alternative states could be identified. The HDX shows that modification of Savinase with inhibitor has very little impact on the stability of hydrogen bonds and solvent accessibility of the backbone. The most pronounced structural heterogeneities detected in the diffraction data are limited to alternative side-chain rotamers and a short peptide segment that has an alternative main-chain conformation in the crystal at cryo conditions. Collectively, our data show that there is very little structural heterogeneity in the resting state of Savinase and hence that Savinase does not rely on conformational selection to drive the catalytic process.

Overview of the SAMPL6 host-guest binding affinity prediction challenge.

Accurately predicting the binding affinities of small organic molecules to biological macromolecules can greatly accelerate drug discovery by reducing the number of compounds that must be synthesized to realize desired potency and selectivity goals. Unfortunately, the process of assessing the accuracy of current computational approaches to affinity prediction against binding data to biological macromolecules is frustrated by several challenges, such as slow conformational dynamics, multiple titratable groups, and the lack of high-quality blinded datasets. Over the last several SAMPL blind challenge exercises, host-guest systems have emerged as a practical and effective way to circumvent these challenges in assessing the predictive performance of current-generation quantitative modeling tools, while still providing systems capable of possessing tight binding affinities. Here, we present an overview of the SAMPL6 host-guest binding affinity prediction challenge, which featured three supramolecular hosts: octa-acid (OA), the closely related tetra-endo-methyl-octa-acid (TEMOA), and cucurbit[8]uril (CB8), along with 21 small organic guest molecules. A total of 119 entries were received from ten participating groups employing a variety of methods that spanned from electronic structure and movable type calculations in implicit solvent to alchemical and potential of mean force strategies using empirical force fields with explicit solvent models. While empirical models tended to obtain better performance than first-principle methods, it was not possible to identify a single approach that consistently provided superior results across all host-guest systems and statistical metrics. Moreover, the accuracy of the methodologies generally displayed a substantial dependence on the system considered, emphasizing the need for host diversity in blind evaluations. Several entries exploited previous experimental measurements of similar host-guest systems in an effort to improve their physical-based predictions via some manner of rudimentary machine learning; while this strategy succeeded in reducing systematic errors, it did not correspond to an improvement in statistical correlation. Comparison to previous rounds of the host-guest binding free energy challenge highlights an overall improvement in the correlation obtained by the affinity predictions for OA and TEMOA systems, but a surprising lack of improvement regarding root mean square error over the past several challenge rounds. The data suggests that further refinement of force field parameters, as well as improved treatment of chemical effects (e.g., buffer salt conditions, protonation states), may be required to further enhance predictive accuracy.

QM/MM Geometry Optimization on Extensive Free-Energy Surfaces for Examination of Enzymatic Reactions and Design of Novel Functional Properties of Proteins.

Many remarkable molecular functions of proteins use their characteristic global and slow conformational dynamics through coupling of local chemical states in reaction centers with global conformational changes of proteins. To theoretically examine the functional processes of proteins in atomic detail, a methodology of quantum mechanical/molecular mechanical (QM/MM) free-energy geometry optimization is introduced. In the methodology, a geometry optimization of a local reaction center is performed with a quantum mechanical calculation on a free-energy surface constructed with conformational samples of the surrounding protein environment obtained by a molecular dynamics simulation with a molecular mechanics force field. Geometry optimizations on extensive free-energy surfaces by a QM/MM reweighting free-energy self-consistent field method designed to be variationally consistent and computationally efficient have enabled examinations of the multiscale molecular coupling of local chemical states with global protein conformational changes in functional processes and analysis and design of protein mutants with novel functional properties.

Adding a Structural Context to the Deprotometalation and Trans-Metal Trapping Chemistry of Phenyl-Substituted Benzotriazole.

Organometallic bases are becoming increasingly complex, because mixing components can lead to bases superior to single-component bases. To better understand this superiority, it is useful to study metalated intermediate structures prior to quenching. This study is on 1-phenyl-1H-benzotriazole, which was previously deprotonated by an in situ ZnCl2 ⋅TMEDA/LiTMP (TMEDA=N,N,N',N'-tetramethylethylenediamine; TMP=2,2,6,6-tetramethylpiperidide) mixture and then iodinated. Herein, reaction with LiTMP exposes the deficiency of the single-component base as the crystalline product obtained was [{4-R-1-(2-lithiophenyl)-1H-benzotriazole⋅3THF}2 ], [R=2-C6 H4 (Ph)NLi], in which ring opening of benzotriazole and N2 extrusion had occurred. Supporting lithiation by adding iBu2 Al(TMP) induces trans-metal trapping, in which C-Li bonds transform into C-Al bonds to stabilise the metalated intermediate. X-ray diffraction studies revealed homodimeric [(4-R'-1-phenyl-1H-benzotriazole)2 ], [R'=(iBu)2 Al(μ-TMP)Li], and its heterodimeric isomer [(4-R'-1-phenyl-1H-benzotriazole){2-R'-1-phenyl-1H-benzotriazole}], whose structure and slow conformational dynamics were probed by solution NMR spectroscopy.

Computationally efficient prediction of area per lipid

Area per lipid (APL) is an important property of biological and artificial membranes. Newly constructed bilayers are characterized by their APL and newly elaborated force fields must reproduce APL. Computer simulations of APL are very expensive due to slow conformational dynamics. The simulated dynamics increases exponentially with respect to temperature. APL dependence on temperature is linear over an entire temperature range. I provide numerical evidence that thermal expansion coefficient of a lipid bilayer can be computed at elevated temperatures and extrapolated to the temperature of interest. Thus, sampling times to predict accurate APL are reduced by a factor of ∼10.

Linking Intrinsic Disorder to Allosteric Regulation in the NMDA Receptor

Ionotropic glutamate receptors mediate excitatory synaptic transmission in the central nervous system. Ligand binding to the extracellular domain triggers opening of a transmembrane ion channel domain allowing sodium influx. The NMDA-sensitive receptor subtype also permits calcium influx, which plays important roles in regulating synaptic strength. With the evolution of vertebrates, the NMDA receptor isoform GluN2B acquired a large cytoplasmic domain (CTD) that enables allosteric modulation by Src kinase. GluN2B-containing receptors are inhibited by extracellular zinc but the block can be reversed by phosphorylation of two tyrosine residues in the distal CTD. The mechanism whereby the CTD affects gating of the ion channel is unknown. The CTD is predicted to contain large stretches of intrinsic disorder. A central palmitoylation motif splits the CTD into two domains. Using ensemble biophysical and biochemical methods, we confirmed that the distal cytoplasmic domain (CTD2) is a disordered globule. Single molecule FRET revealed stochastic transitions that appear to arise from slow conformational dynamics on the second timescale. We found that Src phosphorylation caused a uniform expansion of the polypeptide. We were able to induce compaction in CTD2 by removing proline residues, which retarded the Src-induced expansion. Scanning mutagenesis also revealed a role for proline in aggregation as well as CTD2 binding to small molecules, like the fluorophore ANS, and a PDZ domain involved in synaptic targeting of the receptor. Thus, proline appears to modulate the attractive potential of short linear motifs within regions of intrinsic disorder. In the context of the native NMDA receptor, proline depletion near the Src phosphorylation sites uncoupled allosteric regulation of the channel by zinc and eliminated the zinc block. This suggests that allosteric modulation is linked to the conformational dynamics within the disordered cytoplasmic domain.

Role of conformational entropy in the activity and regulation of the catalytic subunit of protein kinase A

Protein kinase A (PKA) is the archetypical phosphokinase, sharing a catalytic core with the entire protein kinase superfamily. In eukaryotes, the ubiquitous location of PKA makes it one of the most important cellular signaling molecules, involved in a myriad of events. The catalytic subunit of PKA (PKA-C) is one of the most studied enzymes and was the first kinase to be crystallized; however, the effects of ligand binding, post-translational modifications and mutations on the activity of the kinase have been difficult to understand with only structural data. Here, we review our latest NMR studies on PKA-C, the results of which underscore the role of fast and slow conformational dynamics in the activation and inhibition of the kinase.

A Comparative CEST NMR Study of Slow Conformational Dynamics of Small GTPases Complexed with GTP and GTP Analogues

The Ras superfamily of small GTPases are important intracellular signaling molecules, the functions of which are determined by the binding of guanosine nucleotides (GTP = guanosine triphosphate and GDP = guanosine diphosphate). [1] The GTP-bound (“active”) states of these enzymes are capable of interacting with specific downstream effector proteins, thus eliciting a wide range of cellular responses. [2, 3] Mutations that reduce the rate of GTP hydrolysis and thus increase the lifetime of the active GTP-bound state are frequently oncogenic and contribute to the development and metastasis of human cancers. [4] Elegant 31 P NMR studies of GTP-bound Ras showed that the enzyme interconverts between two states, a minor conformer termed state 1 and a major species designated state 2. [5–9] Similar conformational dynamics have been observed in other Ras family GTPases as well. [10–12] State 2 is generally regarded as the conformation competent for binding effector proteins, whereas state 1 exhibits significantly reduced affinity for these molecules. [5–7, 13, 14] Stabilization of the low-affinity state 1 was

A Comparative CEST NMR Study of Slow Conformational Dynamics of Small GTPases Complexed with GTP and GTP Analogues

A Comparative CEST NMR Study of Slow Conformational Dynamics of Small GTPases Complexed with GTP and GTP Analogues

Precise Monitoring of Chemical Changes through Localization Analysis of Dynamic Spectra (LADS)

We present a method for monitoring subtle (sub-wavenumber) dynamics within time-varying spectra. Peak fitting is performed for large numbers of spectra in a series, allowing for monitoring time evolutions of peak positions with high precision and confidence. Sub-wavenumber peak shifts due to physical or chemical changes in the sample can be monitored and their temporal evolution characterized. In surface-enhanced Raman scattering experiments, we were able to distinguish between slow photo-damage and fast conformational change dynamics. Fluctuations in peak positions of Raman spectra recorded from a single yeast cell indicated that no significant irreversible photo-damage occurred, but these fluctuations suggest changes in the trapping conditions or biochemical changes associated with the cellular machinery in the cell. The technique is particularly suitable for applications where dynamics of spectra are of interest.

Correction to Slow Conformational Dynamics in the Cystoviral RNA-Directed RNA Polymerase P2: Influence of Substrate Nucleotides and Template RNA

Correction to Slow Conformational Dynamics in the Cystoviral RNA-Directed RNA Polymerase P2: Influence of Substrate Nucleotides and Template RNA