Hydration and dewetting near graphite-CH(3) and graphite-COOH plates

The dynamics of water near the nanoscale hydrophobic (graphite-CH(3)) and hydrophilic (graphite-COOH) plates has been studied in detail with molecular dynamics simulations in this paper. It is shown that these designed surfaces (by growing a layer of methyl or carboxyl groups on top of graphite) can have a significant impact on the neighboring water dynamics, with the hydrophilic carboxyl surface having even more profound effects. The water hydrogen bond lifetime is much longer near both types of surfaces than that in the bulk, while on the other hand the water diffusion constant is much smaller than that in the bulk. The difference in the diffusion constant can be as large as a factor of 8 and the difference in the hydrogen bond lifetime can be as large as a factor of 2, depending on the distance from the surface. Furthermore, the water molecules in the first solvation shell of surface atoms show a strong bias in hydroxyl group orientation near the surface, confirming some of the previous findings. Finally, the possible water dewetting transition between two graphite-CH(3) plates and the effect of the strength of the solute-solvent attractions on the water drying transition are investigated. The relationship among the dewetting transition critical distance, van der Waals potential well depth, and water contact angle on the graphite-CH(3) surface is also analyzed on the basis of a simple macroscopic theory, which can be used to predict the dewetting transition critical distance.     

Aromatic-aromatic interactions in proteins: beyond the dimer

Aromatic residues are key widespread elements of protein structures and have been shown to be important for structure stability, folding, protein-protein recognition, and ligand binding. The interactions of pairs of aromatic residues (aromatic dimers) have been extensively studied in protein structures. Isolated aromatic molecules tend to form higher order clusters, like trimers, tetramers, and pentamers, that adopt particular well-defined structures. Taking this into account, we have surveyed protein structures deposited in the Protein Data Bank in order to find clusters of aromatic residues in proteins larger than dimers and characterized them. Our results show that larger clusters are found in one of every two unique proteins crystallized so far, that the clusters are built adopting the same trimer motifs found for benzene clusters in vacuum, and that they are clearly nonlocal brining primary structure distant sites together. We extensively analyze the trimers and tetramers conformations and found two main cluster types: a symmetric cluster and an extended ladder. Finally, using calmodulin as a test case, we show aromatic clsuters possible role in folding and protein-protein interactions. All together, our study highlights the relevance of aromatic clusters beyond the dimer in protein function, stability, and ligand recognition.     

Gas-phase H/D exchange and collision cross sections of hemoglobin monomers, dimers, and tetramers

The conformations of gas-phase ions of hemoglobin, and its dimer and monomer subunits have been studied with H/D exchange and cross section measurements. During the H/D exchange measurements, tetramers undergo slow dissociation to dimers, and dimers to monomers, but this did not prevent drawing conclusions about the relative exchange levels of monomers, dimers, and tetramers. Assembly of the monomers into tetramers, hexamers, and octamers causes the monomers to exchange a greater fraction of their hydrogens. Dimer ions, however, exchange a lower fraction of their hydrogens than monomers or tetramers. Solvation of tetramers affects the exchange kinetics. Solvation molecules do not appear to exchange, and solvation lowers the overall exchange level of the tetramers. Cross section measurements show that monomer ions in low charge states, and tetramer ions have compact structures, comparable in size to the native conformations in solution. Dimers have remarkably compact structures, considerably smaller than the native conformation in solution and smaller than might be expected from the monomer or tetramer cross sections. This is consistent with the relatively low level of exchange of the dimers.     

Synthesis and Chemistry of Bicyclo[4.1.0]hept-1,6-ene

Bicyclo[4.1.0]hept-1,6-ene has been generated by elimination of 1-chloro-2-(trimethysilyl)bicyclo[4.1.0]heptane in the gas phase over solid fluoride at 25 degrees C. The cyclopropene dimerizes by a rapid ene reaction forming two diastereomeric cyclopropenes. In tetrahydrofuran or chloroform the ene dimers couple to form a single crystalline triene tetramer, whereas a mixture of tricyclohexane tetramers is formed when the neat dimers are allowed to warm to room temperature. Oxidation by dimethyldioxirane or dioxygen gives carbonyl products. Quantum mechanical calculations yielded an increase in strain of approximately 17 kcal/mol over that for 1,2-dimethylcyclopropene. The potential enegy barrier to flexing (folding) along the fused double bond of bicyclo[4.1.0]hept-1,6-ene is only approximately 1 kcal/mol at the highest level of theory investigated.     

Computational study on the structural diversity of amyloid Beta Peptide (abeta(10-35)) oligomers

We studied the oligomerization of Alzheimer amyloid beta peptide (Abeta) using a replica exchange molecular dynamics (REMD) simulation. The simulation was performed with Abeta(10-35) dimers, trimers, and tetramers. Extensive REMD simulations illustrated several possible oligomer conformations. As the size of the oligomer increased from a dimer to a tetramer, the number of possible configurations was reduced. We identified all the possible conformations for each oligomer and characterized their temperature dependence. It was found that the detailed structures of the oligomers, which may act as folding intermediates, are highly sensitive to the parameters of the simulation environment such as temperature and concentration. Structural diversities of Abeta oligomers suggest multiple pathways of the aggregation process.     

On the manifestation of hydrophobicity at the nanoscale

The manifestation of hydrophobicity at the nanoscale has been shown to depend on the topology of the solute. Using various nanoscopic hydrophobic plates, molecular dynamics simulation has been employed to explore the hydration and dewetting at the nanoscale. The topology of the solute regulates the behavior of nanoconfined water, resulting in any of the wet, dry, and intermittent wet-dry intersolute states. The present result reconciles apparently contrasting literature reports on how water behaves at extended hydrophobic surfaces and sheds light on the mechanism of dewetting.     

Ab initio folding of helix bundle proteins using molecular dynamics simulations

We have demonstrated that ab initio fast folding simulations at 400 K using a GB implicit solvent model with an all-atom based force field can describe the spontaneous formation of nativelike structures for the 36-residue villin headpiece and the 46-residue fragment B of Staphylococcal protein A. An implicit solvent model combined with high-temperature MD makes it possible to perform direct folding simulations of small- to medium-sized proteins by reducing the computational requirements tremendously. In the early stage of folding of the villin headpiece and protein A, initial hydrophobic collapse and rapid formation of helices were found to play important roles. For protein A, the third helix forms first in the early stage of folding and exhibits higher stability. The free energy profiles calculated from the folding simulations suggested that both of the helix-bundle proteins show a two-state thermodynamic behavior and protein A exhibits rather broad native basins.     

Position‐Dependent Effects of Fluorinated Amino Acids on the Hydrophobic Core Formation of a Heterodimeric Coiled Coil

Systematic model investigations of the molecular interactions of fluorinated amino acids within native protein environments substantially improve our understanding of the unique properties of these building blocks. A rationally designed heterodimeric coiled coil peptide (VPE/VPK) and nine variants containing amino acids with variable fluorine content in either position a16 or d19 within the hydrophobic core were synthesized and used to evaluate the impact of fluorinated amino acid substitutions within different hydrophobic protein microenvironments. The structural and thermodynamic stability of the dimers were examined by applying both experimental (CD spectroscopy, FRET, and analytical ultracentrifugation) and theoretical (MD simulations and MM‐PBSA free energy calculations) methods. The coiled coil environment imposes position‐dependent conformations onto the fluorinated side chains and thus affects their packing and relative orientation towards their native interaction partners. We find evidence that such packing effects exert a significant influence on the contribution of fluorine‐induced polarity to coiled coil folding.     

Enhanced Wang Landau sampling of adsorbed protein conformations

Using computer simulations to model the folding of proteins into their native states is computationally expensive due to the extraordinarily low degeneracy of the ground state. In this paper, we develop an efficient way to sample these folded conformations using Wang Landau sampling coupled with the configurational bias method (which uses an unphysical "temperature" that lies between the collapse and folding transition temperatures of the protein). This method speeds up the folding process by roughly an order of magnitude over existing algorithms for the sequences studied. We apply this method to study the adsorption of intrinsically disordered hydrophobic polar protein fragments on a hydrophobic surface. We find that these fragments, which are unstructured in the bulk, acquire secondary structure upon adsorption onto a strong hydrophobic surface. Apparently, the presence of a hydrophobic surface allows these random coil fragments to fold by providing hydrophobic contacts that were lost in protein fragmentation.     

Dewetting of polystyrene thin films on poly(ethylene glycol)-modified surfaces as a simple approach for patterning proteins

A simple technique for patterning proteins utilizing dewetted polystyrene (PS) droplets is demonstrated. A polystyrene thin film was spin coated on a poly(ethylene glycol) (PEG) silane-modified surface. As the PS film dewets from the surface, upon annealing, to form droplets, the PEG-silane-modified surface is exposed, which retains its capability to resist protein adsorption, and the PS droplets allow the selective adsorption of proteins. In contrast to the undewetted flat PS film, the droplet surface had a greater amount of adsorbed proteins. Atomic force microscopy scans reveal that the roughness of the droplet surface is higher, and a multilayer of proteins results on the droplet surface. Moreover, micro- and nanoscale droplet patterns can easily be achieved by tuning the thickness of PS thin films. Because dewetting approaches for generating ordered dewetting droplets have been successfully generated by others, those approaches could be easily combined with this technique to fabricate ordered protein patterns.     

Role of water in mediating the assembly of Alzheimer amyloid-beta Abeta16-22 protofilaments

The role of water in promoting the formation of protofilaments (the basic building blocks of amyloid fibrils) is investigated using fully atomic molecular dynamics simulations. Our model protofilament consists of two parallel beta-sheets of Alzheimer Amyloid-beta 16-22 peptides (Ac-K(16)-L(17)-V(18)-F(19)-F(20)-A(21)-E(22)-NH2). Each sheet presents a distinct hydrophobic and hydrophilic face and together self-assemble to a stable protofilament with a core consisting of purely hydrophobic residues (L(17), F(19), A(21)), with the two charged residues (K(16), E(22)) pointing to the solvent. Our simulations reveal a subtle interplay between a water mediated assembly and one driven by favorable energetic interactions between specific residues forming the interior of the protofilament. A dewetting transition, in which water expulsion precedes hydrophobic collapse, is observed for some, but not all molecular dynamics trajectories. In the trajectories in which no dewetting is observed, water expulsion and hydrophobic collapse occur simultaneously, with protofilament assembly driven by direct interactions between the hydrophobic side chains of the peptides (particularly between F-F residues). For those same trajectories, a small increase in the temperature of the simulation (on the order of 20 K) or a modest reduction in the peptide-water van der Waals attraction (on the order of 10%) is sufficient to induce a dewetting transition, suggesting that the existence of a dewetting transition in simulation might be sensitive to the details of the force field parametrization.     

The native ensemble and folding of a protein molten-globule: functional consequence of downhill folding

The continually emerging functional significance of intrinsic disorder and conformational flexibility in proteins has challenged the long-standing dogma of a well-defined structure contributing to a specific function. Molten-globular states, a class of proteins with significant secondary-structure but a fluid hydrophobic core, is one such example. They have however been difficult to characterize due to the complexity of experimental data and lack of computational avenues. Here, we dissect the folding mechanism of the α-helical molten-globular protein NCBD from three fundamentally different approaches: statistical-mechanical variable barrier model, C(α)-based Gō-model and explicit water all-atom molecular dynamics (MD) simulations. We find that NCBD displays the characteristics of a one-state globally downhill folder but is significantly destabilized. Using simulation techniques, we generate a highly constrained but a heterogeneous native ensemble of the molten-globule for the first time that is consistent with experimental data including small angle X-ray scattering (SAXS), circular dichroism (CD), and nuclear magnetic resonance (NMR). The resulting native ensemble populates conformations reported in other bound-forms providing direct evidence to the mechanism of conformational selection for binding multiple partners in this domain. Importantly, our simulations reveal a connection between downhill folding and large conformational flexibility in this domain that has been evolutionarily selected and functionally exploited resulting in large binding promiscuity. Finally, the multimodel approach we employ here serves as a powerful methodology to study mechanisms and suggests that the thermodynamic features of molten-globules fall within the array of folding mechanisms available to small single-domain proteins.     

Fast folding of a four-helical bundle protein

The FK506-FKBP12 binding-domain of the kinase FRAP (FRB) forms a classic up-down four-helical bundle. The folding pathway of this protein has been investigated using a combination of equilibrium and kinetic studies. The native state of the protein is stable with respect to the unfolded state by some 7 kcal mol(-1) at pH 6.0, 10 degrees C. A kinetic analysis of unfolding and refolding rate constants as a function of chemical denaturant concentration suggests that an intermediate state may be populated during folding at low concentrations of denaturant. The presence of this intermediate state is confirmed by refolding experiments performed in the presence of the hydrophobic dye 8-anilinonaphthalene-1 sulfonate (ANS). ANS binds to the partially folded intermediate state populated during the folding of FRB and undergoes a large change in fluorescence that can be detected using stopped-flow techniques. Analysis of the kinetic data suggests that the intermediate state is compact and it may even be a misfolded species that has to partially unfold before it can reach the transition state. Folding and unfolding rate constants in water are approximately 150-200 s(-1) and 0.005-0.06 s(-1), respectively, at neutral pH and 10 degrees C. The folding of FRB is somewhat slower than for other all-helical proteins, probably as a consequence of the formation of a metastable intermediate state. The folding rate constant in the absence of any populated intermediate can be estimated to be 8800 s(-1). Despite the presence of an intermediate state, which effectively slows folding, the protein still folds rapidly with a half-life of 5 ms at 10 degrees C. The dependence of the rate constants on denaturant concentration indicates that the transition state for folding is compact with some 80% of the surface area exposed in the unfolded state buried in the transition state. Data presented for FRB is compared with kinetic data obtained for other all-helical proteins.     

Mechanism of acylation of lithium phenylacetylide with a Weinreb amide

Additions of lithium phenylacetylide to a Weinreb amide are described. Dimeric lithium acetylide reacts via a monosolvated monomer-based transition structure. The robust tetrahedral intermediate forms sequentially a C(1) 2:2 mixed tetramer with the excess lithium acetylide and a 1:3 (alkoxide-rich) mixed tetramer. The stabilities of the mixed tetramers are consistent with a pronounced autoinhibition.     

Conformational selection or induced fit for Brinker and DNA recognition

Brinker is the key target protein of the Drosophila Decapentaplegic morphogen signalling pathway. Brinker is widely expressed and can bind with DNA. NMR spectra suggest that apo-Brinker is intrinsically unstructured and undergoes a folding transition upon DNA-binding. However, the coupled mechanism of binding and folding is poorly understood. Here, we performed molecular dynamics (MD) simulations for both bound and apo-Brinker to study the mechanism. Room-temperature MD simulations suggest that Brinker becomes more rigid and stable upon DNA-binding. Kinetic analysis of high-temperature MD simulations shows that both bound and apo-Brinker unfold via a two-state process. The time scale of tertiary unfolding is significantly different between bound and apo-Brinker. The predicted Φ-values suggest that there are more residues with native-like transition state ensembles (TSEs) for bound Brinker than for apo-Brinker. The average RMSD differences between bound and apo-Brinker and Kolmogorov-Smirnov (KS) test analysis illustrate that Brinker folding upon DNA-binding might obey induced-fit mechanism based on MD simulations. These methods can be used for the research of other biomolecular folding upon ligand-binding.     

Chemical denaturants inhibit the onset of dewetting

The mechanism by which the aqueous cosolvents guanidinium chloride and urea denature proteins is a matter of controversy. Here, we use all-atom molecular dynamics simulations to study the effect of both denaturants on the dewetting of water confined between nanoseparated hydrophobic plates. It is found that the denaturants inhibit the onset of dewetting, so that it occurs at shorter interplate distances than in pure water. Our results support a role for urea and guanidinium in assisting in the solvation of nonpolar surfaces, thereby weakening hydrophobic effects known to be important for protein stability.     

二维紧密蛋白质链折叠序列的特性研究

采用完全计数法,研究了二维紧密蛋白质链在不同HP序列时的构象性质,特别是具有唯一基态能量的折叠序列的性质.对于具有N个单体的紧密蛋白质链,发现有一定比例的序列为折叠序列.在这些折叠序列中,疏水基团(H)的数目比亲水基团(P)多20%,并同200种真实蛋白质分子的疏水基团和亲水基团的结果进行了比较.对于不同的折叠序列,根据序列中其疏水基团的数目,把具有相同疏水基团数目的序列归在同一类,发现这样的序列在总的序列中的相对含量满足高斯分布.同时还对序列中H(或者P)团族大小及其分别进行了研究,发现折叠序列与无规随机序列不同.还研究了不同折叠序列在不同链长时的比热情况,发现其相转变温度TC主要与链长有关,与折叠序列无关.     

Sterically and guest-controlled self-assembly of calix[4]arene derivatives

In solvents such as chloroform or benzene, tetraurea calix[4]arenes 1 form dimeric capsules in which one solvent molecule is usually included as guest. To explore the structural requirements for the formation of such hydrogen-bonded dimers we replaced one p-tolylurea residue by a simple acetamide function. The resulting calix[4]arene 2 a, substituted at its wide rim with one acetamide and three p-tolylurea functions, assumes a C(1)-symmetrical conformation in apolar solvents as shown by (1)H NMR, which is not compatible with the usual capsule. In the crystalline state, four molecules of 2 a, adopting a pinched cone conformation, assemble into a quasi S(4)-symmetrical tetramer stabilized by a cyclic array of 24 NH.O==C hydrogen bonds and four NH.pi interactions. Four acetamide groups are hydrogen-bonded to each other and pack tightly in the center of the assembly. All polar residues are buried inside the tetramer, the surface of which is lipophilic. Extensive NMR studies revealed similar structures in apolar solvents such as [D]chloroform or [D(6)]benzene for calixacetamides 2 a-c. The formation of these tetramers in solution is critically dependent on the size of the amide fragment, so that propionamide 2 d, butyramide 2 e, and p-tolylamide 2 f form only ill-defined aggregates. This is caused by steric crowding inside the tetrameric assembly. The tetramers persist during molecular dynamics simulations, and the optimized average structure of the MD run is similar to that found in the crystalline state. Theoretical studies revealed that cooperation of hydrogen bonds with multiple NH.pi, C--H.pi, and pi.pi attractions make the tetramer more stable than the capsular dimer with the solvent as guest. In the presence of tetraethylammonium salts, however, compounds 2 a-e form dimeric capsular assemblies, each incorporating a single ammonium cation. Only one of two possible regioisomeric dimers is formed, in which both acetamide groups are surrounded by two urea residues. These examples give striking evidence of how self-assembly in solution can be strongly dependent on subtle structural factors and of how the formation of dimeric capsules can be induced by the presence of an appropriate guest.     

Does the most stable formic acid tetramer have pi stacking or C-H...O interactions?

Density functional theory (DFT), Moller-Plesset (MP) perturbation theory, and coupled-cluster calculations are used to examine low-energy minima on the potential energy surface of the formic acid tetramer (HCOOH)(4). The potential energy surface is rather flat with respect to rotation of one of the dimers, relative to the other dimer in an aligned stack, about the axis passing through the inversion centers of the dimers. Our best calculations suggest that an aligned pi-pi stack of two dimers is very likely to be the global minimum but there are two other pi-pi stacks within 0.5 kcal /mol. Moreover, a fourth pi-pi stack, a planar association of two dimers held together by C-H...O interactions, and a bowl structure all lie within 1 kcal /mol of the lowest-energy structure.     

Structure, energetics, and reactions of alkali tetramers

Electronic structure calculations have been carried out for all possible alkali tetramers that can be formed from X(2) + X(2) → X(2)X(2), X(2) + Y(2) → X(2)Y(2), and XY + XY → X(2)Y(2) alkali dimer association reactions. Vibrationally stable rhombic (D(2h)) and planar (C(s)) structures are found for all possible tetramers formed from the alkali metals, Li to Cs. All tetramer formation reactions (from ground state singlet homonuclear or heteronuclear dimers) are found to be exothermic with binding energies ranging from 6282 cm(-1) for Li(2)Li(2) to 1985 cm(-1) for Cs(2)Cs(2). Extensive calculations, carried out at long-range for several reactant pairs, indicate that there are barrier-less pathways for the formation of tetramers from dimer association reactions. At low temperatures, direct formation of tetramers is unlikely, owing to the large exothermicity associated with these association reactions, but atom exchange reactions (X(2) + Y(2) ? XY + XY) are possible for some species.