Pre-existing soft modes of motion uniquely defined by native contact topology facilitate ligand binding to proteins

Modeling protein flexibility constitutes a major challenge in accurate prediction of protein-ligand and protein-protein interactions in docking simulations. The lack of a reliable method for predicting the conformational changes relevant to substrate binding prevents the productive application of computational docking to proteins that undergo large structural rearrangements. Here, we examine how coarse-grained normal mode analysis has been advantageously applied to modeling protein flexibility associated with ligand binding. First, we highlight recent studies that have shown that there is a close agreement between the large-scale collective motions of proteins predicted by elastic network models and the structural changes experimentally observed upon ligand binding. Then, we discuss studies that have exploited the predicted soft modes in docking simulations. Two general strategies are noted: pregeneration of conformational ensembles that are then utilized as input for standard fixed-backbone docking and protein structure deformation along normal modes concurrent to docking. These studies show that the structural changes apparently "induced" upon ligand binding occur selectively along the soft modes accessible to the protein prior to ligand binding. They further suggest that proteins offer suitable means of accommodating/facilitating the recognition and binding of their ligand, presumably acquired by evolutionary selection of the suitable three-dimensional structure.     

Evidence of protein collective motions on the picosecond timescale

We investigate the presence of structural collective motions on a picosecond timescale for the heme protein, cytochrome c, as a function of oxidation and hydration, using terahertz (THz) time domain spectroscopy and molecular dynamics simulations. The THz response dramatically increases with oxidation, with the largest increase for lowest hydrations, and highest frequencies. For both oxidation states the THz response rapidly increases with hydration saturating above ∼25% (g H₂O/g protein). Quasiharmonic vibrational modes and dipole-dipole correlation functions were calculated from molecular dynamics trajectories. The collective mode density of states alone reproduces the measured hydration dependence, providing strong evidence of the existence of these motions. The large oxidation dependence is reproduced only by the dipole-dipole correlation function, indicating the contrast arises from diffusive motions consistent with structural changes occurring in the vicinity of buried internal water molecules. This source for the observed oxidation dependence is consistent with the lack of an oxidation dependence in nuclear resonant vibrational spectroscopy measurements.     

Discovering conformational sub-states relevant to protein function

Background

Internal motions enable proteins to explore a range of conformations, even in the vicinity of native state. The role of conformational fluctuations in the designated function of a protein is widely debated. Emerging evidence suggests that sub-groups within the range of conformations (or sub-states) contain properties that may be functionally relevant. However, low populations in these sub-states and the transient nature of conformational transitions between these sub-states present significant challenges for their identification and characterization.

Methods and Findings

To overcome these challenges we have developed a new computational technique, quasi-anharmonic analysis (QAA). QAA utilizes higher-order statistics of protein motions to identify sub-states in the conformational landscape. Further, the focus on anharmonicity allows identification of conformational fluctuations that enable transitions between sub-states. QAA applied to equilibrium simulations of human ubiquitin and T4 lysozyme reveals functionally relevant sub-states and protein motions involved in molecular recognition. In combination with a reaction pathway sampling method, QAA characterizes conformational sub-states associated with cis/trans peptidyl-prolyl isomerization catalyzed by the enzyme cyclophilin A. In these three proteins, QAA allows identification of conformational sub-states, with critical structural and dynamical features relevant to protein function.

Conclusions

Overall, QAA provides a novel framework to intuitively understand the biophysical basis of conformational diversity and its relevance to protein function.     

A Physical Picture of Protein Dynamics and Conformational Changes

A physical model is reviewed which explains different aspects of protein dynamics consistently. At low temperatures, the molecules are frozen in conformational substates. Their average energy is 3/2RT. Solid-state vibrations occur on a time scale of femtoseconds to nanoseconds. Above a characteristic temperature, often called the dynamical transition temperature, slow modes of motions can be observed occurring on a time scale between about 140 and 1 ns. These motions are overdamped, quasidiffusive, and involve collective motions of segments of the size of an α-helix. Molecules performing these types of motion are in the “flexible state”. This state is reached by thermal activation. It is shown that these motions are essential for conformational relaxation. Based on this picture, a new approach is proposed to understand conformational changes. It connects structural fluctuations and conformational transitions.     

Coarse-Grained Simulation of Full-Length Integrin Activation

Integrin conformational dynamics are critical to their receptor and signaling functions in many cellular processes, including spreading, adhesion, and migration. However, assessing integrin conformations is both experimentally and computationally challenging because of limitations in resolution and dynamic sampling. Thus, structural changes that underlie transitions between conformations are largely unknown. Here, focusing on integrin αvβ₃, we developed a modified form of the coarse-grained heterogeneous elastic network model (hENM), which allows sampling conformations at the onset of activation by formally separating local fluctuations from global motions. Both local fluctuations and global motions are extracted from all-atom molecular dynamics simulations of the full-length αvβ₃ bent integrin conformer, but whereas the former are incorporated in the hENM as effective harmonic interactions between groups of residues, the latter emerge by systematically identifying and treating weak interactions between long-distance domains with flexible and anharmonic connections. The new hENM model allows integrins and single-point mutant integrins to explore various conformational states, including the initiation of separation between α- and β-subunit cytoplasmic regions, headpiece extension, and legs opening.     

Molecular dynamics simulations of peptides and proteins with amplified collective motions

We present a novel method that uses the collective modes obtained with a coarse-grained model/anisotropic network model to guide the atomic-level simulations. Based on this model, local collective modes can be calculated according to a single configuration in the conformational space of the protein. In the molecular dynamics simulations, the motions along the slowest few modes are coupled to a higher temperature by the weak coupling method to amplify the collective motions. This amplified-collective-motion (ACM) method is applied to two test systems. One is an S-peptide analog. We realized the refolding of the denatured peptide in eight simulations out of 10 using the method. The other system is bacteriophage T4 lysozyme. Much more extensive domain motions between the N-terminal and C-terminal domain of T4 lysozyme are observed in the ACM simulation compared to a conventional simulation. The ACM method allows for extensive sampling in conformational space while still restricting the sampled configurations within low energy areas. The method can be applied in both explicit and implicit solvent simulations, and may be further applied to important biological problems, such as long timescale functional motions, protein folding/unfolding, and structure prediction.     

A Link between Hinge-Bending Domain Motions and the Temperature Dependence of Catalysis in 3-Isopropylmalate Dehydrogenase

Enzyme function depends on specific conformational motions. We show that the temperature dependence of enzyme kinetic parameters can provide insight into these functionally relevant motions. While investigating the catalytic properties of IPMDH from Escherichia coli, we found that its catalytic efficiency (k_cat/K_M,IPM) for the substrate IPM has an unusual temperature dependence, showing a local minimum at ∼35°C. In search of an explanation, we measured the individual constants k_cat and K_M,IPM as a function of temperature, and found that the van 't Hoff plot of K_M,IPM shows sigmoid-like transition in the 20-40°C temperature range. By means of various measurements including hydrogen-deuterium exchange and fluorescence resonance energy transfer, we showed that the conformational fluctuations, including hinge-bending domain motions increase more steeply with temperatures >30°C. The thermodynamic parameters of ligand binding determined by isothermal titration calorimetry as a function of temperature were found to be strongly correlated to the conformational fluctuations of the enzyme. Because the binding of IPM is associated with a hinge-bending domain closure, the more intense hinge-bending fluctuations at higher temperatures increasingly interfere with IPM binding, thereby abruptly increasing its dissociation constant and leading to the observed unusual temperature dependence of the catalytic efficiency.     

An NMR and MD study of complexes of bacteriophage lambda lysozyme with tetra- and hexa-N-acetylchitohexaose

The X-ray structure of lysozyme from bacteriophage lambda (λ lysozyme) in complex with the inhibitor hexa-N-acetylchitohexaose (NAG6) (PDB: 3D3D) has been reported previously showing sugar units from two molecules of NAG6 bound in the active site. One NAG6 is bound with four sugar units in the ABCD sites and the other with two sugar units in the E′F′ sites potentially representing the cleavage reaction products; each NAG6 cross links two neighboring λ lysozyme molecules. Here we use NMR and MD simulations to study the interaction of λ lysozyme with the inhibitors NAG4 and NAG6 in solution. This allows us to study the interactions within the complex prior to cleavage of the polysaccharide. ¹H^N and ¹⁵N chemical shifts of λ lysozyme resonances were followed during NAG4/NAG6 titrations. The chemical shift changes were similar in the two titrations, consistent with sugars binding to the cleft between the upper and lower domains; the NMR data show no evidence for simultaneous binding of a NAG6 to two λ lysozyme molecules. Six 150 ns MD simulations of λ lysozyme in complex with NAG4 or NAG6 were performed starting from different conformations. The simulations with both NAG4 and NAG6 show stable binding of sugars across the D/E active site providing low energy models for the enzyme-inhibitor complexes. The MD simulations identify different binding subsites for the 5th and 6th sugars consistent with the NMR data. The structural information gained from the NMR experiments and MD simulations have been used to model the enzyme-peptidoglycan complex.     

Conformational Transitions in Yeast Chorismate Mutase Important for Allosteric Regulation as Identified by Nuclear Magnetic Resonance Spectroscopy

Proteins fluctuate between different conformations in solution, and these conformational fluctuations can be important for protein function and allosteric regulation. The chorismate mutase from Saccharomyces cerevisiae (ScCM), a key enzyme in the biosynthesis of aromatic amino acids, is allosterically activated and inhibited by tryptophan and tyrosine, respectively. It was initially proposed that in the absence of effector, ScCM fluctuates between activated R and inhibited T conformations according to the Monod-Wyman-Changeux (MWC) model, although a more complex regulation pattern was later suggested by mutagenesis and kinetic data. Here we used NMR relaxation dispersion experiments to understand the conformational fluctuations on the microsecond-to-millisecond timescale that occur in ScCM. In the absence of allosteric effectors, ScCM did not exclusively exchange between T and R conformations, suggesting that the two-state MWC model is insufficient to explain conformational dynamics. Addition of tyrosine led to the quenching of much of the motion on this timescale, while new motions were identified in the presence of tryptophan. These new motions are consistent with conformational fluctuations into an alternative conformation that may be important for enzyme activity.     

Ligand-induced conformational changes in a thermophilic ribose-binding protein

Background

Members of the periplasmic binding protein (PBP) superfamily are involved in transport and signaling processes in both prokaryotes and eukaryotes. Biological responses are typically mediated by ligand-induced conformational changes in which the binding event is coupled to a hinge-bending motion that brings together two domains in a closed form. In all PBP-mediated biological processes, downstream partners recognize the closed form of the protein. This motion has also been exploited in protein engineering experiments to construct biosensors that transduce ligand binding to a variety of physical signals. Understanding the mechanistic details of PBP conformational changes, both global (hinge bending, twisting, shear movements) and local (rotamer changes, backbone motion), therefore is not only important for understanding their biological function but also for protein engineering experiments.

Results

Here we present biochemical characterization and crystal structure determination of the periplasmic ribose-binding protein (RBP) from the hyperthermophile Thermotoga maritima in its ribose-bound and unliganded state. The T. maritima RBP (tmRBP) has 39% sequence identity and is considerably more resistant to thermal denaturation ( ^app T _m value is 108°C) than the mesophilic Escherichia coli homolog (ecRBP) ( ^app T _m value is 56°C). Polar ligand interactions and ligand-induced global conformational changes are conserved among ecRBP and tmRBP; however local structural rearrangements involving side-chain motions in the ligand-binding site are not conserved.

Conclusion

Although the large-scale ligand-induced changes are mediated through similar regions, and are produced by similar backbone movements in tmRBP and ecRBP, the small-scale ligand-induced structural rearrangements differentiate the mesophile and thermophile. This suggests there are mechanistic differences in the manner by which these two proteins bind their ligands and are an example of how two structurally similar proteins utilize different mechanisms to form a ligand-bound state.     

Energy landscape of a native protein: Jumping-among-minima model

We have investigated energy landscape of human lysozyme in its native state by using principal component analysis and a model, jumping-among-minima (JAM) model. These analyses are applied to 1 nsec molecular dynamics trajectory of the protein in water. An assumption embodied in the JAM model allows us to divide protein motions into intra-substate and inter-substate motions. By examining intra-substate motions, it is shown that energy surfaces of individual conformational substates are nearly harmonic and mutually similar. As a result of principal component analysis and JAM model analysis, protein motions are shown to consist of three types of collective modes, multiply hierarchical modes, singly hierarchical modes, and harmonic modes. Multiply hierarchical modes, the number of which accounts only for 0.5% of all modes, dominate contributions to total mean-square atomic fluctuation. Inter-substate motions are observed only in a small-dimensional subspace spanned by the axes of multiplyhierarchical and singly hierarchical modes. Inter-substate motions have two notable time components: faster component seen within 200 psec and slower component. The former involves transitions among the conformational substates of the low-level hierarchy, whereas the latter involves transitions of the higher level substates observed along the first four multiply hierarchical modes. We also discuss dependence of the subspace, which contains conformational substates, on time duration of simulation. Proteins 33:496–517, 1998. © 1998 Wiley-Liss, Inc.     

Crystal Structures of Bacillus cereus NCTU2 Chitinase Complexes with Chitooligomers Reveal Novel Substrate Binding for Catalysis: A CHITINASE WITHOUT CHITIN BINDING AND INSERTION DOMAINS*

Chitinases hydrolyze chitin, an insoluble linear polymer of N-acetyl-d-glucosamine (NAG)_n, into nutrient sources. Bacillus cereus NCTU2 chitinase (ChiNCTU2) predominantly produces chitobioses and belongs to glycoside hydrolase family 18. The crystal structure of wild-type ChiNCTU2 comprises only a catalytic domain, unlike other chitinases that are equipped with additional chitin binding and insertion domains to bind substrates into the active site. Lacking chitin binding and chitin insertion domains, ChiNCTU2 utilizes two dynamic loops (Gly-67—Thr-69 and Ile-106–Val-112) to interact with (NAG)_n, generating novel substrate binding and distortion for catalysis. Gln-109 is crucial for direct binding with substrates, leading to conformational changes of two loops with a maximum shift of ∼4.6 Å along the binding cleft. The structures of E145Q, E145Q/Y227F, and E145G/Y227F mutants complexed with (NAG)_n reveal (NAG)₂, (NAG)₂, and (NAG)₄ in the active site, respectively, implying various stages of reaction: before hydrolysis, E145G/Y227F with (NAG)₄; in an intermediate state, E145Q/Y227F with a boat-form NAG at the −1 subsite, −1-(NAG); after hydrolysis, E145Q with a chair form −1-(NAG). Several residues were confirmed to play catalytic roles: Glu-145 in cleavage of the glycosidic bond between −1-(NAG) and +1-(NAG); Tyr-227 in the conformational change of −1-(NAG); Asp-143 and Gln-225 in stabilizing the conformation of −1-(NAG). Additionally, Glu-190 acts in the process of product release, and Tyr-193 coordinates with water for catalysis. Residues Asp-143, E145Q, Glu-190, and Tyr-193 exhibit multiple conformations for functions. The inhibitors zinc ions and cyclo-(l-His-l-Pro) are located at various positions and confirm the catalytic-site topology. Together with kinetics analyses of related mutants, the structures of ChiNCTU2 and its mutant complexes with (NAG)_n provide new insights into its substrate binding and the mechanistic action.     

Conformational energy calculations of enzyme-substrate complexes of lysozyme. I. Energy minimization of monosaccharide and oligosaccharide inhibitors and substrates of lysozyme

Conformational energy calculations were performed on monosaccharide and oligosaccharide inhibitors and substrates of lysozyme to examine the preferred conformations of these molecules. A grid-search method was used to locate all of the low-energy conformational regions for N-acetyl-β-D -glycosamine (NAG), and energy minimization was then carried out in each of these regions. Three stable positions for the N-acetyl group have ben located, in two of which the plane of the amide unit is normal to the mean plane of the pyranosyl ring. Nine local energy minima were located for the —CH₂OH group. The positions of the two vicinal cis —OH groups are determined predominantly by interactions with either the —CH₂OH or the N-acetyl group. The most stable conformations of β-N-acetylmuramic acid (NAM) were determined from the study of the low-energy conformations of NAG. In the two stable orientations for the D -lactic acid side chain, the O—C—C′ plane (C′ being the carbon atom of the terminal carboxyl group) was found to be normal to the mean plane of the pyranosyl ring. The low-energy positions for the COOH group of NAM are determined mainly by interactions with neighboring groups. The conformational preferences of the α-anomers of NAG and NAM were also explored. The calculated conformation of the N-acetyl group for α-NAG was quite close to that determined by X-ray analysis. Two of the three lowest energy conformations of α-NAM are similar to the corresponding conformations of the β-anomer. A third low-energy structure, which has a hydrogen bond from the NH of the N-acetyl group to the C?O of the lactic acid group, corresponds very closely to the X-ray structure of this molecule. The preferred conformations of the disaccharides NAG–NAG, NAM–NAG and NAG–NAM were also investigated. Two preferred orientations of the reducing pyranosyl ring relative to the nonreducing ring were found for all of these disaccharides, both of which are close to the extended conformation. In one of these conformations, a hydrogen bond can form between the OH group attached to C₃ of the reducing sugar and the ring oxygen of the preceding residue. Each conformation can be stabilized further by a hydrogen bond between the CH₂OH (donor) of residue i + 1 and the C?O of residue i (acceptor). The interactions that determine conformations for all oligosaccharides containing both NAG and NAM are shown to be exclusively intraresidue and nearest neighbor interactions, so that it is possible to predict all stable conformations of oligosaccharides containing NAG and NAM in any sequence.     

Study of the phase transition in lysozyme crystals by Raman spectroscopy

BackgroundRecently, it has been revealed that tetragonal lysozyme crystals show a phase transition at 307 K upon heating. The underlying mechanisms of the phase transition are still not fully understood. Here we focus on the study of high-frequency vibrational modes arising from the protein and their temperature evolution in the vicinity of T_ph as well as on the detailed study of crystalline water dynamics near T_ph.MethodsRaman experiments have been performed at temperatures 295–323 K including T_ph. The low-frequency modes and the modes of fingerprint region, CH- and OH-stretching regions have been analyzed.Results and conclusionsIn spite of the absence of noticeable rearrangements in protein structure, the high-frequency vibrational modes of lysozyme located in the fingerprint region have been found to exhibit the features of critical dynamics near T_ph. Pronounced changes in the dynamics of α-helixes and Tyr residues exposed on the protein surface point to the important role of H-bond rearrangements at the phase transition. Additionally the study of temperature evolution of OH-stretching modes has shown an increase in distortions of tertahedral H-bond network of crystalline water above T_ph. These changes in water dynamics could play a crucial role in the mechanisms of the phase transition.General significanceThe present results shed light on the mechanisms of the phase transition in lysozyme crystals.     

Modeling correlated main-chain motions in proteins for flexible molecular recognition

We describe a new method for modeling protein and ligand main-chain flexibility, and show its ability to model flexible molecular recognition. The goal is to sample the full conformational space, including large-scale motions that typically cannot be reached in molecular dynamics simulations due to the computational intensity, as well as conformations that have not been observed yet by crystallography or NMR. A secondary goal is to assess the degree of flexibility consistent with protein-ligand recognition. Flexibility analysis of the target protein is performed using the graph-theoretic algorithm FIRST, which also identifies coupled networks of covalent and noncovalent bonds within the protein. The available conformations of the flexible regions are then explored with ROCK by random-walk sampling of the rotatable bonds. ROCK explores correlated motions by only sampling dihedral angles that preserve the coupled bond networks in the protein and generates conformers with good stereochemistry, without using a computationally expensive potential function. A representative set of the conformational ensemble generated this way can be used as targets for docking with SLIDE, which handles the flexibility of protein and ligand side-chains. The realism of this protein main-chain conformational sampling is assessed by comparison with time-resolved NMR studies of cyclophilin A motions. ROCK is also effective for modeling the flexibility of large cyclic and polycyclic ligands, as demonstrated for cyclosporin and zearalenol. The use of this combined approach to perform docking with main-chain flexibility is illustrated for the cyclophilin A-cyclosporin complex and the estrogen receptor in complex with zearalenol, while addressing the question of how much flexibility is allowed without hindering molecular recognition.     

Nuclear magnetic resonance‐based determination of dioxygen binding sites in protein cavities

The location and ligand accessibility of internal cavities in cysteine‐free wild‐type T4 lysozyme was investigated using O₂ gas‐pressure NMR spectroscopy and molecular dynamics (MD) simulation. Upon increasing the concentration of dissolved O₂ in solvent to 8.9 mM, O₂‐induced paramagnetic relaxation enhancements (PREs) to the backbone amide and side chain methyl protons were observed, specifically around two cavities in the C‐terminal domain. To determine the number of O₂ binding sites and their atomic coordinates from the 1/r⁶ distance dependence of the PREs, we established an analytical procedure using Akaike's Information Criterion, in combination with a grid‐search. Two O₂‐accessible sites were identified in internal cavities: One site was consistent with the xenon‐binding site in the protein in crystal, and the other site was established to be a novel ligand‐binding site. MD simulations performed at 10 and 100 mM O₂ revealed dioxygen ingress and egress as well as rotational and translational motions of O₂ in the cavities. It is therefore suggested that conformational fluctuations within the ground‐state ensemble transiently develop channels for O₂ association with the internal protein cavities.     

Coherent Neutron Scattering and Collective Dynamics in the Protein,GFP

Collective dynamics are considered to be one of the major properties of soft materials, including biological macromolecules. We present coherent neutron scattering studies of the low-frequency vibrations, the so-called boson peak, in fully deuterated green fluorescent protein (GFP). Our analysis revealed unexpectedly low coherence of the atomic motions in GFP. This result implies a low amount of in-phase collective motion of the secondary structural units contributing to the boson peak vibrations and fast conformational fluctuations on the picosecond timescale. These observations are in contrast to earlier studies of polymers and glass-forming systems, and suggest that random or out-of-phase motions of the β-strands contribute greater than two-thirds of the intensity to the low-frequency vibrational spectra of GFP.     

Coherent Neutron Scattering and Collective Dynamics in the Protein,GFP

Collective dynamics are considered to be one of the major properties of soft materials, including biological macromolecules. We present coherent neutron scattering studies of the low-frequency vibrations, the so-called boson peak, in fully deuterated green fluorescent protein (GFP). Our analysis revealed unexpectedly low coherence of the atomic motions in GFP. This result implies a low amount of in-phase collective motion of the secondary structural units contributing to the boson peak vibrations and fast conformational fluctuations on the picosecond timescale. These observations are in contrast to earlier studies of polymers and glass-forming systems, and suggest that random or out-of-phase motions of the β-strands contribute greater than two-thirds of the intensity to the low-frequency vibrational spectra of GFP.