A molecular mechanics and semiempirical conformational analysis of the herbicide diuron inhibitor of photosystem II

The potential energy surface (PES) for the herbicide diuron (DCMU), a photosystem II inhibitor, has been extensively investigated using the quantum-mechanical semiempirical molecular orbital methods AM1 and PM3 and molecular mechanics method. A detailed conformational search has been carried out which revealed the occurrence of four genuine minimum energy structures. The relative stability of the conformers and rotational barriers to conformational interconversion were evaluated using distinct theoretical approaches. The results showed that thetrans form of the diuron molecule is more stable than thecis form in all methods, and so it may possibly be the biologically active isomer.     

Cross examination of the conformational spaces of a set of peptide chains: Study of oligopeptidase action

A conformational search was carried out for five opioid peptide homologues and for angiotensin II. Density of states versus energy plots were obtained for each peptide, and the occurrence of common main-chain conformations was investigated by searching homologies between strings of four, five, and six contiguous main-chain amino acid residues rotamers. The results were compared to rates of hydrolysis by endooligopeptidase (EOP) 24.15, known for its specificity for substrate conformations. A catalytic assay of the hydrolysis of angiotensin II was also performed. The two best substrates of EOP 24.15 were found to share unique main-chain conformations and the two worst substrates of EOP 24.15 were found to be nonstructurally homologous to each other and the remaining peptide chains. The conformational search is compared to previous experimental and theoretical results. © 1996 John Wiley & Sons, Inc.     

Adsorption of homopolypeptides on gold investigated using atomistic molecular dynamics

We investigate the role of dynamics on adsorption of peptides to gold surfaces using all-atom molecular dynamics simulations in explicit solvent. We choose six homopolypeptides [Ala(10), Ser(10), Thr(10), Arg(10), Lys(10), and Gln(10)], for which experimental surface coverages are not correlated with amino acid level affinities for gold, with the idea that dynamic properties may also play a role. To assess dynamics we determine both conformational movement and flexibility of the peptide within a given conformation. Low conformational movement indicates stability of a given conformation and leads to less adsorption than homopolypeptides with faster conformational movement. Likewise, low flexibility within a given conformation also leads to less adsorption. Neither amino acid affinities nor dynamic considerations alone predict surface coverage; rather both quantities must be considered in peptide adsorption to gold surfaces.     

Modeling of peptides containing D-amino acids: implications on cyclization

Cyclic peptides are therapeutically attractive due to their high bioavailability, potential selectivity, and scaffold novelty. Furthermore, the presence of D-residues induces conformational preferences not followed by peptides consisting of naturally abundant L-residues. Therefore, comprehending how amino acids induce turns in peptides, subsequently facilitating cyclization, is significant in peptide design. Here, we performed 20-ns explicit-solvent molecular dynamics simulations for three diastereomeric peptides with stereochemistries: LLLLL, LLLDL, and LDLDL. Experimentally LLLLL and LDLDL readily cyclize, whereas LLLDL cyclizes in low yield. Simulations at 310 K produced conformations with inter-terminal hydrogen bonds that correlated qualitatively with the experimental cyclization trend. Energies obtained for representative structures from quantum chemical (B3LYP/PCM/cc-pVTZ//HF/6-31G*) calculations predicted pseudo-cyclic and extended conformations as the most stable for LLLLL and LLLDL, respectively, in agreement with the experimental data. In contrast, the most stable conformer predicted for peptide LDLDL was not a pseudo-cyclic structure. Moreover, D-residues preferred the experimentally less populated α_L rotamers even when simulations were performed at a higher temperature and with strategically selected starting conformations. Energies calculated with molecular mechanics were consistent only with peptide LLLLL. Thus, the conformational preferences obtained for the all L-amino acid peptide were in agreement with the experimental observations. Moreover, refinement of the force field is expected to provide far-reaching conformational sampling of peptides containing D-residues to further develop force field-based conformational-searching methods. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Anion influence on extraction and transport of amino acid methyl esters by functionalised calix[4]arene

The liquid–liquid extraction of a series of amino acid methyl esters has been carried out with functionalised calix[4]arene (5,11,17,23-tetrakis(N-methylpiperazino)-25,26,27,28-tetrahydroxycalix[4]arene) from an aqueous phase into a chloroform phase as ion pairs in the presence of picrate ion or tropaeolin 00 as counter ion in order to study the molecular recognition properties of this receptor. The active transport assisted by pH gradient of amino acids as ion pairs through liquid membrane employing the functionalised calix[4]arene as carrier has been investigated. The results showed that the receptor exhibits good extractability towards amino acids and it can also act as carrier through liquid membrane aiming to the separation of amino acids. It was highlighted that the anion nature used as counter ion, the structure of calix[4]arene, and the structure of amino acids are responsible for the experimental results obtained. High yields in both amino acids extraction and transport were obtained for picrate ion used as counter ion.     

Enantioselective binding of α-pinene and of some cyclohexanetriol derivatives by cyclodextrin hosts: A molecular modeling study

We have used molecular modeling to investigate the enantioselective separation of the monoterpene α-pinene on permethylated β-cyclodextrin and on α-cyclodextrin and the enantioselective separation of three cyclohexanetriol derivatives on permethylated β-cyclodextrin. Using the Consistent Valence Force Field (CVFF) from Insight/Discover, we have carried out systematic rigid-body docking grid searches on each of the optical antipodes of the organic guest molecules interacting with the cyclodextrins, followed by minimizations of the low-energy docked structures. A statistical mechanical analysis of the minimized energies yields data that agree in four out of five cases with the experimental elution order of enantiomers. The computed energies of the rigid-body docking before minimizations do not agree with the experimental results, suggesting that a conformational induced fit of the cyclodextrins upon binding of the organic guests may be involved in the mechanism of the chiral recognition. © 1996 by John Wiley & Sons, Inc.     

Synthesis and Conformational Analysis of Hybrid α/β‐Dipeptides Incorporating S‐Glycosyl‐β2,2‐Amino Acids

We synthesized and carried out the conformational analysis of several hybrid dipeptides consisting of an α‐amino acid attached to a quaternary glyco‐β‐amino acid. In particular, we combined a S‐glycosylated β^2,2‐amino acid and two different types of α‐amino acid, namely, aliphatic (alanine) and aromatic (phenylalanine and tryptophan) in the sequence of hybrid α/β‐dipeptides. The key step in the synthesis involved the ring‐opening reaction of a chiral cyclic sulfamidate, inserted in the peptidic sequence, with a sulfur‐containing nucleophile by using 1‐thio‐β‐D ‐glucopyranose derivatives. This reaction of glycosylation occurred with inversion of configuration at the quaternary center. The conformational behavior in aqueous solution of the peptide backbone and the glycosidic linkage for all synthesized hybrid glycopeptides was analyzed by using a protocol that combined NMR experiments and molecular dynamics with time‐averaged restraints (MD‐tar). Interestingly, the presence of the sulfur heteroatom at the quaternary center of the β‐amino acid induced θ torsional angles close to 180° (anti). Notably, this value changed to 60° (gauche) when the peptidic sequence displayed aromatic α‐amino acids due to the presence of CH–π interactions between the phenyl or indole ring and the methyl groups of the β‐amino acid unit.     

Use of constrained synthetic amino acids in beta-helix proteins for conformational control

A highly constrained amino acid has been introduced in the turn region of a beta-helix to increase the conformational stability of the native fold for nanotechnological purposes. The influence of this specific amino acid replacement in the final organization of beta-helix motifs has been evaluated by combining ab initio first-principles calculations on model systems and molecular dynamics simulations of entire peptide segments. The former methodology, which has been applied to a sequence containing three amino acids, has been used to develop adjusted templates. Calculations indicated that 1-amino-2,2-diphenylcyclopropanecarboxylic acid, a constrained cyclopropane analogue of phenylalanine, exhibits a strong tendency to form and promote folded conformations. On the other hand, molecular dynamics simulations are employed to probe the ability of such a synthetic amino acid to enhance the conformational stability of the beta-helix motif, which is the first requirement for further protein nanoengineering. A highly regular segment from a naturally occurring beta-helix protein was selected as a potential nanoconstruct module. Simulations of wild type and mutated segments revealed that the ability of the phenylalanine analogue to nucleate turn conformations enhances the conformational stability of the beta-helix motif in isolated peptide segments.     

Comparative analysis of the conformational profile of substance P using simulated annealing and molecular dynamics

The present study describes an extensive conformational search of substance P using two different computational methods. On the one hand, the peptide was studied using the iterative simulated annealing, and on the other, molecular dynamics simulations at 300 and 400 K. With the former method, the peptide was studied in vacuo with a dielectric constant of 80, whereas using the latter study the peptide was studied in a box of TIP3P water molecules. Analysis of the results obtained using both methodologies was carried out using an in-house methodology using a cluster analysis method based on information theory. Comparison of the two sampling methodologies and the different environment used in the calculations is also analyzed. Finally, the conformational motifs that are characteristic of substance P in a hydrophilic environment are presented and compared with the experimental results available in the literature.     

Repositioning of Ligands That Target the Spike Glycoprotein as Potential Drugs for SARS-CoV-2 in an In Silico Study

The worldwide health emergency of the SARS-CoV-2 pandemic and the absence of a specific treatment for this new coronavirus have led to the use of computational strategies (drug repositioning) to search for treatments. The aim of this work is to identify FDA (Food and Drug Administration)-approved drugs with the potential for binding to the spike structural glycoprotein at the hinge site, receptor binding motif (RBM), and fusion peptide (FP) using molecular docking simulations. Drugs that bind to amino acids are crucial for conformational changes, receptor recognition, and fusion of the viral membrane with the cell membrane. The results revealed some drugs that bind to hinge site amino acids (varenicline, or steroids such as betamethasone while other drugs bind to crucial amino acids in the RBM (naldemedine, atovaquone, cefotetan) or FP (azilsartan, maraviroc, and difluprednate); saquinavir binds both the RBM and the FP. Therefore, these drugs could inhibit spike glycoprotein and prevent viral entry as possible anti-COVID-19 drugs. Several drugs are in clinical studies; by focusing on other pharmacological agents (candesartan, atovaquone, losartan, maviroc and ritonavir) in this work we propose an additional target: the spike glycoprotein. These results can impact the proposed use of treatments that inhibit the first steps of the virus replication cycle.     

Modeling of folding and unfolding mechanisms in alanine-based alpha-helical polypeptides

alpha-Helix formation is known to be opposed by the entropy loss due to the folding and favored by the energy of molecular interactions. However, the underlying mechanism of these factors is still being discussed. Here we have used the experimental and calculation data for short alanine-based peptides embedded in water to model the mechanism of helix folding and unfolding and to calculate microscopically the free energy factors of alanine in the frame of helix coil conformational integrals. Classical helix-coil transition theories take into account the interactions in a peptide chain only if the i, i + 3 peptide bond participates in hydrogen bonding. But quantum mechanical calculations showed that interactions of the i, i + 2 peptide bond play an important role in helix folding too. We also included the short-range repulsive interactions due to molecular steric clashes and the end effects due to polar/hydrogen-bonding interactions at the N and C termini. The helix and coil regions of peptide conformational space were defined using an experimental steric criterion for hydrogen bonding. Arginine helix propensity was discussed and estimated. Monte Carlo numerical simulations of thermodynamics and kinetics for the 21 amino acid alpha-helical polypeptide Ac-A5(AAARA)3A-NMe were carried out and found to be in an agreement with the experimental results.     

Conformational analysis by NMR and molecular modelling of the 41–62 hydrophilic region of HIV-1 encoded virus protein U (Vpu). Effect of the phosphorylation on sites 52 and 56

The peptide of 22 amino acid residues, Vpu_P^41–62, phosphorylated at the two sites Ser⁵² and Ser⁵⁶ has been implicated in the degradation of CD4 receptor molecules, an important stage of the pathways to human immunodeficiency virus type 1 pathology (HIV-1). In order to assess the structural influence of phosphorylation, a conformational analysis by NMR and molecular simulation have been carried out for the phosphorylated Vpu_P^41–62, and non-phosphorylated Vpu^41–62 in both H₂O (at pH 3.5 and 7.2) and a 1:1 mixture of H₂O and trifluoroethanol. Analysis of the short-, medium-range NOE connectivities and of the secondary chemical shifts indicated that the peptide segment (42–49) shows a less well-defined helix propensity. The 50–62 segment forms a loop with the phosphate group pointing away, a short β-strand and a flexible extended ‘tail’ of residues 60–62. Differences in this molecular region 50-62 suggest that conformational changes of Vpu_P, play a potential role in Vpu_P-induced degradation of CD4 molecules.     

Conformationally Constrained Peptidomimetics as Inhibitors of the Protein Arginine Methyl Transferases

Protein arginine N‐methyl transferases (PRMTs) belong to a family of enzymes that modulate the epigenetic code through modifications of histones. In the present study, peptides emerging from a phage display screening were modified in the search for PRMT inhibitors through substitution with non‐proteinogenic amino acids, N‐alkylation of the peptide backbone, and incorporation of constrained dipeptide mimics. One of the modified peptides ( 23 ) showed an increased inhibitory activity towards several PRMTs in the low μm range and the conformational preference of this peptide was investigated and compared with the original hit using circular dichroism and NMR spectroscopy. Introducing two constrained tryptophan residue mimics (l ‐Aia) spaced by a single amino acid was found to induce a unique turn structure stabilized by a hydrogen bond and aromatic π‐stacking interaction between the two l ‐Aia residues.     

Restricted dead‐end elimination: Protein redesign with a bounded number of residue mutations

Dead‐end elimination (DEE) has emerged as a powerful structure‐based, conformational search technique enabling computational protein redesign. Given a protein with n mutable residues, the DEE criteria guide the search toward identifying the sequence of amino acids with the global minimum energy conformation (GMEC). This approach does not restrict the number of permitted mutations and allows the identified GMEC to differ from the original sequence in up to n residues. In practice, redesigns containing a large number of mutations are often problematic when taken into the wet‐lab for creation via site‐directed mutagenesis. The large number of point mutations required for the redesigns makes the process difficult, and increases the risk of major unpredicted and undesirable conformational changes. Preselecting a limited subset of mutable residues is not a satisfactory solution because it is unclear how to select this set before the search has been performed. Therefore, the ideal approach is what we define as the κ‐restricted redesign problem in which any κ of the n residues are allowed to mutate. We introduce restricted dead‐end elimination (rDEE) as a solution of choice to efficiently identify the GMEC of the restricted redesign (the κGMEC). Whereas existing approaches require n‐choose‐κ individual runs to identify the κGMEC, the rDEE criteria can perform the redesign in a single search. We derive a number of extensions to rDEE and present a restricted form of the A* conformation search. We also demonstrate a 10‐fold speed‐up of rDEE over traditional DEE approaches on three different experimental systems. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Validation of the GROMOS 54A7 Force Field Regarding Mixed α/β‐Peptide Molecules

A molecular‐dynamics (MD) simulation study of two heptapeptides containing α‐ and β‐amino acid residues is presented. According to NMR experiments, the two peptides differ in dominant fold when solvated in MeOH: peptide 3 adopts predominantly β‐hairpin‐like conformations, while peptide 8 adopts a 14/15‐helical fold. The MD simulations largely reproduce the experimental data. Application of NOE atom? atom distance restraining improves the agreement with experimental data, but reduces the conformational sampling. Peptide 3 shows a variety of conformations, while still agreeing with the NOE and ³J‐coupling data, whereas the conformational ensemble of peptide 8 is dominated by one helical conformation. The results confirm the suitability of the GROMOS 54A7 force field for simulation or structure refinement of mixed α/β‐peptides in MeOH.     

Helix Formation and Folding in γ‐Peptides and Their Vinylogues

A complete overview of all possible periodic structures with characteristic H‐bonding patterns is provided for oligomers composed of γ‐amino acids (γ‐peptides) and their vinylogues by a systematic conformational search on hexamer model compounds employing ab initio MO theory at various levels of approximation (HF/6‐31G*, DFT/B3LYP/6‐31G*, SCRF/HF/6‐31G*, PCM//HF/6‐31G*). A wide variety of structures with definite backbone conformations and H‐bonds formed in forward and backward directions along the sequence was found in this class of foldamers. All formally conceivable H‐bonded pseudocycles between 7‐ and 24‐membered rings are predicted in the periodic hexamer structures, which are mostly helices. The backbone elongation in comparison to α‐ and β‐peptides allows several possibilities to realize identical H‐bonding patterns. In good agreement with experimental data, helical structures with 14‐ and 9‐membered pseudocycles are most stable. It is shown that the introduction of an (E)‐double bond into the backbone of the γ‐amino acid constituents, which leads to vinylogous γ‐amino acids, supports the folding into helices with larger H‐bonded pseudocycles in the resulting vinylogous γ‐peptides. Due to the considerable potential for secondary‐structure formation, γ‐peptides and their vinylogues might be useful tools in peptide and protein design and even in material sciences.     

Rapid Acid Hydrolysis of Albumin in a Microwave Oven

Microwave-oven-assisted acid hydrolysis of albumin in 7.5 M HCl was carried out in only 12 min in a pressurized Teflon reactor using a power of 390 W. The effect of experimental conditions on the hydrolysis of albumin were evaluated and the results obtained were compared with those found by conventional hydrolysis, using in both cases the ninhydrin reaction for the spectrophotometric determination of free amino acids. A chromatographic study showed that microwave oven heating did not change the nature of the amino acids obtained from the hydrolysis of the protein.     

Conformational study of Met-enkephalin in its zwitterionic form

An extensive exploration of Met-enkephalin in its zwitterionic form has been carried out, in order to characterize the different low-energy conformational domains accessible to this pentapeptide. The study builds on previous studies carried out in our lab for the neutral molecule, which provided the initial geometries from which the conformational space of the charged molecule could be scanned. The initial conformations were subjected to a series of high- and lowtemperature molecular dynamics simulations. Snapshots along each trajectory were taken, minimized, and used as starting points in further MD trajectories until no lower-energy conformers could be characterized. The CHARMm force field was used throughout the study for this purpose. The same search strategy was used in these studies simulating two different environmental conditions, a distance-dependent dielectric ? = r and a high constant dielectric ? = 80. In the low dielectric environment, the formation of the salt bridge dominates the structure. In the high dielectric environment, the screening of the electrostatic interactions results in weaker intramolecular interactions. In both cases, the Gly²–Gly³ β-turn-type structures are preferred over the Gly³–Phe⁴ turns, in marked contrast to what is found for the neutral molecule. The lowest-energy structures from both environmental conditions were reoptimized in the presence of a cluster of explicit water molecules. Reoptimization of the structures considering explicit water structures did not result in significant conformational changes for the structures characterized with the ε = r or ε = 80 environments.     

Synthesis of Nature‐Inspired Medium‐Sized Fused Heterocycles from Amino Acids

Herein, we describe the synthesis of molecular scaffolds consisting of medium‐sized fused heterocycles using amino acids, which are some of the most useful building blocks used by nature as well as chemists to create structural diversity. The acyclic precursors were assembled by using traditional Merrifield solid‐phase peptide synthesis, and cyclization was carried out through acid‐mediated tandem endocyclic N‐acyliminium ion formation, followed by nucleophilic addition with internal nucleophiles. The synthesis of molecular scaffolds consisting of seven‐, eight‐, and nine‐membered rings proceeded with full stereocontrol of the newly generated stereogenic center in most cases.