The role of conformational flexibility on protein supercharging in native electrospray ionization

Effects of covalent intramolecular bonds, either native disulfide bridges or chemical crosslinks, on ESI supercharging of proteins from aqueous solutions were investigated. Chemically modifying cytochrome c with up to seven crosslinks or ubiquitin with up to two crosslinks did not affect the average or maximum charge states of these proteins in the absence of m-nitrobenzyl alcohol (m-NBA), but the extent of supercharging induced by m-NBA increased with decreasing numbers of crosslinks. For the model random coil polypeptide reduced/alkylated RNase A, a decrease in charging with increasing m-NBA concentration attributable to reduced surface tension of the ESI droplet was observed, whereas native RNase A electrosprayed from these same solutions exhibited enhanced charging. The inverse relationship between the extent of supercharging and the number of intramolecular crosslinks for folded proteins, as well as the absence of supercharging for proteins that are random coils in aqueous solution, indicate that conformational restrictions induced by the crosslinks reduce the extent of supercharging. These results provide additional evidence that protein and protein complex supercharging from aqueous solution is primarily due to partial or significant unfolding that occurs as a result of chemical and/or thermal denaturation induced by the supercharging reagent late in the ESI droplet lifetime.     

The valence electronic structure and conformational flexibility of epichlorohydrin

The electronic structure of epichlorohydrin is investigated in the whole valence region by a combined experimental and theoretical study. The issue of controversial assignments of the molecular electronic structure is here addressed. Photoelectron spectra (PES) and Threshold Photoelectron spectra (TPES) of room temperature molecules in the gas phase are recorded. Geometries and energies of the stable conformers due to internal rotation of the C-C-C-Cl dihedral angle, gauche-II (g-II), gauche-I (g-I), and cis, are calculated, and the effect of the conformational flexibility on the photoionization energetics is studied by DFT and 2h-1p Configuration Interaction (CI) methods. Strong breakdown of the Koopmans Theorem (KT) is obtained for the four outermost ionizations, which are further investigated by higher level ab initio calculations. The full assignment of the spectrum is put on a firm basis by the combination of experimental and theoretical results. The orbital composition from correlated calculations is found closer to the DFT orbitals, which are then used to analyze the electronic structure of the molecule. The Highest Occupied Molecular Orbital (HOMO) and HOMO--2 are n(O)/n(Cl) mixed orbitals. The nature of each valence MO is generally preserved in all the conformers, although the magnitude of the n(O)/n(Cl) mixing in HOMO and HOMO--2 varies to some extent with the C-C-C-Cl dihedral angle. The low energy part of the HOMO PE band is predicted to be substantially affected by the conformational flexibility, as experimentally observed in the spectra. The rest of the spectrum is described in terms of the dominant conformer g-II, and a good agreement between experiment and theory is found. The inner-valence PE spectrum is characterized by satellite structures, due to electron correlation effects, which are interpreted by means of 2h-1p CI calculations.     

Aromaticity and conformational flexibility of five-membered monoheterocycles: pyrrole-like and thiophene-like structures

Aromaticity and conformational flexibility of the series of five-membered monoheterocycles with group 14–16 heteroatoms, having one or two lone pairs, were studied with ab initio methods using NICS, ASE and I ₅ indices. For non-planar molecules like phosphole, aromaticity of their planar transition states was also studied, and a special modification of ASE index was proposed to that end. It was found that the presence of two lone pairs is generally preferable for aromaticity of all heterocycles except CPD and silolyl dianions. Heterocycles with group 16 heteroatoms have consistently lower aromaticity compared to other groups. A lot of structures should be classified as moderate aromatic and non-aromatic. Energies of out-of-plane deformation do not correlate with other indices studied, but reveal the same qualitative trends. Generally, aromaticity of five-membered monoheterocycles depends strongly on both heteroatom type and number of lone pairs on it.     

Three-pulse photon echo peak shift spectroscopy as a probe of flexibility and conformational heterogeneity in protein folding

We investigate the equilibrium unfolding of Zn-cytochrome c in guanidine hydrochloride by three-pulse photon echo peak shift (3PEPS) spectroscopy. Unexpectedly, the measurements reveal that inhomogeneous broadening of the sample at the midpoint of the denaturation is larger than that of either native or unfolded states. To interpret this finding, we present simulations of the peak shift for both two-state and three-state unfolding models. Both the denaturant concentration dependence of the asymptotic peak shift (APS) and the wavelength dependence of the APS at the midpoint of the denaturation are different for the two models. Our data are consistent with two-state unfolding.     

Toward functional carboxylate-bridged diiron protein mimics: achieving structural stability and conformational flexibility using a macrocylic ligand framework

A dinucleating macrocycle, H(2)PIM, containing phenoxylimine metal-binding units has been prepared. Reaction of H(2)PIM with [Fe(2)(Mes)(4)] (Mes = 2,4,6-trimethylphenyl) and sterically hindered carboxylic acids, Ph(3)CCO(2)H or Ar(Tol)CO(2)H (2,6-bis(p-tolyl)benzoic acid), afforded complexes [Fe(2)(PIM)(Ph(3)CCO(2))(2)] (1) and [Fe(2)(PIM)(Ar(Tol)CO(2))(2)] (2), respectively. X-ray diffraction studies revealed that these diiron(II) complexes closely mimic the active site structures of the hydroxylase components of bacterial multicomponent monooxygenases (BMMs), particularly the syn disposition of the nitrogen donor atoms and the bridging μ-η(1)η(2) and μ-η(1)η(1) modes of the carboxylate ligands at the diiron(II) centers. Cyclic voltammograms of 1 and 2 displayed quasi-reversible redox couples at +16 and +108 mV vs ferrocene/ferrocenium, respectively. Treatment of 2 with silver perchlorate afforded a silver(I)/iron(III) heterodimetallic complex, [Fe(2)(μ-OH)(2)(ClO(4))(2)(PIM)(Ar(Tol)CO(2))Ag] (3), which was structurally and spectroscopically characterized. Complexes 1 and 2 both react rapidly with dioxygen. Oxygenation of 1 afforded a (μ-hydroxo)diiron(III) complex [Fe(2)(μ-OH)(PIM)(Ph(3)CCO(2))(3)] (4), a hexa(μ-hydroxo)tetrairon(III) complex [Fe(4)(μ-OH)(6)(PIM)(2)(Ph(3)CCO(2))(2)] (5), and an unidentified iron(III) species. Oxygenation of 2 exclusively formed di(carboxylato)diiron(III) compounds, a testimony to the role of the macrocylic ligand in preserving the dinuclear iron center under oxidizing conditions. X-ray crystallographic and (57)Fe M?ssbauer spectroscopic investigations indicated that 2 reacts with dioxygen to give a mixture of (μ-oxo)diiron(III) [Fe(2)(μ-O)(PIM)(Ar(Tol)CO(2))(2)] (6) and di(μ-hydroxo)diiron(III) [Fe(2)(μ-OH)(2)(PIM)(Ar(Tol)CO(2))(2)] (7) units in the same crystal lattice. Compounds 6 and 7 spontaneously convert to a tetrairon(III) complex, [Fe(4)(μ-OH)(6)(PIM)(2)(Ar(Tol)CO(2))(2)] (8), when treated with excess H(2)O.     

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.     

Drug detection based on the conformational changes of calmodulin and the fluorescence of its enhanced green fluorescent protein fusion partner

A homogeneous phase protein-based assay for the high throughput screening of drugs was developed using enhanced green fluorescent protein (EGFP) as the reporter. For that, a fusion protein between calmodulin (CaM) and EGFP was constructed in order to monitor the conformational changes induced in CaM upon binding to tricyclic anti-depressant drugs. In the presence of Ca²⁺, CaM undergoes a conformational change exposing a hydrophobic pocket that interacts with CaM-binding proteins, peptides, and drugs. Further, the conformational changes induced in CaM upon binding to Ca²⁺ and the target analyte drug, leads to a change in the microenvironment of EGFP concomitant with a change in its fluorescence intensity. The observed change in fluorescence intensity of EGFP can be correlated to the concentration of the analyte present in the sample. Further, the response of CaM–EGFP fusion protein in the presence of Ca²⁺ to increasing concentrations of phenothiazines and structurally related tricyclic anti-depressants was investigated. Dose-response curves for various tricyclic anti-depressants were prepared. Moreover, this assay can serve as a model system for other homogeneous binding assays for pharmaceuticals employing genetically fused binding proteins with reporter proteins and may find applications in the high throughput screening of tricyclic anti-depressants.     

Mass spectrometry techniques for detection of ligand-dependent changes in the conformational flexibility of cellular retinol-binding protein type I localized by hydrogen/deuterium exchange

Hydrogen/deuterium exchange, measured by electrospray ionization with orthogonal quadrupole time-of-flight mass spectrometry (ESI-Q-TOFMS) and by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOFMS), was used as a means to probe and map differences in conformational flexibility between the ligand-free and ligand-bound forms of cellular retinol-binding protein type I. Labelled fragments were obtained by digestion of the protein with pepsin. The differences in space-resolved time courses of deuterium incorporation identified regions that exhibit a remarkably higher degree of flexibility in the apo-protein than in the holo-protein. These segments encompass residues that are thought, on the basis of structural homology of the retinol carrier with other members of the intracellular lipid-binding proteins family, to belong to the dynamic portal through which all-trans retinol can access its high-affinity, solvent-shielded, binding site.     

Reducing the conformational flexibility of carbohydrates: locking the 6-hydroxyl group by cyclopropanes

The 6-hydroxyl group of hexopyranosides was stereochemically locked by the spiroannelation of a cyclopropane unit at C-5. The corresponding glucose and mannose derivatives were prepared and their behaviour in glycosidation reactions was studied.     

A molecular dynamics investigation of CDK8/CycC and ligand binding: conformational flexibility and implication in drug discovery

Abnormal activity of cyclin-dependent kinase 8 (CDK8) along with its partner protein cyclin C (CycC) is a common feature of many diseases including colorectal cancer. Using molecular dynamics (MD) simulations, this study determined the dynamics of the CDK8-CycC system and we obtained detailed breakdowns of binding energy contributions for four type-I and five type-II CDK8 inhibitors. We revealed system motions and conformational changes that will affect ligand binding, confirmed the essentialness of CycC for inclusion in future computational studies, and provide guidance in development of CDK8 binders. We employed unbiased all-atom MD simulations for 500 ns on twelve CDK8-CycC systems, including apoproteins and protein–ligand complexes, then performed principal component analysis (PCA) and measured the RMSF of key regions to identify protein dynamics. Binding pocket volume analysis identified conformational changes that accompany ligand binding. Next, H-bond analysis, residue-wise interaction calculations, and MM/PBSA were performed to characterize protein–ligand interactions and find the binding energy. We discovered that CycC is vital for maintaining a proper conformation of CDK8 to facilitate ligand binding and that the system exhibits motion that should be carefully considered in future computational work. Surprisingly, we found that motion of the activation loop did not affect ligand binding. Type-I and type-II ligand binding is driven by van der Waals interactions, but electrostatic energy and entropic penalties affect type-II binding as well. Binding of both ligand types affects protein flexibility. Based on this we provide suggestions for development of tighter-binding CDK8 inhibitors and offer insight that can aid future computational studies.     

The conformational flexibility of nucleic acid bases paired in gas phase: a Car-Parrinello molecular dynamics study

The flexibilities of pyrimidine and imidazole rings in the paired nucleobases are investigated using Car-Parrinello molecular dynamics simulation in gas phase. The pairing influence on the stiffness of rings is analyzed based on the molecular structure of the nucleobases and constraints caused by pairing. We prove that the flexibilities of pyrimidine rings in isolated state have subtle correlation with the degree of aromaticity of the rings. The pairings in nucleic base pairs cause the rings to be more rigid for G, T, and U but more flexible for A and the same for C.     

Dynamic asymmetric organocatalysis: cooperative effects of weak interactions and conformational flexibility in asymmetric organocatalysts

This review describes recent advances by our research group in the design and applications of conformationally flexible guanidine/bisthiourea organocatalysts to facilitate unique bond-forming processes, including retro-free, enantiodivergent, entropy-controlled, and cycle-specific organocatalysis. Special emphasis is placed on the key requirements to attain high selectivity and on functional switching by exploiting the structural flexibility of our organocatalysts.     

Molecular characterization of the 30-AA N-terminal mineral interaction domain of the biomineralization protein AP7

The AP7 protein is one of several mollusk shell proteins which are responsible for aragonite polymorph formation and stabilization within the nacre layer of the Pacific red abalone, H. rufescens. Previously, we demonstrated that the 30-AA N-terminal domain of AP7, denoted as AP7-1, exists as an unfolded sequence and possesses the capability of inhibiting calcium carbonate crystal growth in vitro via growth step frustration or interruption. However, very little is known with regard to the interactive capabilities of this sequence with Ca(II) and with calcium carbonates. Using multidisciplinary techniques, we determine that the AP7-1 polypeptide interacts with Ca(II) ions at the -DD- sequence clusters, yet retains its unfolded, conformationally labile structure in the presence of Ca(II) ions. Further, NMR experiments reveal that the extended structured sequence blocks, -GNGM-, -SVRTQG-, and -ISYL, exhibit motional, chemical exchange, and/or backbone geometry perturbations in response to Ca(II) interactions with AP7-1. Solid-state NMR magic angle spinning studies verify that during the course of in vitro calcium carbonate crystal growth, AP7-1 becomes bound to calcite fragments and cannot be entirely displaced from the mineral fragments using competitive Ca(II) washing. Finally, using a scrambled sequence version of the AP7-1 polypeptide, we observe that sequence scrambling does not adversely affect the crystal growth inhibitory activity of AP7-1, suggesting that the amino acid composition of AP7-1 may be more critical to growth step inhibition than the linear ordering of amino acids.     

Protein dynamics and function: Making new strides with an old warhorse,the alkaline conformational transition of cytochrome c

Protein dynamics is intimately linked to function. In metalloproteins, dynamics are often coupled to redox activity, ligand binding and enzyme function. We provide a concise overview of the field and then focus on the use of the alkaline conformer of cytochrome c as a model system to probe the factors that control the conformational dynamics of proteins in general and metalloproteins in particular. We consider the effects of ligands on metal-mediated dynamics, the interplay between intrinsic metal-ligand dynamics and barriers imposed by the protein scaffold itself, and the effects of local and overall protein stability on dynamics. Discussed within are the collected results from equilibrium thermodynamic methods, pH jump kinetics and conformationally gated redox reactions between small inorganic reagents and metalloproteins used as a means to probe conformational switching in metalloproteins.     

Tunneling and conformational flexibility play critical roles in the isomerization mechanism of vitamin D

The thermal isomerization reaction converting previtamin D to vitamin D is an intramolecular [1,7]-sigmatropic hydrogen shift with antarafacial stereochemistry. We have studied the dynamics of this reaction by means of the variational transition-state theory with multidimensional corrections for tunneling in both gas-phase and n-hexane environments. Two issues that may have important effects on the dynamics were analyzed in depth, i.e., the conformations of previtamin D and the quantum effects associated with the hydrogen-transfer reaction. Of the large number of conformers of previtamin D that were located, there are 16 that have the right disposition to react. The transition-state structures associated with these reaction paths are very close in energy, so all of them should be taken into account for an accurate calculation of both the thermal rate constants and the kinetic isotope effects. This issue is particularly important because the contribution of each of the reaction paths to the total thermal rate constant is quite sensitive to the environment. The dynamics results confirm that tunneling plays an important role and that model systems that were considered previously to study the hydrogen shift reaction cannot mimic the complexity introduced by the flexibility of the rings of previtamin D. Finally, the characterization of the conformers of both previtamin D and vitamin D allowed the calculation of the thermal equilibrium constants of the isomerization process.     

The N-terminal domain of apolipoprotein B-100: structural characterization by homology modeling

Background  

Apolipoprotein B-100 (apo B-100) stands as one of the largest proteins in humans. Its large size of 4536 amino acids hampers the production of X-ray diffraction quality crystals and hinders in-solution NMR analysis, and thus necessitates a domain-based approach for the structural characterization of the multi-domain full-length apo B.