Exploring the Energy Landscapes of Cyclic Tetrapeptides with Discrete Path Sampling

Exploiting Sequence-Dependent Rotamer Information in Global Optimization of Proteins

Rotamers, namely amino acid side chain conformations common to many different peptides, can be compiled into libraries. These rotamer libraries are used in protein modeling, where the limited conformational space occupied by amino acid side chains is exploited. Here, we construct a sequence-dependent rotamer library from simulations of all possible tripeptides, which provides rotameric states dependent on adjacent amino acids. We observe significant sensitivity of rotamer populations to sequence and find that the library is successful in locating side chain conformations present in crystal structures. The library is designed for applications with basin-hopping global optimization, where we use it to propose moves in conformational space. The addition of rotamer moves significantly increases the efficiency of protein structure prediction within this framework, and we determine parameters to optimize efficiency.

Crossing of the Cystic Barriers of Toxoplasma gondii by the Fluorescent Coumarin Tetra-Cyclopeptide

FR235222 is a natural tetra-cyclopeptide with a strong inhibition effect on histone deacetylases, effective on mammalian cells as well as on intracellular apicomplexan parasites, such as Toxoplasma gondii, in the tachyzoite and bradyzoite stages. This molecule is characterized by two parts: the zinc-binding group, responsible for the binding to the histone deacetylase, and the cyclic tetrapeptide moiety, which plays a crucial role in cell permeability. Recently, we have shown that the cyclic tetrapeptide coupled with a fluorescent diethyl-amino-coumarin was able to maintain properties of cellular penetration on human cells. Here, we show that this property can be extended to the crossing of the Toxoplasma gondii cystic cell wall and the cell membrane of the parasite in its bradyzoite form, while maintaining a high efficacy as a histone deacetylase inhibitor. The investigation by molecular modeling allows a better understanding of the penetration mechanism.

Side-Chain Polarity Modulates the Intrinsic Conformational Landscape of Model Dipeptides

The intrinsic conformational preferences of small peptides may provide additional insight into the thermodynamics and kinetics of protein folding. In this study, we explore the underlying energy landscapes of two model peptides, namely, Ac-Ala-NH₂ and Ac-Ser-NH₂, using geometry-optimization-based tools developed within the context of energy landscape theory. We analyze not only how side-chain polarity influences the structural preferences of the dipeptides, but also other emergent properties of the landscape, including heat capacity profiles, and kinetics of conformational rearrangements. The contrasting topographies of the free energy landscape agree with recent results from Fourier transform microwave spectroscopy experiments, where Ac-Ala-NH₂ was found to exist as a mixture of two conformers, while Ac-Ser-NH₂ remained structurally locked, despite exhibiting an apparently rich conformational landscape.

Synthesis of cyclotetrapeptide analogues of c-PLAI and evaluation of their antimicrobial properties

Antimicrobial peptides (AMPs) are interesting compounds owing to their ability to kill several pathogens. In order to identify new AMPs, c-PLAI analogues were synthesized and evaluated together with their linear precursors for their antimicrobial properties against two Gram-positive bacteria (Staphylococcus aureus and Bacillus cereus), two Gram-negative bacteria (Escherichia coli and Klebsiella pneumoniae), and two fungal strains (Candida albicans and Trichophyton mentagrophytes). The new c-PLAI analogues were prepared through a combination of solid- and solution-phase syntheses, as previously employed for the synthesis of c-PLAI. The antimicrobial activity tests showed that the synthetic parent peptide c-PLAI was inactive or weakly active towards the bioindicators employed in the assay. The tests also indicated that cyclic c-PLAI analogues possessed enhanced antimicrobial properties against most of the bacteria and fungi tested. Furthermore, this study revealed that analogues containing cationic lysine residues displayed the highest activity towards most bioindicators. A combination of lysine and aromatic residues yielded analogues with broad-spectrum antimicrobial properties.

Strategies for Fine-Tuning the Conformations of Cyclic Peptides

Cyclic peptides are promising scaffolds for drug development, attributable in part to their increased conformational order compared to linear peptides. However, when optimizing the target-binding or pharmacokinetic properties of cyclic peptides, it is frequently necessary to "fine-tune" their conformations, e.g., by imposing greater rigidity, by subtly altering certain side chain vectors, or by adjusting the global shape of the macrocycle. This review systematically examines the various types of structural modifications that can be made to cyclic peptides in order to achieve such conformational control.

Bioinformatics and Biosimulations as Toolbox for Peptides and Peptidomimetics Design: Where Are We?

Peptides and peptidomimetics are strongly re-emerging as amenable candidates in the development of therapeutic strategies against a plethora of pathologies. In particular, these molecules are extremely suitable to treat diseases in which a major role is played by protein-protein interactions (PPIs). Unlike small organic compounds, peptides display both a high degree of specificity avoiding secondary off-targets effects and a relatively low degree of toxicity. Further advantages are provided by the possibility to easily conjugate peptides to functionalized nanoparticles, so improving their delivery and cellular uptake. In many cases, such molecules need to assume a specific three-dimensional conformation that resembles the bioactive one of the endogenous ligand. To this end, chemical modifications are introduced in the polypeptide chain to constrain it in a well-defined conformation, and to improve the drug-like properties. In this context, a successful strategy for peptide/peptidomimetics design and optimization is to combine different computational approaches ranging from structural bioinformatics to atomistic simulations. Here, we review the computational tools for peptide design, highlighting their main features and differences, and discuss selected protocols, among the large number of methods available, used to assess and improve the stability of the functional folding of the peptides. Finally, we introduce the simulation techniques employed to predict the binding affinity of the designed peptides for their targets.

Single-root networks for describing the potential energy surface of Lennard-Jones clusters

Potential energy surface (PES) holds the key in understanding a number of atomic clusters or molecular phenomena. However, due to the high dimension and incredible complexity of PES, only indirect methods can be used to characterize a PES of a given system in general. In this paper, a branched dynamic lattice searching method was developed to travel the PES, which was described in detail by a single-root network (SRN). The advantage of SRN is that it reflects the topological relation between different conformations and highlights the size of each structure energy trap. On the basis of SRN, to demonstrate how to transform one conformation to another, the transition path that connects two local minima in the PES was constructed. Herein, we take Lennard-Jones (LJ) clusters at the sizes of 38, 55, and 75 as examples. It is found that the PES of these three clusters have many local funnels and each local funnel represents one morphology. If a morphology is located more frequently, it will lie in a larger local funnel. Besides, certain steps of the transition path were generated successfully, such as changing from icosahedral to truncated octahedral of the LJ₃₈-cluster. Though we do not exhibit all the parts of the PES or all transition paths, this method indeed works well in the local area and can be used more widely.

Anti-proliferative potential of cyclotetrapeptides from Bacillus velezensis RA5401 and their molecular docking on G-Protein-Coupled Receptors

Elucidation of bioactive chemical compounds from rhizobacteria is highly utilized in pharmaceuticals and naturopathy, due to their health benefits to human and plants. In current study, four cyclopeptides along with one phenyl amide were isolated from the ethyl acetate extract of Bacillus velezensis sp. RA5401. Their structures were determined and characterized as cycle (L-prolyl-L-leucyl)₂ (1), cyclo (L-prolyl-l-valine)₂ (2), cycle (L-phenylanalyl-L-propyl)₂ (3), cyclo (D-pro-L-tyr-L-pro-L-tyr)₂ (4) and N-(2-phenylethyl)acetamide (5) on the basis of electron spray ionization mass spectrometry (ESI-MS), nuclear magnetic resonance (NMR) techniques and comparison with the literature data. The five compounds have been isolated for the first time from this species. The effect of various concentrations of these compounds on the proliferation of MDA-MB-231 breast cancer cells was examined. It was found that 1 and 2 induced concentration-independent anti-proliferative effects, while 3, 4 and 5 inhibited cancer cell proliferation in a concentration-dependent manner. Furthermore, to determine the suitable binding targets of these compounds within cancer cell line, detailed target prediction and comparative molecular-docking studies were performed. The compounds 1 and 2 hit intracellular anti-cancer targets of proteases family, while compounds 3, 4 and 5 interacted with different membrane receptors of G-Protein-Coupled Receptors (GPCRs). In conclusion, the Bacillus velezensis RA5401 can be an ideal strain to produce anti-proliferative constituents at industrial scale.

Designing Well-Structured Cyclic Pentapeptides Based on Sequence-Structure Relationships

Cyclic peptides are a promising class of molecules for unique applications. Unfortunately, cyclic peptide design is severely limited by the difficulty in predicting the conformations they will adopt in solution. In this work, we use explicit-solvent molecular dynamics simulations to design well-structured cyclic peptides by studying their sequence-structure relationships. Critical to our approach is an enhanced sampling method that exploits the essential transitional motions of cyclic peptides to efficiently sample their conformational space. We simulated a range of cyclic pentapeptides from all-glycine to a library of cyclo-(X1X2AAA) peptides to map their conformational space and determine cooperative effects of neighboring residues. By combining the results from all cyclo-(X1X2AAA) peptides, we developed a scoring function to predict the structural preferences for X1-X2 residues within cyclic pentapeptides. Using this scoring function, we designed a cyclic pentapeptide, cyclo-(GNSRV), predicted to be well structured in aqueous solution. Subsequent circular dichroism and NMR spectroscopy revealed that this cyclic pentapeptide is indeed well structured in water, with a nuclear Overhauser effect and J-coupling values consistent with the predicted structure.

Mapping the sequence-structure relationships of simple cyclic hexapeptides

Cyclic peptides are promising protein-protein interaction modulators with high binding affinities and specificities, as well as enhanced stabilities and oral availabilities over linear analogs. Despite their relatively small size and cyclic architecture, it is currently difficult to predict the favored conformation(s) of most classes of cyclic peptides. An improved understanding of the sequence-structure relationships for cyclic peptides will offer an avenue for the rational design of cyclic peptides as possible therapeutics. In this work, we systematically explored the sequence-structure relationships for two cyclic hexapeptide systems using molecular dynamics simulation techniques. Starting with an all-glycine cyclic hexapeptide, cyclo-G₆, we systematically replaced glycine residues with alanines and characterized the structural ensembles of different variants. The same process was repeated with valines to investigate the effects of larger side chains. An analysis of the origin of structure preferences was performed using thermodynamics decomposition and several general observations are reported.

Toward accurately modeling N-methylated cyclic peptides

Cyclic peptides have unique properties and can target protein surfaces specifically and potently. N-Methylation provides a promising way to further optimize the pharmacokinetic and structural profiles of cyclic peptides. The capability to accurately model structures adopted by N-methylated cyclic peptides would facilitate rational design of this interesting and useful class of molecules. We apply molecular dynamics simulations with advanced enhanced sampling methods to efficiently characterize the structural ensembles of N-methylated cyclic peptides, while simultaneously evaluating the overall performance of several simulation force fields. We find that one of the residue-specific force fields, RSFF2, is able to recapitulate experimental structures of the N-methylated cyclic peptide benchmarks tested here when the correct amide isomers are used as initial configurations and enforced during the simulations. Thus, using our simulation approach, it is possible to accurately and efficiently predict the structures of N-methylated cyclic peptides if sufficient information is available to determine the correct amide cis/trans configuration. Moreover, our results suggest that, upon further optimization of RSFF2 to more reliably predict cis/trans isomers, molecular dynamics simulations will be able to de novo predict N-methylated cyclic peptides in the near future, strongly motivating such continued optimization.

Energy Landscape and Pathways for Transitions between Watson-Crick and Hoogsteen Base Pairing in DNA

The recent discovery that Hoogsteen (HG) base pairs are widespread in DNA across diverse sequences and positional contexts could have important implications for understanding DNA replication and DNA-protein recognition. While evidence is emerging that the Hoogsteen conformation could be a thermodynamically accessible conformation of the DNA duplex and provide a means to expand its functionality, relatively little is known about the molecular mechanism underlying the Watson-Crick (WC) to HG transition. In this Perspective, we describe pathways and kinetics for this transition at an atomic level of detail, using the energy landscape perspective. We show that competition between the duplex conformations results in a double funnel landscape, which explains some recent experimental observations. The interconversion pathways feature a number of intermediates, with a variable number of WC and HG base pairs. The relatively slow kinetics, with possible deviations from two-state behavior, suggest that this conformational switch is likely to be a challenging target for both simulation and experiment.

Synthesis and characterization of water-soluble macrocyclic peptides stabilizing protein α-turn

Macrocyclic peptides mimic tight “non-classical” α-turn type II-α_LS found in proteins, as shown by spectroscopic and computational analysis of their equilibrating conformations.

Crystal and NMR Structures of a Peptidomimetic β-Turn That Provides Facile Synthesis of 13-Membered Cyclic Tetrapeptides

Herein we report the unique conformations adopted by linear and cyclic tetrapeptides (CTPs) containing 2-aminobenzoic acid (2-Abz) in solution and as single crystals. The crystal structure of the linear tetrapeptide H₂ N-d-Leu-d-Phe-2-Abz-d-Ala-COOH (1) reveals a novel planar peptidomimetic β-turn stabilized by three hydrogen bonds and is in agreement with its NMR structure in solution. While CTPs are often synthetically inaccessible or cyclize in poor yield, both 1 and its N-Me-d-Phe analogue (2) adopt pseudo-cyclic frameworks enabling near quantitative conversion to the corresponding CTPs 3 and 4. The crystal structure of the N-methylated peptide (4) is the first reported for a CTP containing 2-Abz and reveals a distinctly planar 13-membered ring, which is also evident in solution. The N-methylation of d-Phe results in a peptide bond inversion compared to the conformation of 3 in solution.

Far-infrared amide IV-VI spectroscopy of isolated 2- and 4-Methylacetanilide

Delocalized molecular vibrations in the far-infrared and THz ranges are highly sensitive to the molecular structure, as well as to intra- and inter-molecular interactions. Thus, spectroscopic studies of biomolecular structures can greatly benefit from an extension of the conventional mid-infrared to the far-infrared wavelength range. In this work, the conformer-specific gas-phase far-infrared spectra of two aromatic molecules containing the peptide -CO-NH- link, namely, 2- and 4-Methylacetanilide, are investigated. The planar conformations with trans configuration of the peptide link have only been observed in the supersonic-jet expansion. The corresponding far-infrared signatures associated with the vibrations of the peptide -CO-NH- moiety, the so-called amide IV-VI bands, have been assigned and compared with the results of density functional theory frequency calculations based on the anharmonic vibrational second-order perturbation theory approach. The analysis of the experimental and theoretical data shows that the amide IV-VI bands are highly diagnostic for the geometry of the peptide moiety and the molecular backbone. They are also strongly blue-shifted upon formation of the NH⋯O-C hydrogen bonding, which is, for example, responsible for the formation of secondary protein structures. Furthermore, the amide IV-VI bands are also diagnostic for the cis configuration of the peptide link, which can be present in cyclic peptides. The experimental gas-phase data presented in this work can assist the vibrational assignment of similar biologically important systems, either isolated or in natural environments.

The reactivity and conformational control of cyclic tetrapeptides derived from aziridine-containing amino acids

Among the smallest of the macrocyclic peptides, 12- and 13-membered cyclic tetrapeptides are particularly noteworthy because they exhibit a broad spectrum of biological activities due to their innate capacity to mimic β-turns in proteins. In this report, we demonstrate that aziridine-containing cyclic tetrapeptides offer a platform to interrogate the conformational properties of tetrapeptides. We show that aziridine ring-opening of 12-membered cyclic tetrapeptides yields exclusively 13-membered α3β macrocycles, regardless of peptide sequence, nucleophile, aziridine β-carbon substitution, or stereochemistry. NMR and computational studies on two related aziridine-containing cyclic tetrapeptides revealed that the amide conformations of their N-acyl aziridines are similar, and are likely the determinant of the observed ring-opening regioselectivity. Interestingly, some of the resulting 13-membered α3β macrocycles were found to be conformationally heterogeneous. This study on the reactivity and conformational control of aziridine-containing cyclic tetrapeptides provides useful insight on the design and development of macrocyclic therapeutics.

Computational methods to design cyclic peptides

Cyclic peptides (CPs) are promising modulators of protein-protein interactions (PPIs), but their application remains challenging. It is currently difficult to predict the structures and bioavailability of CPs. The ability to design CPs using computer modeling would greatly facilitate the development of CPs as potent PPI modulators for fundamental studies and as potential therapeutics. Herein, we describe computational methods to generate CP libraries for virtual screening, as well as current efforts to accurately predict the conformations adopted by CPs. These advances are making it possible to envision robust computational design of active CPs. However, unique properties of CPs pose significant challenges associated with sampling CP conformational space and accurately describing CP energetics. These major obstacles to structure prediction likely must be solved before robust design of active CPs can be reliably achieved.

Conformational Energy Landscape of the Ritonavir Molecule

Conformational polymorphism of ritonavir, a well-known pharmaceutical drug, is intricately linked to its efficacy in the treatment of acquired immunodeficiency syndrome (AIDS). Polymorphic transition from the crystalline form I to form II leads to the loss of bioactivity. The constituent ritonavir molecules adopt a trans configuration about the carbamate torsion angle in the form I crystal, and a cis configuration in the form II crystal. Investigating the energetics and mechanistic features of conformational transitions at the single molecule level is a key step toward decoding the complex features of the solid state polymorphism. In this work, we employ the energy landscape framework to investigate the conformational transitions of an isolated ritonavir molecule. The landscape is explored using discrete path sampling (DPS) and visualized in terms of disconnectivity graphs. We identify two distinct funnels corresponding to the two molecular forms that are identified by crystallography. The two regions can be reliably distinguished using the carbamate torsion angle, and the corresponding interconversion rates are predicted to follow Arrhenius behavior. The results provide mechanistic insight into pathways for cis ↔ trans interconversion at the molecular level and may also help in elucidating the polymorphic transitions in the crystal state.

Flexibility Coexists with Shape-Persistence in Cyanostar Macrocycles

Shape-persistent macrocycles are attractive functional targets for synthesis, molecular recognition, and hierarchical self-assembly. Such macrocycles are noncollapsible and geometrically well-defined, and they are traditionally characterized by having repeat units and low conformational flexibility. Here, we find it necessary to refine these ideas in the face of highly flexible yet shape-persistent macrocycles. A molecule is shape-persistent if it has a small change in shape when perturbed by external stimuli (e.g., heat, light, and redox chemistry). In support of this idea, we provide the first examination of the relationships between a macrocycle's shape persistence, its conformational space, and the resulting functions. We do this with a star-shaped macrocycle called cyanostar that is flexible as well as being shape-persistent. We employed molecular dynamics (MD), density functional theory (DFT), and NMR experiments. Considering a thermal bath as a stimulus, we found a single macrocycle has 332 accessible conformers with olefins undergoing rapid interconversion by up-down and in-out motions on short time scales (0.2 ns). These many interconverting conformations classify single cyanostars as flexible. To determine and confirm that cyanostars are shape-persistent, we show that they have a high 87% shape similarity across these conformations. To further test the idea, we use the binding of diglyme to the single macrocycle as guest-induced stimulation. This guest has almost no effect on the conformational space. However, formation of a 2:1 sandwich complex involving two macrocycles enhances rigidity and dramatically shifts the conformer distribution toward perfect bowls. Overall, the present study expands the scope of shape-persistent macrocycles to include flexible macrocycles if, and only if, their conformers have similar shapes.

Energy Landscapes and Global Optimization of Self-Assembling Cyclic Peptides

Self-assembled cyclic peptide nanotubes have attracted much attention because of their antimicrobial properties. Here, we present calculations on the formation of cyclic peptide dimers using basin-hopping and discrete path sampling. We present an analysis of the basin-hopping move sets that most efficiently explore the conformations of cyclic peptides. Group rotation moves, in which sections of the ring are rotated as a rigid body, are the most effective for cyclic peptides containing up to 20 residues. For cyclic peptide dimers, we find that a combination of group rotation intramolecular moves and rigid body intermolecular moves performs well. Discrete path sampling calculations on the cyclic peptide dimers show significant differences in the dimerization of hexa- and octapeptides.

Protein structure prediction: assembly of secondary structure elements by basin-hopping

The prediction of protein tertiary structure from primary structure remains a challenging task. One possible approach to this problem is the application of basin-hopping global optimization combined with an all-atom force field. In this work, the efficiency of basin-hopping is improved by introducing an approach that derives tertiary structures from the secondary structure assignments of individual residues. This approach is termed secondary-to-tertiary basin-hopping and benchmarked for three miniproteins: trpzip, trp-cage and ER-10. For each of the three miniproteins, the secondary-to-tertiary basin-hopping approach successfully and reliably predicts their three-dimensional structure. When it is applied to larger proteins, correctly folded structures are obtained. It can be concluded that the assembly of secondary structure elements using basin-hopping is a promising tool for de novo protein structure prediction.

Visualizing energy landscapes with metric disconnectivity graphs

The visualization of multidimensional energy landscapes is important, providing insight into the kinetics and thermodynamics of a system, as well the range of structures a system can adopt. It is, however, highly nontrivial, with the number of dimensions required for a faithful reproduction of the landscape far higher than can be represented in two or three dimensions. Metric disconnectivity graphs provide a possible solution, incorporating the landscape connectivity information present in disconnectivity graphs with structural information in the form of a metric. In this study, we present a new software package, PyConnect, which is capable of producing both disconnectivity graphs and metric disconnectivity graphs in two or three dimensions. We present as a test case the analysis of the 69-bead BLN coarse-grained model protein and show that, by choosing appropriate order parameters, metric disconnectivity graphs can resolve correlations between structural features on the energy landscape with the landscapes energetic and kinetic properties.

Molecular dynamics of the proline switch and its role in Crk signaling

The Crk adaptor proteins play a central role as a molecular timer for the formation of protein complexes including various growth and differentiation factors. The loss of regulation of Crk results in many kinds of cancers. A self-regulatory mechanism for Crk was recently proposed, which involves domain–domain rearrangement. It is initiated by a cis–trans isomerization of a specific proline residue (Pro238 in chicken Crk II) and can be accelerated by Cyclophilin A. To understand how the proline switch controls the autoinhibition at the molecular level, we performed large-scale molecular dynamics and metadynamics simulations in the context of short peptides and multidomain constructs of chicken Crk II. We found that the equilibrium and kinetic properties of the macrostates are regulated not only by the local environments of specified prolines but also by the global organization of multiple domains. We observe the two macrostates (cis closed/autoinhibited and trans open/uninhibited) consistent with NMR experiments and predict barriers. We also propose an intermediate state, the trans closed state, which interestingly was reported to be a prevalent state in human Crk II. The existence of this macrostate suggests that the rate of switching off the autoinhibition by Cyp A may be limited by the relaxation rate of this intermediate state.