Folding of the villin headpiece subdomain from random structures. Analysis of the charge distribution as a function of pH

The structure of the 36 residue villin headpiece subdomain is investigated with the electrostatically driven Monte Carlo method. The ECEPP/3 (Empirical Conformational Energy Program for Peptides) force field, plus two different continuum solvation models, were used to describe the conformational energy of the chain with both blocked and unblocked N and C termini. A statistical analysis of an ensemble of ab initio generated conformations was carried out, based on a comparison with a set of ten native-like structures derived from published experimental data, by using rigid geometry and NMR-derived constraints obtained at pH 3.7. The ten native-like structures satisfy the NMR-derived constraints. The whole ensemble of conformations of the terminally unblocked villin headpiece sub-domain, generated by using ECEPP/3 with a continuum solvation model, were subsequently evaluated at pH 3.7 with a potential function that includes ECEPP/3 combined with a fast multigrid boundary element method. At pH 3.7, the lowest-energy conformation found during the conformational search satisfies approximately 70% of both the distance and the dihedral-angle constraints, and possesses the characteristic packing of three phenylalanine residues that constitute the main part of the hydrophobic core of the molecule. On the other hand, computations at pH 3.7 and pH 7.0 for the ten native-like structures satisfying the NMR-derived constraints indicate a substantial change in the charge distribution for each type of amino acid residue with the change in pH. The results of this study provide a basis to understand the effect of the interactions, such as hydrophobicity, charge-charge interaction and solvent polarization, on the stability of this small alpha-helical protein.     

Effect of subdomain interactions on methyl group dynamics in the hydrophobic core of villin headpiece protein

Thermostable villin headpiece protein (HP67) consists of the N‐terminal subdomain (residues 10–41) and the autonomously folding C‐terminal subdomain (residues 42–76) which pack against each other to form a structure with a unified hydrophobic core. The X‐ray structures of the isolated C‐terminal subdomain (HP36) and its counterpart in HP67 are very similar for the hydrophobic core residues. However, fine rearrangements of the free energy landscape are expected to occur because of the interactions between the two subdomains. We detect and characterize these changes by comparing the µs‐ms time scale dynamics of the methyl‐bearing side chains in isolated HP36 and in HP67. Specifically, we probe three hydrophobic side chains at the interface of the two subdomains (L42, V50, and L75) as well as at two residues far from the interface (L61 and L69). Solid‐state deuteron NMR techniques are combined with computational modeling for the detailed characterization of motional modes in terms of their kinetic and thermodynamic parameters. The effect of interdomain interactions on side chain dynamics is seen for all residues but L75. Thus, changes in dynamics because of subdomain interactions are not confined to the site of perturbation. One of the main results is a two‐ to threefold increase in the value of the activation energies for the rotameric mode of motions in HP67 compared with HP36. Detailed analysis of configurational entropies and heat capacities complement the kinetic view of the degree of the disorder in the folded state.     

On the unyielding hydrophobic core of villin headpiece

Villin headpiece (HP67) is a small, autonomously-folding domain that has become a model system for understanding the fundamental tenets governing protein folding. In this communication, we explore the role that Leu61 plays in the structure and stability of the construct. Deletion of Leu61 results in a completely unfolded protein that cannot be expressed in Escherichia coli. Omission of only the aliphatic leucine side chain (HP67 L61G) perturbed neither the backbone conformation nor the orientation of local hydrophobic side chains. As a result, a large, solvent-exposed hydrophobic pocket, a negative replica of the leucine side-chain, was created on the surface. The loss of the hydrophobic interface between leucine 61 and the hydrophobic pocket destabilized the construct by ~3.3 kcal/mol. Insertion of a single glycine residue immediately before Leu61 (HP67 L61[GL]) was also highly destabilizing and had the effect of altering the backbone conformation (α-helix to π-helix) in order to precisely preserve the wild-type position and conformation of all hydrophobic residues, including Leu61. In addition to demonstrating that the hydrophobic side-chain of Leu61 is critically important for the stability of villin headpiece, our results are consistent with the notion that the precise interactions present within the hydrophobic core, rather than the hydrogen bonds that define the secondary structure, specify a protein's fold.     

Folding and stability of helical bundle proteins from coarse‐grained models

We develop a coarse‐grained model where solvent is considered implicitly, electrostatics are included as short‐range interactions, and side‐chains are coarse‐grained to a single bead. The model depends on three main parameters: hydrophobic, electrostatic, and side‐chain hydrogen bond strength. The parameters are determined by considering three level of approximations and characterizing the folding for three selected proteins (training set). Nine additional proteins (containing up to 126 residues) as well as mutated versions (test set) are folded with the given parameters. In all folding simulations, the initial state is a random coil configuration. Besides the native state, some proteins fold into an additional state differing in the topology (structure of the helical bundle). We discuss the stability of the native states, and compare the dynamics of our model to all atom molecular dynamics simulations as well as some general properties on the interactions governing folding dynamics. Proteins 2013; 81:1200–1211. © 2013 Wiley Periodicals, Inc.     

Determination of nucleic acid by [tetra‐(3‐methoxy‐4‐hydroxyphenyl)]–Tb3+ porphyrin as the fluorescence spectral probe in bis(2‐ethylhexyl)sulfosuccinate sodium salt micelle system

A new system for the determination of nucleic acid by rare earth metallic porphyrin of [tetra‐(3‐methoxy‐4‐hydroxyphenyl)]–Tb³⁺ [T(3‐MO‐4HP)–Tb³⁺] porphyrin as fluorescence spectral probe has been developed in this paper. Nucleic acid can enhance the fluorescence intensity of the T(3‐MO‐4HP)–Tb³⁺ porphyrin in the presence of bis(2‐ethylhexyl)sulfosuccinate sodium salt(AOT) micelle. In pH 8.00 Tris–HCl buffer solution, under optimum conditions, the enhanced fluorescence intensity is in proportion to the concentration of nucleic acids in the range of 0.05–3.00 µg mL^?1 for calf thymus DNA (ct DNA) and 0.03–4.80 µg mL^?1 for fish sperm DNA(fs DNA). Their detection limits are 0.03 and 0.01 µg mL^?1, respectively. In addition, the binding interaction mechanism between T(3‐MO‐4HP)–Tb³⁺ porphyrin and ct DNA is also investigated by resonance scattering and fluorescence spectra. The maximum binding number is calculated by molar ratio method. The new system can be used for the determination of nucleic acid in pig liver, yielding satisfactory results. Copyright © 2009 John Wiley & Sons, Ltd.     

Synthesis of 5‐(4′‐carboxyphenyl)‐10,15,20‐tris‐(4 pyridyl)‐porphyrin and its peptidyl phosphonate derivatives

The synthesis and characterization of three new 4‐pyridyl porphyrin‐peptidyl‐phosphonate compounds, containing a diphenyl 3‐pyridylmethyl‐phosphonate moiety, is described in this article. Nitrogen atoms in the pyridine rings of the obtained compounds were alkylated using methyl iodide, to give additional three, water soluble derivatives of these peptidyl‐porphyrin conjugates. All the synthesized compounds could serve as potential photosensitizers for the photodynamic therapy (PDT) method of tumor therapy and displayed activity as inhibitors of aminopeptidase N. Copyright © 2009 European Peptide Society and John Wiley & Sons, Ltd.     

Reconciling the solution and X-ray structures of the villin headpiece helical subdomain: molecular dynamics simulations and double mutant cycles reveal a stabilizing cation-pi interaction

The 36-residue helical subdomain of the villin headpiece, HP36, is one of the smallest cooperatively folded proteins, folding on the microsecond time scale. The domain is an extraordinarily popular model system for both experimental and computational studies of protein folding. The structure of HP36 has been determined using X-ray crystallography and NMR spectroscopy, with the resulting structures exhibiting differences in helix packing, van der Waals contacts, and hydrogen bonding. It is important to determine the solution structure of HP36 with as much accuracy as possible since this structure is widely used as a reference for simulations and experiments. We complement the existing data by using all-atom molecular dynamics simulations with explicit solvent to evaluate which of the experimental models is the better representation of HP36 in solution. After simulation for 50 ns initiated with the NMR structure, we observed that the protein spontaneously adopts structures with a backbone conformation, core packing, and C-capping motif on the third helix that are more consistent with the crystal structure. We also examined hydrogen bonding and side chain packing interactions between D44 and R55 and between F47 and R55, respectively, which were observed in the crystal structure but not in the NMR-based solution structure. Simulations showed large fluctuations in the distance between D44 and R55, while the distance between F47 and R55 remained stable, suggesting the formation of a cation-pi interaction between those residues. Experimental double mutant cycles confirmed that the F47-R55 pair has a larger energetic coupling than the D44-R55 interaction. Overall, these combined experimental and computational studies show that the X-ray crystal structure is the better reference structure for HP36 in solution at neutral pH. Our analysis also shows how detailed molecular dynamics simulations combined with experimental validation can help bridge the gap between NMR and crystallographic methods.     

Relation between TLR4/NF‐κB signaling pathway activation by 27‐hydroxycholesterol and 4‐hydroxynonenal,and atherosclerotic plaque instability

It is now thought that atherosclerosis, although due to increased plasma lipids, is mainly the consequence of a complicated inflammatory process, with immune responses at the different stages of plaque development. Increasing evidence points to a significant role of Toll‐like receptor 4 (TLR4), a key player in innate immunity, in the pathogenesis of atherosclerosis. This study aimed to determine the effects on TLR4 activation of two reactive oxidized lipids carried by oxidized low‐density lipoproteins, the oxysterol 27‐hydroxycholesterol (27‐OH) and the aldehyde 4‐hydroxynonenal (HNE), both of which accumulate in atherosclerotic plaques and play a key role in the pathogenesis of atherosclerosis. Secondarily, it examined their potential involvement in mediating inflammation and extracellular matrix degradation, the hallmarks of high‐risk atherosclerotic unstable plaques. In human promonocytic U937 cells, both 27‐OH and HNE were found to enhance cell release of IL‐8, IL‐1β, and TNF‐α and to upregulate matrix metalloproteinase‐9 (MMP‐9) via TLR4/NF‐κB‐dependent pathway; these actions may sustain the inflammatory response and matrix degradation that lead to atherosclerotic plaque instability and to their rupture. Using specific antibodies, it was also demonstrated that these inflammatory cytokines increase MMP‐9 upregulation, thus enhancing the release of this matrix‐degrading enzyme by macrophage cells and contributing to plaque instability. These innovative results suggest that, by accumulating in atherosclerotic plaques, the two oxidized lipids may contribute to plaque instability and rupture. They appear to do so by sustaining the release of inflammatory molecules and MMP‐9 by inflammatory and immune cells, for example, macrophages, through activation of TLR4 and its NF‐κB downstream signaling.     

Exploring the effect of silicene monolayer on the structure and function of villin headpiece and amyloid fibrils by molecular dynamics simulations

With the development of various nanomaterial expected to be used in biomedical fields, it is more important to evaluate and understand their potential effects on biological system. In this work, two proteins with different structure, Villin Headpiece (HP35) with α‐helix structure and protofibrils Aβ1‐42 with five β‐strand chains, were selected and their interactions with silicene were studied by means of molecular dynamics (MD) simulation to reveal the potential effect of silicene on the structure and function of biomolecules. The obtained results indicated that silicene could rapidly attract HP35 and Aβ1‐42 fibrils onto the surface to form a stable binding. The adsorption strength was moderate and no significant structural distortion of HP35 and Aβ1‐42 fibrils was observed. Moreover, the strength of calculated the H‐bonds in neighbor chain of Aβ1‐42 fibrils indicated that the mild interactions between silicene and fibrils could regularize the structure of Aβ1‐42 fibrils and stabilize the interactions between five chains of fibrils protein, which might enhance the aggregation of Aβ1‐42 fibrils. This study provides a new insight for understanding the interaction between nanomaterials and biomolecules and moves forward the development of silicene into biomedical fields.     

Enantiopure 4‐oxazolin‐2‐ones and 4‐methylene‐2‐oxazolidinones as chiral building blocks in a divergent asymmetric synthesis of heterocycles

Enantiopure 3‐((R)‐ and 3‐((S)‐1‐phenylethyl)‐4‐oxazoline‐2‐ones were evaluated as chiral building blocks for the divergent construction of heterocycles with stereogenic quaternary centers. The N‐(R)‐ or N‐(S)‐1‐phenylethyl group of these compounds proved to be an efficient chiral auxiliary for the asymmetric induction of the 4‐ and 5‐positions of the 4‐oxazolin‐2‐one ring through thermal and MW‐promoted nucleophilic conjugated addition to Michael acceptors and alkyl halides. The resulting adducts were transformed via a cascade process into fused six‐membered carbo‐ and heterocycles. The structure of the reaction products depended on the electrophiles and reaction conditions used. Alternative isomeric 4‐methylene‐2‐oxazolidinones served as chiral precursors for a versatile and divergent approach to highly substituted cyclic carbamates. DFT quantum calculations showed that the formation of bicyclic pyranyl compounds was generated by a diastereoselective concerted hetero‐Diels‐Alder cycloaddition.     

Folding thermodynamics of β‐hairpins studied by replica‐exchange molecular dynamics simulations

We study the differences in folding stability of β‐hairpin peptides, including GB1 hairpin and a point mutant GB1 K10G, as well as tryptophan zippers (TrpZips): TrpZip1, TrpZip2, TrpZip3‐1, and TrpZip4. By performing replica‐exchange molecular dynamics simulations with Amber03* force field (a modified version of Amber ff03) in explicit solvent, we observe ab initio folding of all the peptides except TrpZip3‐1, which is experimentally known to be the least stable among the peptides studied here. By calculating the free energies of unfolding of the peptides at room temperature and folding midpoint temperatures for thermal unfolding of peptides, we find that TrpZip4 and GB1 K10G peptides are the most stable β‐hairpins followed by TrpZip1, GB1, and TrpZip2 in the given order. Hence, the proposed K10G mutation of GB1 peptide results in enhanced stability compared to wild‐type GB1. An important goal of our study is to test whether simulations with Amber 03* model can reproduce experimentally predicted folding stability differences between these peptides. While the stabilities of GB1 and TrpZip1 yield close agreement with experiment, TrpZip2 is found to be less stable than predicted by experiment. However, as heterogenous folding of TrpZip2 may yield divergent thermodynamic parameters by different spectroscopic methods, mismatching of results with previous experimental values are not conclusive of model shortcomings. For most of the cases, molecular simulations with Amber03* can successfully reproduce experimentally known differences between the mutated peptides, further highlighting the predictive capabilities of current state‐of‐the‐art all‐atom protein force fields. Proteins 2015; 83:1307–1315. © 2015 Wiley Periodicals, Inc.     

An improved synthesis of (2S, 4S)‐ and (2S, 4R)‐2‐amino‐4‐methyldecanoic acids: assignment of the stereochemistry of culicinins

An improved synthesis of (2S, 4S)‐ and (2S, 4R)‐2‐amino‐4‐methyldecanoic acids was accomplished using a glutamate derivative as starting material and Evans' asymmetric alkylation as the decisive step. The NMR data of the two diastereomers were measured and compared with those of the natural product. As a result, the stereochemistry of this novel amino acid unit in culicinins was assigned as (2S, 4R). Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Protein instability during HIC: Evidence of unfolding reversibility,and apparent adsorption strength of disulfide bond‐reduced α‐lactalbumin variants

A two‐conformation, four‐state model has been proposed to describe protein adsorption and unfolding behavior on hydrophobic interaction chromatography (HIC) resins. In this work, we build upon previous study and application of a four‐state model to the effect of salt concentration on the adsorption and unfolding of the model protein α‐lactalbumin in HIC. Contributions to the apparent adsorption strength of the wild‐type protein from native and unfolded conformations, obtained using a deuterium labeling technique, reveal the free energy change and kinetics of unfolding on the resin, and demonstrate that surface unfolding is reversible. Additionally, variants of α‐lactalbumin in which one of the disulfide bonds is reduced were synthesized to examine the effects of conformational stability on apparent retention. Below the melting temperatures of the wild‐type protein and variants, reduction of a single disulfide bond significantly increases the apparent adsorption strength (～6–8 kJ/mol) due to increased instability of the protein. Finally, the four‐state model is used to accurately predict the apparent adsorption strength of a disulfide bond‐reduced variant. Biotechnol. Bioeng. 2009;102: 1416–1427. © 2008 Wiley Periodicals, Inc.     

Temperature‐dependent study reveals that dynamics of hydrophobic residues plays an important functional role in the mitochondrial Tim9–Tim10 complex

Protein–protein interaction is a fundamental process in all major biological processes. The hexameric Tim9–Tim10 (translocase of inner membrane) complex of the mitochondrial intermembrane space plays an essential chaperone‐like role during import of mitochondrial membrane proteins. However, little is known about the functional mechanism of the complex because the interaction is weak and transient. This study investigates how electrostatic and hydrophobic interactions affect the conformation and function of the complex at physiological temperatures, using both experimental and computational methods. The results suggest that, first, different complex conformational states exist at equilibrium, and the major difference between these states is the degree of hydrophobic interactions. Second, the conformational change mimics the biological activity of the complex as measured by substrate binding at the same temperatures. Finally, molecular dynamics simulation and detailed energy decomposition analysis provided supporting evidence at the atomic level for the presence of an excited state of the complex, the formation of which is largely driven by the disruption of hydrophobic interactions. Taken together, this study indicates that the dynamics of the hydrophobic residues plays an important role in regulating the function of the Tim9–Tim10 complex. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

NMR studies on thermal stability of α‐helix conformation of melittin in pure ethanol and ethanol–water mixture solvents

Thermal stability of the α‐helix conformation of melittin in pure ethanol and ethanol–water mixture solvents has been investigated by using NMR spectroscopy. With increase in water concentration of the mixture solvents (from 0 wt% to ~71.5 wt%) as well as temperature (from room temperature to 60 °C), the intramolecular hydrogen bonds formed in melittin are destabilized and the α‐helix is partially uncoiled. Further, the hydrogen bonds are found to be more thermally stable in pure ethanol than in pure methanol, suggesting that their stability is enhanced with increase in the size of the alkyl groups of alcohol molecules. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Highly enantioselective Michael addition of pyrazolin‐5‐ones to nitroolefins catalyzed by cinchona alkaloid derived 4‐methylbenzoylthioureas

Cinchona alkaloid‐derived 4‐methyl/nitro benzoylthioureas were synthesized, which smoothly catalyzed the asymmetric Michael addition of pyrazolin‐5‐ones to nitroolefins. The results showed that electronic effects of substituents in the benzene ring of benzoylthioureas have subtle influences on their catalytic abilities and electron donating methyl group is favored than electron withdrawing nitro group. Preliminary Hartree‐Fock calculations revealed that in the catalytic cycle, hydrogen bond energies of the complex formed with 4‐methylbenzoylthioureas are about 0.19 to 1.56 kcal/mol higher than with the corresponding 4‐nitrobenzoylthioureas. 4‐Methylbenzoylthioureas were identified as the most effective catalysts that promoted asymmetric Michael addition of pyrazolin‐5‐ones to nitroolefins to give the S‐ or R‐products with high enantioselectivities.     

Validated spectrofluorimetric method for the determination of carbamazepine in pharmaceutical dosage forms after reaction with 4‐chloro‐7‐‐nitrobenzo‐2‐oxa‐1,3‐diazole (NBD‐Cl)

A sensitive and simple spectrofluorimetric method has been developed and validated for the determination of the anti‐epileptic drug carbamazepine (CBZ) in its dosage forms. The method was based on a nucleophilic substitution reaction of CBZ with 4‐chloro‐7‐nitrobenzo‐2‐ oxa‐1,3‐diazole (NBD‐Cl) in borate buffer (pH 9) to form a highly fluorescent derivative that was measured at 530 nm after excitation at 460 nm. Factors affecting the formation of the reaction product were studied and optimized, and the reaction mechanism was postulated. The fluorescence–concentration plot is rectilinear over the range of 0.6–8 µg/mL with limit of detection of 0.06 µg/mL and limit of quantitation of 0.19 µg/mL. The method was applied to the analysis of commercial tablets and the results were in good agreement with those obtained using the reference method. Validation of the analytical procedures was evaluated according to ICH guidelines. Copyright © 2015 John Wiley & Sons, Ltd.     

β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

Regulation of the 5‐HT3A Receptor‐Mediated Current by Alkyl 4‐Hydroxybenzoates Isolated from the Seeds of Nelumbo nucifera

Four known alkyl 4‐hydroxybenzoates, i.e., methyl 4‐hydroxybenzoate ( 1 ), ethyl 4‐hydroxybenzoate ( 2 ), propyl 4‐hydroxybenzoate ( 3 ), and butyl 4‐hydroxybenzoate ( 4 ), were isolated from the seeds of Nelumbo nucifera Gaertner (Nymphaeaceae) for the first time. The structures of the isolates were identified by 1D‐ and 2D‐NMR spectroscopy and comparison with published values. The compounds were evaluated for their effects on the 5‐HT‐stimulated inward current (I_5‐HT) mediated by the human 5‐HT₃A receptors expressed in Xenopus oocytes. Compounds 1 and 2 enhanced the I_5‐HT, but 4 reduced it. These results indicate that 4 is an inhibitor of the 5‐HT₃A receptors expressed in Xenopus oocytes.