Effects of ligand binding upon flexibility of proteins

Enzyme dynamics: Looking beyond a single structure

Conventional understanding of how enzymes function strongly emphasizes the role of structure. However, increasing evidence clearly indicates that enzymes do not remain fixed or operate exclusively in or close to their native structure. Different parts of the enzyme (from individual residues to full domains) undergo concerted motions on a wide range of time-scales, including that of the catalyzed reaction. Information obtained on these internal motions and conformational fluctuations has so far uncovered and explained many aspects of enzyme mechanisms, which could not have been understood from a single structure alone. Although there is wide interest in understanding enzyme dynamics and its role in catalysis, several challenges remain. In addition to technical difficulties, the vast majority of investigations are performed in dilute aqueous solutions, where conditions are significantly different than the cellular milieu where a large number of enzymes operate. In this review, we discuss recent developments, several challenges as well as opportunities related to this topic. The benefits of considering dynamics as an integral part of the enzyme function can also enable new means of biocatalysis, engineering enzymes for industrial and medicinal applications.

Predicted dynamical couplings of protein residues characterize catalysis, transport and allostery

MOTIVATION

Protein function is intrinsically linked to native dynamics, but the systematic characterization of functionally relevant dynamics remains elusive besides specific examples. Here we exhaustively characterize three types of dynamical couplings between protein residues: co-directionality (moving along collinear directions), coordination (small fluctuations of the interatomic distance) and deformation (the extent by which perturbations applied at one residue modify the local structure of the other one), which we analytically compute through the torsional network model.

RESULTS

We find that ligand binding sites are characterized by large within-site coordination and co-directionality, much larger than expected for generic sets of residues with equivalent sequence distances. In addition, catalytic sites are characterized by high coordination couplings with other residues in the protein, supporting the view that the overall protein structure facilitates the catalytic dynamics. The binding sites of allosteric effectors are characterized by comparably smaller coordination and higher within-site deformation than other ligands, which supports their dynamic nature. Allosteric inhibitors are coupled to the active site more frequently through deformation than through coordination, while the contrary holds for activators. We characterize the dynamical couplings of the sodium-dependent Leucine transporter protein (LeuT). The couplings between and within sites progress consistently along the transport cycle, providing a mechanistic description of the coupling between the uptake and release of ions and substrate, and they highlight qualitative differences between the wild-type and a mutant for which chloride is necessary for transport.

AVAILABILITY AND IMPLEMENTATION

The program tnm is freely available at https://github.com/ugobas/tnm.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Distinctive communication networks in inactive states of β2 -adrenergic receptor: Mutual information and entropy transfer analysis

Mutual information and entropy transfer analysis employed on two inactive states of human beta-2 adrenergic receptor (β₂ -AR) unraveled distinct communication pathways. Previously, a so-called "highly" inactive state of the receptor was observed during 1.5 microsecond long molecular dynamics simulation where the largest intracellular loop (ICL3) was swiftly packed onto the G-protein binding cavity, becoming entirely inaccessible. Mutual information quantifying the degree of correspondence between backbone-C_α fluctuations was mostly shared between intra- and extra-cellular loop regions in the original inactive state, but shifted to entirely different regions in this latest inactive state. Interestingly, the largest amount of mutual information was always shared among the mobile regions. Irrespective of the conformational state, polar residues always contributed more to mutual information than hydrophobic residues, and also the number of polar-polar residue pairs shared the highest degree of mutual information compared to those incorporating hydrophobic residues. Entropy transfer, quantifying the correspondence between backbone-C_α fluctuations at different timesteps, revealed a distinctive pathway directed from the extracellular site toward intracellular portions in this recently exposed inactive state for which the direction of information flow was the reverse of that observed in the original inactive state where the mobile ICL3 and its intracellular surroundings drove the future fluctuations of extracellular regions.

Intrinsic Dynamics and Causality in Correlated Motions Unraveled in Two Distinct Inactive States of Human β2-Adrenergic Receptor

The alternative inactive state of the human β₂-adrenergic receptor originally exposed in molecular dynamics simulations was investigated using various analysis tools to evaluate causality between correlated residue-pair fluctuations and suggest allosteric communication pathways. A major conformational shift observed in the third intracellular loop (ICL3) displayed a novel inactive state, featuring an inaccessible G protein binding site blocked by ICL3 and an expanded orthosteric ligand binding site. Residue-based mean-square fluctuation and stiffness calculations revealed a significant mobility decrease in ICL3, which induced a mobility increase in the remaining loop regions. This indicates conformational entropy loss in one mobile region being compensated by residual intermolecular motions in other mobile regions. Moreover, the extent of significantly correlated motions decreased, and correlations that once existed between transmembrane helices shifted toward regions with increased mobility. Conditional time-delayed cross-correlation analysis identified distinct driver-follower relationship profiles. Prior to its packing, freely moving ICL3 was markedly driven by transmembrane helix-8 whereas once packed, ICL3 controlled future fluctuations of nearby helices. Moreover, two transmembrane helices, (H5 and H6), started to control future fluctuations of a remote site, the extracellular loop, ECL2. This clearly suggests that allosteric coupling between extra- and intracellular parts intensified, in agreement with the receptor's well recognized feature, which is the inverse proportionality between activity and the degree of coupling.

Intrinsic K-Ras dynamics: A novel molecular dynamics data analysis method shows causality between residue pair motions

K-Ras is the most frequently mutated oncogene in human cancers, but there are still no drugs that directly target it in the clinic. Recent studies utilizing dynamics information show promising results for selectively targeting mutant K-Ras. However, despite extensive characterization, the mechanisms by which K-Ras residue fluctuations transfer allosteric regulatory information remain unknown. Understanding the direction of information flow can provide new mechanistic insights for K-Ras targeting. Here, we present a novel approach –conditional time-delayed correlations (CTC) – using the motions of all residue pairs of a protein to predict directionality in the allosteric regulation of the protein fluctuations. Analyzing nucleotide-dependent intrinsic K-Ras motions with the new approach yields predictions that agree with the literature, showing that GTP-binding stabilizes K-Ras motions and leads to residue correlations with relatively long characteristic decay times. Furthermore, our study is the first to identify driver-follower relationships in correlated motions of K-Ras residue pairs, revealing the direction of information flow during allosteric modulation of its nucleotide-dependent intrinsic activity: active K-Ras Switch-II region motions drive Switch-I region motions, while α-helix-3L7 motions control both. Our results provide novel insights for strategies that directly target mutant K-Ras.

A new concept for molecular engineering of artificial enzymes: a multiscale simulation

We propose a new concept for the design of artificial enzymes from synthetic protein-like copolymers and non-natural functional monomers which in terms of their affinity for water can be divided into two categories: hydrophobic and hydrophilic. Hydrophilic monomers comprise catalytically active groups similar to those in the corresponding amino acid residues. A key ingredient of our approach is that the target globular conformation of protein-like, core-shell morphology with multiple catalytic groups appears spontaneously in the course of controlled radical polymerization in a selective solvent. As a proof of concept, we construct a fully synthetic analog of serine hydrolase, e.g.α-chymotrypsin, using the conformation-dependent sequence design approach and multiscale simulation that combines the methods of "mesoscale chemistry" and atomistic molecular dynamics (MD). A 100 ns GPU-accelerated MD simulation of the designed polymer-supported catalyst in the aqueous environment provides valuable information on the structural organization of this system that has been synthesized in our Lab.