Nitriles at Silica Interfaces Resemble Supported Lipid Bilayers.

Nitriles are important solvents not just for bulk reactions but also for interfacial processes such as separations, heterogeneous catalysis, and electrochemistry. Although nitriles have a polar end and a lipophilic end, the cyano group is not hydrophilic enough for these substances to be thought of as prototypical amphiphiles. This picture is now changing, as research is revealing that at a silica surface nitriles can organize into structures that, in many ways, resemble lipid bilayers. This unexpected organization may be a key component of unique interfacial behavior of nitriles that make them the solvents of choice for so many applications. The first hints of this lipid-bilayer-like (LBL) organization of nitriles at silica interfaces came from optical Kerr effect (OKE) experiments on liquid acetonitrile confined in the pores of sol-gel glasses. The orientational dynamics revealed by OKE spectroscopy suggested that the confined liquid is composed of a relatively immobile sublayer of molecules that accept hydrogen bonds from the surface silanol groups and an interdigitated, antiparallel layer that is capable of exchanging into the centers of the pores. This picture of acetonitrile has been borne out by molecular dynamics simulations and vibrational sum-frequency generation (VSFG) experiments. Remarkably, these simulations further indicate that the LBL organization is repeated with increasing disorder at least 20 Å into the liquid from a flat silica surface. Simulations and VSFG and OKE experiments indicate that extending the alkyl chain to an ethyl group leads to the formation of even more tightly packed LBL organization featuring entangled alkyl tails. When the alkyl portion of the molecule is a bulky t-butyl group, packing constraints prevent well-ordered LBL organization of the liquid. In each case, the surface-induced organization of the liquid is reflected in its interfacial dynamics. Acetonitrile/water mixtures are favored solvent systems for separations technologies such as hydrophilic interaction chromatography. Simulations had suggested that although a monolayer of water partitions to the silica surface in such mixtures, acetonitrile tends to associate with this monolayer. VSFG experiments reveal that, even at high water mole fractions, patches of well-ordered acetonitrile bilayers remain at the silica surface. Due to its ability to donate and accept hydrogen bonds, methanol also partitions to a silica surface in acetonitrile/methanol mixtures and can serve to take the place of acetonitrile in the sublayer closest to the surface. These studies reveal that liquid nitriles can exhibit an unexpected wealth of new organizational and dynamic behaviors at silica surfaces, and presumably at the surfaces of other chemically important materials as well. This behavior cannot be predicted from the bulk organization of these liquids. Our new understanding of the interfacial behavior of these liquids will have important implications for optimizing a wide range of chemical processes in nitrile solvents.

How osmolytes influence hydrophobic polymer conformations: A unified view from experiment and theory.

It is currently the consensus belief that protective osmolytes such as trimethylamine N-oxide (TMAO) favor protein folding by being excluded from the vicinity of a protein, whereas denaturing osmolytes such as urea lead to protein unfolding by strongly binding to the surface. Despite there being consensus on how TMAO and urea affect proteins as a whole, very little is known as to their effects on the individual mechanisms responsible for protein structure formation, especially hydrophobic association. In the present study, we use single-molecule atomic force microscopy and molecular dynamics simulations to investigate the effects of TMAO and urea on the unfolding of the hydrophobic homopolymer polystyrene. Incorporated with interfacial energy measurements, our results show that TMAO and urea act on polystyrene as a protectant and a denaturant, respectively, while complying with Tanford-Wyman preferential binding theory. We provide a molecular explanation suggesting that TMAO molecules have a greater thermodynamic binding affinity with the collapsed conformation of polystyrene than with the extended conformation, while the reverse is true for urea molecules. Results presented here from both experiment and simulation are in line with earlier predictions on a model Lennard-Jones polymer while also demonstrating the distinction in the mechanism of osmolyte action between protein and hydrophobic polymer. This marks, to our knowledge, the first experimental observation of TMAO-induced hydrophobic collapse in a ternary aqueous system.

Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field.

Designing tight-binding ligands is a primary objective of small-molecule drug discovery. Over the past few decades, free-energy calculations have benefited from improved force fields and sampling algorithms, as well as the advent of low-cost parallel computing. However, it has proven to be challenging to reliably achieve the level of accuracy that would be needed to guide lead optimization (∼5× in binding affinity) for a wide range of ligands and protein targets. Not surprisingly, widespread commercial application of free-energy simulations has been limited due to the lack of large-scale validation coupled with the technical challenges traditionally associated with running these types of calculations. Here, we report an approach that achieves an unprecedented level of accuracy across a broad range of target classes and ligands, with retrospective results encompassing 200 ligands and a wide variety of chemical perturbations, many of which involve significant changes in ligand chemical structures. In addition, we have applied the method in prospective drug discovery projects and found a significant improvement in the quality of the compounds synthesized that have been predicted to be potent. Compounds predicted to be potent by this approach have a substantial reduction in false positives relative to compounds synthesized on the basis of other computational or medicinal chemistry approaches. Furthermore, the results are consistent with those obtained from our retrospective studies, demonstrating the robustness and broad range of applicability of this approach, which can be used to drive decisions in lead optimization.

How force unfolding differs from chemical denaturation.

Single-molecule force spectroscopies are remarkable tools for studying protein folding and unfolding, but force unfolding explores protein configurations that are potentially very different from the ones traditionally explored in chemical or thermal denaturation. Understanding these differences is crucial because such configurations serve as starting points of folding studies, and thus can affect both the folding mechanism and the kinetics. Here we provide a detailed comparison of both chemically induced and force-induced unfolded state ensembles of ubiquitin based on extensive, all-atom simulations of the protein either extended by force or denatured by urea. As expected, the respective unfolded states are very different on a macromolecular scale, being fully extended under force with no contacts and partially extended in urea with many nonnative contacts. The amount of residual secondary structure also differs: A significant population of α-helices is found in chemically denatured configurations but such helices are absent under force, except at the lowest applied force of 30 pN where short helices form transiently. We see that typical-size helices are unstable above this force, and ß-sheets cannot form. More surprisingly, we observe striking differences in the backbone dihedral angle distributions for the protein unfolded under force and the one unfolded by denaturant. A simple model based on the dialanine peptide is shown to not only provide an explanation for these striking differences but also illustrates how the force dependence of the protein dihedral angle distributions give rise to the worm-like chain behavior of the chain upon force.

Rate limit of protein elastic response is tether dependent.

The elastic restoring force of tissues must be able to operate over the very wide range of loading rates experienced by living organisms. It is surprising that even the fastest events involving animal muscle tissues do not surpass a few hundred hertz. We propose that this limit is set in part by the elastic dynamics of tethered proteins extending and relaxing under a changing load. Here we study the elastic dynamics of tethered proteins using a fast force spectrometer with sub-millisecond time resolution, combined with Brownian and Molecular Dynamics simulations. We show that the act of tethering a polypeptide to an object, an inseparable part of protein elasticity in vivo and in experimental setups, greatly reduces the attempt frequency with which the protein samples its free energy. Indeed, our data shows that a tethered polypeptide can traverse its free-energy landscape with a surprisingly low effective diffusion coefficient D(eff) ~ 1,200 nm(2)/s. By contrast, our Molecular Dynamics simulations show that diffusion of an isolated protein under force occurs at D(eff) ~ 10(8) nm(2)/s. This discrepancy is attributed to the drag force caused by the tethering object. From the physiological time scales of tissue elasticity, we calculate that tethered elastic proteins equilibrate in vivo with D(eff) ~ 10(4)-10(6) nm(2)/s which is two to four orders magnitude smaller than the values measured for untethered proteins in bulk.

Probing static disorder in Arrhenius kinetics by single-molecule force spectroscopy.

The widely used Arrhenius equation describes the kinetics of simple two-state reactions, with the implicit assumption of a single transition state with a well-defined activation energy barrier DeltaE, as the rate-limiting step. However, it has become increasingly clear that the saddle point of the free-energy surface in most reactions is populated by ensembles of conformations, leading to nonexponential kinetics. Here we present a theory that generalizes the Arrhenius equation to include static disorder of conformational degrees of freedom as a function of an external perturbation to fully account for a diverse set of transition states. The effect of a perturbation on static disorder is best examined at the single-molecule level. Here we use force-clamp spectroscopy to study the nonexponential kinetics of single ubiquitin proteins unfolding under force. We find that the measured variance in DeltaE shows both force-dependent and independent components, where the force-dependent component scales with F(2), in excellent agreement with our theory. Our study illustrates a novel adaptation of the classical Arrhenius equation that accounts for the microscopic origins of nonexponential kinetics, which are essential in understanding the rapidly growing body of single-molecule data.

Diversity of chemical mechanisms in thioredoxin catalysis revealed by single-molecule force spectroscopy.

Thioredoxins (Trxs) are oxidoreductase enzymes, present in all organisms, that catalyze the reduction of disulfide bonds in proteins. By applying a calibrated force to a substrate disulfide, the chemical mechanisms of Trx catalysis can be examined in detail at the single-molecule level. Here we use single-molecule force-clamp spectroscopy to explore the chemical evolution of Trx catalysis by probing the chemistry of eight different Trx enzymes. All Trxs show a characteristic Michaelis-Menten mechanism that is detected when the disulfide bond is stretched at low forces, but at high forces, two different chemical behaviors distinguish bacterial-origin from eukaryotic-origin Trxs. Eukaryotic-origin Trxs reduce disulfide bonds through a single-electron transfer reaction (SET), whereas bacterial-origin Trxs show both nucleophilic substitution (S(N)2) and SET reactions. A computational analysis of Trx structures identifies the evolution of the binding groove as an important factor controlling the chemistry of Trx catalysis.

Dewetting and hydrophobic interaction in physical and biological systems.

Hydrophobicity manifests itself differently on large and small length scales. This review focuses on large-length-scale hydrophobicity, particularly on dewetting at single hydrophobic surfaces and drying in regions bounded on two or more sides by hydrophobic surfaces. We review applicable theories, simulations, and experiments pertaining to large-scale hydrophobicity in physical and biomolecular systems and clarify some of the critical issues pertaining to this subject. Given space constraints, we cannot review all the significant and interesting work in this active field.

Role of the active-site solvent in the thermodynamics of factor Xa ligand binding.

Understanding the underlying physics of the binding of small-molecule ligands to protein active sites is a key objective of computational chemistry and biology. It is widely believed that displacement of water molecules from the active site by the ligand is a principal (if not the dominant) source of binding free energy. Although continuum theories of hydration are routinely used to describe the contributions of the solvent to the binding affinity of the complex, it is still an unsettled question as to whether or not these continuum solvation theories describe the underlying molecular physics with sufficient accuracy to reliably rank the binding affinities of a set of ligands for a given protein. Here we develop a novel, computationally efficient descriptor of the contribution of the solvent to the binding free energy of a small molecule and its associated receptor that captures the effects of the ligand displacing the solvent from the protein active site with atomic detail. This descriptor quantitatively predicts (R(2) = 0.81) the binding free energy differences between congeneric ligand pairs for the test system factor Xa, elucidates physical properties of the active-site solvent that appear to be missing in most continuum theories of hydration, and identifies several features of the hydration of the factor Xa active site relevant to the structure-activity relationship of its inhibitors.

A Systematic Comparison of Pairwise and Many-Body Silica Potentials.

The role of many-body effects in modeling silica was investigated using self-consistent force matching. Both pairwise and polarizable classical force fields were developed systematically from ab initio density functional theory force calculations, allowing for a direct comparison of the role of polarization in silica. It was observed that the pairwise potential performed remarkably well at reproducing the basic silica tetrahedral structure. However, the Si-O-Si angle that links the silica tetrahedra showed small but distinct differences with the polarizable potential, a result of the inability of the pairwise potential to properly account for variations in the polarization of the oxygens. Furthermore, the transferability of the polarizable potential was investigated and suggests that additional forces may be necessary to more completely describe silica annealing.

Nanoscale dewetting transition in protein complex folding.

In a previous study, a surprising drying transition was observed to take place inside the nanoscale hydrophobic channel in the tetramer of the protein melittin. The goal of this paper is to determine if there are other protein complexes capable of displaying a dewetting transition during their final stage of folding. We searched the entire protein data bank (PDB) for all possible candidates, including protein tetramers, dimers, and two-domain proteins, and then performed the molecular dynamics (MD) simulations on the top candidates identified by a simple hydrophobic scoring function based on aligned hydrophobic surface areas. Our large scale MD simulations found several more proteins, including three tetramers, six dimers, and two two-domain proteins, which display a nanoscale dewetting transition in their final stage of folding. Even though the scoring function alone is not sufficient (i.e., a high score is necessary but not sufficient) in identifying the dewetting candidates, it does provide useful insights into the features of complex interfaces needed for dewetting. All top candidates have two features in common: (1) large aligned (matched) hydrophobic areas between two corresponding surfaces, and (2) large connected hydrophobic areas on the same surface. We have also studied the effect on dewetting of different water models and different treatments of the long-range electrostatic interactions (cutoff vs PME), and found the dewetting phenomena is fairly robust. This work presents a few proteins other than melittin tetramer for further experimental studies of the role of dewetting in the end stages of protein folding.

Signatures of hydrophobic collapse in extended proteins captured with force spectroscopy.

We unfold and extend single proteins at a high force and then linearly relax the force to probe their collapse mechanisms. We observe a large variability in the extent of their recoil. Although chain entropy makes a small contribution, we show that the observed variability results from hydrophobic interactions with randomly varying magnitude from protein to protein. This collapse mechanism is common to highly extended proteins, including nonfolding elastomeric proteins like PEVK from titin. Our observations explain the puzzling differences between the folding behavior of highly extended proteins, from those folding after chemical or thermal denaturation. Probing the collapse of highly extended proteins with force spectroscopy allows separation of the different driving forces in protein folding.

An experimental and theoretical study of the formation of nanostructures of self-assembled cyanuric acid through hydrogen bond networks on graphite.

The self-assembly of cyanuric acid into ordered nanostructures on a crystalline substrate, highly ordered pyrolytic graphite (HOPG), has been investigated at low temperature under ultrahigh vacuum (UHV) conditions by means of scanning tunneling microscopy in conjunction with theoretical simulations. Many domains with different self-assembly patterns were observed. One such domain represents the formation of an open 2D rosette (cyclic) structure from the self-assembly process, the first observation of this type of structure for pure cyanuric acid on a graphite substrate. Each self-assembled domain exhibits characteristic superstructures formed through different hydrogen bond networks at low coverage and low deposition rate. Experimental observation of coexistent, two-dimensional crystalline structures with distinct hydrogen bond patterns is supported by energy minimizations and molecular dynamics calculations, which show multiple stable structures for this molecule when self-assembled on graphite.

Destruction of long-range interactions by a single mutation in lysozyme.

We propose a mechanism, based on a > or =10-micros molecular dynamics simulation, for the surprising misfolding of hen egg-white lysozyme caused by a single mutation (W62G). Our simulations of the wild-type and mutant lysozymes in 8 M urea solution at biological temperature (with both pH 2 and 7) reveal that the mutant structure is much less stable than that of the wild type, with the mutant showing larger fluctuations and less native-like contacts. Analysis of local contacts reveals that the Trp-62 residue is the key to a cooperative long-range interaction within the wild type, where it acts like a bridge between two neighboring basic residues. Thus, a native-like cluster or nucleation site can form near these residues in the wild type but not in the mutant. The time evolution of the secondary structure also exhibits a quicker loss of the beta-sheets in the mutant than in the wild type, whereas some of the alpha-helices persist during the entire simulation in both the wild type and the mutant in 8 M urea (even though the tertiary structures are basically all gone). These findings, while supporting the general conclusions of a recent experimental study by Dobson and coworkers [Klein-Seetharam J, Oikama M, Grimshaw SB, Wirmer J, Duchardt E, Ueda T, Imoto T, Smith LJ, Dobson CM, Schwalbe H (2002) Science 295:1719-1722], provide a detailed but different molecular picture of the misfolding mechanism.

Motifs for molecular recognition exploiting hydrophobic enclosure in protein-ligand binding.

The thermodynamic properties and phase behavior of water in confined regions can vary significantly from that observed in the bulk. This is particularly true for systems in which the confinement is on the molecular-length scale. In this study, we use molecular dynamics simulations and a powerful solvent analysis technique based on inhomogenous solvation theory to investigate the properties of water molecules that solvate the confined regions of protein active sites. Our simulations and analysis indicate that the solvation of protein active sites that are characterized by hydrophobic enclosure and correlated hydrogen bonds induce atypical entropic and enthalpic penalties of hydration. These penalties apparently stabilize the protein-ligand complex with respect to the independently solvated ligand and protein, which leads to enhanced binding affinities. Our analysis elucidates several challenging cases, including the super affinity of the streptavidin-biotin system.

Frustrated ostwald ripening in self-assembled monolayers of cruciform pi-systems.

This study details a scanning tunneling microscopy investigation into the mechanism of chiral grain growth in highly ordered, self-assembled monolayer films composed of cruciform pi-systems. Although the molecules themselves are achiral, when they adsorb from solution onto graphite, they adopt a gear-like conformation that, by virtue of the surface, is chiral. These handed subunits arrange themselves into enantiomeric two-dimensional domains. The unique finding from this study is that Ostwald ripening is frustrated between domain boundaries that are of opposite chirality because direct interconversion between the chiral units on the surface is energetically inhibited.

Why is the partial molar volume of CO2 so small when dissolved in a room temperature ionic liquid? Structure and dynamics of CO2 dissolved in [Bmim+] [PF6(-)].

When supercritical CO2 is dissolved in an ionic liquid, its partial molar volume is much smaller than that observed in most other solvents. In this article we explore in atomistic detail and explain in an intuitive way the peculiar volumetric behavior experimentally observed when supercritical CO2 is dissolved in 1-butyl-3-methylimidazolium hexafluorophosphate ([Bmim+] [PF6(-)]). We also provide physical insight into the structure and dynamics occurring across the boundary of the CO2 ionic liquid interface. We find that the liquid structure of [Bmim+] [PF6(-)] in the presence of CO2 is nearly identical to that in the neat ionic liquid (IL) even at fairly large mole fractions of CO2. Our simulations indicate, in agreement with experiments, that partial miscibilities of one fluid into the other are very unsymmetrical, CO2 being highly soluble in the ionic liquid phase while the ionic liquid is highly insoluble in the CO2 phase. We interpret our results in terms of the size and shape of spontaneously forming cavities in the ionic liquid phase, and we propose that CO2 occupies extremely well-defined locations in the IL. Even though our accurate prediction of cavity sizes in the neat IL indicates that these cavities are small compared with the van der Waals radius of a single carbon or oxygen atom, CO2 appears to occupy a space that was for the most part a priori "empty".

Integrated Modeling Program, Applied Chemical Theory (IMPACT).

We provide an overview of the IMPACT molecular mechanics program with an emphasis on recent developments and a description of its current functionality. With respect to core molecular mechanics technologies we include a status report for the fixed charge and polarizable force fields that can be used with the program and illustrate how the force fields, when used together with new atom typing and parameter assignment modules, have greatly expanded the coverage of organic compounds and medicinally relevant ligands. As we discuss in this review, explicit solvent simulations have been used to guide our design of implicit solvent models based on the generalized Born framework and a novel nonpolar estimator that have recently been incorporated into the program. With IMPACT it is possible to use several different advanced conformational sampling algorithms based on combining features of molecular dynamics and Monte Carlo simulations. The program includes two specialized molecular mechanics modules: Glide, a high-throughput docking program, and QSite, a mixed quantum mechanics/molecular mechanics module. These modules employ the IMPACT infrastructure as a starting point for the construction of the protein model and assignment of molecular mechanics parameters, but have then been developed to meet specialized objectives with respect to sampling and the energy function.

Importance of accurate charges in molecular docking: quantum mechanical/molecular mechanical (QM/MM) approach.

The extent to which accuracy of electric charges plays a role in protein-ligand docking is investigated through development of a docking algorithm, which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations. In this algorithm, fixed charges of ligands obtained from force field parameterization are replaced by QM/MM calculations in the protein environment, treating only the ligands as the quantum region. The algorithm is tested on a set of 40 cocrystallized structures taken from the Protein Data Bank (PDB) and provides strong evidence that use of nonfixed charges is important. An algorithm, dubbed "Survival of the Fittest" (SOF) algorithm, is implemented to incorporate QM/MM charge calculations without any prior knowledge of native structures of the complexes. Using an iterative protocol, this algorithm is able in many cases to converge to a nativelike structure in systems where redocking of the ligand using a standard fixed charge force field exhibits nontrivial errors. The results demonstrate that polarization effects can play a significant role in determining the structures of protein-ligand complexes, and provide a promising start towards the development of more accurate docking methods for lead optimization applications.

Ultra-high vacuum scanning tunneling microscopy and theoretical studies of 1-halohexane monolayers on graphite.

A simple model system for the 2D self-assembly of functionalized organic molecules on surfaces was examined in a concerted experimental and theoretical effort. Monolayers of 1-halohexanes were formed through vapor deposition onto graphite surfaces in ultrahigh vacuum. Low-temperature scanning tunneling microscopy allowed the molecular conformation, orientation, and monolayer crystallographic parameters to be determined. Essentially identical noncommensurate monolayer structures were found for all 1-halohexanes, with differences in image contrast ascribed mainly to electronic factors. Energy minimizations and molecular dynamics simulations reproduced structural parameters of 1-bromohexane monolayers quantitatively. An analysis of interactions driving the self-assembly process revealed the crucial role played by small but anisotropic electrostatic forces associated with the halogen substituent. While alkyl chain dispersion interactions drive the formation of a close-packed adsorbate monolayer, electrostatic headgroup forces are found to compete successfully in the control of both the angle between lamella and backbone axes and the angle between surface and backbone planes. This competition is consistent with energetic tradeoffs apparent in adsorption energies measured in earlier temperature-programmed desorption studies. In accordance with the higher degree of disorder observed in scanning tunneling microscopy images of 1-fluorohexane, theoretical simulations show that electrostatic forces associated with the fluorine substituent are sufficiently strong to upset the delicate balance of interactions required for the formation of an ordered monolayer. The detailed dissection of the driving forces for self-assembly of these simple model systems is expected to aid in the understanding of the more complex self-assembly processes taking place in the presence of solvent.