Conformational properties of a model alanyl dipeptide and of alanine-derived oligopeptides: Effects of solvation in water and in organic solvents—A combined SIBFA/continuum reaction field, ab initio self-consistent field, and density functional theory investigation

Role of intramolecular hydrogen bonds in promoting electron flow through amino acid and oligopeptide conjugates

Elucidating the factors that control charge transfer rates in relatively flexible conjugates is of importance for understanding energy flows in biology as well as assisting the design and construction of electronic devices. Here, we report ultrafast electron transfer (ET) and hole transfer (HT) between a corrole (Cor) donor linked to a perylene-diimide (PDI) acceptor by a tetrameric alanine (Ala)₄ Selective photoexcitation of the donor and acceptor triggers subpicosecond and picosecond ET and HT. Replacement of the (Ala)₄ linker with either a single alanine or phenylalanine does not substantially affect the ET and HT kinetics. We infer that electronic coupling in these reactions is not mediated by tetrapeptide backbone nor by direct donor-acceptor interactions. Employing a combination of NMR, circular dichroism, and computational studies, we show that intramolecular hydrogen bonding brings the donor and the acceptor into proximity in a "scorpion-shaped" molecular architecture, thereby accounting for the unusually high ET and HT rates. Photoinduced charge transfer relies on a (Cor)NH^…O=C-NH^…O=C(PDI) electronic-coupling pathway involving two pivotal hydrogen bonds and a central amide group as a mediator. Our work provides guidelines for construction of effective donor-acceptor assemblies linked by long flexible bridges as well as insights into structural motifs for mediating ET and HT in proteins.

Conformational Interconversions of Amino Acid Derivatives

Exhaustive conformational interconversions including transition structure analyses of N-acetyl-l-glycine-N-methylamide as well as its alanine, serine, and cysteine analogues have been investigated at the MP2/6-31G** level, yielding a total of 142 transition states. Improved estimates of relative energies were obtained by separately extrapolating the Hartree-Fock and MP2 energies to the basis set limit and adding the difference between CCSD(T) and MP2 results with the cc-pVDZ basis set to the extrapolated MP2 results. The performance of eight empirical force fields (AMBER94, AMBER14SB, MM2, MM3, MMFFs, CHARMM22_CMAP, OPLS_2005, and AMOEBAPRO13) in reproducing ab initio energies of transition states was tested. Our results indicate that commonly used class I force fields employing a fixed partial charge model for the electrostatic interaction provide mean errors in the ∼10 kJ/mol range for energies of conformational transition states for amino acid conformers. Modern reparametrized versions, such as CHARMM22_CMAP, and polarizable force fields, such as AMOEBAPRO13, have slightly lower mean errors, but maximal errors are still in the 35 kJ/mol range. There are differences between the force fields in their ability for reproducing conformational transitions classified according to backbone/side-chain or regions in the Ramachandran angles, but the data set is likely too small to draw any general conclusions. Errors in conformational interconversion barriers by ∼10 kJ/mol suggest that the commonly used force field may bias certain types of transitions by several orders of magnitude in rate and thus lead to incorrect dynamics in simulations. It is therefore suggested that information for conformational transition states should be included in parametrizations of new force fields.

Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy

We present an overview of the SIBFA polarizable molecular mechanics procedure, which is formulated and calibrated on the basis of quantum chemistry (QC). It embodies nonclassical effects such as electrostatic penetration, exchange-polarization, and charge transfer. We address the issues of anisotropy, nonadditivity, and transferability by performing parallel QC computations on multimolecular complexes. These encompass multiply H-bonded complexes and polycoordinated complexes of divalent cations. Recent applications to the docking of inhibitors to Zn-metalloproteins are presented next, namely metallo-beta-lactamase, phosphomannoisomerase, and the nucleocapsid of the HIV-1 retrovirus. Finally, toward third-generation intermolecular potentials based on density fitting, we present the development of a novel methodology, the Gaussian electrostatic model (GEM), which relies on ab initio-derived fragment electron densities to compute the components of the total interaction energy. As GEM offers the possibility of a continuous electrostatic model going from distributed multipoles to densities, it allows an inclusion of short-range quantum effects in the molecular mechanics energies. The perspectives of an integrated SIBFA/GEM/QM procedure are discussed.

Intramolecular interaction energies in model alanine and glycine tetrapeptides. Evaluation of anisotropy, polarization, and correlation effects. A parallel ab initio HF/MP2, DFT, and polarizable molecular mechanics study

An extension of the SIBFA polarizable molecular mechanics procedure to flexible oligopeptides is reported. The procedure is evaluated by computing the relative conformational energies, deltaE(conf), of the alanine tetrapeptide in 10 representative conformations, which were originally derived by Beachy et al. (J Am Chem Soc 1997, 119, 5908) to benchmark molecular mechanics procedures with respect to ab initio computations. In the present study, a particular emphasis is on the separable nature of the components of the energy and the particular impact of the polarization energy component on deltaE(conf). We perform comparisons with respect to single-point HF, DFT, LMP2, and MP2 computations done at the SIBFA-derived energy minima. Such comparisons are made first for the 10 conformers derived from phi/psi torsional angle energy-minimization (the rigid rotor approach), and, in a second step, after allowing additional relaxation of the C(alpha) centered valence angles. In both series of energy-minimization, the SIBFA deltaE(conf) compared best with the LMP2 results using the 6-311G** basis set, the rms being 1.3 kcal/mol. In the absence of the polarization component, the rms is 3.5 kcal/mol. In both series of minimizations, the magnitudes of deltaE(conf), computed as differences with respect to the most stable conformer taken as energy zero, decrease along the series: HF > DFT > LMP2 > SIBFA > MP2, indicative of increasing stabilization of the most highly folded conformers.

Binding of D- and L-captopril inhibitors to metallo-beta-lactamase studied by polarizable molecular mechanics and quantum mechanics

The bacterial Zn2+ metallo-beta-lactamase from B. fragilis is a zinc-enzyme with two potential metal ion binding sites. It cleaves the lactam ring of antibiotics, thus contributing to the acquired resistance of bacteria against antibiotics. The present study bears on the binuclear form of the enzyme. We compare several possible binding modes of captopril, a mercaptocarboxamide inhibitor of several zinc-metalloenzymes. Two diastereoisomers of captopril were considered, with either a D- or an L-proline residue. We have used the polarizable molecular mechanics procedure SIBFA (Sum of Interactions Between Fragments ab initio computed). Two beta-lactamase models were considered, encompassing 104 and 188 residues, respectively. The energy balances included the inter and intramolecular interaction energies as well as the contribution from solvation computed using a continuum reaction field procedure. The thiolate ion of the inhibitor is binding to both metal ions, expelling the bridging solvent molecule from the uncomplexed enzyme. Different competing binding modes of captopril were considered, either where the inhibitor binds in a monodentate mode to the zinc cations only with its thiolate ion, or in bidentate modes involving additional zinc binding by its carboxylate or ketone carbonyl groups. The additional coordination by the inhibitor's carboxylate or carbonyl group always occurs at the zinc ion, which is bound by a histidine, a cysteine, and an aspartate side chain. For both diastereomers, the energy balances favor monodentate binding of captopril via S-. The preference over bidentate binding is small. The interaction energies were recomputed in model sites restricted to captopril, the Zn2+ cations, and their coordinating end side chains from beta-lactamase (98 atoms). The interaction energies and their ranking among competing arrangements were consistent with those computed by ab initio HF and DFT procedures.

The crystal state conformation of Aib-rich segments of peptaibol antibiotics

Ac-(Aib-Ala)3-OH (a protected segment of the peptaibols gliodeliquescin and paracelsin), Z-Leu-Aib-Val-Aib-Gly-OtBu (a segment of [Leu]7-gliodeliquescin), Z-Val-Aib-Aib-Gln-OtBu (a common segment of alamethicin, paracelsin, and hypelcin), and Ac-Aib-Pro-(Aib-Ala)2-OMe and Z-Aib-Pro-(Aib-Ala)2-OMe, which represent differently N(alpha)-protected 1-6 segments of alamethicin and hypelcin, have been synthesized by solution methods. The crystal-state conformations of these five Aib-containing peptides have been determined by X-ray diffraction analysis. We have confirmed that the 3(10)-helical structure is preferentially adopted by Aib-rich short peptides. An experimentally unambiguous proof for the 3(10)-->alpha-helix conversion has been provided by the two differently N-blocked -Aib-Pro-(Aib-Ala)2-OMe hexapeptides. The beta-bend ribbon conformation, commonly observed in the (Aib-Pro)n sequential oligopeptides, is not found in the -Aib-Pro-Aib-Ala-Aib-Ala-sequence. As expected on the basis of the L-configuration of the C(alpha)-monoalkylated residues, a right-handed helix screw sense was found in all peptides investigated.