Can a synthetic thread act as an electrochemically switchable molecular device?

Here we report that at room temperature in acetonitrile after the reduction of the naphthalimide-site, a synthetic molecular thread undergoes a complete conformational change which makes possible an efficient conversion of chemical energy into mechanical work; such results point out the ability of the thread to act as a molecular device under electrochemical control.     

Dynamic Interconversion between Solomon Link and Trapezoidal Metallacycle Ensembles Accompanying Conformational Change of the Linker

A novel template-free Cp*Rh-based molecular Solomon link has been established through selection of the flexible ligand L as a linker and the half-sandwich rhodium(III) dinuclear fragment B1 as a rigid capping unit. Furthermore, we demonstrate that the self-assembly of the Solomon link based on the flexible ligand is both solvent- and concentration-dependent: the Solomon link is formed in concentrated methanolic solutions, whereas formation of a dinuclear trapezoidal rectangle is favored at low concentrations or in acetonitrile or DMF solutions. Remarkably, alteration of the solvent or concentration can promote a unique and dynamic interconversion between the two molecular species, accompanying conformational change of the ligand. The synthetic outcomes are supported by single-crystal X-ray diffraction analysis.     

ZIBgridfree: efficient conformational analysis by partition-of-unity coupling

Obtaining a sufficient sampling of conformational space is a common problem in molecular simulation. We present the implementation of an umbrella-like adaptive sampling approach based on function-based meshless discretization of conformational space that is compatible with state of the art molecular dynamics code and that integrates an eigenvector-based clustering approach for conformational analysis and the computation of inter-conformational transition rates. The approach is applied to three example systems, namely $n$ -pentane, alanine dipeptide, and a small synthetic host-guest system, the latter two including explicitly modeled solvent.     

Solvent dependence on conformational transition, dipole moment, and molecular geometry of 1,2-dichloroethane: insight from Car-Parrinello molecular dynamics calculations

We have investigated the molecular geometry and dipole moment distribution for the major conformational states of 1,2-dichloroethane (DCE) in three different solvents under ambient conditions using the Car-Parrinello mixed quantum mechanics/molecular mechanics method. The solvents studied were water, DCE, and chloroform. Within the time scale investigated, we find a conformational equilibrium existing between the gauche and trans forms of DCE in all three solvents. In the chloroform solvent, the conformational transition occurs more frequently than in water solvent and in liquid DCE (i.e., DCE solute in DCE solvent). The population of gauche conformer is more in the case of water solvent, while the trans conformer dominates in chloroform solvent. We report a bimodal nature of the dipole moment distribution for DCE in all three solute-solvents studied, where the peaks are attributed to the trans and gauche conformational states. The dipole moments of both of the conformational states increase with increasing polarity of the solvent. Also, with an increase in solvent polarity, an increase in the C-Cl bond length and magnitude of atomic charges in DCE has been observed. The increase in atomic charges of DCE is almost twice when the solvent is changed from chloroform to water.     

Impact of the solvent on the conformational isomerism of calix[4]arenes: a study based on continuum solvation models

The influence of solvation on the conformational isomerism of calix[4]arene and p-tert-butylcalix[4]arene has been investigated by using the continuum model reported by Miertus, Scrocco, and Tomasi (MST). The quantum mechanical (QM) and semiclassical (SC) formalisms of the MST model have been considered for two different solvents (chloroform and water). The suitability of the QM-MST and SC-MST methods has been examined by comparison with previous results derived from classical molecular dynamics (MD) simulations with explicit solvent molecules. The application of the continuum model to the solute configurations generated by using in vacuo classical MD simulations provides a fast strategy to evaluate the effects of the solvent on the conformational preferences of calixarenes. These encouraging results allow us to propose the use of continuum models to solutes with complex molecular structures, which are traditionally studied by MD simulations.     

Infrared study of intercomponent interactions in a switchable hydrogen-bonded rotaxane

The macrocycle in rotaxane 1 is preferentially hydrogen bonded to the succinamide station in the neutral form, but can be moved to the naphthalimide station by one-electron reduction of the latter. The hydrogen bonding between the amide NH groups of the macrocycle and the C double bond O groups in the binding stations in the thread was studied with IR spectroscopy in different solvents in both states. In addition, the solvent effect on the vibrational frequencies was analyzed; a correlation with the solvent acceptor number (AN) was observed. The conformational switching upon reduction could be detected by monitoring the hydrogen-bond-induced shifts of the nu(CO) frequencies of the C double bond O groups of the succinamide and the reduced naphthalimide stations. The macrocycle was found to shield the encapsulated station from the solvent: wavenumbers of nu(CO) bands of the C double bond O groups residing inside the macrocycle cavity remain unaffected by the solvent polarity.     

Ab initio studies of solvent effect on molecular conformation and vibrational spectra of diacetamide

The conformational equilibria and vibrational spectra of diacetamide have been investigated by ab initio molecular orbital studies using the basis sets 6-31g(d,p) and 6-31++g(d,p) at Hartree-Fock and MP2 levels. The vibrational spectra of diacetamide have been satisfactorily interpreted taking into consideration the agreement between the calculated harmonic vibrational frequencies, infrared and Raman band intensities and shifts in deuterated molecules with those observed. The solvent effects were investigated by the self-consistent reaction field (SCRF) theory. The effect of solvent on the conformational equilibria and vibrational spectra is discussed. The calculated changes in the geometry and vibrational spectra on going from the gas phase to the solvent medium are in accord with the increasing weight of the dipolar resonance structure of the amide group in more polar solvents.     

A study of hydration effects on the conformational aspects of gaba mediators

The conformation of numerous chemical compounds is strongly influenced by solvents. Knowledge of their structure in solution is necessary, especially for a discussion of the biological and pharmacological activity of the molecules. The neurotransmitters and their agonists and antagonists are known to be flexible molecules that interact with quite distinct receptors. The conformational properties of several GABA (γ-aminobutyric acid) mediators have been studied by the MNDO technique. The optimized geometries of the molecules have been used to study the solvent effects on their conformational properties considering the supermolecule approach for their first hydration shell. A conjectural pharmacophoric pattern for several GABA inhibitors has already been suggested from the molecular electrostatic potentials (MEP ) of several molecules using a localized bond orbital technique. In the present work, MEP calculations have been carried out considering a solvent effect on the MNDO optimized geometries to investigate any deviation from the earlier results.     

Conformational analysis of siloxane-based enzyme-mimic precursors[1]

The conformational flexibility of six hybrid organodisiloxane oligomers were studied using the Low Mode-Monte Carlo conformational search method with the MM2* force field and the Generalized Born/Surface Area continuum solvent model for water. These systems have enzyme-like properties as synthetic acyltransferases and contain aminopyridine groups in various states of protonation. An ensemble of low energy structures was generated and used to investigate the dependence of molecular shape and flexibility on protonation, which plays an important role in catalyst solubility and self-association. The results as measured by the number of unique conformations, end-to-end or longest intramolecular distance and radius of gyration of the conformational point cloud indicate that the number of protonated pyridines plays a significant role in the overall molecular shape. A similar study was also carried out on various POSS-substitutive organodisiloxane oligomers.     

Conformational behavior of comb-like polyelectrolytes in selective solvent: computer simulation study

In this work, we present molecular dynamics simulations of comb-like polyelectrolytes in selective solvent. The studied polymers have a neutral backbone and polyelectrolyte side chains. The solvent is poor for the backbone and the theta solvent for the side chains. The polymers are modeled on a coarse-grained level with implicit solvent. The simulations show that the comb-like polyelectrolytes tend to form intramolecular self-organized structures of the pearl necklace type. This type of conformational behavior has been predicted by Borisov and Zhulina (Borisov, O. V.; Zhulina, E. B. Macromolecules 2005, 38, 2506) for neutral comb-like copolymers in selective solvent. The present study shows that comb-like polyelectrolytes in selective solvent exhibit the same type of behavior; however, it can be controlled by one additional parameter, the degree of dissociation of the grafts. The local conformational characteristics are studied using the ensemble-averaged bond angle cosines as functions of monomer position in the chain, which reveal structural details invisible by other means.     

Conformational equilibria of 7-benzyl-2-iodo-9-oxa-7-azabicyclo[4.3.0]nonan-8-one in solution. Correlations between conformational distribution and solvent solvatochromic parameters

A quantitative study of the variation of the conformational equilibria of 7-benzyl-2-iodo-9-oxa-7-azabicyclo[4.3.0]nonan-8-one 1 in 10 solvents has been carried out. The experimental composition in each solvent has been obtained from experimental NMR vicinal H-H coupling constants together with molecular modeling. The solvent properties, particularly polarity and hydrogen bonding ability, were described according to Kamlet and Taft using experimental parameters. Very good linear relationships were obtained between the equilibrium constants of each single conformational equilibrium and the polarity and hydrogen bonding parameters of the solvent. These linear relationships allow an accurate prediction of the conformational composition in any solvent as well as a thorough understanding of the influence of each separate parameter on the conformational equilibrium composition.     

Analyzing the robustness of the MM/PBSA free energy calculation method: application to DNA conformational transitions

The ability to predict and characterize free energy differences associated with conformational equilibria or the binding of biomolecules is vital to understanding the molecular basis of many important biological functions. As biological studies focus on larger molecular complexes and properties of the genome, proteome, and interactome, the development and characterization of efficient methods for calculating free energy becomes increasingly essential. The aim of this study is to examine the robustness of the end-point free energy method termed the molecular mechanics Poisson-Boltzmann solvent accessible surface area (MM/PBSA) method. Specifically, applications of MM/PBSA to the conformational equilibria of nucleic acid (NA) systems are explored. This is achieved by comparing A to B form DNA conformational free energy differences calculated using MM/PBSA with corresponding free energy differences determined with a more rigorous and time-consuming umbrella sampling algorithm. In addition, the robustness of NA MM/PBSA calculations is also evaluated in terms of the sensitivity towards the choice of force field and the choice of solvent model used during conformational sampling. MM/PBSA calculations of the free energy difference between A-form and B-form DNA are shown to be in very close agreement with the PMF result determined using an umbrella sampling approach. Further, it is found that the MM/PBSA conformational free energy differences were also in agreement using either the CHARMM or AMBER force field. The influence of ionic strength on conformational stability was particularly insensitive to the choice of force field. Finally, it is also shown that the use of a generalized Born implicit solvent during conformational sampling results in free energy estimates that deviate slightly from those obtained using explicitly solvated MD simulations in these NA systems.     

Cages, baskets, ladders, and tubes: conformational studies of polyhedral oligomeric silsesquioxanes

The conformational flexibility of a series of cage, basket, ladder, and tube polyhedral oligomeric silsesquioxanes (POSS) has been examined using the Low Mode:Monte Carlo conformational search method in conjunction with the MM3/GBSA(CHCl3) surface. An ensemble of low energy structures was generated and used to explore the molecular shape and flexibility of each system. The results indicate that, except for the ladder molecule, the incompletely condensed systems that are studied are relatively rigid. Even in cases where the molecule is able to adopt numerous low energy conformations, the overall shape remains cage-like and the conformations differ only by small angles or substituent orientations. The ladder molecule is the most flexible and this ensemble clusters into two families: one that is cage-like and the other that is more open and ladder-like. The conformational flexibilities in the gas and solvent phases, as approximated using the GBSA continuum solvent model, are very similar.     

Calculating nuclear magnetic resonance chemical shifts in solvated systems

The nuclear magnetic resonance (NMR) chemical shift is extremely sensitive to molecular geometry, hydrogen bonding, solvent, temperature, pH, and concentration. Calculated magnetic shielding constants, converted to chemical shifts, can be valuable aids in NMR peak assignment and can also give detailed information about molecular geometry and intermolecular effects. Calculating chemical shifts in solution is complicated by the need to include solvent effects and conformational averaging. Here, we review the current state of NMR chemical shift calculations in solution, beginning with an introduction to the theory of calculating magnetic shielding in general, then covering methods for inclusion of solvent effects and conformational averaging, and finally discussing examples of applications using calculated chemical shifts to gain detailed structural information.     

Molecular recognition: improved binding of biotin derivatives with synthetic receptors

NMR titrations and Monte Carlo conformational searches have been used to study the molecular recognition features of five urea derivatives with two synthetic hosts. We have improved the binding constant (Kb) values for all the studied guests and measured the largest binding constant of a complex involving a biotin derivative (biotin methyl ester) bound to a synthetic host by means of several interaction points and not only through the urea moiety.     

Investigations into conformational transitions and solvation structure of a 7-piperidino-5,9-methanobenzo[8] annulene in water

Solvation shell structure of a 7-piperidino-5,9-methanobenzo[8] annulene (PMA) in water has been investigated in ambient conditions using both molecular dynamics (MD) and Car-Parrinello molecular dynamics (CPMD) calculations. From the MD calculations, we find that this molecule exists in three major conformational states out of which two are in twist-boat forms and one in chair form. Due to the limited time scale accessible in CPMD simulations, we have studied all the three conformational states separately using CPMD. The molecular geometry, electronic charge distribution and solvation structure for all three forms are investigated. The stability order of the chair and twist-boat conformations in water solvent has been reversed when compared to the gaseous phase results and in the case of polar aprotic solvents (J. Org. Chem., 1999, 61, 5979). From the radial distribution function, we find that the solvent density around the chair form is significantly lower, which has to be directly related to the smaller solvent accessible area for this conformation and this is in complete agreement with earlier reports. Among the findings are that the solvation shell structure around the nitrogen atom in the chair form of PMA is considerably different from the open conformational forms or the twist-boat forms. The dipole moment for the closed form is found to be significantly larger when compared to the twist-boat forms.     

Contraction process of an electroactive actuator based on a one microsecond atomistic molecular dynamics simulation

The contraction process of an electroactive actuator constituted by calix[4]arene units and quaterthiophene segments has been investigated at the microscopic level by using atomistic molecular dynamics simulations in dichloromethane solution using explicit solvent molecules. Results derived from a 1 mus trajectory of the oxidized and deprotonated actuator indicate that the contraction occurs through a non-concerted mechanism in which each actuating units present in the system behave independently. The efficiency of the contraction process can be reduced by the presence of secondary conformational transitions in the calix[4]arene scaffolds. Accordingly, the drastic reduction of the molecular length expected during the contraction process can be limited by such transitions, which involve the rotational isomerism of a phenolate ring. However, such type of conformational transitions does not compromise the actuator power due to its intrinsic capacity to adopt compact molecular arrangements. On the other hand, the rate of the contraction process is influenced by the presence of solvent molecules, which have been found to reduce it by a factor of about 1000.     

De novo prediction of the structures of M. tuberculosis membrane proteins

The structures of four integral membrane proteins from the Mycobacterium tuberculosis (TB) gene, Rv2433c, Rv1861, Rv1616, and Rv3069, have been de novo predicted by combining a generalized Born implicit solvent/membrane model with replica exchange molecular dynamics simulations to sample the conformational space of each protein.     

Mosaic energy landscapes of liquids and the control of protein conformational dynamics by glass-forming solvents

Using recent advances in the Random First-Order Transition (RFOT) Theory of glass-forming liquids, we explain how the molecular motions of a glass-forming solvent distort the protein's boundary and slave some of the protein's conformational motions. Both the length and time scales of the solvent imposed constraints are provided by the RFOT theory. Comparison of the protein relaxation rate to that of the solvent provides an explicit lower bound on the size of the conformational space explored by the protein relaxation. Experimental measurements of slaving of myoglobin motions indicate that a major fraction of functionally important motions have significant entropic barriers.     

IR spectroscopy on jet-cooled isolated two-station rotaxanes

High-resolution IR spectroscopy has been employed to study isolated, switchable [2]rotaxanes. IR absorption spectra of two-station rotaxanes, their separate thread, and macrocycle components, as well as those of the individual stations incorporated into the thread, have been measured in the 1800-1000 cm(-1) region. These spectra have been fully analyzed, aided by quantum chemical predictions of the IR spectra. From these analyses, a comprehensive picture emerges of the conformational structure and binding interactions between the mechanically interlocked components of the rotaxane.