What induces pocket openings on protein surface patches involved in protein–protein interactions?

We previously showed for the proteins BCL-X_L, IL-2, and MDM2 that transient pockets at their protein–protein binding interfaces can be identified by applying the PASS algorithm to molecular dynamics (MD) snapshots. We now investigated which aspects of the natural conformational dynamics of proteins induce the formation of such pockets. The pocket detection protocol was applied to three different conformational ensembles for the same proteins that were extracted from MD simulations of the inhibitor bound crystal conformation in water and the free crystal/NMR structure in water and in methanol. Additional MD simulations studied the impact of backbone mobility. The more efficient CONCOORD or normal mode analysis (NMA) techniques gave significantly smaller pockets than MD simulations, whereas tCONCOORD generated pockets comparable to those observed in MD simulations for two of the three systems. Our findings emphasize the influence of solvent polarity and backbone rearrangements on the formation of pockets on protein surfaces and should be helpful in future generation of transient pockets as putative ligand binding sites at protein–protein interfaces. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Effects of conformational flexibility of alkyl chains of cations on diffusion of ions in ionic liquids

Molecular dynamics simulations of ionic liquids [1-alkyl-3-methylimidazolium (alkyl = ethyl, butyl and hexyl), N-butylpyridinium, N-butyl-N,N,N-trimethylammonium and N-butyl-N-methylpyrrolidinium cations combined with the (CF(3)SO(2))(2)N(-) (TFSA) anion] show that the conformational flexibility of the alkyl chains in the cations is one of the important factors determining the diffusion of ions. Artificial constraint imposed on the internal rotation of alkyl chains significantly decreases the self-diffusion coefficients of cations and anions. The internal rotation of the C-N bond connecting the alkyl chain and the aromatic ring has large effects on the diffusion of ions in imidazolium and pyridinium based ionic liquids. The calculated self-diffusion coefficients of cations and anions decrease 20-40% by imposing the torsional constraint of the C-N bond. On the other hand the torsional constraint of the C-N bond does not largely change the diffusion of ions in the quaternary alkyl ammonium based ionic liquids. The conformational flexibility of the terminal C-C-C-C bond of the alkyl chains has large effects on the diffusion of ions in the quaternary alkyl ammonium based ionic liquids. The influence of the electrostatic interactions and the high density of ionic liquids on the diffusion of ions were studied. The electrostatic interactions have the paramount importance on the slow diffusion of ions in ionic liquids, while the high density of ionic liquids is also responsible for the slow diffusion. The electrostatic interactions and the high density of ionic liquids enhance the effects of the torsional constraint on the diffusion of ions, which suggests that the charge-ordering structure and small free volume originated in the strong electrostatic interactions are the causes of the significant effects of the conformational flexibility on the diffusion of ions in ionic liquids.     

Structural analysis of the conformational flexibility of tris(pyrazolyl)borate ligands and their analogues

Database analysis and molecular mechanics were used to determine the conformational flexibility of tridentate scorpionate ligands. The tris(pyrazolyl)methane and tris(pyrazolyl)borate ligands act like molecular vises, opening their tripodal structure for larger metals and closing around smaller metal ions. Tris(3-tert-butylpyrazolyl)methane has significant preference for larger metal ions than its unsubstituted parent compound. Tris(pyrazolyl)methanes and tris(pyrazolyl)borates have similar conformational flexibilities. Placing sterically hindered groups on the central carbon or boron has only a minor effect on the geometry of the tris(pyrazolyl)methanes and tris(pyrazolyl)borates. However, it does influence the flexibility of the ligands, particularly when they have to open far from their ideal geometry, which commonly occurs.     

Chemical modulation of peptoids: synthesis and conformational studies on partially constrained derivatives

The high conformational flexibility of peptoids can generate problems in biomolecular selectivity as a result of undesired off-target interactions. This drawback can be counterbalanced by restricting the original flexibility to a certain extent, thus leading to new peptidomimetics. By starting from the structure of an active peptoid as an apoptosis inhibitor, we designed two families of peptidomimetics that bear either 7-substituted perhydro-1,4-diazepine-2,5-dione 2 or 3-substituted 1,4-piperazine-2,5-dione 3 moieties. We report an efficient, solid-phase-based synthesis for both peptidomimetic families 2 and 3 from a common intermediate. An NMR spectroscopic study of 2a,b and 3a,b showed two species in solution in different solvents that interconvert slowly on the NMR timescale. The cis/trans isomerization around the exocyclic tertiary amide bond is responsible for this conformational behavior. The cis isomers are more favored in nonpolar environments, and this preference is higher for the six-membered-ring derivative 3a,b. We propose that the hydrogen-bonding pattern could play an important role in the cis/trans equilibrium process. These hydrogen bonds were characterized in solution, in the solid state (i.e., by using X-ray studies), and by molecular modeling of simplified systems. A comparative study of a model peptoid 10 containing the isolated tertiary amide bond under study outlined the importance of the heterocyclic moiety for the prevalence of the cis configuration in 2a and 3a. The kinetics of the cis/trans interconversion in 2a, 3a, and 10 was also studied by variable-temperature NMR spectroscopic analysis. The full line-shape analysis of the NMR spectra of 10 revealed negligible entropic contribution to the energetic barrier in this conformational process. A theoretical analysis of 10 supported the results observed by NMR spectroscopic analysis. Overall, these results are relevant for the study of the peptidomimetic/biological-target interactions.     

Conformational studies of the robust 2-Cys peroxiredoxin Salmonella typhimurium AhpC by solution phase hydrogen/deuterium (H/D) exchange monitored by electrospray ionization mass spectrometry

This is the first comprehensive HX-MS study of a "robust" 2-Cys peroxiredoxin (Prx), namely Salmonella typhimurium AhpC (StAhpC). Prx proteins control intracellular peroxide levels and are abundant antioxidant proteins in eukaryotes, archaea and bacteria. Crystal structural analyses and structure/activity studies of several bacterial and mammalian 2-Cys Prxs have revealed that the activity of 2-Cys Prxs is regulated by redox-dependent oligmerization and a sensitivity of the active site cysteine residue to overoxidation. The propensity to overoxidation is linked to the conformational flexibility of the peroxidatic active site loop. The HX-MS results emphasize the modulation of the conformational motility of the active site loop by disulfide formation. To obtain information on the conformational impact of decamer formation on the active site loop motility, mutants with Thr77 substituted by Ile, a decamer-disrupting mutation or by Val, a decamer-stabilizing mutation, were studied. For the isoleucine mutant, enhanced mobility was observed for regions encompassing the α4 helix located in the dimer-dimer interface and regions surrounding the peroxidatic loop. In contrast, the T77V mutation resulted in an increase in conformational stability in most regions of the protein except for the active site loop and the region encompassing the resolving cysteine.     

Conformational flexibility in open chain molecules: chiral and achiral alkanes

Conformational partition functions of chiral and achiral alkanes have been computed by using a continuum approach (instead of rotational isomeric state approximations). The accessible conformational space per bond depends upon the structure of the compound and is only in the range of 5–13% of the maximum accessible range. In order to partly overcome the intrinsic ambiguity of the term “conformational flexibility”, the distinction between number flexibility (a measure of the number of accessible energy minima) and space flexibility (a measure of the total allotted space) is proposed. Further, the conformational versatility of each bond of a molecule is evaluated in terms of the a priori probability density function of that bond, and it is shown that the use of this function permits a comparison of the relative conformational flexibilities of the individual bonds, which is particularly useful for molecules having more than two rotation angles (where the conventional energy maps cannot be used). Optical rotations are calculated for a series of chiral alkanes by combining the continuum approach for conformational analysis and a recent optical activity calculation scheme. Contributions of single bonds to the molar optical rotation are evaluated and discussed. The influence of temperature upon conformational and chiral properties is evaluated.     

Molecular dynamics simulation study on zwitterionic structure to maintain the normal conformations of Glutathione

Molecular dynamics simulations were applied to normal conformational Glutathione (GSH) and GSH over zwitterionic and hydrophobic surfaces respectively. Conformational analysis of GSH during the simulation time on RMSD, conformational flexibility and dihedral distribution were performed. The re- sults showed that zwitterionic structure maintains the normal conformations of GSH to a better extent, which should be a first good proof of the hypothesis of "maintain of normal structure".     

Structural comparison of endomorphin-2 and its conformationally restricted analog

In the present study, the effect of a conformational constraint introduced into the endomorphin-2 (Tyr-Pro-Phe-Phe-NH₂, EM-2) structure was studied using computational analysis and radioligand binding assay. EM-2 was modified by connecting nitrogen atoms of both phenylalanine residues by a methylene bridge. The obtained analog did not bind to the μ- or δ-opioid receptors in the in vitro studies. The computational analysis of this analog showed twisted, type IV turns and the absence of canonical β-turns typical for the EM-2 structure, which can be explained by the lack of hydrogen bonds involving Phe⁴. Our results show that the introduction of chemical constraint in the EM-2 structure has a significant effect on opioid receptor affinity and in vitro bioactivity.     

Minimizing the Entropy Penalty for Ligand Binding: Lessons from the Molecular Recognition of the Histo Blood‐Group Antigens by Human Galectin‐3

Ligand conformational entropy plays an important role in carbohydrate recognition events. Glycans are characterized by intrinsic flexibility around the glycosidic linkages, thus in most cases, loss of conformational entropy of the sugar upon complex formation strongly affects the entropy of the binding process. By employing a multidisciplinary approach combining structural, conformational, binding energy, and kinetic information, we investigated the role of conformational entropy in the recognition of the histo blood‐group antigens A and B by human galectin‐3, a lectin of biomedical interest. We show that these rigid natural antigens are pre‐organized ligands for hGal‐3, and that restriction of the conformational flexibility by the branched fucose (Fuc) residue modulates the thermodynamics and kinetics of the binding process. These results highlight the importance of glycan flexibility and provide inspiration for the design of high‐affinity ligands as antagonists for lectins.     

Asymmetric catalysis with 7-ring chelate diphosphines: DIOP,BINAP and conformational mobility

The conformations of two classes of 7-ring diphosphine metal chelate have been analysed, based on DIOP- or BINAP-type structures. A combination of X-ray analysis, DFT calculations, data analysis based on the CDS structure database and solution NMR studies has been employed. The conformational flexibility of DIOP-type structures has been endorsed, and the scope of BINAP complex flexibility defined. BINAP complexes possess an intrinsic conformational mobility. Analysis of the metal-adjacent torsion angles C–P–M–P′ and C′–P′–M–P provides a useful probe for the ligand–metal environment, and may be more generally useful.     

蒙特卡洛模拟退火与距离限制相结合的方法——在多肽溶液构象分析中的应用

Distance geometry and molecular dynamics are currently employed in determining molecular structures with interatomic distances from NMR NOESY experiment. Because of the flexibility of peptide, distances obtained from NMR are usually not sufficient to confine its structure. Both distance geometry and molecular dynamics will bias in the conformational space at this circumstance. Constraint Monte Carlo simulated annealing was established to solve this problem. Distance constraints were included into the ECEPP/2 force field by introducing a harmonic energy term. Conformational analysis of a pentapeptide with eight interatomic distances from NMR was carried out as a test. By comparison of the 100 conformers obtained from constraint simulated annealing and the 100 conformers from distance geometry calculation, it was found that constraint simulated annealing can cover the outcomes of distance geometry and at the same time give more con-formers fitting to the experimental data. The result shows that constraint Monte-Carlo simulated annealing is more valid in constructing peptide structures from NMR distances than currently employed methods when no sufficient distances from NMR are available.     

Primary and secondary structure dependence of peptide flexibility assessed by fluorescence-based measurement of end-to-end collision rates

The intrachain fluorescence quenching of the fluorophore 2,3-diazabicyclo[2.2.2]oct-2-ene (DBO) is measured in short peptide fragments, namely the two strands and the turn of the N-terminal beta-hairpin of ubiquitin. The investigated peptides adopt a random-coil conformation in aqueous solution according to CD and NMR experiments. The combination of quenchers with different quenching efficiencies, namely tryptophan and tyrosine, allows the extrapolation of the rate constants for end-to-end collision rates as well as the dissociation of the end-to-end encounter complex. The measured activation energies for fluorescence quenching demonstrate that the end-to-end collision process in peptides is partially controlled by internal friction within the backbone, while measurements in solvents of different viscosities (H2O, D2O, and 7.0 M guanidinium chloride) suggest that solvent friction is an additional important factor in determining the collision rate. The extrapolated end-to-end collision rates, which are only slightly larger than the experimental rates for the DBO/Trp probe/quencher system, provide a measure of the conformational flexibility of the peptide backbone. The chain flexibility is found to be strongly dependent on the type of secondary structure that the peptides represent. The collision rates for peptides derived from the beta-strand motifs (ca. 1 x 10(7) s(-1)) are ca. 4 times slower than that derived from the beta-turn. The results provide further support for the hypothesis that chain flexibility is an important factor in the preorganization of protein fragments during protein folding. Mutations to the beta-turn peptide show that subtle sequence changes strongly affect the flexibility of peptides as well. The protonation and charge status of the peptides, however, are shown to have no significant effect on the flexibility of the investigated peptides. The meaning and definition of end-to-end collision rates in the context of protein folding are critically discussed.     

A simple sheathless CE-MS interface with a sub-micrometer electrical contact fracture for sensitive analysis of peptide and protein samples

Online coupling of capillary electrophoresis (CE) to electrospray ionization mass spectrometry (MS) has shown considerable potential, however, technical challenges have limited its use. In this study, we have developed a simple and sensitive sheathless CE-MS interface based on the novel concept of forming a sub-micrometer fracture directly in the capillary. The simple interface design allowed the generation of a stable ESI spray capable of ionization at low nanoliter flow-rates (45–90 nL/min) for high sensitivity MS analysis of challenging samples like those containing proteins and peptides. By analysis of a model peptide (leucine enkephalin), a limit of detection (LOD) of 0.045 pmol/μL (corresponding to 67 attomol in a sample volume of ∼15 nL) was obtained. The merit of the CE-MS approach was demonstrated by analysis of bovine serum albumin (BSA) tryptic peptides. A well-resolved separation profile was achieved and comparable sequence coverage was obtained by the CE-MS method (73%) compared to a representative UPLC-MS method (77%). The CE-MS interface was subsequently used to analyse a more complex sample of pharmaceutically relevant human proteins including insulin, tissue factor and α-synuclein. Efficient separation and protein ESI mass spectra of adequate quality could be achieved using only a small amount of sample (30 fmol). In addition, analysis of ubiquitin samples under both native and denatured conditions, indicate that the CE-MS setup can facilitate native MS applications to probe the conformational properties of proteins. Thus, the described CE-MS setup should be useful for a wide range of high-sensitivity applications in protein research.     

Why does the silica-binding protein “Si-tag” bind strongly to silica surfaces? Implications of conformational adaptation of the intrinsically disordered polypeptide to solid surfaces

We recently reported that the bacterial 50S ribosomal protein L2 binds strongly to silica surfaces even in the presence of high salt concentrations, detergents, and denaturants such as 8 M urea. We designated L2 as Si-tag, a fusion tag for immobilizing functional proteins on silica materials. Here we discuss the remarkable properties of the Si-tag polypeptide in order to understand the mechanism underlying this binding. Experimental and theoretical studies have shown that the 60-aa N-terminal region and the 71-aa C-terminal region, both of which are rich in positively charged residues, lack a well-defined three-dimensional structure under physiological conditions. This lack of a stable tertiary structure suggests that Si-tag belongs to a family of intrinsically disordered (ID) proteins that exist as dynamic ensembles of rapidly fluctuating structures in aqueous solution. Because of its inherent flexibility, Si-tag could form a large intermolecular interface and optimize its structure for surface interactions by conformational adaptation at the binding interface. Such conformational adaptation occurring concomitantly with binding is common to many ID proteins and is called "coupled folding and binding". Through this conformational adaptation, Si-tag could optimize the interactions between its positively charged side chains and ionized surface silanol groups and between its apolar side chains and hydrophobic surface siloxane sites. The cumulative contribution of these contacts would significantly strengthen the binding of Si-tag, resulting in strong, virtually irreversible binding. Our study suggests that flexible ID proteins have tremendous potential for connecting biomolecules to inorganic materials.     

Use of essential and molecular dynamics to study gammaB-crystallin unfolding after non-enzymic post-translational modifications

Essential and Molecular Dynamics (ED/MD) have been used to model the conformational changes of a protein implicated in a conformational disease--cataract, the largest cause of blindness in the world-after non-enzymic post-translational modification. Cyanate modification did not significantly alter flexibility, while the Schiff's base adduct produced a more flexible N-terminal domain, and intra-secondary structure regions, than either the cyanate adduct or the native structure. Glycation also increased linker flexibility and disrupted the charge network. A number of post-translational adducts showed structural disruption around Cys15 and increased linker flexibility; this may be important in subsequent protein aggregation. Our modelling results are in accord with experimental evidence, and show that ED/MD is a useful tool in modelling conformational changes in proteins implicated in disease processes.     

Solvent induced polymorphism in supramolecular 1,3,5-benzenetribenzoic acid monolayers

This work presents a scanning tunneling microscopy (STM) based study of benzenetribenzoic acid (BTB) monolayer structures at the liquid-solid interface. On graphite(0001) the tailored molecules self-assemble into 2D supramolecular host systems, suitable for the incorporation of other nanoscopic objects. Two crystallographically different BTB structures were found-both hydrogen bonded networks. A specific structure was deliberately selected by solvent identity. One of the BTB polymorphs is a 6-fold chicken-wire structure with circular, approximately 2.8 nm wide cavities. The other structure exhibits an oblique unit cell and a different hydrogen bonding pattern. The large cavity size of the chicken-wire structure was made possible through comparatively strong 2-fold hydrogen bonds between carboxylic groups. In addition, the low conformational flexibility of BTB was supportive to combat the tendency for dense packing.     

X-ray reflectivity study of ultrathin liquid films of diphenylsiloxane-dimethylsiloxane copolymers

Using X-ray reflectivity, we observe drastic differences in the interfacial structure and molecular ordering of diphenylsiloxane-dimethylsiloxane copolymer thin films deposited on hydroxylated versus H-terminated (etched) silicon wafers. We find that substrate type and comonomer ratio determine the conformational arrangements in these liquid films. High-energy bonding between the substrate and the molecules and an increase in rigidity of the molecules due to replacement of methyl groups by phenyl groups leads to a specific molecular ordering at the liquid/solid interface and pronounced density oscillations in this region. The observed structural reorganizations are explained by the interplay and the established balance between the chain flexibility and the polymer-substrate interactions.     

Nematic-amorphous polymer interfaces in the presence of a compatibilizer

We introduce and apply a variant of a dynamic self-consistent field simulation in two dimensions to predict the structure of interfaces between a nematic and an amorphous polymer compatibilized by a diblock copolymer. First, we investigate the effect of the nematic order on the polymer polymer interface without compatibilizer. Then we include the compatibilizer and consider two interfacial setups previously used in experiments, i.e., the bilayer setup and the trilayer setup. In the bilayer setup the diblock copolymer is mixed into the amorphous homopolymer and migrates to the interface in the course of the simulation forming a layered structure. We compare the amount of copolymer at the interface for initial concentrations of the copolymer below and above the critical micelle concentration. In the trilayer setup the initial thickness of the diblock copolymer is varied. The resulting interfacial morphology evolves in the competition between the lamellar structure induced by the interface and a micellar structure, which is intrinsic to the copolymer.     

A CoMFA analysis with conformational propensity: An attempt to analyze the SAR of a set of molecules with different conformational flexibility using a 3D-QSAR method

CoMFA analysis, a widely used 3D-QSAR method, has limitations to handle a set of SAR data containing diverse conformational flexibility since it does not explicitly include the conformational entropic effects into the analysis. Here, we present an attempt to incorporate the conformational entropy effects of a molecule into a 3D-QSAR analysis. Our attempt is based on the assumption that the conformational entropic loss of a ligand upon making a ligand-receptor complex is small if the ligand in an unbound state has a conformational propensity to adopt an active conformation in a complex state. For a QSAR analysis, this assumption was interpreted as follows: a potent ligand should have a higher conformational propensity to adopt an `active-conformation'-like structure in an unbound state than an inactive one. The conformational propensity value was defined as the populational ratio, N<sub>active</sub>/N<sub>stable</sub>, of the number of energetically stable conformers, N<sub>stable</sub>, to the number of `active-conformation'-like structures, N_active. The latter number was calculated by counting the number of conformers that satisfied the structural parameters deduced from the active conformation. A set of SAR data of imidazoleglycerol phosphate dehydratase inhibitors containing 20 molecules with different conformational flexibility was used as a training set for developing a 3D structure-activity relationship by a CoMFA analysis with the conformational propensity value. This resulted in a cross-validated squared correlation coefficient of the CoMFA model with the conformational propensity value (R ² _cross = 0.640) higher than that of the standard CoMFA model (R ² _cross = 0.431). Then we evaluated the quality of the CoMFA models by predicting the inhibitory activity for a new molecule.