Charge conductivity in peptides: dynamic simulations of a bifunctional model supporting experimental data

Stitching together a nm thick peptide-based semiconductor sheet using UV light

Langmuir monolayer allows for a two-dimensional nano-scale organization of amphiphilic molecules. We have adapted this technique to measure lateral and transverse conductivity in confined peptide nanosheets for the first time. We reported that two retro-isomers amphipathic peptides form stable monolayers showing a semiconductor-like behavior. Both peptides exhibit the same hydrophobicity and surface stability. They differ in the lateral conductivity and current-voltage due to the asymmetric peptide bond backbone orientation at the interface. Both peptides contain several tyrosines allowing the lateral crosslinking in neighboring molecules induced by UVB. UVB-light induces changes in the lateral conductivity and current-voltage behavior as well as monolayer heterogeneity monitored by Brewster Angle Microscopy. The semiconductor properties depend on the peptide bond backbone orientation and tyrosine crosslinking. Our results indicate that one may design extended nano-sheets with particular electric properties, reminiscent of semiconductors. We propose to exploit such properties for biosensing and neural interfaces.

Molecular Dynamics Simulations of Dna and Its Complexes

This article describes how classical molecular simulation methods are being used to gain a molecular-level understanding of the interaction mechanisms responsible for DNA–ligand recognition, and that govern the response of DNA to ligand binding. Case studies using a variety of different ligands—including small pharmaceutical drugs, proteins and lipids—are used to illustrate the power of modern molecular dynamics simulation methods for understanding how we may control the function and structure of DNA.

Computational studies on ground and excited state charge transfer properties of peptidomimetics

Chemical modifications at various peptide positions result in peptidomimetics with unique physical and chemical properties that can be used for a range of applications. Among many peptidomimetics, ureidopeptides are interesting due to their ability to act as donor-bridge-acceptor systems through which charge transfer occurs in one direction and can be triggered by an electrochemical pulse without perturbing the nuclear position. In this regard, some UP mimetics with different chromophoric units are studied in this work to understand their role using DFT based methods. Computational results and natural charge analysis provide evidence for the extensive contribution of the substituents to the excitation and hole migration dynamics. Further, the results show that the UP backbone preserves its uni-directional charge transfer phenomenon from the ureido to carboxylate terminal irrespective of the terminal groups and position. However, the substituent affects the excitation energies and the time scales of the hole migration. Among the substituents studied here, fluorine migrates to the hole within a shorter time scale while phenyl groups take longer.

Peptides as Bio-inspired Molecular Electronic Materials

Understanding the electronic properties of single peptides is not only of fundamental importance to biology, but it is also pivotal to the realization of bio-inspired molecular electronic materials. Natural proteins have evolved to promote electron transfer in many crucial biological processes. However, their complex conformational nature inhibits a thorough investigation, so in order to study electron transfer in proteins, simple peptide models containing redox active moieties present as ideal candidates. Here we highlight the importance of secondary structure characteristic to proteins/peptides, and its relevance to electron transfer. The proposed mechanisms responsible for such transfer are discussed, as are details of the electrochemical techniques used to investigate their electronic properties. Several factors that have been shown to influence electron transfer in peptides are also considered. Finally, a comprehensive experimental and theoretical study demonstrates that the electron transfer kinetics of peptides can be successfully fine tuned through manipulation of chemical composition and backbone rigidity. The methods used to characterize the conformation of all peptides synthesized throughout the study are outlined, along with the various approaches used to further constrain the peptides into their geometric conformations. The aforementioned sheds light on the potential of peptides to one day play an important role in the fledgling field of molecular electronics.

A controllable mechanistic transition of charge transfer in helical peptides: from hopping to superexchange

A controllable mechanistic transition of charge transfer in helical peptides is demonstrated as a direct result of side-bridge gating.

Computational Studies on Structural, Excitation, and Charge-Transfer Properties of Ureidopeptidomimetics

Peptides with ureido group enclosing backbones are considered peptidomimetics and are known for their higher stabilities, biocompatibilities, antibiotic, inhibitor, and charge-transduction activities. These peptidomimetics have some unique applications, which are quite different from those of natural peptides. Hence, it is imperative to appreciate their properties at a microscopic level. In this regard, this work outlines, in detail, the charge transfer (CT) properties, hole-migration dynamics, and electronic structures of various experimentally comprehended ureidopeptidomimetic models using density functional theory (DFT). Time-dependent DFT and complete active space self-consistent field computations on basic models provide the necessary evidence for the viability of CT from the end enfolding the ureido group to the other end with a carboxylate entity. This donor-to-acceptor CT has been reflected in excitation studies, in which the higher intensity band corresponds to CT from the π orbital of the ureido group to the π* orbital of the carboxylate entity. Further, hole-migration studies have shown that charge can evolve from the ureido end, whereas the hole generated at the carboxylate end does not migrate. However, hole migration has been reported to occur from both ends (amino and carboxylate ends) in glycine oligopeptides, and our studies show that the ability to transfer and migrate charge can be tuned by modifying the donor and acceptor functional groups in both the neutral and cationic charge states. We have analyzed the possibility of hole migration following ionization using DFT-based wave-packet propagation and found its occurrence on a ∼2-5 fs time scale, which reflects the charge-transduction ability of peptidomimetics.

A model for ultra-fast charge transport in membrane proteins

Isolated proteins have recently been observed to transport charge and reactivity over very long distances with extraordinary rates and near perfect efficiencies in spite of their site. This is not the case if the peptide is in water, where the efficiency of charge hopping to the next site is reduced to approximately 2%. Here, water is not an ideal solvent for charge transport. The issue at hand is how to explain such enormous charge transfer quenching in water compared to another typical medium, namely lipid. We performed molecular dynamics simulations to computationally substantiate the novel long-distance charge transfer yield of the polypeptides in lipids. This is characterized by the charge transfer persistent-distance decay constant and not by the rate, which is seldom, if ever, measured and hence not directly addressed here. This model can encompass an extremely wide range of yields over very long distances in peptides in various media. The calculations here demonstrate the good charge transport efficiency in lipids in contrast to the poor efficiency in water. The protein charge transport also exhibits a very strong anisotropic effect in lipids. The peptide secondary structure effect of charge transfer in membranes is analyzed in contrast to that in water. These results suggest that this model can be useful for the prediction of charge transfer efficiency in various environments of interest and indicate that the charge transfer is highly efficient in membrane proteins.

Influence of ionization on the conformational preferences of peptide models. Ramachandran surfaces of N-formyl-glycine amide and N-formyl-alanine amide radical cations

Ramachandran maps of neutral and ionized HCO-Gly-NH2 and HCO-Ala-NH2 peptide models have been built at the B3LYP/6-31++G(d,p) level of calculation. Direct optimizations using B3LYP and the recently developed MPWB1K functional have also been carried out, as well as single-point calculations at the CCSD(T) level of theory with the 6-311++G(2df,2p) basis set. Results indicate that for both peptide models ionization can cause drastic changes in the shape of the PES in such a way that highly disallowed regions in neutral PES become low-energy regions in the radical cation surface. The structures localized in such regions, epsilonL+* and epsilonD+* are highly stabilized due to the formation of 2-centre-3-electron interactions between the two carbonyl oxygens. Inclusion of solvent effects by the conductor-like polarizable continuum model (CPCM) shows that the solute-solvent interaction energy plays an important role in determining the stability order.

Probing Ultrafast Purely Electronic Charge Migration in Small Peptides

A pump–probe experiment that can examine a pure charge migration on a time scale short compared to the onset of nuclear motion is discussed. The mass spectrometric studies of Schlag et al. suggest that short peptide terminated by an aromatic amino acid are particularly suitable test compounds. The pump pulse needs to ionize the molecule on a time scale short compared to the period of the electronic motion, typically sub-fs. However, ionization occurs preferentially when the electrical field of the light is maximal so that the duration of the pulse envelope can be somewhat longer. Detection by photoelectron spectrometry of the peptide cation, to produce a dication, is shown to be able to probe the electronic rearrangement.

Charge transport in proteins probed by resonant photoemission

The degrees of charge localization in the highest occupied and lowest unoccupied molecular orbitals (HOMO and LUMO) of the bacterial surface layer protein of Bacillus sphaericus NCTC 9602 were studied by resonant photoemission. In agreement with a charge transport hopping mechanism that involves torsional motions of the peptide backbone, the lifetime of electrons excited into the LUMO was found to be approximately 100 fs.

Electron transfer through a self-assembled monolayer of a double-helix peptide with linking the terminals by ferrocene

A unique molecular structure, a double-helix peptide, was self-assembled on gold, and the electron transfer through the monolayer was studied. The double-helix peptide consists of two 9mer 3(10)-helical peptide chains having a disulfide group at each N terminal and being linked by a ferrocene dicarboxylic acid between the C terminals. Each helical peptide chain has three naphthyl groups in a linear arrangement along the helix. The monolayer properties and the electron transfer from the ferrocene unit to gold were studied with reference peptides with a similar double helix but without naphthyl groups, a single helix with a dicarboxylic ferrocene unit, and a single helix with a monocarboxylic ferrocene unit. It was demonstrated that the naphthyl groups on the side chains had no effect on electron transfer, and the electron-transfer rate in the double-helix monolayer was not promoted, despite the two electron pathways in the molecule. We propose that in the double-helix monolayer, molecular motions are suppressed, possibly by its rigid structure tethered by the two linkers on gold to cancel out acceleration effects of the 2-fold electron pathways and the ferrocene substitution number. The factors that affect the electron-transfer reaction across the helical peptide SAMs are discussed in depth.

Distal charge transport in peptides

Biological systems often transport charges and reactive processes over substantial distances. Traditional models of chemical kinetics generally do not describe such extreme distal processes. In this Review, an atomistic model for a distal transport of information, which was specifically developed for peptides, is considered. Chemical reactivity is taken as the result of distal effects based on two-step bifunctional kinetics involving unique, very rapid motional properties of peptides in the subpicosecond regime. The bifunctional model suggests highly efficient transport of charge and reactivity in an isolated peptide over a substantial distance; conversely, a very low efficiency in a water environment was found. The model suggests ultrafast transport of charge and reactivity over substantial molecular distances in a peptide environment. Many such domains can be active in a protein.

Charge Transfer Study through the Determination of the Ionization Energies of Tetrapeptides X3-Tyr, X = Gly, Ala, or Leu. Influence of the Inclusion of One Glycine in Alanine and Leucine Containing Peptides

The energies of the fundamental and several excited states of tetrapeptide radical cations were determined at the outer valence Green's function (OVGF) level, at three geometries corresponding to the lowest energy conformations: two for the neutral and one for the cation. The conformations were optimized at the density functional theory level within the B3LYP framework. It was found that, from a purely energetic point of view, a charge initially created on the tyrosine chromophore could migrate without any geometrical change and without further activation once the excited electronic state of the ionized chromophore was formed. This migration could reach the NH(2) terminus for the neutral conformations but should stop at the adjacent peptide link for the cation conformation. These results stress the probable influence of the electronic coupling between the states rather than the existence of a barrier on the charge pathway to explain the difference between the peptides in the charge-transfer process leading to the loss of an iminium [NH(2)=CHR](+) cation. The dissociation energy of the asymptote related to the formation of this NH(2) terminus iminium cation was calculated for few species and it appears that the excess energy available for dissociation is significant when starting from the lowest energy conformations of the neutral or the cation, provided that the charge transfer is effective. It was also found that the amino acids did not conserve their energetic properties and their zero order energy levels turned to a complete new energetic scheme corresponding to the conformation of the peptide.

Charge Transfer in Polypeptides: Effect of Secondary Structures on Charge-Transfer Integral and Site Energies

We have theoretically studied the charge transfer in glycine polypeptide using quantum mechanical models based on a tight-binding Hamiltonian approach. The charge-transfer integrals and site energies involved in the transport of positive charge through the peptide bond in glycine polypeptide have been calculated. The charge-transfer integrals and site energies have been calculated directly from the matrix elements of the Kohn-Sham Hamiltonian defined in terms of the molecular orbitals of the individual fragments of the glycine polypeptide. In addition to this, we have calculated the rate of charge transfer between a neighboring amino acid subgroup through the Marcus rate equation. These calculations have been performed for the different secondary structures of the glycine model peptide such as linear, alpha-helix, 3(10)-helix, and antiparallel beta-sheet by varying the dihedral angles omega, varphi, and psi along the Calpha-carbon of amino acid subgroup. Present theoretical results confirm that the charge transfer through the peptide bond is strongly affected by the conformations of the oligopeptide.

Demonstration of a new biosensing concept for immunodiagnostic applications based on change in surface conductance of antibodies after biomolecular interactions

We report an important observation that the surface conductivity of antibody layer immobilized on polylysine-coated glass substrate decreases upon the formation of complex with their specific antigens. This change in conductivity has been observed for both monoclonal and polyclonal antibodies. The conductance of monoclonal mouse IgG immobilized on polylysine-coated glass substrate changed from 1.02x10(-8) ohm(-1) to 1.41x10(-11) ohm(-1) at 10 V when complex is formed due to the specific biomolecular interactions with rabbit anti-mouse IgG F(ab')(2). Similar behavior was observed when the same set up was tested in two clinical assays: (1) anti-Leishmania antigen polyclonal antibodies taken from Kala Azar positive patient serum interacting with Leishmania promastigote antigen, and (2) anti-p21 polyclonal antibodies interacting with p21 antigen. The proposed concept can represent a new immunodiagnostic technique and may have wide ranging applications in biosensors and nanobiotechnology too.

An electronic time scale in chemistry

Ultrafast, subfemtosecond charge migration in small peptides is discussed on the basis of computational studies and compared with the selective bond dissociation after ionization as observed by Schlag and Weinkauf. The reported relaxation could be probed in real time if the removal of an electron could be achieved on the attosecond time scale. Then the mean field seen by an electron would be changing rapidly enough to initiate the migration. Tyrosine-terminated tetrapeptides have a particularly fast charge migration where in <1 fs the charge arrives at the other end. A femtosecond pulse can be used to observe the somewhat slower relaxation induced by correlation between electrons of different spins. A slower relaxation also is indicated when removing a deeper-lying valence electron. When a chromophoric amino acid is at one end of the peptide, the charge can migrate all along the peptide backbone up to the N end, but site-selective ionization is probably easier to detect for tryptophan than for tyrosine.

Electron transfer processes in DNA: mechanisms, biological relevance and applications in DNA analytics

In principle, DNA-mediated charge transfer processes can be categorized as oxidative hole transfer and reductive electron transfer. With respect to the routes of DNA damage most of the past research has been focused on the investigation of oxidative hole transfer or transport. On the other hand, the transport or transfer of excess electrons has a large potential for biomedical applications, mainly for DNA chip technology.

Femtosecond Dynamics after Ionization: 2-Phenylethyl-N,N-dimethylamine as a Model System for Nonresonant Downhill Charge Transfer in Peptides

The cation of 2-phenylethyl-N,N-dimethylamine (PENNA) offers two local sites for the charge: the amine group and 0.7 eV higher in energy the phenyl chromophore. In this paper, we investigate the dynamics of the charge transfer (CT) from the phenyl to the amine site. We present a femtosecond resonant two-color photoionization spectrum which shows that the femtosecond pump laser pulse is resonant in the phenyl chromophore. As shown previously with resonant wavelengths the aromatic phenyl chromophore can be then selectively ionized. Because the state "charge in the phenyl chromophore" is the first excited state in the PENNA cation, it can relax to the lower-energetic state "charge in the amine site". To follow this CT dynamics, femtosecond probe photoabsorption of green light (vis) is used. The vis light is absorbed by the charged phenyl chromophore, but not by the neutral phenyl and the neutral or cationic amine group. Thus, the absorption of vis photons of the probe laser pulse is switched off by the CT process. For detection of the resonant absorption of two or more vis photons in the cation the intensity of a fragmentation channel is monitored which opens only at high internal energy. The CT dynamics in PENNA cations has a time constant of 80 +/- 28 fs and is therefore not a purely electronic process. Because of its structural similarity to phenylalanine, PENNA is a model system for a downhill charge transfer in peptide cations.

Protein dynamics and electron transfer: electronic decoherence and non-Condon effects

We compute the autocorrelation function of the donor-acceptor tunneling matrix element <T(DA)(t)T(DA)(0)> for six Ru-azurin derivatives. Comparison of this decay time to the decay time of the time-dependent Franck-Condon factor {computed by Rossky and coworkers [Lockwood, D. M., Cheng, Y.-K. & Rossky, P. J. (2001) Chem. Phys. Lett. 345, 159-165]} reveals the extent to which non-Condon effects influence the electron-transfer rate. <T(DA)(t)T(DA)(0)> is studied as a function of donor-acceptor distance, tunneling pathway structure, tunneling energy, and temperature to explore the structural and dynamical origins of non-Condon effects. For azurin, the correlation function is remarkably insensitive to tunneling pathway structure. The decay time is only slightly shorter than it is for solvent-mediated electron transfer in small organic molecules and originates, largely, from fluctuations of valence angles rather than bond lengths.

Evidence Against the Hopping Mechanism as an Important Electron Transfer Pathway for Conformationally Constrained Oligopeptides

The rate constant of intramolecular electron transfer through oligopeptides based on the alpha-aminoisobutyric acid residue was determined as a function of the peptide length and found to depend weakly on the donor-acceptor separation. By measuring the electron-transfer activation energy and estimating the energy gap between donor and bridge, we were able to discard the electron hopping mechanism.

Understanding Electron Transfer across Negatively-Charged Aib Oligopeptides

The physicochemical effects modulating the conformational behavior and the rate of intramolecular dissociative electron transfer in phthalimide-Aibn-peroxide peptides (n = 0-3) have been studied by an integrated density functional/continuum solvent model. We found that three different orientations of the phthalimide ring are possible, labeled Phihel, PhiC7, and PhipII. In the condensed phase, they are very close in energy when the system is neutral and short. When the peptide chain length increases and the system is negatively charged, Phihel becomes instead the most stable conformer. Our calculations confirm that the 3(10)-helix is the most stable secondary structure for the peptide bridge. However, upon charge injection in the phthalimide end of the phthalimide-Aib3-peroxide, the peptide bridge can adopt an alpha-helix conformation as well. The study of the dependence of the frontier orbitals on the length and on the conformation of the peptide bridge (in agreement with experimental indications) suggests that for n = 3 the process could be influenced by a 3(10) --> alpha-helix conformational transition of the peptide chain.

Dynamic nature of the intramolecular electronic coupling mediated by a solvent molecule: a computational study

We present a combined Molecular Dynamics/Quantum Chemical study of the solvent-mediated electronic coupling between an electron donor and acceptor in a C-clamp molecule. We characterize the coupling fluctuations due to the solvent motion for different solvents (acetonitrile, benzene, 1,3-diisopropyl-benzene) for the charge separation and the charge recombination processes. The time scale for solvent-induced coupling fluctuation is approximately 0.1 ps. The effect of these fluctuations on the observed rate is discussed using a recently developed theoretical model. We show that, while the microscopic charge transfer process is very complicated and its computational modeling very subtle, the macroscopic phenomenology can be captured by the standard models. Analyzing the contribution to the coupling given by different solvent orbitals, we find that many solvent orbitals mediate the electron transfer and that paths through different solvent orbitals can interfere constructively or destructively. A relatively small subset of substrate-solvent configurations dominate contributions to solvent-mediated coupling. This subset of configurations is related to the electronic structure of the C-clamp molecule.

Formation of molecular radical cations of enkephalin derivatives via collision-induced dissociation of electrospray-generated copper (II) complex ions of amines and peptides

Fragmentation of some electrospray-generated complex ions, [63CuII(amine)M].2+, where M is an enkephalin derivative, produces the radical cation of the peptide, M.+. This ion has only been observed when M contains a tyrosyl or tryptophanyl residue plus a basic residue, typically arginyl or lysyl. A typical viable amine is diethylenetriamine. Collision-induced dissociation (CID) of the M.+ ion yields a prominent [M - 106].+ product ion for tyrosine-containing peptides, and a prominent [M - 129].+ ion for a tryptophan-containing peptide. These fragment ions are formed as a result of elimination of the tyrosyl and tryptophanyl side chains. Dissociation of these ions, in turn, produces second generation product ions, many of which are typically absent in the fragmentation of protonated peptide ions. Structures for some of these unusual ions are proposed.