Computing magnetic anisotropy constants of single molecule magnets

We present here a theoretical approach to compute the molecular magnetic anisotropy parameters, D (M) and E (M) for single molecule magnets in any given spin eigenstate of exchange spin Hamiltonian. We first describe a hybrid constant M (S) valence bond (VB) technique of solving spin Hamiltonians employing full spatial and spin symmetry adaptation and we illustrate this technique by solving the exchange Hamiltonian of the Cu6Fe8 system. Treating the anisotropy Hamiltonian as perturbation, we compute the D (M)and E(M) values for various eigenstates of the exchange Hamiltonian. Since, the dipolar contribution to the magnetic anisotropy is negligibly small, we calculate the molecular anisotropy from the single-ion anisotropies of the metal centers. We have studied the variation of D (M) and E(M) by rotating the single-ion anisotropies in the case of Mn12Ac and Fe-8 SMMs in ground and few low-lying excited states of the exchange Hamiltonian. In both the systems, we find that the molecular anisotropy changes drastically when the single-ion anisotropies are rotated. While in Mn12Ac SMM D (M) values depend strongly on the spin of the eigenstate, it is almost independent of the spin of the eigenstate in Fe-8 SMM. We also find that the D (M)value is almost insensitive to the orientation of the anisotropy of the core Mn(IV) ions. The dependence of D (M) on the energy gap between the ground and the excited states in both the systems has also been studied by using different sets of exchange constants.

Single-molecule thermodynamics: the heat distribution function of a charged particle in a static magnetic field

Interest in the applicability of fluctuation theorems to the thermodynamics of single molecules in external potentials has recently led to calculations of the work and total entropy distributions of Brownian oscillators in static and time-dependent electromagnetic fields. These calculations, which are based on solutions to a Smoluchowski equation, are not easily extended to a consideration of the other thermodynamic quantity of interest in such systems-the heat exchanges of the particle alone-because of the nonlinear dependence of the heat on a particle's stochastic trajectory. In this paper, we show that a path integral approach provides an exact expression for the distribution of the heat fluctuations of a charged Brownian oscillator in a static magnetic field. This approach is an extension of a similar path integral approach applied earlier by our group to the calculation of the heat distribution function of a trapped Brownian particle, which was found, in the limit of long times, to be consistent with experimental data on the thermal interactions of single micron-sized colloids in a viscous solvent.

Single-Molecule DNA Analysis Reveals That Yeast Hop1 Protein Promotes DNA Folding and Synapsis: Implications for Condensation of Meiotic Chromosomes

During meiosis, long-range interaction between homologous chromosomes is thought to be crucial for homology recognition, exchange of DNA strands, and production of normal haploid gametes. However, little is known about the identity of the proteins involved and the actual molecular mechanism(s) by which chromosomes recognize and recombine with their appropriate homologous partners. Single-molecule analyses have the potential to provide insights into our understanding of this fascinating and long-standing question. Using atomic force microscopy and magnetic tweezers techniques, we discovered that Hop1 protein, a key structural component of Saccharomyces cerevisiae synaptonemal complex, exhibits the ability to bridge noncontiguous DNA segments into intramolecular stem-loop structures in which the DNA segments appear to be fully synapsed within the filamentous protein stems. Additional evidence suggests that Hop1 folds DNA into rigid protein DNA filaments and higher-order nucleoprotein structures. Importantly, Hop1 promotes robust intra- and intermolecular synapsis between double-stranded DNA molecules, suggesting that juxtaposition of DNA sequences may assist in strand exchange between homologues by recombination-associated proteins. Finally, the evidence from ensemble experiments is consistent with the notion that Hop1 causes rigidification of DNA molecules. These results provide the first direct evidence for long-range protein-mediated DNA DNA synapsis, independent of crossover recombination, which is presumed to occur during meiotic recombination.

Coherent (photon) vs incoherent (current) detection of multidimensional optical signals from single molecules in open junctions

The nonlinear optical response of a current-carrying single molecule coupled to two metal leads and driven by a sequence of impulsive optical pulses with controllable phases and time delays is calculated. Coherent (stimulated, heterodyne) detection of photons and incoherent detection of the optically induced current are compared. Using a diagrammatic Liouville space superoperator formalism, the signals are recast in terms of molecular correlation functions which are then expanded in the many-body molecular states. Two dimensional signals in benzene-1,4-dithiol molecule show cross peaks involving charged states. The correlation between optical and charge current signal is also observed. (C) 2015 AIP Publishing LLC.

Reversible Aptamer-Au Plasmon Rulers for Secreted Single Molecules

Plasmon rulers, consisting of pairs of gold nanoparticles, allow single-molecule analysis without photobleaching or blinking; however, current plasmon rulers are irreversible, restricting detection to only single events. Here, we present a reversible plasmon ruler, comprised of coupled gold nanoparticles linked by a single aptamer, capable of binding individual secreted molecules with high specificity. We show that the binding of target secreted molecules to the reversible plasmon ruler is characterized by single-molecule sensitivity, high specificity, and reversibility. Such reversible plasmon rulers should enable dynamic and adaptive live-cell measurement of secreted single molecules in their local microenvironment.

1-Anilino-8-naphthalene sulfonate (ANS) binding to proteins revisited: Investigation of dye binding to molten globule states by electrospray ionization mass spectrometry

The binding of 1-anilino-8-naphthalene-sulfonic acid to globular proteins at acidic pH has been investigated by electrospray ionization mass spectrometry (ESIMS). Mass spectra of apomyoglobin recorded in the pH range 2−7 establish that maximal ANS binding is observed at pH 4.0. As many as seven distinct species may be observed in the gas phase which correspond to protein molecules containing one to six molecules of bound ANS. At neutral pH only a single molecule of ANS is bound. In the case of cytochrome c, maximal binding is observed at pH 4.0, with five molecules being bound. Binding is suppressed at neutral pH. In both cases ESIMS demonstrates maximal ANS binding at pH values where the proteins have been reported to exist in molten globule states. ANS binding is not observed for lysozyme, which has a tightly folded structure over the entire pH range. Reduction of disulfide bonds in lysozyme leads to the detection of ANS-bound species at neutral pH. Binding is suppressed at low pH due to complete unfolding of the reduced protein. The results suggest that ESIMS may provide a convenient method of probing the stoichiometry and distribution of dye complexes with molten protein globules

Stretching single polymer chains of donor-acceptor foldamers: toward the quantitative study on the extent of folding

Single-molecule force spectroscopy has proven to be an efficient tool for the quantitative characterization of flexible foldamers on the single-molecule level in this study. The extent of folding has been estimated quantitatively for the first time to the best of our knowledge, which is crucial for a better understanding of the ``folding-process'' on single-molecule level. Therefore, this study may provide a guidance to regulate folding for realizing rational control over the functions of bulk materials.

N-Urethane-Protected Amino Alkyl Isothiocyanates: Synthesis, Isolation, Characterization, and Application to the Synthesis of Thioureidopeptides

Synthetically useful N-Fmoc amino-alkyl isothiocyanates have been described, starting from protected amino acids. These compounds have been synthesized in excellent yields by thiocarbonylation of the monoprotected 1,2-diamines with CS2/TEA/p-TsCl, isolated as stable solids, and completely characterized. The procedure has been extended to the synthesis of amino alkyl isothiocyanates from Boc- and Z-protected amino acids as well. The utility of these isothiocyanates for peptidomimetics synthesis has been demonstrated by employing them in the preparation of a series of dithioureidopeptide esters. Boc-Gly-OH and Boc-Phe-OH derived isothiocyanates 9a and 9c have been obtained as single crystals and their structures solved through X-ray diffraction. They belong to the orthorhombic crystal system, and have a single molecule in the asymmetric unit (Z′ = 1). 9a crystallizes in the centrosymmetric space group Pbca, while 9c crystallizes in the noncentrosymmetric space group P212121.

Nonspecifically bound proteins spin while diffusing along DNA

It is known that DNA-binding proteins can slide along the DNA helix while searching for specific binding sites, but their path of motion remains obscure. Do these proteins undergo simple one-dimensional (1D) translational diffusion, or do they rotate to maintain a specific orientation with respect to the DNA helix? We measured 1D diffusion constants as a function of protein size while maintaining the DNA-protein interface. Using bootstrap analysis of single-molecule diffusion data, we compared the results to theoretical predictions for pure translational motion and rotation-coupled sliding along the DNA. The data indicate that DNA-binding proteins undergo rotation-coupled sliding along the DNA helix and can be described by a model of diffusion along the DNA helix on a rugged free-energy landscape. A similar analysis including the 1D diffusion constants of eight proteins of varying size shows that rotation-coupled sliding is a general phenomenon. The average free-energy barrier for sliding along the DNA was 1.1 +/- 0.2 k(B)T. Such small barriers facilitate rapid search for binding sites.

Computational systems approach for drug target discovery

Importance of the field: The shift in focus from ligand based design approaches to target based discovery over the last two to three decades has been a major milestone in drug discovery research. Currently, it is witnessing another major paradigm shift by leaning towards the holistic systems based approaches rather the reductionist single molecule based methods. The effect of this new trend is likely to be felt strongly in terms of new strategies for therapeutic intervention, new targets individually and in combinations, and design of specific and safer drugs. Computational modeling and simulation form important constituents of new-age biology because they are essential to comprehend the large-scale data generated by high-throughput experiments and to generate hypotheses, which are typically iterated with experimental validation. Areas covered in this review: This review focuses on the repertoire of systems-level computational approaches currently available for target identification. The review starts with a discussion on levels of abstraction of biological systems and describes different modeling methodologies that are available for this purpose. The review then focuses on how such modeling and simulations can be applied for drug target discovery. Finally, it discusses methods for studying other important issues such as understanding targetability, identifying target combinations and predicting drug resistance, and considering them during the target identification stage itself. What the reader will gain: The reader will get an account of the various approaches for target discovery and the need for systems approaches, followed by an overview of the different modeling and simulation approaches that have been developed. An idea of the promise and limitations of the various approaches and perspectives for future development will also be obtained. Take home message: Systems thinking has now come of age enabling a `bird's eye view' of the biological systems under study, at the same time allowing us to `zoom in', where necessary, for a detailed description of individual components. A number of different methods available for computational modeling and simulation of biological systems can be used effectively for drug target discovery.

Experimental and Theoretical Charge Density Analysis of Polymorphic Structures: The Case of Coumarin 314 Dye

Experimental charge density distributions in two known conformational polymorphs (orange and yellow) of coumarin 314 dye are analyzed based on multipole modeling of X-ray diffraction data collected at 100 K. The experimental results are compared with the charge densities derived from multipole modeling of theoretical structure factors obtained from periodic quantum calculation with density functional theory (DFT) method and B3LYP/6-31G(d,p) level of theory. The presence of disorder at the carbonyl oxygen atom of ethoxycarbonyl group in the yellow form, which was not identified earlier, is addressed here. The investigationof intermolecular interactions, based on Hirshfeld surface analysis and topological properties via quantum theory of atoms in molecule and total electrostatic interaction energies, revealed significant differences between the polymorphs. The differences of electrostatic nature in these two polymorphic forms were unveiled via construction of three-dimensional deformation electrostatic potential maps plotted over the molecular surfaces. The lattice energies evaluated from ab initio calculations on the two polymorphic forms indicate that the yellow form is likely to be the most favorable thermodynamically. The dipole moments derived from experimental and theoretical charge densities and also from Lorentz tensor approach are compared with the single-molecule dipole moments. In each case, the differences of dipole moments between the polymorphs are identified.

Photophysics of conjugated polymers: interplay between Forster energy migration and defect concentration in shaping a photochemical funnel in PPV

Recent single molecule experiments have suggested the existence of a photochemical funnel in the photophysics of conjugated polymers, like poly[2-methoxy-5-(2'-ethylhexyl)oxy-1,4-phenylenevinylene] (MEH-PPV). The funnel is believed to be a consequence of the presence of conformational or chemical defects along the polymer chain and efficient non-radiative energy transfer among different chromophore segments. Here we address the effect of the excitation energy dynamics on the photophysics of PPV. The PPV chain is modeled as a polymer with the length distribution of chromophores given either by a Gaussian or by a Poisson distribution. We observe that the Poisson distribution of the segment lengths explains the photophysics of PPV better than the Gaussian distribution. A recently proposed version of an extended particle-in-a-box' model is used to calculate the exciton energies and the transition dipole moments of the chromophores, and a master equation to describe the excitation energy transfer among different chromophores. The rate of energy transfer is assumed to be given here, as a first approximation, by the well-known Forster expression. The observed excitation population dynamics confirms the photochemical funneling of excitation energy from shorter to longer chromophores of the polymer chain. The time scale of spectral shift and energy transfer for our model polymer, with realistic values of optical parameters, is in the range of 200-300 ps. We find that the excitation energy may not always migrate towards the longest chromophore segments in the polymer chain as the efficiency of energy transfer between chromophores depends on the separation distance between the two and their relative orientation.

Spatial confinement of exciton transfer and the role of conformational order in organic nanoparticles

Organic nanoparticles consisting of single conjugated polymer chains were investigated as a function of degree of conjugation by means of single-molecule spectroscopy. The degree of conjugation was synthetically controlled. For highly conjugated chains, singlet excitons are efficiently funneled over nanometer distances to a small number of sites. In contrast, chains with less conjugation and a high number of saturated bonds do not exhibit energy funneling due to a highly disordered conformation.

Ligand-modulated Parallel Mechanical Unfolding Pathways of Maltose-binding Proteins

Protein folding and unfolding are complex phenomena, and it is accepted that multidomain proteins generally follow multiple pathways. Maltose-binding protein (MBP) is a large (a two-domain, 370-amino acid residue) bacterial periplasmic protein involved in maltose uptake. Despite the large size, it has been shown to exhibit an apparent two-state equilibrium unfolding in bulk experiments. Single-molecule studies can uncover rare events that are masked by averaging in bulk studies. Here, we use single-molecule force spectroscopy to study the mechanical unfolding pathways of MBP and its precursor protein (preMBP) in the presence and absence of ligands. Our results show that MBP exhibits kinetic partitioning on mechanical stretching and unfolds via two parallel pathways: one of them involves a mechanically stable intermediate (path I) whereas the other is devoid of it (path II). The apoMBP unfolds via path I in 62% of the mechanical unfolding events, and the remaining 38% follow path II. In the case of maltose-bound MBP, the protein unfolds via the intermediate in 79% of the cases, the remaining 21% via path II. Similarly, on binding to maltotriose, a ligand whose binding strength with the polyprotein is similar to that of maltose, the occurrence of the intermediate is comparable (82% via path I) with that of maltose. The precursor protein preMBP also shows a similar behavior upon mechanical unfolding. The percentages of molecules unfolding via path I are 53% in the apo form and 68% and 72% upon binding to maltose and maltotriose, respectively, for preMBP. These observations demonstrate that ligand binding can modulate the mechanical unfolding pathways of proteins by a kinetic partitioning mechanism. This could be a general mechanism in the unfolding of other large two-domain ligand-binding proteins of the bacterial periplasmic space.

A Distributed Approach to Interference Cancellation

In this paper, we propose and analyze a novel idea of performing interference cancellation (IC) in a distributed/cooperative manner, with a motivation to provide multiuser detection (MUD) benefit to nodes that have only a single user detection capability. In the proposed distributed interference cancellation (DIC) scheme, during phase-1 of transmission, an MUD capable cooperating relay node estimates all the sender nodes' bits through multistage interference cancellation. These estimated bits are then sent by the relay node on orthogonal tones in phase-2 of transmission. The destination nodes receive these bit estimates and use them for interference estimation/cancellation, thus achieving IC benefit in a distributed manner. For this DIC scheme, we analytically derive an exact expression for the bit error rate (BER) in a basic five-node network (two source-destination node pairs and a cooperating relay node) on AWGN channels. Analytical BER results are shown to match with simulation results. For more general system scenarios, including more than two source-destination pairs and fading channels without and with space-time coding, we present simulation results to establish the potential for improved performance in the proposed distributed approach to interference cancellation. We also present a linear version of the proposed DIC.