Biomedical subjects

Bruce J Berne

Publications and source records attributed to Bruce J Berne.

11 recordsLinked to original sources

PubMed2007

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.

Amino Acid Motifs↗

PubMed2006

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.

Journal Article↗

PubMed2005

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".

Journal Article↗

PubMed2005

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.

Algorithms↗

PubMed2005

Self-assembly of small polycyclic aromatic hydrocarbons on graphite: a combined scanning tunneling microscopy and theoretical approach.

Self-assembled monolayers of chrysene and indene on graphite have been observed and characterized individually with scanning tunneling microscopy (STM) at 80 K under low-temperature, ultrahigh vacuum conditions. These molecules are small, polycyclic aromatic hydrocarbons (PAHs) containing no alkyl chains or functional groups that are known to promote two-dimensional self-assembly. Energy minimization and molecular dynamics simulations performed for small groups of the molecules physisorbed on graphite provide insight into the monolayer structure and forces that drive the self-assembly. The adsorption energy for a single chrysene molecule on a model graphite substrate is calculated to be 32 kcal/mol, while that for indene is 17 kcal/mol. Two distinct monolayer structures have been observed for chrysene, corresponding to high- and low-density assemblies. High-resolution STM images taken of chrysene with different bias polarities reveal distinct nodal structure that is characteristic of the molecular electronic state(s) mediating the tunneling process. Density functional theory calculations are utilized in the assignment of the observed electronic states and possible tunneling mechanism. These results are discussed within the context of PAH and soot particle formation, because both chrysene and indene are known reaction products from the combustion of small hydrocarbons. They are also of fundamental interest in the fields of nanotechnology and molecular electronics.

Journal Article↗

PubMed2005

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.

Journal Article↗

PubMed2005

Drying and hydrophobic collapse of paraffin plates.

We perform molecular dynamics simulations of the hydrophobic collapse of two paraffin plates to examine how the collapse is mediated by realistic paraffin-water attractive van der Waals forces. We explore several aspects of the drying transition between the plates, including the critical separation for drying and the critical size of the vapor bubble required for the nucleation of the drying event. We also investigate the kinetics of hydrophobic collapse and find that the hydrophobic collapse occurs in about 100 ps. We compare these results with the simulations with the plate-water van der Waals attractions turned off and with recent results on the hydrophobic collapse of multidomain proteins. Last, we discuss the relationship among the dewetting transition critical distance, van der Waals potential well depth, and water contact angle on solute surface using a simple macroscopic theory.

Chemistry, Physical↗

PubMed2005

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.

Computer Simulation↗

PubMed2004

Hydrophobic collapse in multidomain protein folding.

We performed molecular dynamics simulations of the collapse of a two-domain protein, the BphC enzyme, into a globular structure to examine how water molecules mediate hydrophobic collapse of proteins. In the interdomain region, liquid water persists with a density 10 to 15% lower than in the bulk, even at small domain separations. Water depletion and hydrophobic collapse occur on a nanosecond time scale, which is two orders of magnitude slower than that found in the collapse of idealized paraffin-like plates. When the electrostatic protein-water forces are turned off, a dewetting transition occurs in the interdomain region and the collapse speeds up by more than an order of magnitude. When attractive van der Waals forces are turned off as well, the dewetting in the interdomain region is more profound, and the collapse is even faster.

Computer Simulation↗

PubMed2002

Can a continuum solvent model reproduce the free energy landscape of a beta -hairpin folding in water?

The folding free energy landscape of the C-terminal beta-hairpin of protein G is explored using the surface-generalized Born (SGB) implicit solvent model, and the results are compared with the landscape from an earlier study with explicit solvent model. The OPLSAA force field is used for the beta-hairpin in both implicit and explicit solvent simulations, and the conformational space sampling is carried out with a highly parallel replica-exchange method. Surprisingly, we find from exhaustive conformation space sampling that the free energy landscape from the implicit solvent model is quite different from that of the explicit solvent model. In the implicit solvent model some nonnative states are heavily overweighted, and more importantly, the lowest free energy state is no longer the native beta-strand structure. An overly strong salt-bridge effect between charged residues (E42, D46, D47, E56, and K50) is found to be responsible for this behavior in the implicit solvent model. Despite this, we find that the OPLSAA/SGB energies of all the nonnative structures are higher than that of the native structure; thus the OPLSAA/SGB energy is still a good scoring function for structure prediction for this beta-hairpin. Furthermore, the beta-hairpin population at 282 K is found to be less than 40% from the implicit solvent model, which is much smaller than the 72% from the explicit solvent model and approximately equal 80% from experiment. On the other hand, both implicit and explicit solvent simulations with the OPLSAA force field exhibit no meaningful helical content during the folding process, which is in contrast to some very recent studies using other force fields.

Animals↗

PubMed2002

The calculation of transport properties in quantum liquids using the maximum entropy numerical analytic continuation method: application to liquid para-hydrogen.

We present a method based on augmenting an exact relation between a frequency-dependent diffusion constant and the imaginary time velocity autocorrelation function, combined with the maximum entropy numerical analytic continuation approach to study transport properties in quantum liquids. The method is applied to the case of liquid para-hydrogen at two thermodynamic state points: a liquid near the triple point and a high-temperature liquid. Good agreement for the self-diffusion constant and for the real-time velocity autocorrelation function is obtained in comparison to experimental measurements and other theoretical predictions. Improvement of the methodology and future applications are discussed.

Diffusion↗