Mesoscopic analysis of conformational and entropic contributions to nonspecific adsorption of HP copolymer chains using dynamic Monte Carlo simulations

Dynamic Monte Carlo simulations of short linear HP-type copolymers exhibiting proteinlike characteristics are used to investigate both chain dynamics and changes in chain conformational entropy and their contributions to the energetics of adsorption onto a solid-liquid interface. The dMC results show that the conformations and energies of adsorbed chains are highly degenerate. The ensemble-averaged energy of the adsorbed state is dependent on temperature, chain sequence, native-state stability, and sorbent surface geometry and hydrophobicity. Mesoscopic thermodynamic analyses reveal that, although increased chain conformational entropy contributes to the driving force for adsorption in certain cases, many conditions exist where the change in conformational entropy is either negligible or unfavorable due to constraints imposed by the need to form a large and specific number of favorable intra- and intermolecular contacts and by the impenetrable nature of the sorbent surface. Step-number-averaged energy trajectories, based on sampling of a large number of energy trajectories and thus conformational states at each step number, suggest that the search for a global energy minimum is gradual, so that adsorption is first reversible but becomes apparently irreversible with longer exposure to the sorbent. These results appear to be connected to the conformational adaptability of the chain both on the surface and in solution, and an adsorption model taking chain conformational dynamics into account is proposed.     

Mesoscopic dynamic Monte Carlo simulations of the adsorption of proteinlike HP chains within laterally constricted spaces

Two-dimensional dynamic Monte Carlo simulations are applied to the protein-like HP chain model to investigate the influence of lateral confinement of the adsorbed chain on adsorption thermodynamics and the ensemble of accessible chain conformations. The structure of the model makes it possible to enumerate all possible chain conformations and thereby define with precision the relation between adsorption thermodynamics and changes in accessible chain conformations resulting from the adsorption process. Lateral confinement of the adsorbed chain is shown to dramatically reduce the number of accessible energy states and unique chain conformations such that, under certain conditions, adsorption is predicted to actually stabilize the chain against denaturation. Lateral confinement preferentially eliminates expanded conformations of the adsorbed chain, shifting the equilibrium from the unfolded state toward the native state. As a result, the conformational entropy of the adsorbed chain is predicted to be lower than that of the chain free in solution. The protein-like HP chain responds to an increase in the hydrophobicity of the sorbent surface by strongly favoring those conformations that minimize the overall internal energy of the system. As a result, adsorption severely destabilizes the native-state conformation. The ability of our simulation results to provide insights into underlying mechanisms for nonspecific protein adsorption is illustrated through qualitative comparison with activity data for hen egg-white lysozyme adsorbed on silica at different surface concentrations.     

Determination of the adsorption free energy for peptide-surface interactions by SPR spectroscopy

To understand and predict protein adsorption behavior, we must first understand the fundamental interactions between the functional groups presented by the amino acid residues making up a protein and the functional groups presented by the surface. Limited quantitative information is available, however, on these types of submolecular interactions. The objective of this study was therefore to develop a reliable method to determine the standard state adsorption free energy (delta Go ads) of amino acid residue-surface interactions using surface plasma resonance (SPR) spectroscopy. Two problems are commonly encountered when using SPR for peptide adsorption studies: the need to account for "bulk-shift" effects and the influence of peptide-peptide interactions at the surface. Bulk-shift effects represent the contribution of the bulk solute concentration to the SPR response that occurs in addition to the response due to adsorption. Peptide-peptide interactions, which are assumed to be zero for Langmuir adsorption, can greatly skew the isotherm shape and result in erroneous calculated values of delta Go ads. To address these issues, we have developed a new approach for the determination of delta Go ads using SPR that is based on the chemical potential. In this article, we present the development of this new approach and its application for the calculation of delta Go ads for a set of peptide-surface systems where the peptide has a host-guest amino acid sequence of TGTG-X-GTGT (where G and T are glycine and threonine residues and X represents a variable residue) and the surface consists of alkanethiol self-assembled monolayers (SAMs) with methyl (CH 3) and hydroxyl (OH) functionality. This new approach enables bulk-shift effects to be directly determined from the raw SPR versus peptide concentration data plots and the influence of peptide-peptide interaction effects to be minimized, thus providing a very straightforward and accurate method for the determination of delta Go ads for peptide adsorption. Further studies are underway to characterize delta Go ads for a large library of peptide-SAM combinations.     

The influence of intramolecular hydrogen bonds on the adsorption properties of aromatic alcohols and thiols

The thermodynamic characteristics of adsorption equilibria of primary aromatic alcohols and thiols on graphitized carbon black at 300 K were determined by the molecular-statistical method. The influence of intramolecular effects under the action of the force field of the sorbent on the conformation of the molecules studied was considered in comparison with n-alkylbenzenes having similar structures. An increase in the chain length by one-CH₂ group was shown to influence the formation of intramolecular H-bonds in alcohol and thiol molecules. In adsorption, a considerable fraction of molecules assumed configurations close to planar. An exception was the nonplanar conformation of the 2-phenylethanol molecule stabilized by the intramolecular H-bond.     

Photochemistry and photophysics of organomonosilane and oligosilanes: Updating their studies on conformation and intramolecular interactions

This review article deals with recent studies on photochemistry and photophysics of organosilicon monomer to oligomers (here, chain length of Si is mostly less than 10). Arylmonosilanes have been caught attention as the fluorescent and charge-transporting materials. Increasing in silicon chain length, an importance of conformation for silicon single chain in organooligosilane is significant. The conformation affects the spectroscopic, photochemical, and photophysical characteristics. This review is mainly focused on the properties of arylsilanes and α,ω-diaryloligosilanes in which three types of intramolecular interactions, the σ–σ interaction between silicon atoms in a chain unit, the π–π interaction between two aromatic groups, and the σ–π interaction between an aromatic group and a silicon chain unit are expected. As a result, intramolecular excimer/exciplex formation, charge-transfer (CT) complex formation, excitation energy or photoinduced electron transfers were observed by the steady-state and the time-resolved fluorescence (TR-FL) spectroscopy. For the better understanding, the recent conformational studies of single oligosilane chain are also reviewed.     

Bond dissocation and conformational energetics of tetrasulfur: a quantum Monte Carlo study

Variational Monte Carlo (VMC) and fixed-node diffusion Monte Carlo (DMC) calculations are performed for S4. The effect of single- and multireference trial functions, as well as choice of orbitals, is investigated for its effect on the quality of the Monte Carlo estimates. Estimates of symmetric (two S2 molecules) and asymmetric (S atom and S3 molecule) bond dissociation are reported. The conformational change of S4 from C2v to D2h defines a double-well potential and is also estimated. Multireference DMC with natural orbitals (DMC/NO) estimates the energy of the conformational change as 1.20(20) kcal/mol; the dissociation of the long S-S single bond is estimated at 21.1(1.3) kcal/mol, and the asymmetric bond energy is estimated as 53.2(2.4) kcal/mol. An estimate of the total atomization energy using multireference DMC/NO gives a value of 219.5(2.2) kcal/mol. The relative quality of result and implications for simplified trial function design are discussed.     

Conformational and spectroscopic analysis of the tyrosyl radical dipeptide analogue in the gas phase and in aqueous solution by a density functional/continuum solvent model

The conformational and spectroscopic properties of the tyrosyl radical dipeptide analogue (T(R)DA) are investigated both in gas phase and in aqueous solution by means of density functional calculations. Electronic interactions between backbone and side chain, determining the relative stability of the different energy minimums, depend on the electronic state of the phenoxy substituent. As a consequence, (i) the conformational behavior of T(R)DA is quite different from that of the tyrosine dipeptide analogue, and (ii) the energy required for the homolytic breaking of the OH bond depends on the adopted conformation. The calculated hyperfine coupling constants are in good agreement with the available experimental results. Side-chain-backbone interactions cause an asymmetrization of the magnetic properties of the phenoxy ring and deviations from McConnell relationship. Solvent effects, taken into account by means of a combined discrete/continuum model, significantly affect both the conformational and the magnetic behavior of T(R)DA.     

Spectroscopic examination of chain conformation and bonding in poly(phenylene oxide)–polystyrene blends

Ultraviolet and Fourier-transform infrared spectroscopy are used to characterize the chain conformation of 2,6-dimethyl poly(phenylene oxide) (2MPPO) in compatible 2MPPO-polystyrene (PS) blends. Blending with PS increases the intramolecular energy of the 2MPPO. Atactic PS induces a greater distortion from the minimum energy conformation of 2MPPO than isotactic PS. The dispersion interaction between the phenyl ring of PS and the phenylene ring of 2MPPO is found to be quite strong. These large intermolecular forces are thought to supercede the increase in 2MPPO conformational energy and to be responsible for the blend compatibility.     

Underpotential deposition of hydrogen on benzene-modified Pt(111) in aqueous H2SO4

The Pt(111) electrode is modified by an overlayer of C6H6 (ads) upon its cycling in the 0.05-0.80 V range in aq H2SO4 + 1 mM C6H6. The C6H6 (ads) overlayer significantly changes the underpotential-deposited H (H(UPD)) and anion adsorption, and cyclic-voltammetry (CV) profiles show a sharp cathodic peak and an asymmetric anodic one in the 0.05-0.80 V potential range. The C6H6 (ads) layer blocks the (bi)sulfate adsorption but facilitates the adsorption of one monolayer of H(UPD). Cycling of the benzene-modified Pt(111) in benzene-free aq 0.05 H2SO4 from 0.05 to 0.80 V results in a partial desorption of C6H6 (ads) and in a partial recovery of the CV profile characteristic of an unmodified Pt(111). The peak potential of the cathodic and anodic feature is independent of the scan rate, s (10 < or = s < or = 100 mV s(-1)), and the peak current density increases linearly with an increase of the scan rate. Temperature variation modifies the peak potential and current density but does not affect the charge density of the cathodic or anodic feature. Temperature-dependent studies allow us to determine the thermodynamic state function for the H(UPD) adsorption and desorption. Delta G degrees(ads)(H(UPD))assumes values from -4 to -12 kJ mol(-1), while has values from 9 to 14 kJ mol(-1). The values of delta Delta G degrees (delta Delta G degrees = delat Delta G degrees(ads) + delta Delta D degrees(des)) decrease almost linearly from 6 kJ mol(-1) at theta(H(UPD) --> 0 to 0 kJ mol(-1) at theta(H(UPD) --> 1. The nonzero values of delta Delta G degrees testify that the adsorbing and desorbing H(UPD) adatoms interact with an energetically different substrate. The lateral interactions changed from repulsive (omega = 29 kJ mol(-1) at theta(H(UPD) --> 0) to attractive (omega = -28 kJ mol(-1) at theta(H(UPD) --> 1) as the H(UPD) coverage increases. The values of delta S degrees(ads)(H(UPD)) increase from 19 to 56 J K(-1) mol(-1), while those of delta S degrees(des)(H(UPD)) decrease from 45 to -30 J K(-1) mol(-1) with an increase of H(UPD) coverage. The values of delta H degrees(des)(H(UPD)) and delta H degrees(des)(H(UPD)) vary from 0 to 27 kJ mol(-1). The Pt(111)-H(UPD) surface bond energy at the benzene-modified Pt(111) electrode falls in the 191-218 kJ mol(-1) range and is weaker than in the case of the unmodified Pt(111) electrode in the same electrolyte.     

Adsorption of organic molecules on a porous polymer surface modified with the supramolecular structure of melamine–cyanuric acid

The adsorption of organic molecules on the surface of a porous polymeric sorbent modified with a mixed cyanuric acid–melamine supramolecular structure is studied. The parameters of thermodynamic adsorption are considered and the contributions from intermolecular interactions to the Helmholtz energy of adsorption are assessed. Analysis of the molar changes in internal energy and adsorption entropy shows that the supramolecular structure formed on the surface could not exhibit dimension effects, indicating there were no cavities. The contributions from nonspecific interactions to the Helmholtz energy of adsorption generally fall, while those of specific interactions increase, indicating an increase in the polarity of the sorbent surface.

     

A robust force field based method for calculating conformational energies of charged drug-like molecules

The binding affinity of a drug-like molecule depends among other things on the availability of the bioactive conformation. If the bioactive conformation has a significantly higher energy than the global minimum energy conformation, then the molecule is unlikely to bind to its target. Determination of the global minimum energy conformation and calculation of conformational penalties of binding is a prerequisite for prediction of reliable binding affinities. Here, we present a simple and computationally efficient procedure to estimate the global energy minimum for a wide variety of structurally diverse molecules, including polar and charged compounds. Identifying global energy minimum conformations of such compounds with force field methods is problematic due to the exaggeration of intramolecular electrostatic interactions. We demonstrate that the global energy minimum conformations of zwitterionic compounds generated by conformational analysis with modified electrostatics are good approximations of the conformational distributions predicted by experimental data and with molecular dynamics performed in explicit solvent. Finally the method is used to calculate conformational penalties for zwitterionic GluA2 agonists and to filter false positives from a docking study.     

Interfacial properties of poly(maleic acid-alt-1-alkene) disodium salts at water/hydrocarbon interfaces

The interfacial properties of poly(maleic acid-alt-1-alkene) disodium salts at hydrocarbon/water interfaces are determined. In all the studied systems, the interfacial tension decreases markedly with the polyelectrolyte concentration as the side-chain length increases. The results of the standard free energy of adsorption, DeltaG(ads)(0), are a linear function of the number of carbon atoms in the polyelectrolyte side chain. The contribution to DeltaG(ads)(0) per mol of methylene group varies from -0.64 to -0.52 kJ/mol for the n-octane/water to n-dodecane/water interfaces. DeltaG(ads)(0) data also reveal that the adsorption process is mainly determined by adsorption efficiency. Comparatively, the adsorption effectiveness seems to play a less important role. The theoretical interaction energies calculated for the insertion of one hydrocarbon molecule into the space formed by two neighboring polyelectrolyte side chains are in good agreement with the experimental results. The latter results are consistent with van der Waals-type interactions between the hydrocarbon molecules and the polyelectrolyte side chains.     

Conformational analysis of TMC114, a novel HIV-1 protease inhibitor

TMC114, a potent novel HIV-1 protease inhibitor, remains active against a broad spectrum of mutant viruses. In order to bind to a variety of mutants, the compound needs to make strong, preferably backbone, interactions and have enough conformational flexibility to adapt to the changing geometry of the active site. The conformational analysis of TMC114 in the gas phase yielded 43 conformers in which five types of intramolecular H-bond interactions could be observed. All 43 conformers were subject to both rigid and flexible ligand docking in the wild-type and a triple mutant (L63P/V82T/I84V) of HIV-1 protease. The largest binding energy was calculated for the conformations that are close to the conformation observed in the X-ray complexes of TMC114 and HIV-1 protease.     

Molecular modeling of polymers: I. Correct and efficient enumeration of intrachain conformational energetics

Polymer conformational analyses can require being able to model the intramolecular energetics of a very long (infinite) chain employing calculations carried out on a relatively short chain sequence. A method to meet this need, based upon symmetry considerations and molecular mechanics energetics, has been developed. Given N equivalent degrees of freedom in a linear polymer chain, N unique molecular groups are determined within the chain. A molecular unit is defined as a group of atoms containing backbone rotational degrees of conformational freedom on each of its ends. The interaction of these N molecular groups, each with a finite number of nearest neighbors, properly describe the intramolecular energetics of a long (infinite) polymer chain. Thus, conformational energetics arising from arbitrarily distant neighbor interactions can be included in the estimation of statistical and thermodynamic properties of a linear polymeric system. This approach is called the polymer reduced interaction matrix method (PRIMM) and the results of applying it to isotactic polystyrene (I-PS) are presented by way of example.     

The molecular and crystal structure of tert-butyl Nalpha-tert-butoxycarbonyl-L-(S-trityl)cysteinate and the conformation-stabilizing function of weak intermolecular bonding

The title compound, C31H37NO4S [systematic name: (R)-tert-butyl-2-[(tert-butoxycarbonyl)amino]-3-(tritylsulfanyl)propanoate] is an L-cysteine derivative with three functions: NH2, COOH and SH, blocked by protecting groups tert-butoxycarbonyl, tert-butyl and trityl, respectively. The main chain of the molecule adopts the extended, nearly all-trans C5 conformation with the intramolecular N-H...O=C hydrogen bond. The urethane group is not involved in any intermolecular hydrogen bonding. Only weak intermolecular hydrogen bonds and hydrophobic contacts are observed in the crystal structure. These are C-H...O hydrogen bonds and CH/pi interactions with donor...acceptor distances, C...O ca. 3.5 A and C...C ca. 3.7 A, respectively. The first type of interaction links phenyl H-atoms and carbonyl groups. The second type of interaction is formed between a methyl group of the tert-butyl fragment and a trityl phenyl ring. The resulting molecular conformation in the crystal is very close to an ab initio minimum energy conformer of the isolated molecule. The extended C5 conformation of the main peptide chain is the same and there is slight discrepancy in the disposition of trityl phenyl rings. Their small dislocation creates the possibility of forming the entire network above of extensive, specific, weak intermolecular interactions; these constrain the molecule and permit it to retain the minimum energy C5 conformation of its main chain in the solid state. In contrast, in n-hexane solution, where such specific interactions cannot occur, only a small population of the molecules adopts the extended C5 conformation.     

Conformational aspects of glutathione tripeptide: electron density topological & natural bond orbital analyses

Glutathione tripeptide (γ-glutamyl-cysteinyl-glycine) is a flexible molecule and its conformational energy landscape is strongly influenced by forming intramolecular hydrogen bond, its charge and the environment. This study employs DFT-B3LYP method with the 6-31+G (d,p) basis set to carry out conformational analysis of neutral, zwitterionic, cationic, and anionic forms of glutathione. In analyzing the structural characteristics of these structures, intramolecular hydrogen bonds were identified and characterized in details by topological parameters such as electron density ρ(r) and Laplacian of electron density $ \nabla^{2} $ ρ(r) from Bader’s atom in molecules theory. Charge transfer energies based on natural bond orbital analysis are also considered to interpret these intramolecular hydrogen bonds. Our results show that these hydrogen bonds are partially electrostatic and partially covalent in nature, in which the covalent contribution increases as the stabilization energy of hydrogen bond increases. Furthermore, the back bone and side chain (Ramachandran map) orientations of various ionic forms of glutathione have been studied and conformation of each constitution of glutathione tripeptide (i.e., Glu, Cys, and Gly moieties) was determined. In most species side chain conformation were found to be hindered gauche–gauche orientation by intramolecular hydrogen bonds.     

Chiral Photoresponsive Tetrathiazoles That Provide Snapshots of Folding States

Herein, we designed chiral photoresponsive tetra(2‐phenylthiazole)s, which induce a diastereoselective 6π‐electrocyclization reaction in a helically folded structure to freeze the conformational interconversions. The folding conformation with one helical turn of tetra(2‐phenylthiazole)s was supported by multiple intramolecular noncovalent interactions including vicinal S???N interheteroatom interactions and CH–π and π–π stacking interactions between nonadjacent units, as found in X‐ray crystal structures as well as quantum chemical calculations. The introduction of a chiral group at both ends of tetra(2‐phenylthiazole) dictates the preferential folding into a one‐handed helix conformation by the simultaneous operation of S???O and multiple CH–π interactions that involve the chiral end groups. Since the tetra(2‐phenylthiazole)s possess two equivalent photoreactive 6π‐electron systems and the folded conformation is suitable for photoinduced electrocyclization reaction, they undergo a photocyclization reaction in a stereoselective manner to memorize the chirality of the helix in a resulting diastereomeric closed form.     

A hydrogen-bonding-modulated molecular rotor: environmental effect in the conformational stability of peptidomimetic macrocyclic cyclophanes

The conformational behavior of 16- to 18-membered ring peptidomimetic p-cyclophanes 1a,b-3a,b has been studied by NMR. The cycles bearing 16 and 17 atoms showed a dynamic process within the NMR time scale, produced by the rotation of the aromatic p-diphenylene moiety with respect to the macrocyclic main plane. The temperature dependence of 1H NMR spectra has been studied in order to get activation parameters of the energetic barrier for the process (VT-NMR and line shape analysis). The rate of the movement clearly depends on the macrocyclic ring size and the nature of the peptidomimetic side chain. Entropic and enthalpic contributions to the free energy of activation are discussed. The rotation of the aromatic ring is closely related to the intramolecular hydrogen bonding pattern, as suggested by temperature factors of NH chemical shifts (DeltadeltaNH/DeltaT) and molecular modeling. The interconnected roles of the solvation and the intramolecular H-bonds have been established by measurements (VT-NMR and DeltadeltaNH/DeltaT) in environments of different polarities and H-bonding abilities. We concluded that the conformational stability of the systems directly depends on the stability of the intramolecular H-bonding pattern. We finally showed how one of these peptidomimetics behaves as a methanol-dependent artificial molecular rotor. In this simple molecular device, the well-defined molecular rotation is tuned by the competition between intramolecular hydrogen bonds and interactions with the solvent.     

Design, synthesis, and computational studies of inhibitors of Bcl-XL

One of the primary objectives in the design of protein inhibitors is to shape the three-dimensional structures of small molecules to be complementary to the binding site of a target protein. In the course of our efforts to discover potent inhibitors of Bcl-2 family proteins, we found a unique folded conformation adopted by tethered aromatic groups in the ligand that significantly enhanced binding affinity to Bcl-XL. This finding led us to design compounds that were biased by nonbonding interactions present in a urea tether to adopt this bioactive, folded motif. To characterize the key interactions that induce the desired conformational bias, a series of substituted N,N'-diarylureas were prepared and analyzed using X-ray crystallography and quantum mechanical calculations. Stabilizing pi-stacking interactions and destabilizing steric interactions were predicted to work in concert in two of the substitution patterns to promote the bioactive conformation as a global energy minimum and result in a high target binding affinity. Conversely, intramolecular hydrogen bonding present in the third substitution motif promotes a less active, extended conformer as the energetically favored geometry. These findings were corroborated when the inhibition constant of binding to Bcl-XL was determined for fully elaborated analogues bearing these structural motifs. Finally, we obtained the NMR solution structure of the disubstituted N,N'-diarylurea bound to Bcl-XL demonstrating the folded conformation of the urea motif engaged in extensive pi-interactions with the protein.     

Molecular dynamics simulation study on controlling the adsorption behavior of polyethylene by fine tuning the surface nanodecoration of graphite

Molecular dynamics simulations are applied to study the adsorption of polyethylene with different chain lengths on patterned graphite surfaces that contain nanoscale protrusions. The influence of the nanostructure on the strong attractive interaction inherently in the hydrophobic polyethylene and hydrophobic graphite system is investigated by modifying the top surface area and the height and the shape of the protrusions. The results are analyzed in terms of the chain configuration, the adsorption energy, the global orientational order parameter, and the normalized surface-chain contacting pair number in the first adsorption layer. When the size of the protrusion increases, the adsorption energy, the order parameter, and the normalized surface-chain contacting pair number decrease at a fixed chain length. When the size of the protrusion is fixed, the average adsorption energy per monomer and the order parameter decrease with increasing chain length because of the stronger intramolecular interactions between the monomers. Changing the protrusion shape in a suitable way will effectively reduce the strong surface-chain interaction.