The effect of vicinal polar and charged groups on hydrophobic hydration

Existence of Quantum Pharmacology in Sartans: Evidence in Isolated Rabbit Iliac Arteries

Quantum pharmacology introduces theoretical models to describe the possibility of ultra-high dilutions to produce biological effects, which may help to explain the placebo effect observed in hypertensive clinical trials. To determine this within physiology and to evaluate novel ARBs, we tested the ability of known angiotensin II receptor blockers (ARBs) (candesartan and telmisartan) used to treat hypertension and other cardiovascular diseases, as well as novel ARBs (benzimidazole-N-biphenyl tetrazole (ACC519T), benzimidazole-bis-N,N'-biphenyl tetrazole (ACC519T(2)) and 4-butyl-N,N0-bis[[20-2Htetrazol-5-yl)biphenyl-4-yl]methyl)imidazolium bromide (BV6(K+)2), and nirmatrelvir (the active ingredient in Paxlovid) to modulate vascular contraction in iliac rings from healthy male New Zealand White rabbits in responses to various vasopressors (angiotensin A, angiotensin II and phenylephrine). Additionally, the hemodynamic effect of ACC519T and telmisartan on mean arterial pressure in conscious rabbits was determined, while the ex vivo ability of BV6(K+)2 to activate angiotensin-converting enzyme-2 (ACE2) was also investigated. We show that commercially available and novel ARBs can modulate contraction responses at ultra-high dilutions to different vasopressors. ACC519T produced a dose-dependent reduction in rabbit mean arterial pressure while BV6(K+)2 significantly increased ACE2 metabolism. The ability of ARBs to inhibit contraction responses even at ultra-low concentrations provides evidence of the existence of quantum pharmacology. Furthermore, the ability of ACC519T and BV6(K+)2 to modulate blood pressure and ACE2 activity, respectively, indicates their therapeutic potential against hypertension.

Progress in Simulation Studies of Insulin Structure and Function

Insulin is a peptide hormone known for chiefly regulating glucose level in blood among several other metabolic processes. Insulin remains the most effective drug for treating diabetes mellitus. Insulin is synthesized in the pancreatic β-cells where it exists in a compact hexameric architecture although its biologically active form is monomeric. Insulin exhibits a sequence of conformational variations during the transition from the hexamer state to its biologically-active monomer state. The structural transitions and the mechanism of action of insulin have been investigated using several experimental and computational methods. This review primarily highlights the contributions of molecular dynamics (MD) simulations in elucidating the atomic-level details of conformational dynamics in insulin, where the structure of the hormone has been probed as a monomer, dimer, and hexamer. The effect of solvent, pH, temperature, and pressure have been probed at the microscopic scale. Given the focus of this review on the structure of the hormone, simulation studies involving interactions between the hormone and its receptor are only briefly highlighted, and studies on other related peptides (e.g., insulin-like growth factors) are not discussed. However, the review highlights conformational dynamics underlying the activities of reported insulin analogs and mimetics. The future prospects for computational methods in developing promising synthetic insulin analogs are also briefly highlighted.

Insights into the Sensitivity of Arginine Concentration to Preserve the Folded Form of Insulin Monomer under Thermal Stress

Arginine, although popularly known as aggregation suppressor additive, has been found to quench proteins' structure and function by destabilizing their conformations. Driven by such controversial evidence, in this work we performed a series of atomistic molecular dynamics simulations of insulin monomer, a biologically active hormone protein, in arginine solution of varying concentrations (0.5, 1, and 2 M) at ambient and elevated temperature (400 K) to explore the arginine concentration driven structure-based stability of the protein. Our study reveals that the flexibility of the protein's structure is dependent on the arginine concentration, and among all the used solutions, 2 M arginine, a "neutral crowder" that mimics the cellular environment, can preserve the native folded form of the protein at ambient temperature in an excellent manner. Further, while the protein unfolds at 400 K in pure water, this solution worked satisfactorily to preserve the protein's folded conformation more firmly than the other solutions. The replica-exchange MD of insulin in 2 M arginine solution further supports the fact. In this aspect an important issue in molecular pharmacology is to identify and recognize the physical origin of the stability of a protein, i.e, in this case, how arginine directs the conformational flexibility of the protein and preserves its native folded form. We identified that the exclusion of arginine from the protein surface increases the local structuration of water around the protein, thereby preserving its "biological water" layer, and makes the protein more hydrated at 2 M concentration as compared to the other arginine solutions. Additionally, our microscopic investigation on the interactions of the protein-solvation layer revealed that the structural heterogeneity of the protein surface, arising from the differential physicochemical nature of the amino acid residues, controls the favorable formation of sluggish water-arginine mixed solvation layer at higher arginine concentration that helps the protein to maintain its structural rigidity. Importantly, apart from the protein-solvent hydrogen-bonding interactions, the anion-pi interactions, established between the carboxyl group of arginine and the aromatic amino acid residues of insulin, were recognized to facilitate the protein to maintain its native folded form at the experimental temperatures.

Extent of Surface Force Additivity on Chemically Heterogeneous Substrates at Varied Orientations

Surface interactions between chemically mixed surfaces, as well as those among dissolved biomolecules, comprise distinct contributions from polar and hydrophobic moieties. These contributions are often context dependent. Approximate compliance to the Cassie additivity equation for the wetting free energies on mixed surfaces in water is, however, indicative of similarly additive forces between individual surface elements, suggesting a quadratic interpolation model for total force from the forces between pure surfaces. We use molecular dynamics/umbrella sampling simulations of parallel and nonparallel mixed surfaces with demonstrable Cassie-like behavior to verify how well the total surface force between the heterogeneous, molecularly rough surfaces can be approximated as a combination of forces among the homogeneous ones. When accounting for dissimilar distances of approach between functional groups of different types, our results for graphene surfaces with mixed methyl and nitrile coating show such a superposition to provide a reasonable first order approximation of interactions between the platelets. Deviations from additivity are more prominent in parallel-plate configurations, at high content of hydrophobic groups, and small separations. The inclusion of water polarizability does not visibly alter the observed behavior regardless of platelet orientations. The outcome of this study determines the necessary molecular conditions for observing force additivity that emphasize the context dependence of hydrophobic interaction in the presence of polar groups. This notion provides guidelines for the syntheses of new, chemically heterogeneous materials with tailored function-oriented properties in aqueous media.

Hydrophobic hydration driven self-assembly of curcumin in water: similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that the curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of β-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. We have constructed a free energy landscape of aggregation where the metatsable state has been identified. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with the total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting rather strongly with the hydrophilic groups at the surface of the curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer.

High-dilution effects revisited. 1. Physicochemical aspects

Several lines of evidence suggest that homeopathic high dilutions (HDs) can effectively have a pharmacological action, and so cannot be considered merely placebos. However, until now there has been no unified explanation for these observations within the dominant paradigm of the dose-response effect. Here the possible scenarios for the physicochemical nature of HDs are reviewed. A number of theoretical and experimental approaches, including quantum physics, conductometric and spectroscopic measurements, thermoluminescence, and model simulations investigated the peculiar features of diluted/succussed solutions. The heterogeneous composition of water could be affected by interactive phenomena such as coherence, epitaxy and formation of colloidal nanobubbles containing gaseous inclusions of oxygen, nitrogen, carbon dioxide, silica and, possibly, the original material of the remedy. It is likely that the molecules of active substance act as nucleation centres, amplifying the formation of supramolecular structures and imparting order to the solvent. Three major models for how this happens are currently being investigated: the water clusters or clathrates, the coherent domains postulated by quantum electrodynamics, and the formation of nanoparticles from the original solute plus solvent components. Other theoretical approaches based on quantum entanglement and on fractal-type self-organization of water clusters are more speculative and hypothetical. The problem of the physicochemical nature of HDs is still far from to be clarified but current evidence strongly supports the notion that the structuring of water and its solutes at the nanoscale can play a key role.

Wetting and interfacial properties of water nanodroplets in contact with graphene and monolayer boron-nitride sheets

Born-Oppenheim quantum molecular dynamics (QMD) simulations are performed to investigate wetting, diffusive, and interfacial properties of water nanodroplets in contact with a graphene sheet or a monolayer boron-nitride (BN) sheet. Contact angles of the water nanodroplets on the two sheets are computed for the first time using QMD simulations. Structural and dynamic properties of the water droplets near the graphene or BN sheet are also studied to gain insights into the interfacial interaction between the water droplet and the substrate. QMD simulation results are compared with those from previous classic MD simulations and with the experimental measurements. The QMD simulations show that the graphene sheet yields a contact angle of 87°, while the monolayer BN sheet gives rise to a contact angle of 86°. Hence, like graphene, the monolayer BN sheet is also weakly hydrophobic, even though the BN bonds entail a large local dipole moment. QMD simulations also show that the interfacial water can induce net positive charges on the contacting surface of the graphene and monolayer BN sheets, and such charge induction may affect electronic structure of the contacting graphene in view that graphene is a semimetal. Contact angles of nanodroplets of water in a supercooled state on the graphene are also computed. It is found that under the supercooled condition, water nanodroplets exhibit an appreciably larger contact angle than under the ambient condition.

Probing surface tension additivity on chemically heterogeneous surfaces by a molecular approach

Surface free energy of a chemically heterogeneous surface is often treated as an approximately additive quantity through the Cassie equation [Cassie ABD (1948) Discuss Faraday Soc 3:11-16]. However, deviations from additivity are common, and molecular interpretations are still lacking. We use molecular simulations to measure the microscopic analogue of contact angle, Θ(c), of aqueous nanodrops on heterogeneous synthetic and natural surfaces as a function of surface composition. The synthetic surfaces are layers of graphene functionalized with prototypical nonpolar and polar head group: methyl, amino, and nitrile. We demonstrate positive as well as negative deviations from the linear additivity. We show the deviations reflect the uneven exposure of mixture components to the solvent and the linear relation is recovered if fractions of solvent-accessible surface are used as the measure of composition. As the spatial variations in polarity become of larger amplitude, the linear relation can no longer be obtained. Protein surfaces represent such natural patterned surfaces, also characterized by larger patches and roughness. Our calculations reveal strong deviations from linear additivity on a prototypical surface comprising surface fragments of melittin dimer. The deviations reflect the disproportionately strong influence of isolated polar patches, preferential wetting, and changes in the position of the liquid interface above hydrophobic patches. Because solvent-induced contribution to the free energy of surface association grows as cos Θ(c), deviations of cos Θ(c) from the linear relation directly reflect nonadditive adhesive energies of biosurfaces.

The hydrophobic effect: molecular dynamics simulations of water confined between extended hydrophobic and hydrophilic surfaces

Structural and dynamic properties of water confined between two parallel, extended, either hydrophobic or hydrophilic crystalline surfaces of n-alkane C(36)H(74) or n-alcohol C(35)H(71)OH, are studied by molecular dynamics simulations. Electron density profiles, directly compared with corresponding experimental data from x-ray reflectivity measurements, reveal a uniform weak de-wetting characteristic for the extended hydrophobic surface, while the hydrophilic surface is weakly wetted. These microscopic data are consistent with macroscopic contact angle measurements. Specific water orientation is present at both surfaces. The ordering is characteristically different between the surfaces and of longer range at the hydrophilic surface. Furthermore, the dynamic properties of water are different at the two surfaces and different from the bulk behavior. In particular, at the hydrophobic surface, time-correlation functions reveal that water molecules have characteristic diffusive behavior and orientational ordering due to the lack of hydrogen bonding interactions with the surface. These observations suggest that the altered dynamical properties of water in contact with extended hydrophobic surfaces together with a partial drying of the surfaces are more indicative of the hydrophobic effect than structural ordering, which we suggest to be independent of surface topology.

Effect of Surface Polarity on Water Contact Angle and Interfacial Hydration Structure

We perform molecular dynamics simulations of water in the presence of hydrophobic/hydrophilic walls at T = 300 K and P = 0 GPa. For the hydrophilic walls, we use a hydroxylated silica model introduced in previous simulations [Lee, S. H.; Rossky, P. J. J. Chem. Phys. 1994, 100, 3334. Giovambattista, N.; Rossky, P. J.; Debenedetti, P. G.; Phys. Rev. E 2006, 73, 041604.]. By rescaling the physical partial atomic charges by a parameter 0 <or= k <or= 1, we can continuously transform the hydrophilic walls (hydroxylated silica, k = 1) into hydrophobic apolar surfaces (k = 0). From a physical point of view, k is the normalized magnitude of a surface dipole moment, and thus it quantifies the polarity of the surface. We calculate the contact angle of water for 0 <or= k <or= 1. We find that, at least for the present homogeneous, atomically flat, and defect-free surface model, the magnitude of the surface dipole correlates with the contact angle in a one-to-one correspondence. In particular, we find that polar surfaces with 0 < k <or= kc = 0.4 are macroscopically hydrophobic; that is, the contact angle is larger than 90 degrees . For the cutoff value k = kc, the magnitude of the dipole moment of the polar silica surface unit is 41% that of the water molecule dipole moment. We also study the water orientation distributions next to the walls (a microscopic property). We find that these distributions also correlate with the contact angle in a one-to-one correspondence. Thus, the structure of confined water, the surface polarity, and the contact angle are in a direct correspondence to each other, and therefore, each quantifies the hydrophobicity/hydrophilicity of the surface.

Oxygen accessibility to ribonuclease a: quantitative interpretation of nuclear spin relaxation induced by a freely diffusing paramagnet

The nuclear spin relaxation induced by a freely diffusing paramagnetic center provides a direct measure of intermolecular accessibility. A number of factors are involved in a quantitative interpretation of relaxation data including excluded volume effects, solvation differences, and the details of the electron spin relaxation in the paramagnetic center. In the case where the electron relaxation time is short compared with correlation times describing the electron-nuclear coupling, the nuclear spin relaxation rates may be related to the effective local concentration of the paramagnetic center at different locations about the solute of interest. The local concentrations may in turn be related to differences in the local free energies of interaction between the diffusing paramagnet and the cosolute. We demonstrate this approach for the case of ribonuclease A and deduce surface free energy differences for a large number of protein proton sites. We find that the oxygen accessibility is poorly represented by hard-sphere models such as computed solvent or steric accessibility. There is a distribution of local intermolecular interactions with a width of the order of RT that dominates the report of the intermolecular exploration of the protein by this simple solute.

Structural and dynamical analysis of the hydration of the Alzheimer's beta-amyloid peptide

An analysis of the water molecules in the first solvation shell obtained from the molecular dynamics simulation of the amyloid beta(10-35)NH2 peptide and the amyloid beta(10-35)NH2E22Q "Dutch" mutant peptide is presented. The structure, energetics, and dynamics of water in the hydration shell have been investigated using a variety of measures, including the hydrogen bond network, the water residence times for all the peptide residues, the diffusion constant, experimentally determined HN amide proton exchange, and the transition probabilities for water to move from one residue to another or into the bulk. The results of the study indicate that: (1) the water molecules at the peptide-solvent interface are organized in an ordered structure similar for the two peptide systems but different from that of the bulk, (2) the peptide structure inhibits diffusion perpendicular to the peptide surface by a factor of 3 to 5 relative to diffusion parallel to the peptide surface, which is comparable to diffusion of bulk water, (3) water in the first solvation shell shows dynamical relaxation on fast (1-2 ps) and slow (10-40 ps) time scales, (4) a novel solvent relaxation master equation is shown to capture the details of the fast relaxation of water in the peptide's first solvation shell, (5) the interaction between the peptide and the solvent is stronger in the wild type than in the E22Q mutant peptide, in agreement with earlier results obtained from computer simulations [Massi, F.; Straub, J. E. Biophys J 2001, 81, 697] correlated with the observed enhanced activity of the E22Q mutant peptide.

Hydrophobicity at a Janus interface

Water confined between adjoining hydrophobic and hydrophilic surfaces (a Janus interface) is found to form stable films of nanometer thickness whose responses to shear deformations are extraordinarily noisy. The power spectrum of this noise is quantified. In addition, the frequency dependence of the complex shear modulus is a power law with slope one-half, indicating a distribution of relaxation processes rather than any dominant one. The physical picture emerges that whereas surface energetics encourage water to dewet the hydrophobic side of the interface, the hydrophilic side constrains water to be present, resulting in a flickering, fluctuating complex.

Kinetic evidence for hydrophobically stabilized encounter complexes formed by hydrophobic esters in aqueous solutions containing monohydric alcohols

The pH-independent hydrolysis of four esters, p-methoxyphenyl 2,2-dichloroethanoate (1a), p-methoxyphenyl 2,2-dichloropropanoate (1b), p-methoxyphenyl 2,2-dichlorobutanoate (1c), and p-methoxyphenyl 2,2-dichloropentanoate (1d), in dilute aqueous solution has been studied as a function of the molality of added cosolutes ethanol, 1-propanol, and 1-butanol. The rate constants for the neutral hydrolysis decrease with increasing cosolute concentration. These kinetic medium effects respond to both the hydrophobicity of the ester and of the monohydric alcohol. The observed rate effects were analyzed using both a thermodynamic and a kinetic model. The kinetic model suggests a molecular picture of a hydrophobically stabilized encounter complex, with equilibrium constants K(ec) often smaller than unity, in which the cosolute blocks the reaction center of the hydrolytic ester for attack by water. The formation of these encounter complexes leads to a dominant initial-state stabilization as follows from the thermodynamic model. Decreases in both apparent enthalpies and entropies of activation for these hydrolysis reactions correspond to unfavorable enthalpies and favorable entropies of complexation, which confirms that the encounter complexes are stabilized by hydrophobic interactions.

Probing the origins of increased activity of the E22Q "Dutch" mutant Alzheimer's beta-amyloid peptide

The amyloid peptide congener A beta(10--35)-NH(2) is simulated in an aqueous environment in both the wild type (WT) and E22Q "Dutch" mutant forms. The origin of the noted increase in deposition activity resulting from the Dutch mutation is investigated. Multiple nanosecond time scale molecular dynamics trajectories were performed and analyzed using a variety of measures of the peptide's average structure, hydration, conformational fluctuations, and dynamics. The results of the study support the conclusions that 1) the E22Q mutant and WT peptide are both stable in "collapsed coil" conformations consistent with the WT structure of, J. Struct. Biol. 130:130--141); 2) the E22Q peptide is more flexible in solution, supporting early claims that its equilibrium structural fluctuations are larger than those of the WT peptide; and 3) the local E22Q mutation leads to a change in the first solvation layer in the region of the peptide's "hydrophobic patch," resulting in a more hydrophobic solvation of the mutant peptide. The simulation results support the view that the noted increase in activity due to the Dutch mutation results from an enhancement of the desolvation process that is an essential step in the aggregation of the peptide.

Effects of electric charges on hydrophobic forces. II

We study by molecular-dynamics simulations the effect of electric charges of either sign on hydrophobic interactions and on the dynamics of hydration water, using explicit water and very simplified solutes. Results show that the presence of a charged solute can disrupt the "hydrophobic contact bond" between two apolar solutes nearby, by forcing them towards a different configuration. As a consequence of different structural changes of the solvent caused by charges of opposite sign, the effect is markedly charge-sign-dependent. Analogous weaker effects appear to be induced by the presence of one additional apolar element. The dynamics of hydration water around each solute is also seen to be strongly influenced by the presence of other (charged or uncharged) nearby solutes. Comparison between our results on hydration water dynamics around charged solutes and available experimental data allows sorting out the effects of solute charge sign and size. Our results also offer a plain interpretation of the equivalence of the effects on water structure due to solute ions and to high pressures. These results reflect at a basic paradigmatic level the immensely more complex cases of well-known phenomena such as salting-in and salting-out, and of protein conformational changes caused, e.g., by the arrival of a charged or of an apolar group (phosphorilation or methylation). As it will be discussed, they help in the direction of Delbruck's desirable "progress towards a radical physical explanation" for this class of phenomena.