Computational techniques have the potential to accelerate the design and optimization of nanomaterials for applications such as drug delivery and contaminant removal; however, the success of such techniques requires reliable models of nanomaterial surfaces as well as accurate descriptions of their interactions with relevant solutes. In the present work, we evaluate the ability of selected models of naked and hydroxylated carbon nanotubes to predict adsorption equilibrium constants for about 30 small aromatic compounds with a variety of functional groups. The equilibrium constants determined using molecular dynamics coupled with free-energy calculation techniques are directly compared to those derived from experimental measurements. The calculations are highly predictive of the relative adsorption affinities of the compounds, with excellent correlation (r ≥ 0.9) between calculated and measured values of the logarithm of the adsorption equilibrium constant. Moreover, the agreement in absolute terms is also ...

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) has been recognized as an important activator of certain transient receptor potential (TRP) channels. More specifically, TRPV1 is a pain receptor activated by a wide range of stimuli. However, whether or not PI(4,5)P2 is a TRPV1 agonist remains open to debate. Utilizing a combined approach of mutagenesis and molecular modeling, we identified a PI(4,5)P2 binding site located between the TRP box and the S4-S5 linker. At this site, PI(4,5)P2 interacts with the amino acid residues R575 and R579 in the S4-S5 linker and with K694 in the TRP box. We confirmed that PI(4,5)P2 behaves as a channel agonist and found that R575, R579 and K694 mutations to alanine reduce PI(4,5)P2 binding affinity. Additionally, in silico mutations R575A, R579A and K694A showed that the reduction in binding affinity results from the delocalization of PI(4,5)P2 in the binding pocket. Molecular dynamics simulations indicate that PI(4,5)P2 binding induces conformational rearrangements of the structure formed by S6 and the TRP domain, which cause an opening of the lower TRPV1 channel gate.

The photophysics and photochemistry of rose bengal (RB) and methylene blue (MB) bound to human serum albumin (HSA) have been investigated under a variety of experimental conditions. Distribution of the dyes between the external solvent and the protein has been estimated by physical separation and fluorescence measurements. The main localization of protein-bound dye molecules was estimated by the intrinsic fluorescence quenching, displacement of fluorescent probes bound to specific protein sites, and by docking modelling. All the data indicate that, at low occupation numbers, RB binds strongly to the HSA site I, while MB localizes predominantly in the protein binding site II. This different localization explains the observed differences in the dyes' photochemical behaviour. In particular, the environment provided by site I is less polar and considerably less accessible to oxygen. The localization of RB in site I also leads to an efficient quenching of the intrinsic protein fluorescence (ascribed to the nearby Trp residue) and the generation of intra-protein singlet oxygen, whose behaviour is different to that observed in the external solvent or when it is generated by bound MB.

We investigated two critical aspects of rose Bengal (RB) photosensitized protein cross-linking that may underlie recently developed medical applications. Our studies focused on the binding of RB to collagen by physical interaction and the effect of this binding and certain amino acids on RB photochemistry. Molecular dynamics simulations and free-energy calculation techniques, complemented with isothermal titration calorimetry, provided insight into the binding between RB and a collagen-like peptide (CLP) at the atomic level. Electrostatic interactions dominated, which is consistent with the finding that RB bound equally well to triple helical and single chain collagen. The binding free energy ranged from −5.7 to −3 kcal/mol and was strongest near the positively charged amino groups at the N-terminus and on lysine side chains. At high RB concentration, a maximum of 16 ± 3 bound dye molecules per peptide was found, which is consistent with spectroscopic evidence for aggregated RB bound to collagen or the CLP. Within a tissue-mimetic collagen matrix, RB photobleached rapidly, probably due to electron transfer to certain protein amino acids, as was demonstrated in solutions of free RB and arginine. In the presence of arginine and low oxygen concentrations, a product absorbing at 510 nm formed, presumably due to dehalogenation after electron transfer to RB. In the collagen matrix without arginine, the dye generated singlet oxygen as well as the 510 nm product. These results provide the first evidence of the effects of a tissue-like environment on the photochemical mechanisms of rose Bengal.

Adsorption of organic molecules from aqueous solution to the surface of carbon nanotubes or graphene is an important process in many applications of these materials. Here we use molecular dynamics simulation, supplemented by analytical chemistry, to explore in detail the adsorption thermodynamics of a diverse set of aromatic compounds on graphenic materials, elucidating the effects of the solvent, surface coverage, surface curvature, defects, and functionalization by hydroxy groups. We decompose the adsorption free energies into entropic and enthalpic components and find that different classes of compounds—such as phenols, benzoates, and alkylbenzenes—can easily be distinguished by the relative contributions of entropy and enthalpy to their adsorption free energies. Overall, entropy dominates for the more hydrophobic compounds, while enthalpy plays the greatest role for more hydrophilic compounds. Experiments and independent simulations using two different force field frameworks (CHARMM and Amber) support the robustness of these conclusions. We determine that concave curvature is generally associated with greater adsorption affinity, more favorable enthalpy, and greater contact area, while convex curvature reduces both adsorption enthalpy and contact area. Defects on the graphene surfaces can create concave curvature, resulting in localized binding sites. As the graphene surface becomes covered with aromatic solutes, the affinity for adsorbing an additional solute increases until a complete monolayer is formed, driven by more favorable enthalpy and partially canceled by less favorable entropy. Similarly, hydroxylation of the surface leads to preferential adsorption of the aromatic solutes to remaining regions of bare graphene, resulting in less favorable adsorption entropy, but compensated by an increase in favorable enthalpic interactions.