How protein chemists learned about the hydrophobic factor

Free-energy decomposition of salt effects on the solubilities of small molecules and the role of excluded-volume effects

The roles of cations and anions are different in the perturbation on solvation, and thus, the analyses of the separated contributions from cations and anions are useful to establish molecular pictures of ion-specific effects. In this work, we investigate the effects of cations, anions, and water separately in the solvation of n-alcohols and n-alkanes by free-energy decomposition. By utilising energy-representation theory of solvation, we address the contributions arising from the direct solute-solvent interactions and the excluded-volume effects. It is found that the change in solvation of n-alcohols and n-alkanes upon addition of salt depends primarily on the anion species. The direct interaction between the anion and solute is in agreement with the Setschenow coefficient in terms of the ranking of salting-in and salting-out for n-alkanes, which corresponds to the extent of accumulation of the anion on the solute surface. For each of the n-alcohols and n-alkanes examined, the excluded-volume component in the Setschenow coefficient is well correlated to the (total) Setschenow coefficient when the salt effects are concerned. The ranking of the excluded-volume component in the variation of the salt species is parallel to the water contribution, which is correlated further to the change in the water density upon the addition of the salt.

A thermodynamic analysis of CLC transporter dimerization in lipid bilayers

The CLC-ec1 chloride/proton antiporter is a membrane-embedded homodimer with subunits that can dissociate and associate, but the thermodynamic driving forces favor the assembled dimer at biological densities. Yet, the physical reasons for this stability are confounding as dimerization occurs via the burial of hydrophobic interfaces away from the lipid solvent. For binding of nonpolar surfaces in aqueous solution, the driving force is often attributed to the hydrophobic effect, but this should not apply in the membrane since there is very little water. To investigate this further, we quantified the thermodynamic changes associated with CLC dimerization in membranes by carrying out a van 't Hoff analysis of the temperature dependency of the free energy of dimerization, ΔG°. To ensure that the reaction reached equilibrium at different temperatures, we utilized a Förster resonance energy transfer assay to report on relaxation kinetics of subunit exchange as a function of temperature. Equilibration times were then applied to measure CLC-ec1 dimerization isotherms at different temperatures using the single-molecule subunit-capture photobleaching analysis approach. The results demonstrate that the dimerization free energy of CLC in Escherichia coli-like membranes exhibits a nonlinear temperature dependency corresponding to a large, negative change in heat capacity, a signature of solvent ordering effects such as the hydrophobic effect. Consolidating this with our previous molecular analyses suggests that the nonbilayer defect required to solvate the monomeric state is one source of the observed change in heat capacity and indicates the existence of a generalizable driving force for protein association in membranes.

Rapid prediction and analysis of protein intrinsic disorder

Protein intrinsic disorder is found in all kingdoms of life and is known to underpin numerous physiological and pathological processes. Computational methods play an important role in characterizing and identifying intrinsically disordered proteins and protein regions. Herein, we present a new high-efficiency web-based disorder predictor named Rapid Intrinsic Disorder Analysis Online (RIDAO) that is designed to facilitate the application of protein intrinsic disorder analysis in genome-scale structural bioinformatics and comparative genomics/proteomics. RIDAO integrates six established disorder predictors into a single, unified platform that reproduces the results of individual predictors with near-perfect fidelity. To demonstrate the potential applications, we construct a test set containing more than one million sequences from one hundred organisms comprising over 420 million residues. Using this test set, we compare the efficiency and accessibility (i.e., ease of use) of RIDAO to five well-known and popular disorder predictors, namely: AUCpreD, IUPred3, metapredict V2, flDPnn, and SPOT-Disorder2. We show that RIDAO yields per-residue predictions at a rate two to six orders of magnitude greater than the other predictors and completely processes the test set in under an hour. RIDAO can be accessed free of charge at https://ridao.app.

Biophysical studies of amorphous protein aggregation and in vivo immunogenicity

Amorphous protein aggregates are oligomers that lack specific, high-order structures. Soluble amorphous aggregates are smaller than ~1 µm. Despite their lack of high-order structure, amorphous protein aggregates exhibit specific biophysical properties such as reversibility of formation, density, conformation, and biochemical stability. Our mutational analysis using a Solubility Controlling Peptide (SCP) tag strongly suggests that amorphous aggregation of small globular proteins can significantly increase in vivo immune response and that the magnitude of enhanced immune responses depends on the aggregates' biophysical and biochemical properties. We propose that SCP tags might help develop subunit (protein) adjuvant-free (immunostimulant-free) vaccines by controlling the aggregation propensity of target proteins.

Amyloid β 42 fibril structure based on small-angle scattering

Amyloid fibrils are associated with a number of neurodegenerative diseases, including fibrils of amyloid β42 peptide (Aβ42) in Alzheimer's disease. These fibrils are a source of toxicity to neuronal cells through surface-catalyzed generation of toxic oligomers. Detailed knowledge of the fibril structure may thus facilitate therapeutic development. We use small-angle scattering to provide information on the fibril cross-section dimension and shape for Aβ42 fibrils prepared in aqueous phosphate buffer at pH = 7.4 and pH 8.0 under quiescent conditions at 37 °C from pure recombinant Aβ42 peptide. Fitting the data using a continuum model reveals an elliptical cross-section and a peptide mass-per-unit length compatible with two filaments of two monomers, four monomers per plane. To provide a more detailed atomistic model, the data were fitted using as a starting state a high-resolution structure of the two-monomer arrangement in filaments from solid-state NMR (Protein Data Bank ID 5kk3). First, a twofold symmetric model including residues 11 to 42 of two monomers in the filament was optimized in terms of twist angle and local packing using Rosetta. A two-filament model was then built and optimized through fitting to the scattering data allowing the two N-termini in each filament to take different conformations, with the same conformation in each of the two filaments. This provides an atomistic model of the fibril with twofold rotation symmetry around the fibril axis. Intriguingly, no polydispersity as regards the number of filaments was observed in our system over separate samples, suggesting that the two-filament arrangement represents a free energy minimum for the Aβ42 fibril.

Structure and adsorption behavior of high hydrostatic pressure-treated β-lactoglobulin

HYPOTHESIS

High hydrostatic pressure treatment causes structural changes in interfacial-active β-lactoglobulin (β-lg). We hypothesized that the pressure-induced structural changes affect the intra- and intermolecular interactions which determine the interfacial activity of β-lg. The conducted experimental and numerical investigations could contribute to the mechanistic understanding of the adsorption behavior of proteins in food-related emulsions.

EXPERIMENTS

We treated β-lg in water at pH 7 with high hydrostatic pressures up to 600 MPa for 10 min at 20 °C. The secondary structure was characterized with Fourier-transform infrared spectroscopy (FTIR) and circular dichroism (CD), the surface hydrophobicity and charge with fluorescence-spectroscopy and ζ-potential, and the quaternary structure with membrane-osmometry, analytical ultracentrifugation (AUC) and mass spectrometry (MS). Experimental analyses were supported through molecular dynamic (MD) simulations. The adsorption behavior was investigated with pendant drop analysis.

FINDINGS

MD simulation revealed a pressure-induced molten globule state of β-lg, confirmed by an unfolding of β-sheets with FTIR, a stabilization of α-helices with CD and loss in tertiary structure induced by an increase in surface hydrophobicity. Membrane-osmometry, AUC and MS indicated the formation of non-covalently linked dimers that migrated slower through the water phase, adsorbed more quickly due to hydrophobic interactions with the oil, and lowered the interfacial tension more strongly than reference β-lg.

Bridging Gaussian Density Fluctuations from Microscopic to Macroscopic Volumes: Applications to Non-Polar Solute Hydration Thermodynamics

The hydration of hydrophobic solutes is intimately related to the spontaneous formation of cavities in water through ambient density fluctuations. Information theory-based modeling and simulations have shown that water density fluctuations in small volumes are approximately Gaussian. For limiting cases of microscopic and macroscopic volumes, water density fluctuations are known exactly and are rigorously related to the density and isothermal compressibility of water. Here, we develop a theory—interpolated gaussian fluctuation theory (IGFT)—that builds an analytical bridge to describe water density fluctuations from microscopic to molecular scales. This theory requires no detailed information about the water structure beyond the effective size of a water molecule and quantities that are readily obtained from water’s equation-of-state—namely, the density and compressibility. Using simulations, we show that IGFT provides a good description of density fluctuations near the mean, that is, it characterizes the variance of occupancy fluctuations over all solute sizes. Moreover, when combined with the information theory, IGFT reproduces the well-known signatures of hydrophobic hydration, such as entropy convergence and solubility minima, for atomic-scale solutes smaller than the crossover length scale beyond which the Gaussian assumption breaks down. We further show that near hydrophobic and hydrophilic self-assembled monolayer surfaces in contact with water, the normalized solvent density fluctuations within observation volumes depend similarly on size as observed in the bulk, suggesting the feasibility of a modified version of IGFT for interfacial systems. Our work highlights the utility of a density fluctuation-based approach toward understanding and quantifying the solvation of non-polar solutes in water and the forces that drive them toward surfaces with different hydrophobicities.

In silico investigation of potential inhibitors to main protease and spike protein of SARS-CoV-2 in propolis

Docking analysis of propolis's natural compound was successfully performed against SARS-CoV-2 main protease (Mpro) and spike protein subunit 2 (S2). Initially, the propolis's protein was screened using chromatography analysis and successfully identified 22 compounds in the propolis. Four compounds were further investigated, i.e., neoblavaisoflavone, methylophiopogonone A, 3'-Methoxydaidzin, and genistin. The binding affinity of 3'-Methoxydaidzin was -7.7 kcal/mol, which is similar to nelfinavir (control), while the others were -7.6 kcal/mol. However, we found the key residue of Glu A:166 in the methylophiopogonone A and genistin, even though the predicted binding energy slightly higher than nelfinavir. In contrast, the predicted binding affinity of neoblavaisoflavone, methylophiopogonone A, 3'-Methoxydaidzin, and genistin against S2 were -8.1, -8.2, -8.3, and -8.3 kcal/mol, respectively, which is far below of the control (pravastatin, -7.3 kcal/mol). Instead of conventional hydrogen bonding, the π bonding influenced the binding affinity against S2. The results reveal that this is the first report about methylophiopogonone A, 3'-Methoxydaidzin, and genistin as candidates for anti-viral agents. Those compounds can then be further explored and used as a parent backbone molecule to develop a new supplementation for preventing SARS-CoV-2 infections during COVID-19 outbreaks.

Quantifying biomolecular hydrophobicity: Single molecule force spectroscopy of class II hydrophobins

Hydrophobins are surface-active proteins produced by filamentous fungi. The amphiphilic structure of hydrophobins is very compact, containing a distinct hydrophobic patch on one side of the molecule, locked by four intramolecular disulfide bridges. Hydrophobins form dimers and multimers in solution to shield these hydrophobic patches from water exposure. Multimer formation in solution is dynamic, and hydrophobin monomers can be exchanged between multimers. Unlike class I hydrophobins, class II hydrophobins assemble into highly ordered films at the air-water interface. In order to increase our understanding of the strength and nature of the interaction between hydrophobins, we used atomic force microscopy for single molecule force spectroscopy to explore the molecular interaction forces between class II hydrophobins from Trichoderma reesei under different environmental conditions. A genetically engineered hydrophobin variant, NCys-HFBI, enabled covalent attachment of proteins to the apex of the atomic force microscopy cantilever tip and sample surfaces in controlled orientation with sufficient freedom of movement to measure molecular forces between hydrophobic patches. The measured rupture force between two assembled hydrophobins was ∼31 pN, at a loading rate of 500 pN/s. The results indicated stronger interaction between hydrophobins and hydrophobic surfaces than between two assembling hydrophobin molecules. Furthermore, this interaction was stable under different environmental conditions, which demonstrates the dominance of hydrophobicity in hydrophobin-hydrophobin interactions. This is the first time that interaction forces between hydrophobin molecules, and also between naturally occurring hydrophobic surfaces, have been measured directly at a single-molecule level.

Desiccation-induced cell damage in bacteria and the relevance for inoculant production

Plant growth-promoting bacteria show great potential for use in agriculture although efficient application remains challenging to achieve. Cells often lose viability during inoculant production and application, jeopardizing the efficacy of the inoculant. Since desiccation has been documented to be the primary stress factor affecting the decrease in survival, obtaining xerotolerance in plant growth-promoting bacteria is appealing. The molecular damage that occurs by drying bacteria has been broadly investigated, although a complete view is still lacking due to the complex nature of the process. Mechanic, structural, and metabolic changes that occur as a result of water depletion may potentially afflict lethal damage to membranes, DNA, and proteins. Bacteria respond to these harsh conditions by increasing production of exopolysaccharides, changing composition of the membrane, improving the stability of proteins, reducing oxidative stress, and repairing DNA damage. This review provides insight into the complex nature of desiccation stress in bacteria in order to facilitate strategic choices to improve survival and shelf life of newly developed inoculants. KEY POINTS: Desiccation-induced damage affects most major macromolecules in bacteria. Most bacteria are not xerotolerant despite multiple endogenous adaption mechanisms. Sensitivity to drying severely hampers inoculant quality.

Alternative Hydrophobic Core in Proteins—The Effect of Specific Synergy

Proteins with a high degree of sequence similarity representing different structures provide a key to understand how protein sequence codes for 3D structure. An analysis using the fuzzy oil drop model was carried out on two pairs of proteins with different secondary structures and with high sequence identities. It has been shown that distributions of hydrophobicity for these proteins are approximated well using single 3D Gaussian function. In other words, the similar sequences fold into different 3D structures, however, alternative structures also have symmetric and monocentric hydrophobic cores. It should be noted that a significant change in the helical to beta-structured form in the N-terminal section takes places in the fragment much preceding the location of the mutated regions. It can be concluded that the final structure is the result of a complicated synergy effect in which the whole chain participates simultaneously.

The Amyloid as a Ribbon-Like Micelle in Contrast to Spherical Micelles Represented by Globular Proteins

Selected amyloid structures available in the Protein Data Bank have been subjected to a comparative analysis. Classification is based on the distribution of hydrophobicity in amyloids that differ with respect to sequence, chain length, the distribution of beta folds, protofibril structure, and the arrangement of protofibrils in each superfibril. The study set includes the following amyloids: Aβ (1-42), which is listed as Aβ (15-40) and carries the D23N mutation, and Aβ (11-42) and Aβ (1-40), both of which carry the E22Δ mutation, tau amyloid, and α-synuclein. Based on the fuzzy oil drop model (FOD), we determined that, despite their conformational diversity, all presented amyloids adopt a similar structural pattern that can be described as a ribbon-like micelle. The same model, when applied to globular proteins, results in structures referred to as "globular micelles," emerging as a result of interactions between the proteins' constituent residues and the aqueous solvent. Due to their composition, amyloids are unable to attain entropically favorable globular forms and instead attempt to limit contact between hydrophobic residues and water by producing elongated structures. Such structures typically contain quasi hydrophobic cores that stretch along the fibril's long axis. Similar properties are commonly found in ribbon-like micelles, with alternating bands of high and low hydrophobicity emerging as the fibrils increase in length. Thus, while globular proteins are generally consistent with a 3D Gaussian distribution of hydrophobicity, the distribution instead conforms to a 2D Gaussian distribution in amyloid fibrils.

Quo Vadis Biomolecular NMR Spectroscopy?

In-cell nuclear magnetic resonance (NMR) spectroscopy offers the possibility to study proteins and other biomolecules at atomic resolution directly in cells. As such, it provides compelling means to complement existing tools in cellular structural biology. Given the dominance of electron microscopy (EM)-based methods in current structure determination routines, I share my personal view about the role of biomolecular NMR spectroscopy in the aftermath of the revolution in resolution. Specifically, I focus on spin-off applications that in-cell NMR has helped to develop and how they may provide broader and more generally applicable routes for future NMR investigations. I discuss the use of ‘static’ and time-resolved solution NMR spectroscopy to detect post-translational protein modifications (PTMs) and to investigate structural consequences that occur in their response. I argue that available examples vindicate the need for collective and systematic efforts to determine post-translationally modified protein structures in the future. Furthermore, I explain my reasoning behind a Quinary Structure Assessment (QSA) initiative to interrogate cellular effects on protein dynamics and transient interactions present in physiological environments.

A protein interaction free energy model based on amino acid residue contributions: Assessment of point mutation stability of T4 lysozyme

Here we present a model to estimate the interaction free energy contribution of each amino acid residue of a given protein. Protein interaction energy is described in terms of per-residue interaction factors, μ. Multibody interactions are implicitly captured in μ through the combination of amino acid terms (γ) guided by local conformation indices (σ). The model enables construction of an interaction factor heat map for a protein in a given fold, allows prima facie assessment of the degree of residue-residue interaction, and facilitates a qualitative and quantitative evaluation of protein association properties. The model was used to compute thermal stability of T4 bacteriophage lysozyme mutants across seven sites. Qualitative assessment of mutational effects provides a straightforward rationale regarding whether a particular site primarily perturbs native or non-native states, or both. The presented model was found to be in good agreement with experimental mutational data (R 2 = 0.73) and suggests an approach by which to convert structure space into energy space.

How medicinal chemists learned about log P

Although log P is now recognized to be a key factor that determines the bioactivity of a molecule, the focus of medicinal chemists on hydrophobicity and log P started with the quantitative structure-activity relationships (QSAR) publications of Hansch and Fujita. Their original publication represents a dramatic change of focus to incorporate consideration of log P after a decade of work unsuccessfully attempting to use the Hammett equation to explain the structure-activity relationships of plant growth regulators. QSAR allows one to explore the quantitative relationship between log P and biological activity even when other factors also influence potency. In particular, Hansch's publications of thousands of QSAR equations demonstrate that a relationship of biological activity with log P is indeed a general phenomenon. Hansch's group also provided data and tools that enable others to explore the relationship between log P and the biological activity of compounds of interest.

The Excess Chemical Potential of Water at the Interface with a Protein from End Point Simulations

We use end point simulations to estimate the excess chemical potential of water in the homogeneous liquid and at the interface with a protein in solution. When the pure liquid is taken as the reference, the excess chemical potential of interfacial water is the difference between the solvation free energy of a water molecule at the interface and in the bulk. Using the homogeneous liquid as an example, we show that the solvation free energy for growing a water molecule can be estimated by applying UWHAM to the simulation data generated from the initial and final states (i.e., "the end points") instead of multistate free energy perturbation simulations because of the possible overlaps of the configurations sampled at the end points. Then end point simulations are used to estimate the solvation free energy of water at the interface with a protein in solution. The estimate of the solvation free energy at the interface from two simulations at the end points agrees with the benchmark using 32 states within a 95% confidence interval for most interfacial locations. The ability to accurately estimate the excess chemical potential of water from end point simulations facilitates the statistical thermodynamic analysis of diverse interfacial phenomena. Our focus is on analyzing the excess chemical potential of water at protein receptor binding sites with the goal of using this information to assist in the design of tight binding ligands.

Biophysical studies of protein solubility and amorphous aggregation by systematic mutational analysis and a helical polymerization model

At concentrations above solubility, a protein aggregates, most often into amorphous aggregates, and loses its function. However, unlike amyloidogenic aggregates, which are β-sheeted fibrillar aggregates often related to neurodegenerative diseases, amorphous aggregates, where proteins aggregate/oligomerize without forming specific high-order structures, are rarely the focus of biophysical studies. Hence, protein solubility with respect to amorphous aggregation remains to be fully characterized from a biophysical viewpoint. Here, I briefly describe the structural nature of proteins in amorphous aggregates before discussing systematic mutational analyses that aim to rationalize the contribution of individual amino acids to the solubility of a protein. The discussion is expected to demonstrate that protein solubility, and, accordingly, amorphous aggregation, can be understood using thermodynamic and biophysical rationales similar to those used in the study of protein stability or, more recently, amyloidogenesis. Finally, I will argue that the mathematical formalism of the helical polymerization model (HPM) proposed by Oosawa, Kasai, and Asakura's group can be readily adapted to provide a thermodynamic description of a system containing amorphous aggregates and soluble particles. The HPM and HPM-derived models imply the presence of nuclei or seeds for amorphous aggregates, similar to those hypothesized in crystallogenesis and amyloidogenesis.

Interaction of organic ions with proteins

In this study we have investigated how different proteins interact with big organic ions. Two ions that are similar in size and chemical structure (Ph₄B^- anion and Ph₄As⁺ cation) were studied. The proteins chosen are the two major allergenic proteins of cow's milk, β-lactoglobulin and β-casein, and bovine serum albumin, BSA, as the reference protein. First, a quantitative study to determine the hydrophobic degree of the proteins was performed. Then, electrokinetic and stability measurements on protein-coated polystyrene (PS) microspheres as a function of the tetraphenyl ion concentration were carried out. Our results show that the affinity of the organic ions depends on the hydrophobicity of the interface. Big charge inversions and re-stabilization patterns were observed at very low concentrations of tetraphenyl ions for the most hydrophobic protein studied (with β-casein). Besides, the ionic concentrations needed to destabilize these colloidal systems were roughly one order of magnitude lower for the anion than for the cation. In addition, we studied conformational changes of the adsorbed proteins with a quartz crystal microbalance. Proteins were adsorbed onto hydrophobic flat substrates and then exposed to the tetraphenyl ions. The protein films swelled or collapsed as a function of the accumulation of tetraphenyl ions. Similarly to the electrokinetic/stability studies, the ionic concentration necessary to trigger structural changes of the protein films was one order of magnitude larger for the cation than for the anion. All the results evidence that the accumulation of these organic ions on an interface depends directly on its degree of hydrophobicity. We attribute the different interactions of the anion and the cation with these interfaces to their dissimilar hydration, which makes the anion show a more hydrophobic behaviour than the cation.

Hydrophobic Association in Mixed Urea-TMAO Solutions

The formation of a hydrophobic core is key to the folding and resulting function of most proteins in the cell. In several organisms, as well as in many in vitro experiments, protein folding is modulated by the presence of osmolytes, but the mechanism by which hydrophobic association occurs is not well understood. We present a study of the solvation thermodynamics of hydrophobic self-association in mixed-osmolyte urea-TMAO solutions, with neopentane as a model hydrophobic molecule. Using molecular dynamics simulations and the Kirkwood-Buff theory of solutions, we show that a sensitive balance between the TMAO-water and the TMAO-urea interactions governs the osmolyte-induced changes in hydrophobic association in mixed urea-TMAO solutions. This balance must be correctly incorporated in force-field parametrization because hydrophobic association can be either enhanced or prevented all together by slightly increasing or decreasing the osmolyte-water affinity and osmolyte-osmolyte self-affinity of TMAO molecules.

Reconciling the understanding of 'hydrophobicity' with physics-based models of proteins

The idea that a 'hydrophobic energy' drives protein folding, aggregation, and binding by favoring the sequestration of bulky residues from water into the protein interior is widespread. The solvation free energies (ΔGsolv) of small nonpolar solutes increase with surface area (A), and the free energies of creating macroscopic cavities in water increase linearly with A. These observations seem to imply that there is a hydrophobic component (ΔGhyd) of ΔGsolv that increases linearly with A, and this assumption is widely used in implicit solvent models. However, some explicit-solvent molecular dynamics studies appear to contradict these ideas. For example, one definition (ΔG(LJ)) of ΔGhyd is that it is the free energy of turning on the Lennard-Jones (LJ) interactions between the solute and solvent. However, ΔG(LJ) decreases with A for alanine and glycine peptides. Here we argue that these apparent contradictions can be reconciled by defining ΔGhyd to be a near hard core insertion energy (ΔGrep), as in the partitioning proposed by Weeks, Chandler, and Andersen. However, recent results have shown that ΔGrep is not a simple function of geometric properties of the molecule, such as A and the molecular volume, and that the free energy of turning on the attractive part of the LJ potential cannot be computed from first-order perturbation theory for proteins. The theories that have been developed from these assumptions to predict ΔGhyd are therefore inadequate for proteins.

Molecular modeling of ions at interfaces: exploring similarities to hydrophobic solvation through the lens of induced aqueous interfacial fluctuations

Ion specific effects are ubiquitous in chemistry and biology.

Synthesis, insecticidal activities and structure-activity relationship studies of novel anthranilic diamides containing pyridylpyrazole-4-carboxamide

BACKGROUND

Anthranilic diamide insecticides containing pyridylpyrazole-5-carboxamide are extremely important in modern agriculture. New structurally modified compounds with high insecticidal activity were discovered by designing a series of novel pyridylpyrazole-4-carboxamides (9I to 9IV) and pyridylpyrazole-4-carboxamides (10I to 10IV), the latter designed by the cyclisation of two amides. The structure-activity relationship (SAR) between the two series is of interest.

RESULTS

Two series of novel anthranilic diamides containing pyridylpyrazole-4-carboxamide were synthesised and characterised via melting point, (1)H NMR, (13)C NMR, MS and elemental analyses. The insecticidal activities of these compounds against Plutella xylostella were evaluated. At a concentration of 100 mg L(-1), the compounds with unmodified amide moieties (9I to 9IV) exhibited much better larvicidal activities than the other derivative compounds (10I to 10IV). Most of the compounds 9I to 9IV showed over 90% larvicidal activity at 100 mg L(-1). Furthermore, compounds 9IIIa, 9IIIc, 9IIId and 9IVd displayed significant insecticidal activity at 10 mg L(-1). Density functional theory (DFT) calculation was carried out to provide more information regarding SAR.

CONCLUSION

Thirty-two new anthranlic diamides containing pyridylpyrazole-4-carboxamide were designed and obtained. SAR analysis and DFT calculation results revealed that the amide moiety had a very important effect on bioactivity. This work has provided information that could aid investigations of novel insecticides.

Salting-Out of Methane in the Aqueous Solutions of Urea and Glycine-Betaine

We have studied the hydrophobic association and solvation of methane molecules in aqueous solutions of urea and glycine betaine (GB). We have calculated the potentials of mean force (PMFs) between methane molecules in water, aqueous GB, aqueous urea and aqueous urea-GB mixtures. The PMFs and equilibrium constants indicate that both urea and GB increase the hydrophobic association of methane. Calculation of thermodynamic parameters shows that the association of methane is stabilized by entropy whereas solvation is favored by enthalpy. In the case of the water-urea-GB mixture, both hydrophobic association and solvation are stabilized by entropy. From the investigation of radial distribution functions, running coordination numbers and excess coordination numbers, we infer that both urea and GB are preferentially excluded from methane surface in the mixtures of osmolytes and methane is preferentially solvated by water molecules in all the mixtures. The favorable exclusion of both urea and GB from the methane surface suggests that both urea and GB increase the interaction between methane molecules, i.e., salting-out of methane. We observe that addition of both urea and GB to water enhances local water structure. The calculated values of diffusion constants of water also suggest enhanced water-water interactions in the presence of urea and GB. The calculated free energies of methane in these mixtures show that methane is less soluble in the mixtures of urea and GB than in water. The data on solvation free energies support the observations obtained from the PMFs of methane molecules.

Distance-dependent hydrophobic-hydrophobic contacts in protein folding simulations

Successful prediction of protein folding from an amino acid sequence is a challenge in computational biology. In order to reveal the geometric constraints that drive protein folding, highlight those constraints kept or missed by distinct lattices and for establishing which class of intra- and inter-secondary structure element interactions is the most relevant for the correct folding of proteins, we have calculated inter-alpha carbon distances in a set of 42 crystal structures consisting of mainly helix, sheet or mixed conformations. The inter-alpha carbon distances were also calculated in several lattice "hydrophobic-polar" models built from the same protein set. We found that helix structures are more prone to form "hydrophobic-hydrophobic" contacts than beta-sheet structures. At a distance lower than or equal to 3.8 Å (very short-range interactions), "hydrophobic-hydrophobic" contacts are almost absent in the native structures, while they are frequent in all the analyzed lattice models. At distances in-between 3.8 and 9.5 Å (short-/medium-range interactions), the best performing lattice for reproducing mainly helix structures is the body-centered-cubic lattice. If protein structures contain sheet portions, lattice performances get worse, with few exceptions observed for double-tetrahedral and body-centered-cubic lattices. Finally, we can observe that ab initio protein folding algorithms, i.e. those based on the employment of lattices and Monte Carlo simulated annealings, can be improved simply and effectively by preventing the generation of "hydrophobic-hydrophobic" contacts shorter than 3.8 Å, by monitoring the "hydrophobic-hydrophobic/polar-polar" contact ratio in short-/medium distance ranges and by using preferentially a body-centered-cubic lattice.

Low polarity water, a novel transition species at the polyethylene–water interface

Polyethylene sandwich or a single window cell: dark bars represent polyethylene (PE) windows (double or single) of a cell. Wavy lines are water (W) and low polarity water (LPW). Subtraction of the single window spectrum from the double window spectrum leaves the LPW spectrum as illustrated in the figure.

Specific ion effects on the hydrophobic interaction of benzene self-assembled monolayers

The ions, Ca²⁺, Mg²⁺, Na⁺ and K⁺, decrease the hydrophobic attraction (in this order) between benzene-terminated self assembled monolayers by affecting the creation of bridging capillaries and by charging the surfaces.

Electrostatic contribution from solvent in modulating single-walled carbon nanotube association

We perform all-atom molecular dynamics simulations to compute the potential of mean force (PMF) between two (10,10) single-walled carbon nanotubes solvated in pure nonpolarizable SPC/E and polarizable TIP4P-FQ water, at various temperatures. In general, the reversible work required to bring two nanotubes from a dissociated state (free energy reference) to contact state (free energy minimum) is more favorable and less temperature-dependent in TIP4P-FQ than in SPC/E water models. In contrast, molecular properties and behavior of water such as the spatially-resolved water number density (intertube, intratube, or outer regions), for TIP4P-FQ are more sensitive to temperature than SPC/E. Decomposition of the solvent-induced PMF into different spatial regions suggests that TIP4P-FQ has stronger temperature dependence; the opposing destabilizing/stabilizing contributions from intertube water and more distal water balance each other and suppress the temperature dependence of total association free energy. Further investigation of hydrogen bonding network in intertube water reveals that TIP4P-FQ retains fewer hydrogen bonds than SPC/E, which correlates with the lower water number density in this region. This reduction of hydrogen bonds affects the intertube water dipoles. As the intertube volume decreases, TIP4P-FQ dipole moment approaches the gas phase value; the distribution of dipole magnitude also becomes narrower due to less average polarization/perturbation from other water molecules. Our results imply that the reduction of water under confinement may seem trivial, but underlying effects to structure and free energetics are non-negligible.

A chiral, low-cytotoxic [Ni15]-wheel complex

A new chiral [Ni15] complex with a Schiff-base ligand derived from o-vanillin and L-glutamic acid is presented, emphasizing the properties relevant for biology and materials science. The formation of the complex molecules in solution is confirmed by AFM and dynamic light scattering studies. The compound is weakly antiferromagnetic with considerable admixture of excited states, comprising negligibly interacting [Ni3] units. Studies of the interactions with two cell lines indicate low cytotoxicity.

Forces stabilizing proteins

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site-directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a -CH2- group on folding contributes 1.1±0.5 kcal/mol to protein stability. (2) The burial of non-polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1±0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.

Amino acid distribution rules predict protein fold

In the present article, we provide a brief overview of the main approaches to analysing the sequence-structure relationship of proteins and outline a novel method of structure prediction. The proposed method involves finding a set of rules that describes a correlation between the distribution of residues in a sequence and the essential structural characteristics of a protein structure. The residue distribution rules specify the 'favourable' residues that are required in certain positions of a polypeptide chain in order for it to assume a particular protein fold, and the 'unfavourable' residues incompatible with the given fold. Identification of amino acid distribution rules derives from examination of inter-residue contacts. We describe residue distribution rules for a large group of β-sandwich-like proteins characterized by a specific arrangement of strands in their two β-sheets. It was shown that this method has very high accuracy (approximately 85%). The advantage of the residue rule approach is that it makes possible prediction of protein folding even in polypeptide chains that have very low global sequence similarities, as low as 18%. Another potential benefit is that a better understanding of which residues play essential roles in a given protein fold may facilitate rational protein engineering design.

Patchy supramolecules as versatile tools to probe hydrophobicity in nanoglobular systems

We describe precise supramolecules that enable the evaluation of the effective hydrophobicity of amphiphilic or "patchy" nanoglobular systems. These supramolecules exhibit the lower critical solution temperature phenomenon, which provides a quantitative measure of their effective hydrophobicity. Specifically, two isomeric 8-aryl-2'-deoxyguanosine derivatives with a transposed pair of methylene groups self-assemble into hexadecameric nanoglobular supramolecular G-quadruplexes (SGQs) that show large differences in their transition temperatures as determined by turbidity and differential scanning calorimetry studies. Molecular modeling studies suggested that differential clustering of the hydrophobic patches on the surface is responsible for the striking differences between the two isomeric supramolecules.