Side-chain conformational changes upon Protein-Protein Association

Conformational changes upon protein-protein association are the key element of the binding mechanism. The study presents a systematic large-scale analysis of such conformational changes in the side chains. The results indicate that short and long side chains have different propensities for the conformational changes. Long side chains with three or more dihedral angles are often subject to large conformational transition. Shorter residues with one or two dihedral angles typically undergo local conformational changes not leading to a conformational transition. A relationship between the local readjustments and the equilibrium fluctuations of a side chain around its unbound conformation is suggested. Most of the side chains undergo larger changes in the dihedral angle most distant from the backbone. The frequencies of the core-to-surface interface transitions of six nonpolar residues and Tyr are larger than the frequencies of the opposite surface-to-core transitions. The binding increases both polar and nonpolar interface areas. However, the increase of the nonpolar area is larger for all considered classes of protein complexes, suggesting that the protein association perturbs the unbound interfaces to increase the hydrophobic contribution to the binding free energy. To test modeling approaches to side-chain flexibility in protein docking, conformational changes in the X-ray set were compared with those in the docking decoy sets. The results lead to a better understanding of the conformational changes in proteins and suggest directions for efficient conformational sampling in docking protocols.     

Induced fit or conformational selection? The role of the semi-closed state in the maltose binding protein

A full characterization of the thermodynamic forces underlying ligand-associated conformational changes in proteins is essential for understanding and manipulating diverse biological processes, including transport, signaling, and enzymatic activity. Recent experiments on the maltose binding protein (MBP) have provided valuable data about the different conformational states implicated in the ligand recognition process; however, a complete picture of the accessible pathways and the associated changes in free energy remains elusive. Here we describe results from advanced accelerated molecular dynamics (aMD) simulations, coupled with adaptively biased force (ABF) and thermodynamic integration (TI) free energy methods. The combination of approaches allows us to track the ligand recognition process on the microsecond time scale and provides a detailed characterization of the protein's dynamic and the relative energy of stable states. We find that an induced-fit (IF) mechanism is most likely and that a mechanism involving both a conformational selection (CS) step and an IF step is also possible. The complete recognition process is best viewed as a "Pac Man" type action where the ligand is initially localized to one domain and naturally occurring hinge-bending vibrations in the protein are able to assist the recognition process by increasing the chances of a favorable encounter with side chains on the other domain, leading to a population shift. This interpretation is consistent with experiments and provides new insight into the complex recognition mechanism. The methods employed here are able to describe IF and CS effects and provide formally rigorous means of computing free energy changes. As such, they are superior to conventional MD and flexible docking alone and hold great promise for future development and applications to drug discovery.     

Rotamer libraries and probabilities of transition between rotamers for the side chains in protein-protein binding

Conformational changes in the side chains are essential for protein-protein binding. Rotameric states and unbound- to-bound conformational changes in the surface residues were systematically studied on a representative set of protein complexes. The side-chain conformations were mapped onto dihedral angles space. The variable threshold algorithm was developed to cluster the dihedral angle distributions and to derive rotamers, defined as the most probable conformation in a cluster. Six rotamer libraries were generated: full surface, surface noninterface, and surface interface-each for bound and unbound states. The libraries were used to calculate the probabilities of the rotamer transitions upon binding. The stability of amino acids was quantified based on the transition maps. The noninterface residues' stability was higher than that of the interface. Long side chains with three or four dihedral angles were less stable than the shorter ones. The transitions between the rotamers at the interface occurred more frequently than on the noninterface surface. Most side chains changed conformation within the same rotamer or moved to an adjacent rotamer. The highest percentage of the transitions was observed primarily between the two most occupied rotamers. The probability of the transition between rotamers increased with the decrease of the rotamer stability. The analysis revealed characteristics of the surface side-chain conformational transitions that can be utilized in flexible docking protocols.     

163 Enhanced sampling of peptides and proteins with a new biasing replica exchange method

Classical MD simulations (cMD) are limited by the sampling of relevant states of the peptides. Replica exchange (REMD) methods aim to search the conformational space of proteins more efficiently (reviewed in Ostermeir & Zacharias, 2013). We have developed a Hamiltonian REMD method that takes advantage of an intrinsic property of proteins, the specific Φ ? dihedral angle combinations along the polymer backbone. By employing a coupled two-dimensional biasing potential the energy barriers along the polymer backbone are reduced more effectively than by a previous approach based on a one-D biasing potential (Kannan & Zacharias, 2007). Thus, adjacent amino acids along the polymers backbone can easily switch between favourable regions in the Ramachandran plot. Additionally, energy barriers of rotameric states of amino acid side chains of proteins are also biased in the replica runs. The method improves the sampling of conformational substates of proteins at a modest number of replicas (nine replicas in the standard set-up with one replica running without biasing potential) compared to much larger numbers necessary in the case of standard temperature (T)-REMD simulations. A further improvement is achieved by a dynamical adjustment of the penalty potential levels in the replicas such that high exchange rates and improved mixing of conformations between different replicas are guaranteed. The biasing potential (BP)-REMD method turns out to be suitable to speed up both the folding of spaghetti-like test peptides and the refinement of loop decoy structures. Starting from extended structures, an α-helical oligo-alanine and β-hairpin chignolin and the Trp-cage protein fold more rapidly in near-native structures than in cMD simulations. The BP-REMD simulations not only accelerate the folding process of test proteins but also enlarge the variety of sampled configurations in conformational space. Since flexible parts of the protein can be penalized selectively, this method provides a precise tool to investigate regions of interest of the protein.     

DFTMD studies of β-cellobiose: conformational preference using implicit solvent

Previous DFT in vacuo studies on the conformational preferences for cellobiose showed that upon optimization the φ_H-anti conformations were of lower energy than the syn forms. Upon optimization using an implicit solvation method, COSMO, the syn or observed form was still not predicted to be of lower energy than the φ_H-anti form, even though optimization after addition of several explicit water molecules did show a relative energy difference favoring the syn form. In order to examine the predictive ability of COSMO on this carbohydrate, constant energy dynamics, DFTMD, simulations were carried out on low energy syn and φ_H-anti conformations with and without COSMO included during the dynamics. The resulting analysis confirmed that when COSMO is included in the dynamics, the syn conformations become energetically favored over the φ_H-anti forms suggesting that both solvent and entropy play roles in dictating the solution conformation of cellobiose. Analysis of the dynamic runs includes distributions of selected dihedral angles versus time, conformational transitions, and populations of some quasi-planar, boat, skew forms during the simulations.     

Conformations of blood group H-active oligosaccharides of ovarian cyst mucins

Spectroscopic data and conformational energy calculations are reported for eight oligosaccharides from ovarian cyst mucins and from human milk, the nonreducing terminals of which have fucose (α1 → 2)galactose linked either (β1 → 3) (type I) or (β1 → 4)(type II) to N-acetylglucosamine or in (β1 → 3) linkage to galactosaminitol. The fully assigned proton nmr spectra are reported along with nuclear Overhauser enhancement (NOE) data. Amide proton coupling constants and vacuum-uv CD spectra provide information on the amide plane orientation and amide environment. Our results imply that the fucosidic dihedral angles are similar for all three cases and that the substantial differences in the chemical shifts of the fucosyl protons of type I, type II, and 3-ol chains result from different perturbations by the amide group of the residue to which the β-galactose is linked. Stereopair diagrams of conformational models for both type I and II H chains are presented that are consistent with NOE, coupling constants, conformational energy calculations, and the CD data. While the temperature dependence of the observed NOE of penta- and hexasaccharides indicates that their rotational correlation times are strongly temperature dependent, we conclude that the conformations are essentially independent of temperature.     

Conformational Behavior and Flexibility of Terminally Blocked Cysteine and Cystine

Conformational potential energy hypersurfaces, PES, for the terminally blocked L-Cysteine, L,L-Cystine and D,L-Cystine have been analyzed by means of molecular mechanics in combination with the programs ROSE, CICADA, PANIC and COMBINE. Low energy conformations and conformational transitions, conformational channels, have been located. Global and fragmental flexibility and conformational softness have been calculated for each conformer as well as for the entire molecule. The PES analyses were used for simulation of conformational movement based on Boltzmann probability of the points obtained on the PES. Boltzmann travelling revealed interesting correlated conformational movement where three or even more dihedral angles changed simultaneously. It could be shown that conformational behavior and flexibility were strongly influenced by the absolute configurations of the amino acids in the peptides.     

A method for characterizing transition concertedness from polymer dynamics computer simulations

A statistical method based on classifying the transitions among a set of dihedral angles within an “energy transfer window” is developed, and used to analyze Brownian (BD) and molecular dynamics (MD) simulations of the acyl chains in a lipid bilayer, and MD of neat hexadecane. It is shown for the BD simulation that when a transition of the dihedral angle in the center of the chain occurs, a transition of a particular next nearest neighbor (or angle 2-apart) will follow concertedly with a probability of approximately 0.10 within a lime window of approximately 3 ps. The MD bilayer simulations, which are based on a more flexible model of the hydrocarbon chains, yield corresponding concerted transition probabilities of approximately 0.083 and window sizes of 1–2 ps. An analysis of angles 4-apart yields concerted transition probabilities of 0.03 and 0.04 for the BD and MD bilayer simulations, respectively, and window sizes close to those of the corresponding 2-apart cases. Statistical hypothesis testing very strongly rejects the assertion that these follower transitions are occurring at random. Similar analysis reveals marginal or no evidence of concertedness between 1-apart (nearest neighbor) and between 3-apart dihedral angle transitions. The pattern of concertedness for hexadecane is qualitatively similar to that of the lipid chains, although concertedness is somewhat stronger for the 3-apart transitions and somewhat weaker for those 4-apart. Finally, it is suggested that the diffusion of small solute molecules in membranes is better facilitated by non concerted transitions, which are associated with relatively large displacements of the chains, than by concerted transitions, which do little to change the chain shape. © 1995 John Wiley & Sons, Inc.     

Coarse‐grained simulations of proton‐dependent conformational changes in lactose permease

During lactose/H⁺ symport, the Escherichia coli lactose permease (LacY) undergoes a series of global conformational transitions between inward‐facing (open to cytoplasmic side) and outward‐facing (open to periplasmic side) states. However, the exact local interactions and molecular mechanisms dictating those large‐scale structural changes are not well understood. All‐atom molecular dynamics simulations have been performed to investigate the molecular interactions involved in conformational transitions of LacY, but the simulations can only explore early or partial global structural changes because of the computational limits (< 100 ns). In this work, we implement a hybrid force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid, to investigate the proton‐dependent dynamics and conformational changes of LacY. The effects of the protonation states on two key glutamate residues (Glu325 and Glu269) have been studied. Our results on the salt‐bridge dynamics agreed with all‐atom simulations at early short time period, validating our simulations. From our microsecond simulations, we were able to observe the complete transition from inward‐facing to outward‐facing conformations of LacY. Our results showed that all helices have participated during the global conformational transitions and helical movements of LacY. The inter‐helical distances measured in our simulations were consistent with the double electron‐electron resonance experiments at both cytoplasmic and periplasmic sides. Our simulations indicated that the deprotonation of Glu325 induced the opening of the periplasmics side and partial closure of the cytoplasmic side of LacY, while protonation of the Glu269 caused a stable cross‐domain salt‐bridge (Glu130‐Arg344) and completely closed the cytoplasmic side. Proteins 2016; 84:1067–1074. © 2016 Wiley Periodicals, Inc.     

Variation of the NH-C alpha-H coupling constant with dihedral angle in the NMR spectra of peptides

Proton magnetic resonance data and conformational calculations of a series of model compounds containing a NH-C^αH group substituted as in peptides have been used to generate a proton–proton coupling constant–dihedral angle relation for the peptide unit. For those substances used in which the dihedral angle about the N-C^α bond is not fixed, the angle distribution was calculated from conformational theory. Using eight examples in which the number of theoretical assumptions were least, the best values of the coefficients A, B, and C in the expression J(θ) = Acos²θ + B cosθ + Csin²θ were found by a least-squares procedure to be 7.9, ?1.55, and 1.35, respectively. This relation gives reasonable values for the dihedral angles ? in cyclic oligopeptide structures for which the availability of both NMR data and other structural information allow comparison. When applied to N-acetylamino acid N-methylamides having side chains extending beyond C^β, however, agreement with the calculated conformational distribution was found for Leu, Met, and Trp, but observed values of J were larger than expected for Val, He, Phe, and Tyr, These disagreements are considered to be the result of interactions not yet taken into account in the usual conformational calculations.     

Predicting order of conformational changes during protein conformational transitions using an interpolated elastic network model

The decryption of sequence of structural events during protein conformational transitions is essential to a detailed understanding of molecular functions ofvarious biological nanomachines. Coarse‐grained models have proven useful by allowing highly efficient simulations of protein conformational dynamics. By combining two coarse‐grained elastic network models constructed based on the beginning and end conformations of a transition, we have developed an interpolated elastic network model to generate a transition pathway between the two protein conformations. For validation, we have predicted the order of local and global conformational changes during key ATP‐driven transitions in three important biological nanomachines (myosin, F₁ ATPase and chaperonin GroEL). We have found that the local conformational change associated with the closing of active site precedes the global conformational change leading to mechanical motions. Our finding is in good agreement with the distribution of intermediate experimental structures, and it supports the importance of local motions at active site to drive or gate various conformational transitions underlying the workings of a diverse range of biological nanomachines. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Theoretical probes of conformational fluctuations in S-peptide and RNase A/3′–UMP enzyme product complex

The dynamic properties of the RNase A/3′–UMP enzyme/product complex and the S-peptide of RNase A have been investigated by molecular dynamics simulations using suitable generalization of ideas introduced to probe the energy landscape in structural glasses. We introduce two measures, namely, the kinetic energy fluctuation metric and the force metric, both of which are used to calculate the time needed for sampling the conformation space of the molecules. The calculation of the fluctuation metric requires a single trajectory whereas the force metric is computed using two independent trajectories. The vacuum MD simulations show that for both systems the time required for kinetic energy equipartitioning is surprisingly long even at high temperatures. We show that the force metric is a powerful means of probing the nature and relative importance of conformational substates which determine the dynamics at low temperatures. In particular the time dependence of the non-bonded force metric is used to demonstrate that at low temperatures the system is predominantly localized hi a single cluster of conformational substates. The force metric is used to show that relaxation of long range (in sequence space) interactions must be mediated by a sequence of local dihedral angle transitions. We also argue that the time needed for compact structure formation is intimately related to the time needed for the relaxation of the dihedral angle degrees of freedom. The tame for non-bonded interactions, which drive protein molecules to fold under appropriate conditions, to relax becomes extremely long as the temperature is lowered suggesting that the formation of maximally compact structure hi proteins must be a very slow process. © 1993 Wiley-Liss, Inc.     

Flow-induced structural transition in the beta-switch region of glycoprotein Ib

The impact of fluid flow on structure and dynamics of biomolecules has recently gained much attention. In this article, we present a molecular-dynamics algorithm that serves to generate stable water flow under constant temperature, for the study of flow-induced protein behavior. Flow simulations were performed on the 16-residue β-switch region of platelet glycoprotein Ibα, for which crystal structures of its N-terminal domain alone and in complex with the A1 domain of von Willebrand factor have been solved. Comparison of the two structures reveals a conformational change in this region, which, upon complex formation, switches from an unstructured loop to a β-hairpin. Interaction between glycoprotein Ibα and von Willebrand factor initiates platelet adhesion to injured vessel walls, and the adhesion is enhanced by blood flow. It has been hypothesized that the loop to β-hairpin transition in glycoprotein Ibα is induced by flow before binding to von Willebrand factor. The simulations revealed clearly a flow-induced loop→β-hairpin transition. The transition is dominated by the entropy of the protein, and is seen to occur in two steps, namely a dihedral rotation step followed by a side-group packing step.     

α/310-Helix transitions in α-methylalanine homopeptides: Conformational transition pathway and potential of mean force

The conformational transition between the α- and 3₁₀-helical states of α-methylalanine homopeptides is studied with molecular mechanics. Conformational transition pathways for Ace-(MeA)_n-NMe with n = 7, 9, and 11 are obtained with the algorithms of Elber and co-workers [R. Czerminski & R. Elber (1990) International Journal of Quantum Chemistry, Vol. 24, pp. 167–186; A. Ulitsky & R. Elber (1990) Journal of Chemical Physics, Vol. 92, pp. 1510–1511]. The free energy surface, or potential of mean force, for the conformational transition of Ace-(MeA)₉-NMe is calculated from molecular dynamics simulations, and a method is presented for the decomposition of the free energy surface into the constituent energetic and entropic terms, via the calculation of the required temperature derivatives in situ. For the AMBER/OPLS model employed here, the conformational transition pathways each contain a single 3₁₀-helical-like transition state, and the transition state potential energy relative to the 3₁₀-conformation is 3 kcal/mol, independent of peptide length. Entropic stabilization in the barrier region significantly lowers the activation free energies for the forward and reverse transitions from the estimates of the barrier heights based simply on potential energy alone. © 1994 John Wiley & Sons, Inc.     

Effect of side chains on the conformational energy and rotational strength of the n-pi transition for some alpha-helical poly-alpha-amino acids

Calculations of the dependence of the conformational energy and the rotational strength of the amide n–π* electronic transition (in a series of α-helical polyhel-α- amino acids with different side chains) on conformation have been carried out. The conformational energies were computed by procedures developed in this laboratory; the computation of rotational strengths was carried out by the method of Schellman and Oriel, with a slight modification. Polyamino acids with both nonpolar and polar side chains were considered; in the latter case, it was assumed that the only influence of the polar side chain was on the backbone conformation and on the electrostatic field which perturbs the amide chromophore of the backbone. Only conformations in the range of backbone dihedral angles of the right- and left-handed a-helices were considered, and the assumption of regularity (i.e., uniformity of dihedral angles in every residue) was made. The rotational strength per residue was found to vary markedly with chain length (in oligomers of up to 40 residues long); both the conformational energy per residue and the rotational strength per residue were found to vary significantly with the backbone conformation, which in turn depends on the nature of the side chain. The geometry of the hydrogen bond in the α-helical backbone is the most important factor which influences the dependence of the rotational strength on conformation. The implications of these results, for the interpretation of experimental circular dichroism and optical rotatory dispersion data, are discussed.     

Comparative density functional theory based study of the reactivity of Cu,Ag, and Au nanoparticles and of (111) surfaces toward CO oxidation and NO2 reduction

Development of multi-target drugs is becoming increasingly attractive in the repertoire of protein kinase inhibitors discovery. In this study, we carried out molecular docking, molecular dynamics simulations, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energy calculations, principal component analysis (PCA), and dynamical cross-correlation matrices (DCCM) to dissect the molecular mechanism for the valmerin-19 acting as a dual inhibitor for glycogen synthase kinase 3β (GSK3β) and cyclin-dependent kinase 5 (CDK5). Detailed MM-PBSA calculations revealed that the binding free energies of the valmerin-19 to GSK3β/CDK5 were calculated to be ?12.60?±?2.28 kcal mol^-1 and ?11.85?±?2.54 kcal mol^-1, respectively, indicating that valmerin-19 has the potential to act as a dual inhibitor of GSK3β/CDK5. The analyses of PCA and DCCM results unraveled that binding of the valmerin-19 reduced the conformational dynamics of GSK3β/CDK5 and the valmerin-19 bound to GSK3β/CDK5 might occur mostly through a conformational selection mechanism. This study may be helpful for the future design of novel and potent dual GSK3β/CDK5 inhibitors.     

Investigation of selective binding of inhibitors to PTP1B and TCPTP by accelerated molecular dynamics simulations

Abstract

Protein tyrosine phosphatase 1B (PTP1B), a key negative regulator in insulin signaling pathways, is regarded as a potential target for the treatment of type II diabetes and obesity. However, the mechanism underlying the selectivity of PTP1B inhibitors against T-cell protein tyrosine phosphatase (TCPTP) remains controversial, which is due to the high similarity between PTP1B and TCPTP sequence and the fact that no ligand–protein complex of TCPTP has been established yet. Here, the accelerated molecular dynamics (aMD) method was used to investigate the structural dynamics of PTP1B and TCPTP that are bound by two chemically similar inhibitors with distinct selectivity. The conformational transitions during the “open” to “close” states of four complexes were captured, and free energy profiles of important residue pairs were analyzed in detail. Additional MM-PBSA calculations confirmed that the binding free energies of final states were consistent with the experimental results, and the energetic contributions of important residues were further investigated by alanine scanning mutagenesis. By comparing the four complexes, the different conformational behavior of WPD-loop, R-loop, and the second pTyr binding site induced by inhibitors were featured and found to be crucial for the selectivity of inhibitors. This study provides new mechanistic insights of specific binding of inhibitors to PTP1B and TCPTP, which can be exploited to the further structural-based inhibitor design.

Communicated by Ramaswamy H. Sarma     

Equilibrium folding and unfolding pathways for a model protein

The folding–unfolding process of reduced bovine pancreatic trypsin inhibitor was investigated with an idealized model employing approximate free energies. The protein is regarded to consist of only C^α and C^β atoms. The backbone dihedral angles are the only conformational variables and are permitted to take discrete values at every 10°. Intraresidue energies consist of two terms: an empirical part taken from the observed frequency distributions of (?,ψ) and an additional favorable energy assigned to the native conformation of each residue. Interresidue interactions are simplified by assuming that there is an attractive energy operative only between residue pairs in close contact in the native structure. A total of 230,000 molecular conformations, with no atomic overlaps, ranging from the native state to the denatured state, are randomly generated by changing the sampling bias. Each conformation is classified according to its conformational energy, F; a conformational entropy, S(F) is estimated for each value of F from the number of samples. The dependence of S(F) on energy reveals that the folding–unfolding transition for this idealized model is an “all-or-none” type; this is attributable to the specific long-range interactions. Interresidue contact probabilities, averaged over samples representing various stages of folding, serve to characterize folding intermediates. Most probable equilibrium pathways for the folding–unfolding transition are constructed by connecting conformationally similar intermediates. The specific details obtained for bovine pancreatic trypsin inhibitor are as follows: (1) Folding begins with the appearance of nativelike medium-range contacts at a β-turn and at the α-helix. (2) These grow to include the native pair of interacting β-strands. This state includes intact regular secondary conformations, as well as the interstrand sheet contacts, and corresponds to an activated state with the highest free energy on the pathway. (3) Additional native long-range contacts are completely formed either toward the amino terminus or toward the carboxyl terminus. (4) In a final step, the missing contacts appear. Although these folding pathways for this model are not consistent with experimental reports, it does indicate multiple folding pathways. The method is general and can be applied to any set of calculated conformational energies and furthermore permits investigation of gross folding features.     

Equilibrium transitions between side‐chain conformations in leucine and isoleucine

Despite recent improvements in computational methods for protein design, we still lack a quantitative, predictive understanding of the intrinsic probabilities for amino acids to adopt particular side‐chain conformations. Surprisingly, this question has remained unsettled for many years, in part because of inconsistent results from different experimental approaches. To explicitly determine the relative populations of different side‐chain dihedral angles, we performed all‐atom hard‐sphere Langevin Dynamics simulations of leucine (Leu) and isoleucine (Ile) dipeptide mimetics with stereo‐chemical constraints and repulsive‐only steric interactions between non‐bonded atoms. We determine the relative populations of the different χ₁ and χ₂ dihedral angle combinations as a function of the backbone dihedral angles ? and ψ. We also propose, and test, a mechanism for inter‐conversion between the different side‐chain conformations. Specifically, we discover that some of the transitions between side‐chain dihedral angle combinations are very frequent, whereas others are orders of magnitude less frequent, because they require rare coordinated motions to avoid steric clashes. For example, to transition between different values of χ₂, the Leu side‐chain bond angles κ₁ and κ₂ must increase, whereas to transition in χ₁, the Ile bond angles λ₁ and λ₂ must increase. These results emphasize the importance of computational approaches in stimulating further experimental studies of the conformations of side‐chains in proteins. Moreover, our studies emphasize the power of simple steric models to inform our understanding of protein structure, dynamics, and design. Proteins 2015; 83:1488–1499. © 2015 Wiley Periodicals, Inc.