Polyproline II helix conformation in a proline-rich environment: a theoretical study

Interest centers here on whether a polyproline II helix can propagate through adjacent non-proline residues, and on shedding light on recent experimental observations suggesting the presence of significant PP(II) structure in a short alanine-based peptide with no proline in the sequence. For this purpose, we explored the formation of polyproline II helices in proline-rich peptides with the sequences Ac-(Pro)(3)-X-(Pro)(3)-Gly-Tyr-NH(2), with X = Pro (PPP), Ala (PAP), Gln (PQP), Gly (PGP), and Val (PVP), and Ac-(Pro)(3)-Ala-Ala-(Pro)(3)-Gly-Tyr-NH(2) (PAAP), by using a theoretical approach that includes a solvent effect as well as cis <--> trans isomerization of the peptide groups and puckering conformations of the pyrrolidine ring of the proline residues. Since (13)C chemical shifts have proven to be useful for identifying secondary-structure preferences in proteins and peptides, and because values of the dihedral angles (phi,psi) are the main determinants of their magnitudes, we have, therefore, computed the Boltzmann-averaged (13)C chemical shifts for the guest residues in the PXP peptide (X = Pro, Ala, Gln, Gly, and Val) with a combination of approaches, involving molecular mechanics, statistical mechanics, and quantum mechanics. In addition, an improved procedure was used to carry out the conformational searches and to compute the solvent polarization effects faster and more accurately than in previous work. The current theoretical work and additional experimental evidence show that, in short proline-rich peptides, alanine decreases the polyproline II helix content. In particular, the theoretical evidence accumulated in this work calls into question the proposal that alanine has a strong preference to adopt conformations in the polyproline II region of the Ramachandran map.

Molecular modeling study of dihydrofolate reductase inhibitors. Molecular dynamics simulations, quantum mechanical calculations, and experimental corroboration

A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

A substituted sulfonamide and its Co (II), Cu (II), and Zn (II) complexes as potential antifungal agents

A sulfonamide 1-tosyl-1-H-benzo(d)imidazol-2-amine (TBZA) and three new complexes of Co(II), Cu(II), and Zn(II) have been synthesized. The compounds have been characterized by elemental analyses, FTIR, ¹H, and ¹³C-NMR spectroscopy. The structure of the TBZA, and its Co(II) and Cu(II) complexes, was determined by X-ray diffraction methods. TBZA and its Co(II) complex crystallize in the triclinic P-1 space group, while the Cu(II) complex crystallizes in the monoclinic P2₁/c space group. Antifungal activity was screened against eight pathogenic yeasts: Candida albicans (DMic 972576), Candida krusei (DMic 951705), Candida glabrata (DMic 982882), Candida tropicalis (DMic 982884), Candida dubliniensis (DMic 93695), Candida guilliermondii (DMic 021150), Cryptococcus neoformans (ATCC 24067), and Cryptococcus gattii (ATCC MYA-4561). Results on the inhibition of various human (h) CAs, hCA I, II, IV, VII, IX, and XII, and pathogenic beta and gamma CAs are also reported.

Performance of density functional models to reproduce observed (13)C(alpha) chemical shifts of proteins in solution

The purpose of this work is to test several density functional models (namely, OPBE, O3LYP, OPW91, BPW91, OB98, BPBE, B971, OLYP, PBE1PBE, and B3LYP) to determine their accuracy and speed for computing (13)C(alpha) chemical shifts in proteins. The test is applied to 10 NMR-derived conformations of the 76-residue alpha/beta protein ubiquitin (protein data bank id 1D3Z). With each functional, the (13)C(alpha) shielding was computed for 760 amino acid residues by using a combination of approaches that includes, but is not limited to, treating each amino acid X in the sequence as a terminally blocked tripeptide with the sequence Ac-GXG-NMe in the conformation of the regularized experimental protein structure. As computation of the (13)C(alpha) chemical shifts, not their shielding, is the main goal of this work, a computation of the (13)C(alpha) shielding of the reference, namely, tetramethylsilane, is investigated here and an effective and a computed tetramethylsilane shielding value for each of the functionals is provided. Despite observed small differences among all functionals tested, the results indicate that four of them, namely, OPBE, OPW91, OB98, and OLYP, provide the most accurate functionals with which to reproduce observed (13)C(alpha) chemical shifts of proteins in solution, and are among the faster ones. This study also provides evidence for the applicability of these functionals to proteins of any size or class, and for the validation of our previous results and conclusions, obtained from calculations with the slower B3LYP functional.

Unblocked statistical-coil tetrapeptides and pentapeptides in aqueous solution: a theoretical study

NMR studies of the molecular conformations of peptides and proteins rely on a comparison of the relevant spectral parameters with the corresponding values for so-called statistical-coil polypeptides. For this reason, it is necessary to characterize the experimental ensemble of states populated by statistical-coil peptides. Such a characterization, however, has proven to be both difficult and sensitive to changes in many environmental parameters such as solvent composition, temperature, pH, as well as the neighboring amino acids in the sequence. As a consequence, a series of significant discrepancies has been reported for some experimentally observed parameters, such as chemical shifts, or vicinal coupling constants, (3)J(NHalpha), whose values appear to be incompatible with a statistical-coil ensemble. In this work, we report the results of a molecular mechanics study of a series of unblocked tetra- and pentapeptides under different pH conditions. These calculations were carried out with explicit consideration of both the coupling between the process of proton binding/release and conformation adopted by the molecule at a given pH and the contribution of the conformational entropy to the total free energy. Good agreement was found between the calculated and experimentally determined values of the vicinal coupling constant, (3)J(NHalpha), the alpha-proton chemical shift, and the (13)C(alpha) chemical shift. All the evidence accumulated in these theoretical calculations helps to rationalize some of the unsettled anomalies observed experimentally, and to provide an understanding of the effect of pH and amino acid sequence on the conformational preferences of statistical-coil peptides.

New antitumoral acetogenin 'Guanacone type' derivatives: isolation and bioactivity. Molecular dynamics simulation of diacetyl-guanacone

We describe herein the isolation and semisynthesis of four acetogenin derivatives (1-4) as well as their ability to inhibit the mitochondrial respiratory chain and several tumor cell lines. In addition, four nanoseconds (ns) of MD simulation of compound 4, in a fully hydrated POPC bilayer, is reported.

Structure-antifeedant activity relationship of clerodane diterpenoids. Comparative study with withanolides and azadirachtin

A structure-antifeedant activity relationship (SAR) study of clerodane diterpenoids was carried out. Attention was focused on the feeding-deterrent activities exhibited toward Tenebrio molitor by clerodane diterpenoids and withanolides. Azadirachtin was chosen as a reference compound. SAR studies on the clerodane compounds indicate that the stereoelectronic factors are more important than the hydrophobic aspects as determinants of antifeedant activity. A furan ring in the side chain and a carbonyl alpha,beta-unsaturated (or spiro-epoxide) group appear to be indispensable for the biological response. A conformational study indicate that the optimum interatomic distance between these moieties is a range between 9.5 and 10.5 A. In addition, a similar stereoelectronic response was found among withanolides and azadirachtin. On the basis of these results it is reasonable to imagine a closely related chemical mechanism for these compounds.

Peptide model XXVIII: An exploratory ab initio and density functional study on the side-chain-backbone interaction in N-acetyl-L-cysteine- N-methylamide and N-formyl-L-cysteinamide in their γL-backbone conformations

A conformational and electronic study on the energetically preferred conformations (γL) of N- and C-protected L-cysteine (P-CONH-CH(CH2SH)-CONH-Q, where P and Q may be H or Me) was carried out. After restraining the backbone (BB) conformation to its global minimum (γL or C7eq), all nine possible side-chain (SC) conformations were subjected to geometry optimization at the HF/3–21G and the B3LYP/6–31G(d,p) levels of theory. Seven of the nine side-chain conformers were located on the potential-energy surface. All conformers were subjected to an AIM (atoms in molecules) analysis. This study indicates that three of the seven optimized conformers exhibited either or both SC [Formula: see text] BB- or BB [Formula: see text] SC-type intramolecular hydrogen bonding. Five conformers, however, had distances between a proton and a heteroatom that suggested hydrogen bonding.Key words: L-cysteine diamides, side-chain potential-energy surface, ab initio and DFT geometry optimization, AIM analysis, intramolecular hydrogen bonding.

Unblocked statistical-coil tetrapeptides in aqueous solution: quantum-chemical computation of the carbon-13 NMR chemical shifts

We recently reported a theoretical characterization of representative ensembles of statistical-coil conformations for tetrapeptides with unblocked termini in aqueous solution, at pH 7. The results showed good agreement between the computed Boltzmann-averaged and experimentally-determined values for both the vicinal coupling constants (3)J(NHalpha) and the alpha-proton chemical shifts. Here, we carry out a cluster analysis of the ensembles of conformations generated in that study, and use them to compute the Boltzmann-averaged values of the quantum-chemical (13)C chemical shifts for different amino acids in the unblocked tetrapeptides GGXA (where X stands for Phe, Arg, His, Glu, Ile, Lys, Gln, Tyr, Leu, Thr, Ala, Gly and Val). The values of the (13)C chemical shifts in these thirteen amino acids (for which experimental data are available) were computed by using Density Functional Theory with a 6-311+G(2d,p) basis set. Good agreement is found in terms of both the correlation coefficient (R) and standard deviations of the difference between the computed Bolztmann-averaged and the NMR-determined values for the (13)C chemical shifts. These results suggest that it may be possible to build a reliable theoretically-derived database of chemical shifts for statistical-coil residues. The results of the current study contribute to our understanding of the relations between chemical shifts, dihedral angles and vicinal coupling constants, (3)J(NHalpha). In addition, they can shed light as to how the statistical-coil conformation is related to the conformational preference of more structured states, such as the alpha-helical conformation.

An Exhaustive Conformational Analysis of N-Acetyl-l-cysteine-N-methylamide. Identification of the Complete Set of Interconversion Pathways on the ab Initio and DFT Potential Energy Hypersurface

The full conformational space of N-acetyl-l-cysteine-N-methylamide was explored by ab initio (RHF/ 6-31G(d)) and DFT (B3LYP/6-31G(d)) computations. Multidimensional conformational analysis predicts 81 structures in N-acetyl-l-cysteine-N-methylamide, but only 47 relaxed structures were previously determined at the RHF/3-21G level of theory. These structures were now optimized using RHF/6-31G(d) and B3LYP/6-31G(d) approaches. Seven conformational migrations were observed when recalculated at higher level of theory. Besides these major changes, only smaller conformational shifts were operative for the remaining stationary points. The exploration of the whole conformational space of N-acetyl-l-cysteine-N-methylamide, including the transition-state structures allowing the conformational interconversion among the low-energy forms, was analyzed in this study. Our results offer new insights into the influence of polar side chains on the conformational preferences of peptide structures.

Fast and accurate computation of the 13C chemical shifts for an alanine-rich peptide

The purpose of this work is, first, to present a fast and accurate technique to compute Boltzmann-averaged values of the quantum-chemical 13C chemical shifts for each amino acid in oligopeptides, demonstrated here by an application to the peptide Ac-XXAAAAAAAOO-NH2 (where X denotes diaminobutyric acid, A is alanine, and O is ornithine) [XAO] and, second, to discuss the capability of the 13Calpha and 13Cbeta chemical shifts to distinguish the PP(II) conformation from the alpha-helix and statistical-coil conformations. Use is made of a combination of approaches, summarized as follows: (1) derivation of an ensemble of conformations by using a molecular mechanics technique; (2) use of a clustering procedure to form families and build a reduced set of conformations consisting of the lowest-energy conformations of each family, and (3) computation of the 13C chemical shifts for the lowest-energy conformations of each family, using a quantum-chemical approach that treats a selected residue, or group of residues, with a 6-311+G(2d,p) locally-dense basis set, while the remaining residues in the sequence are treated with a 3-21G basis set. The whole procedure is quite accurate and speeds up the computation of the Boltzmann-averaged values of the 13C-chemical shifts by several orders of magnitude. The present application sheds some light on the conformational preference for alanine and non-alanine residues to occupy the PP(II) helical region of the Ramachandran map.

Position dependence of the 13C chemical shifts of alpha-helical model peptides. Fingerprint of the 20 naturally occurring amino acids

The position dependence of the (13)C chemical shifts was investigated at the density functional level for alpha-helical model peptides represented by the sequence Ac-(Ala)(i)-X-(Ala)(j)-NH(2), where X represents any of the 20 naturally occurring amino acids, with 0 < or = i < or = 8 and i + j = 8. Adoption of the locally dense basis approach for the quantum chemical calculations enabled us to reduce the length of the chemical-shift calculations while maintaining good accuracy of the results. For the 20 naturally occurring amino acids in alpha-helices, there is (1) significant variability of the computed (13)C shielding as a function of both the guest residue (X) and the position along the sequence; for example, at the N terminus, the (13)C(alpha) and (13)C(beta) shieldings exhibit a uniform pattern of variation with respect to both the central or the C-terminal positions; (2) good agreement between computed and observed (13)C(alpha) and (13)C(beta) chemical shifts in the interior of the helix, with correlation coefficients of 0.98 and 0.99, respectively; for (13)C(alpha) chemical shifts, computed in the middle of the helix, only five residues, namely Asn, Asp, Ser, Thr, and Leu, exhibit chemical shifts beyond the observed standard deviation; and (3) better agreement for four of these residues (Asn, Asp, Ser, and Thr) only for the computed values of the (13)C(alpha) chemical shifts at the N terminus. The results indicate that (13)C(beta), but not (13)C(beta), chemical shifts are sensitive enough to reflect the propensities of some amino acids for specific positions within an alpha-helix, relative to the N and C termini of peptides and proteins.

Evaluating the conformational space of the active site of D2 dopamine receptor. Scope and limitations of the standard docking methods

We report here for the first time the potential energy surfaces (PES) of phenyletilamine (PEA) and meta-tyramine (m-OH-PEA) at the D₂ dopamine receptor (D2DR) binding site. PESs not only allow us to observe all the critical points of the surface (minimums, maximums, and transition states), but also to note the ease or difficulty that each local minima have for their conformational inter-conversions and therefore know the conformational flexibility that these ligands have in their active sites. Taking advantage of possessing this valuable information, we analyze how accurate a standard docking study is in these cases. Our results indicate that although we have to be careful in how to carry out this type of study and to consider performing some extra-simulations, docking calculations can be satisfactory. In order to analyze in detail the different molecular interactions that are stabilizing the different ligand-receptor (L-R) complexes, we carried out quantum theory of atoms in molecules (QTAIM) computations and NMR shielding calculations. Although some of these techniques are a bit tedious and require more computational time, our results demonstrate the importance of performing computational simulations using different types of combined techniques (docking/MD/hybrid QM-MM/QTAIM and NMR shielding calculations) in order to obtain more accurate results. Our results allow us to understand in details the molecular interactions stabilizing and destabilizing the different L-R complexes reported here. Thus, the different activities observed for dopamine (DA), m-OH-PEA, and PEA can be clearly explained at molecular level.

N-Sulfonyl-1,2,3,4-tetrahydroisoquinoline Derivatives: Synthesis, Antimicrobial Evaluations, and Theoretical Insights

Microbial contamination remains a significant economic challenge in the food industry, emphasizing the need for innovative antimicrobial solutions. In this study, we synthesized N-sulfonyl-1,2,3,4-tetrahydroisoquinolines (NSTHIQ) derivatives using an environmentally friendly Preyssler heteropolyacid catalyst, obtaining moderate to high yields (35-91 %) under mild conditions. Two derivatives (5 and 6) exhibited significant antifungal properties against various fungal species, including Aspergillus spp, Penicillium spp, and Botrytis cinerea. ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) analysis revealed the absence of hepatic toxicity in all compounds, making derivatives 2, 3, 4, and 5 potential candidates for further development. However, derivatives 6 and 7 exhibited immunotoxicity. In support of our experimental findings, reactivity indices were computed using Density Functional Theory principles, deriving valuable insights into the chemical properties of these derivatives. This study underscores the potential of NSTHIQ compounds as potent antifungal agents, coupled with the importance of employing environmentally friendly catalysts in drug discovery.