Conformations and free energy landscapes of polyproline peptides

Clustering one million molecular structures on GPU within seconds

Structure clustering is a general but time-consuming work in the study of life science. Up to now, most published tools do not support the clustering analysis on graphics processing unit (GPU) with root mean square deviation metric. In this work, we specially write codes to do the work. It supports multiple threads on multiple GPUs. To show the performance, we apply the program to a 33-residue fragment in protein Pin1 WW domain mutant. The dataset contains 1,400,000 snapshots, which are extracted from an enhanced sampling simulation and distribute widely in the conformational space. Various testing results present that our program is quite efficient. Particularly, with two NVIDIA RTX4090 GPUs and single precision data type, the clustering calculation on 1 million snapshots is completed in a few seconds (including the uploading time of data from memory to GPU and neglecting the reading time from hard disk). This is hundreds of times faster than central processing unit. Our program could be a powerful tool for fast extraction of representative states of a molecule among its thousands to millions of candidate structures.

Understanding Polyproline's Unusual Thermoresponsive Properties Using a Polyproline-Based Double Hydrophilic Block Copolymer

Polyproline is a unique thermoresponsive polymer characterized by large thermal and conformational hysteresis. This article employs polyproline-based double hydrophilic block copolymers (PNIPAMn-b-PLPm) to gain insight into polyproline's thermoresponsive mechanism. The amine-terminated poly(N-isopropylacrylamide) (NH2-PNIPAMm) was used as the macroinitiator for ring-opening polymerization of proline-NCA monomers, resulting in various block copolymers (PNIPAMn-b-PLPm) with varying PLP block lengths. Block copolymers' thermal phase transitions were compared with their homopolymer counterparts using turbidimetry, variable-temperature NMR, dynamic light scattering, and circular dichroism spectroscopy. These experiments revealed that regardless of their compositions, all block copolymers exhibited a two-stage collapse (TCP(PLP) > TCP(PNIPAM)) during the heating cycle. In contrast, only one clearing temperature (TCL) was observed during cooling. The observed clearing temperature is closely correlated to the clearing temperature of PNIPAM blocks, suggesting the role of water-soluble PNIPAM blocks in resolving the PLP blocks. Moreover, thermal and conformational hysteresis related to the polyproline block is significantly suppressed in the presence of a PNIPAM block. Linking PNIPAM blocks has two significant effects on PLP segments' thermoresponsive behavior. For example, during the heating cycle, the precollapsed PNIPAM chains (as TCP(PNIPAM) < TCP(PLP)) prevent orderly aggregation within the PLP block. Meanwhile, during the cooling cycle below the clearing temperature of the PNIPAM block, the PNIPAM chains impart water solubility (as TCL(PNIPAM) > TCL(PLP)) to the collapsed PLP chains. Overall, the PNIPAM block imparts water solubility and perturbs PLP chains to form the native aggregate structure, suppressing the hysteresis effect. Accordingly, the large thermal and conformational hysteresis associated with native PLP chains appears to result from a noninterfering aggregation above the critical temperature.

Electrochemical Sensing of Phenylalanine using Polyaniline-Based Molecularly Imprinted Polymers

Polyaniline (PANI)-based molecularly imprinted polymers were investigated for their efficacy in sensing phenylalanine (Phe) when fabricated on both glassy carbon electrode (GCE) and indium tin oxide (ITO) sheets. This study highlights the superior performance of PANI-MIP/ITO over PANI-MIP/GCE for sensing Phe, with clear and distinct redox responses. Molecular computation helps to understand the interaction mechanism between PANI and Phe, where molecular crowding, aggregated clusters, hydrogen bonding, and π-π stacking facilitate stable interactions. We tested the specificity of Phe sensing by PANI-MIP with different amino acids such as cysteine, tryptophan, and tyrosine as well as organic molecules such as ascorbic acid, allantoin, sucrose, and urea, confirming its remarkable electrochemical efficiency. The oxidation response curve yielded a limit of detection of 4.88 μM and a limit of quantification of 16.3 μM, comparable to or better than earlier reported sensors. This work demonstrates the promise of MIP-based electrochemical sensing. It also lays the groundwork for future investigations into optimizing PANI-MIPs with nanocomposites to develop more selective and stable sensors.

Digital Light Processing of Thermoresponsive Hydrogels from Polyproline-Based Star Polypeptides

The first report of star poly(L-proline) crosslinkers is disclosed for digital light processing 3D printing of thermoresponsive hydrogels. Through chain end functionalization of star poly(L-proline)s with methacryloyl groups, access to high-resolution defined 3D hydrogel structures via digital light processing is achieved through photoinitiated free radical polymerization. Changing the poly(L-proline) molecular weight has a direct influence on both thermoresponsiveness and printability, while shape-morphing behavior can be induced thermally.

Proline, a unique amino acid whose polymer, polyproline II helix, and its analogues are involved in many biological processes: a review

Proline is a unique amino acid in that its side-chain is cyclised to the backbone, thus giving proline an exceptional rigidity and a considerably restricted conformational space. Polyproline forms two well-characterized helical structures: a left-handed polyproline helix (PPII) and a right-handed polyproline helix (PPI). Usually, sequences made only of prolyl residues are in PPII conformation, but even sequences not rich in proline but which are rich in glycine, lysine, glutamate, or aspartate have also a tendency to form PPII helices. Currently, the only way to study unambiguously PPII structure in solution is to use spectroscopies based on optical activity such as circular dichroism, vibrational circular dichroism and Raman optical activity. The importance of the PPII structure is emphasized by its ubiquitous presence in different organisms from yeast to human beings where proline-rich motifs and their binding domains are believed to be involved in vital biological processes. Some of the domains that are bound by proline-rich motifs include SH3 domains, WW domains, GYF domains and UEV domains, etc. The PPII structure has been demonstrated to be essential to biological activities such as signal transduction, transcription, cell motility, and immune response.

Proline Analogues

Within the canonical repertoire of the amino acid involved in protein biogenesis, proline plays a unique role as an amino acid presenting a modified backbone rather than a side-chain. Chemical structures that mimic proline but introduce changes into its specific molecular features are defined as proline analogues. This review article summarizes the existing chemical, physicochemical, and biochemical knowledge about this peculiar family of structures. We group proline analogues from the following compounds: substituted prolines, unsaturated and fused structures, ring size homologues, heterocyclic, e.g., pseudoproline, and bridged proline-resembling structures. We overview (1) the occurrence of proline analogues in nature and their chemical synthesis, (2) physicochemical properties including ring conformation and cis/trans amide isomerization, (3) use in commercial drugs such as nirmatrelvir recently approved against COVID-19, (4) peptide and protein synthesis involving proline analogues, (5) specific opportunities created in peptide engineering, and (6) cases of protein engineering with the analogues. The review aims to provide a summary to anyone interested in using proline analogues in systems ranging from specific biochemical setups to complex biological systems.

Deducing the conformational space for an octa-proline helix

A molecular dyad, PY-P8-PER, comprising a proline octamer sandwiched between pyrene and perylene terminals has been synthesized in order to address the dynamics of electronic energy transfer (EET) along the oligo-proline chain. A simple pyrene-based control compound equipped with a bis-proline attachment serves as a reference for spectroscopic studies. The N-H NMR signal at the terminal pyrene allows distinction between cis and trans amides and, although the crystal structure for the control has the trans conformation, temperature-dependent NMR studies provide clear evidence for trans/cis isomerisation in D6-DMSO. Polar solvents tend to stabilise the trans structure for the pyrene amide group, even for longer oligo-proline units. Circular dichroism shows that the proline spacer for PY-P8-PER exists mainly in the all-trans geometry in methanol. Preferential excitation of the pyrene chromophore is possible at wavelengths in the 320-350 nm range and, for the dyad, is followed by efficacious EET to the perylene emitter. The probability for intramolecular EET, obtained from analysis of steady-state spectroscopic data, is ca. 80-90% in solvents of disparate polarity. Comparison with the Förster critical distance suggests the terminals are ca. 18 Å apart. Time-resolved fluorescence spectroscopy, in conjunction with DFT calculations, indicates the dyad exists as a handful of conformers displaying a narrow range of EET rates. Optimisation of a distributive model allows accurate simulation of the EET dynamics in terms of reasonable structures based on isomerisation of certain amide groups.

Structural changes in chromosomes driven by multiple condensin motors during mitosis

We create a computational framework that utilizes loop extrusion (LE) by multiple condensin I/II motors to predict changes in chromosome organization during mitosis. The theory accurately reproduces the experimental contact probability profiles for the mitotic chromosomes in HeLa and DT40 cells. The LE rate is smaller at the start of mitosis and increases as the cells approach metaphase. Condensin II-mediated mean loop size is about six times larger than loops because of condensin I. The loops, which overlap each other, are stapled to a central dynamically changing helical scaffold formed by the motors during the LE process. A polymer physics-based data-driven method that uses the Hi-C contact map as the only input shows that the helix is characterized as random helix perversions (RHPs) in which the handedness changes randomly along the scaffold. The theoretical predictions, which are testable using imaging experiments, do not contain any parameters.

SARS-CoV-2 spike opening dynamics and energetics reveal the individual roles of glycans and their collective impact

The trimeric spike (S) glycoprotein, which protrudes from the SARS-CoV-2 viral envelope, binds to human ACE2, initiated by at least one protomer's receptor binding domain (RBD) switching from a "down" (closed) to an "up" (open) state. Here, we used large-scale molecular dynamics simulations and two-dimensional replica exchange umbrella sampling calculations with more than a thousand windows and an aggregate total of 160 μs of simulation to investigate this transition with and without glycans. We find that the glycosylated spike has a higher barrier to opening and also energetically favors the down state over the up state. Analysis of the S-protein opening pathway reveals that glycans at N165 and N122 interfere with hydrogen bonds between the RBD and the N-terminal domain in the up state, while glycans at N165 and N343 can stabilize both the down and up states. Finally, we estimate how epitope exposure for several known antibodies changes along the opening path. We find that the BD-368-2 antibody's epitope is continuously exposed, explaining its high efficacy.

Simulating Polyproline II-Helix-Rich Peptides with the Latest Kirkwood-Buff Force Field: A Direct Comparison with AMBER and CHARMM

We simulated the dynamics of a set of peptides characterized by ensembles rich in PPII-helical content, to assess the ability of the most recent Kirkwood-Buff force field (KBFF20) to sample this conformational peculiarity. KBFF has been previously shown to capably reproduce experimental dimensions of disordered proteins, while being limited in confidently sampling structured proteins. Further development of the force field bridged this gap. It is however still unknown what are the main differences between KBFF and AMBER/CHARMM force fields. A direct comparison is now possible as both AMBER/CHARMM force fields have been used to sample peptides rich in PPII-helical content. We found that KBFF20 samples' PPII-helical content qualitatively matches both AMBER and CHARMM force fields, with the main difference being the KBFF ability to populate the α_R region of the Ramachandran plot in the set of simulated peptides. Overall, KBFF20 is a well-balanced force field, able to sample the dynamics of both structured and unstructured proteins.

Light-Harvesting Crystals Formed from BODIPY-Proline Biohybrid Conjugates: Antenna Effects and Excitonic Coupling

A boron dipyrromethene (BODIPY) derivative bearing a cis-proline residue at the meso-position crystallizes in the form of platelets with strong (i.e., Φ_F = 0.34) red fluorescence, but the absorption and emission spectra differ markedly from those for dilute solutions. A key building block for the crystal is a pseudo-dimer where hydrogen bonding aligns the proline groups and separates the terminal chromophores by ca. 25 Å. Comparison with a covalently linked bichromophore suggests that one-dimensional (1D) excitonic coupling between the terminals is too small to perturb the optical properties. However, accretion of the pseudo-dimer forms narrow channels possessing a high density of chromophores. The resultant absorption spectrum exhibits strong excitonic splitting, which can be explained quantitatively using the extended dipole approach and allowing for coupling between ca. 30 BODIPY units. Fluorescence, which decays with a lifetime of 2.2 ns, is assigned to a delocalized and (slightly) super-radiant BODIPY dimer situated at the interface and populated via electronic energy transfer from the interior.

An Unbound Proline-Rich Signaling Peptide Frequently Samples Cis Conformations in Gaussian Accelerated Molecular Dynamics Simulations

Disordered proline-rich motifs are common across the proteomes of many species and are often involved in protein-protein interactions. Proline is a unique amino acid due to the covalent bond between the backbone nitrogen and the proline side chain. The resulting five-membered ring allows proline to sample the cis state about its peptide bond, which other residues cannot do as readily. Because proline-rich disordered sequences exist as ensembles that likely include structures with the proline peptide bond in cis, a robust methodology to accurately account for these conformations in the overall ensemble is crucial. Observing the cis conformations of proline in a disordered sequence is challenging both experimentally and computationally. Nitrogen-hydrogen NMR spectroscopy cannot directly observe proline residues, which lack an amide bond, and computational methods struggle to overcome the large kinetic barrier between the cis and trans states, since isomerization usually occurs on the order of seconds. In the current work, Gaussian accelerated molecular dynamics was used to overcome this free energy barrier and simulate proline isomerization in a tetrapeptide (KPTP) and in the 12-residue proline-rich SH3 binding peptide, ArkA. We found that Gaussian accelerated molecular dynamics, when combined with a lowered peptide bond dihedral angle potential energy barrier (15 kcal/mol), allowed sufficient sampling of the proline cis and trans states on a microsecond timescale. All ArkA prolines spend a significant fraction of time in cis, leading to a more compact ensemble with less polyproline II helix structure than an ArkA ensemble with all peptide bonds in trans. The ensemble containing cis prolines also matches more closely to in vitro circular dichroism data than the all-trans ensemble. The ability of the ArkA prolines to isomerize likely affects the peptide's ability to bind its partner SH3 domain, and should be studied further. This is the first molecular dynamics simulation study of proline isomerization in a biologically relevant proline-rich sequence that we know of, and a similar protocol could be applied to study multi-proline isomerization in other proline-containing proteins to improve conformational diversity and agreement with in vitro data.

Unlocking the Stereoselectivity and Substrate Acceptance of Enzymes: Proline-Induced Loop Engineering Test

Protein stability and evolvability influence each other. Although protein dynamics play essential roles in various catalytically important properties, their high flexibility and diversity makes it difficult to incorporate such properties into rational engineering. Therefore, how to unlock the potential evolvability in a user-friendly rational design process remains a challenge. In this endeavor, we describe a method for engineering an enantioselective alcohol dehydrogenase. It enables synthetically important substrate acceptance for 4-chlorophenyl pyridine-2-yl ketone, and perfect stereocontrol of both (S)- and (R)-configured products. Thermodynamic analysis unveiled the subtle interaction between enzyme stability and evolvability, while computational studies provided insights into the origin of selectivity and substrate recognition. Preparative-scale synthesis of the (S)-product (73 % yield; >99 % ee) was performed on a gram-scale. This proof-of-principle study demonstrates that interfaced proline residues can be rationally engineered to unlock evolvability and thus provide access to new biocatalysts with highly improved catalytic performance.

Polar and Helical Isomorphous Crystals of Proline Derivatives: Influence of a Fluorine Atom on the Electric Susceptibility

Two novel nonlinear optical isomorphous crystals of proline derivatives with alkyne functionality have been obtained (Boc-L-ProNH(CH2)2CCH and Boc-cis-4-fluoro-L-ProNH(CH2)2CCH). Both derivatives, which differ only by the substitution of a H atom to a F atom, adopt the same polar and columnar right-handed helix arrangement in the crystalline state. In addition, adjacent polar helical columns all point in the same direction, thus generating a macrodipole and a crystalline system conducive for second harmonic generation (SHG) properties. This isomorphous crystal system constitutes an interesting tool to study the effect of the fluorine atom on the dipole moment and on the first hyperpolarizability. Starting from the PBC optimized geometries of the crystals, the macroscopic second-order nonlinearity, χ(2), of the newly synthesized crystals has been estimated by quantum chemical calculations. These χ(2) responses are of the same order of magnitude as those of inorganic proline derivatives while smaller than those observed in crystals of push–pull π-conjugated molecules.

Graphic Abstract

Infrared harmonic features of collagen models at B3LYP-D3: From amide bands to the THz region

In this paper, we have studied the vibrational spectral features for the collagen triple helix using a dispersion corrected hybrid density functional theory (DFT-D) approach. The protein is simulated by an infinite extended polymer both in the gas phase and in a water micro-solvated environment. We have adopted proline-rich collagen models in line with the high content of proline in natural collagens. Our scaled harmonic vibrational spectra are in very good agreement with the experiments and allow for the peak assignment of the collagen amide I and III bands, supporting or questioning the experimental interpretation by means of vibrational normal modes analysis. Furthermore, we demonstrated that IR spectroscopy in the THz region can detect the small variations inherent to the triple helix helicity (10/3 over 7/2), thus elucidating the packing state of the collagen. So far, identifying the collagen helicity is only possible by means of crystal x-ray diffraction.

Investigating the folding mechanism of the N-terminal domain of ribosomal protein L9

Protein folding is a popular topic in the life science. However, due to the limited sampling ability of experiments and simulations, the general folding mechanism is not yet clear to us. In this work, we study the folding of the N-terminal domain of ribosomal protein L9 (NTL9) in detail by a mixing replica exchange molecular dynamics method. The simulation results are close to previous experimental observations. According to the Markov state model, the folding of the protein follows a nucleation-condensation path. Moreover, after the comparison to its 39-residue β-α-β motif, we find that the helix at the C-terminal has a great influence on the folding process of the intact protein, including the nucleation of the key residues in the transition state ensemble and the packing of the hydrophobic residues in the native state.

Simulating selective binding of a biological template to a nanoscale architecture: a core concept of a clamp-based binding-pocket-favored N-terminal-domain assembly

The biological template and its mutants have vital significance in next generation remediation, electrochemical, photovoltaic, catalytic, sensing and digital memory devices. However, a microscopic model describing the biotemplating process is generally lacking on account of modelling complexity, which has prevented widespread commercial use of biotemplates. Here, we demonstrate M13-biotemplating kinetics in atomic resolution by leveraging large-scale molecular dynamics (MD) simulations. The model reveals the assembly of gold nanoparticles on two experimentally-based M13 phage types using full M13-capsid structural models and with polarizable gold nanoparticles in explicit solvent. Both mechanistic and structural insights into the selective binding affinity of the M13 phage to gold nanoparticles are obtained based on a previously unconsidered clamp-based binding-pocket-favored N-terminal-domain assembly and also on surface-peptide flexibility. These results provide a deeper level of understanding of protein sequence-based affinity and open the route for genetically engineering a wide range of 3D electrodes for high-density low-cost device integration.

A Global Minimization Toolkit for Batch-Fitting and χ2 Cluster Analysis of CW-EPR Spectra

Electron paramagnetic resonance spectroscopy (EPR) is a uniquely powerful technique for characterizing conformational dynamics at specific sites within a broad range of molecular species in water. Computational tools for fitting EPR spectra have enabled dynamics parameters to be determined quantitatively. These tools have dramatically broadened the capabilities of EPR dynamics analysis, however, their implementation can easily lead to overfitting or problems with self-consistency. As a result, dynamics parameters and associated properties become difficult to reliably determine, particularly in the slow-motion regime. Here, we present an EPR analysis strategy and the corresponding computational tool for batch-fitting EPR spectra and cluster analysis of the χ2 landscape in Linux. We call this tool CSCA (Chi-Squared Cluster Analysis). The CSCA tool allows us to determine self-consistent rotational diffusion rates and enables calculations of activation energies of diffusion from Arrhenius plots. We demonstrate CSCA using a model system designed for EPR analysis: a self-assembled nanoribbon with radical electron spin labels positioned at known distances off the surface. We anticipate that the CSCA tool will increase the reproducibility of EPR fitting for the characterization of dynamics in biomolecules and soft matter.

Simulation Study of the Plasticity of k-Turn Motif in Different Environments

The k-turn is a widespread and important motif in RNA. According to the internal hydrogen bond network, it has two stable states, called N1 and N3. The relative stability between the states changes with the environment. It is able to accept different conformations in different environments. This is called the "plasticity" of a molecule. In this work, we study the plasticity of k-turn by the mixing REMD method in explicit solvent. The results are concluded as follows. First, N1 and N3 are almost equally stable when k-turn is in the solvent alone. The molecule is quite flexible as a hinge. However, after binding to different proteins, such as the proteins L7Ae and L24e, k-turn falls into one global minimum. The preferred state could be either N1 or N3. On the contrary, the other nonpreferred state becomes unstable with a weaker binding affinity to the protein. It reveals that RNA-binding protein is able to modulate the representative state of k-turn at equilibrium. This is in agreement with the findings in experiments. Moreover, free energy calculations show that the free energy barrier between the N1 and N3 states of k-turn increases in the complexes. The state-to-state transition is greatly impeded. We also give a deep discussion on the mechanism of the high plasticity of k-turn in different environments.

Allosteric inhibition of human exonuclease1 (hExo1) through a novel extended β-sheet conformation

BACKGROUND

Human Exonuclease1 (hExo1) participates in the resection of DNA double-strand breaks by generating long 3'-single-stranded DNA overhangs, critical for homology-based DNA repair and activation of the ATR-dependent checkpoint. The C-terminal region is essential for modulating the activity of hExo1, containing numerous sites of post-translational modification and binding sites for partner proteins.

METHODS

Analytical Ultracentrifugation (AUC), Dynamic Light Scattering (DLS), Circular Dichroism (CD) spectroscopy and enzymatic assays.

RESULTS

AUC and DLS indicates the C-terminal region has a highly extended structure while CD suggest a tendency to adopt a novel left-handed β-sheet structure, together implying the C-terminus may exhibit a transient fluctuating structure that could play a role in binding partner proteins known to regulate the activity of hExo1. Interaction with 14-3-3 protein has a cooperative inhibitory effect upon DNA resection activity, which indicates an allosteric transition occurs upon binding partner proteins.

CONCLUSIONS

This study has uncovered that hExo1 consist of a folded N-terminal nuclease domain and a highly extended C-terminal region which is known to interact with partner proteins that regulates the activity of hExo1. A positively cooperative mechanism of binding allows for stringent control of hExo1 activity. Such a transition would coordinate the control of hExo1 by hExo1 regulators and hence allow careful coordination of the process of DNA end resection.

SIGNIFICANCE

The assays presented herein could be readily adapted to rapidly identify and characterise the effects of modulators of the interaction between the 14-3-3 proteins and hExo1. It is conceivable that small molecule modulators of 14-3-3 s-hExo1 interaction may serve as effective chemosensitizers for cancer therapy.

Energy Landscapes of Deoxyxylo- and Xylo-Nucleic Acid Octamers

Artificial analogues of the natural nucleic acids have attracted interest as a diverse class of information storage molecules capable of self-replication. In this study, we use the computational potential energy landscape framework to investigate the structural and dynamical properties of xylo- and deoxyxylo-nucleic acids (XyNA and dXyNA), which are derived from their respective RNA and DNA analogues by inversion of a single chiral center in the sugar moiety of the nucleotides. For an octameric XyNA sequence and the analogue dXyNA, we observe facile conformational transitions between a left-handed helix, which is the free energy global minimum, and a ladder-type structure with approximately zero helicity. The competing ensembles are better separated in the dXyNA, making it a more suitable candidate for a molecular switch, whereas the XyNA exhibits additional flexibility. Both energy landscapes exhibit greater frustration than we observe in RNA or DNA, in agreement with the higher degree of optimization expected from the principle of minimal frustration in evolved biomolecules.

Flexible and Extended Linker Domains Support Efficient Targeting of Heh2 to the Inner Nuclear Membrane

The nuclear pore complex (NPC) is embedded in the nuclear envelope and forms the main gateway to the nuclear interior including the inner nuclear membrane (INM). Two INM proteins in yeast are selectively imported. Their sorting signals consist of a nuclear localization signal, separated from the transmembrane domain by a long intrinsically disordered (ID) linker. We used computational models to predict the dynamic conformations of ID linkers and analyzed the INM targeting efficiency of proteins with linker regions with altered Stokes radii and decreased flexibilities. We find that flexibility, Stokes radius, and the frequency at which the linkers are at an extended end-to-end distance larger than 25 nm are good predictors for the targeting of the proteins. The data are consistent with a transport mechanism in which INM targeting of Heh2 is dependent on an ID linker that facilitates the crossing of the approximately 25-nm thick NPC scaffold.

FSATOOL: A useful tool to do the conformational sampling and trajectory analysis work for biomolecules

Reliable conformational sampling and trajectory analysis are always important to the study of the folding or binding mechanisms of biomolecules. Generally, one has to prepare many complicated parameters and follow a lot of steps to obtain the final data. The whole process is too complicated to new users. In this article, we provide a convenient and user-friendly tool that is compatible to AMBER, called fast sampling and analysis tool (FSATOOL). FSATOOL has some useful features. First and the most important, the whole work is extremely simplified into two steps, one is the fast sampling procedure and the other is the trajectory analysis procedure. Second, it contains several powerful sampling methods for the simulation on graphics process unit, including our previous mixing replica exchange molecular dynamics method. The method combines the advantages of the biased and unbiased simulations. Finally, it extracts the dominant transition pathways automatically from the folding network by Markov state model. Users do not need to do the tedious intermediate steps by hand. To illustrate the usage of FSATOOL in practice, we perform one simulation for a RNA hairpin in explicit solvent. All the results are presented. © 2019 Wiley Periodicals, Inc.

Investigation of the cis–trans structures and isomerization of oligoprolines by using Raman spectroscopy and density functional theory calculations: solute–solvent interactions and effects of terminal positively charged amino acid residues

Using low-wavenumber Raman spectroscopy in combination with theoretical calculations via solid-state density functional theory (DFT)-D3, we studied the vibrational structures and interaction with solvent of poly-l-proline and the oligoproline P12 series. The P12 series includes P12, the positively charged amino acid residue (arginine and lysine) N-terminus proline oligomers RP11 and KP11, and the C-terminus P11R and P11K. We assigned the spring-type phonon mode to 74–76 cm⁻¹ bands for the PPI and PPII conformers and the carbonyl group ring-opening mode 122 cm⁻¹ in the PPI conformer of poly-l-proline. Amide I and II were assigned based on the results of mode analysis for O, N, and C atom displacements. The broad band feature of the H-bond transverse mode in the Raman spectra indicates that the positively charged proline oligomers PPII form H-bonds with water in the solid phase, whereas P12 is relatively more hydrophobic. In propanol, the PPI conformer of the P12 series forms less H-bond network with the solvent. The PPII conformer exhibits a distinct Raman band at 310 cm⁻¹, whereas the PPI has bands at 365, 660, and 960 cm⁻¹ with reasonable intensity that can be used to quantitatively determine these two conformational forms. The 365 cm⁻¹ mode comprising the motion of a C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O group turning to the helix axis was used to monitor the isomerization reaction PPI ↔ PPII. In pure propanol, RP11 and KP11 were found to have mostly PPI present, but P11R and P11K preferred PPII. After adding 20% water, the PPI in P11R and P11K was completely converted to PPII, whereas a small fraction of PPI remained in RP11 and KP11. The substituted positively charged amino acid affected the balance of the PPI/PPII population ratio.

The low-wavenumber Raman spectra in combination with theoretical calculations via solid-state density functional theory (DFT)-D3 are displayed. The vibrational structures and interaction with solvent of poly-l-proline and the oligoproline P12 series are identified.

Combining the biased and unbiased sampling strategy into one convenient free energy calculation method

Constructing a free energy landscape for a large molecule is difficult. One has to use either a high temperature or a strong driving force to enhance the sampling on the free energy barriers. In this work, we propose a mixed method that combines these two kinds of acceleration strategies into one simulation. First, it applies an adaptive biasing potential to some replicas of the molecule. These replicas are particularly accelerated in a collective variable space. Second, it places some unbiased and exchangeable replicas at various temperature levels. These replicas generate unbiased sampling data in the canonical ensemble. To improve the sampling efficiency, biased replicas transfer their state variables to the unbiased replicas after equilibrium by Monte Carlo trial moves. In comparison to previous integrated methods, it is more convenient for users. It does not need an initial reference biasing potential to guide the sampling of the molecule. And it is also unnecessary to insert many replicas for the requirement of passing the free energy barriers. The free energy calculation is accomplished in a single stage. It samples the data as fast as a biased simulation and it processes the data as simple as an unbiased simulation. The method provides a minimalist approach to the construction of the free energy landscape. © 2019 Wiley Periodicals, Inc.

Probing the functional conformations of an atypical proline-rich fusion peptide

Simulations confirm a propensity for extended and solvent exposed conformations of the p15 fusion peptide capable of membrane targeting.

The structure of Legionella effector protein LpnE provides insights into its interaction with Oculocerebrorenal syndrome of Lowe (OCRL) protein

Legionella pneumophila is a freshwater bacterium that replicates in predatory amoeba and alveolar macrophage. The ability of L. pneumophila to thrive in eukaryotic host cells is conferred by the Legionella containing vacuole (LCV). Formation and intracellular trafficking of the LCV are governed by an arsenal of effector proteins, many of which are secreted by the Icm/Dot Type 4 Secretion System. One such effector, known as LpnE (L. pneumophila Entry), has been implicated in facilitating bacterial entry into host cells, LCV trafficking, and substrate translocation. LpnE belongs to a subfamily of tetratricopeptide repeat proteins known as Sel1-like repeats (SLRs). All eight of the predicted SLRs in LpnE are required to promote host cell invasion. Herein, we report that LpnE(1-375) localizes to cis-Golgi in HEK293 cells via its signal peptide (aa 1-22). We further verify the interaction of LpnE(73-375) and LpnE(22-375) with Oculocerebrorenal syndrome of Lowe protein (OCRL) residues 10-208, restricting the known interacting residues for both proteins. To further characterize the SLR region of LpnE, we solved the crystal structure of LpnE(73-375) to 1.75Å resolution. This construct comprises all SLRs, which are arranged in a superhelical fold. The α-helices forming the inner concave surface of the LpnE superhelix suggest a potential protein-protein interaction interface. DATABASE: Coordinates and structure factors were deposited in the Protein Data Bank with the accession number 6DEH.

A spin-1 representation for dual-funnel energy landscapes

The interconversion between the left- and right-handed helical folds of a polypeptide defines a dual-funneled free energy landscape. In this context, the funnel minima are connected through a continuum of unfolded conformations, evocative of the classical helix-coil transition. Physical intuition and recent conjectures suggest that this landscape can be mapped by assigning a left- or right-handed helical state to each residue. We explore this possibility using all-atom replica exchange molecular dynamics and an Ising-like model, demonstrating that the energy landscape architecture is at odds with a two-state picture. A three-state model-left, right, and unstructured-can account for most key intermediates during chiral interconversion. Competing folds and excited conformational states still impose limitations on the scope of this approach. However, the improvement is stark: Moving from a two-state to a three-state model decreases the fit error from 1.6 kBT to 0.3 kBT along the left-to-right interconversion pathway.

About Chemical Strategies to Fabricate Cell-Instructive Biointerfaces with Static and Dynamic Complexity

Properly functioning cell-instructive biointerfaces are critical for healthy integration of biomedical devices in the body and serve as decisive tools for the advancement of our understanding of fundamental cell biological phenomena. Studies are reviewed that use covalent chemistries to fabricate cell-instructive biointerfaces. These types of biointerfaces typically result in a static presentation of predefined cell-instructive cues. Chemically defined, but dynamic cell-instructive biointerfaces introduce spatiotemporal control over cell-instructive cues and present another type of biointerface, which promises a more biomimetic way to guide cell behavior. Therefore, strategies that offer control over the lateral sorting of ligands, the availability and molecular structure of bioactive ligands, and strategies that offer the ability to induce physical, chemical and mechanical changes in situ are reviewed. Specific attention is paid to state-of-the-art studies on dynamic, cell-instructive 3D materials. Future work is expected to further deepen our understanding of molecular and cellular biological processes investigating cell-type specific responses and the translational steps toward targeted in vivo applications.

Theoretical and NMR Conformational Studies of β-Proline Oligopeptides With Alternating Chirality of Pyrrolidine Units

Synthetic β-peptides are potential functional mimetics of native α-proteins. A recently developed, novel, synthetic approach provides an effective route to the broad group of β-proline oligomers with alternating patterns of stereogenic centers. Conformation of the pyrrolidine ring, Z/E isomerism of β-peptide bonds, and hindered rotation of the neighboring monomers determine the spatial structure of this group of β-proline oligopeptides. Preferences in their structural organization and corresponding thermodynamic properties are determined by NMR spectroscopy, restrained molecular dynamics and quantum mechanics. The studied β-proline oligopeptides exist in dimethyl sulfoxide solution in a limited number of conformers, with compatible energy of formation and different spatial organization. In the β-proline tetrapeptide with alternating chirality of composing pyrrolidine units, one of three peptide bonds may exist in an E configuration. For the alternating β-proline pentapeptide, the presence of an E configuration for at least of one β-peptide bond is mandatory. In this case, three peptide bonds synchronously change their configurations. Larger polypeptides may only exist in the presence of several E configurations of β-peptide bonds forming a wave-like extended structure.

Probing helical transitions in a DNA duplex

The complex conformational change from B-DNA to Z-DNA requires inversion of helix-handedness. Multiple degrees of freedom are intricately coupled during this transition, and formulating an appropriate reaction coordinate that captures the underlying complexity would be problematic. In this contribution, we adopt an alternative approach, based on the potential energy landscape perspective, to construct a kinetic transition network. Microscopic insight into the B → Z transition is provided in terms of geometrically defined discrete paths consisting of local minima and the transition states that connect them. We find that the inversion of handedness can occur via two competing mechanisms, either involving stretched intermediates, or a B-Z junction, in agreement with previous predictions. The organisation of the free energy landscape further suggests that this process is likely to be slow under physiological conditions. Our results represent a key step towards decoding the more intriguing features of the B → Z transition, such as the role of ionic strength and negative supercoiling in reshaping the landscape.

Constructing a multidimensional free energy surface like a spider weaving a web

Complete free energy surface in the collective variable space provides important information of the reaction mechanisms of the molecules. But, sufficient sampling in the collective variable space is not easy. The space expands quickly with the number of the collective variables. To solve the problem, many methods utilize artificial biasing potentials to flatten out the original free energy surface of the molecule in the simulation. Their performances are sensitive to the definitions of the biasing potentials. Fast-growing biasing potential accelerates the sampling speed but decreases the accuracy of the free energy result. Slow-growing biasing potential gives an optimized result but needs more simulation time. In this article, we propose an alternative method. It adds the biasing potential to a representative point of the molecule in the collective variable space to improve the conformational sampling. And the free energy surface is calculated from the free energy gradient in the constrained simulation, not given by the negative of the biasing potential as previous methods. So the presented method does not require the biasing potential to remove all the barriers and basins on the free energy surface exactly. Practical applications show that the method in this work is able to produce the accurate free energy surfaces for different molecules in a short time period. The free energy errors are small in the cases of various biasing potentials. © 2017 Wiley Periodicals, Inc.

Structure and Dynamics of DNA and RNA Double Helices of CAG and GAC Trinucleotide Repeats

CAG trinucleotide repeats are known to cause 10 late-onset progressive neurodegenerative disorders as the repeats expand beyond a threshold, whereas GAC repeats are associated with skeletal dysplasias and expand from the normal five to a maximum of seven repeats. The TR secondary structure is believed to play a role in CAG expansions. We have carried out free energy and molecular dynamics studies to determine the preferred conformations of the A-A noncanonical pairs in (CAG)_n and (GAC)_n trinucleotide repeats (n = 1, 4) and the consequent changes in the overall structure of the RNA and DNA duplexes. We find that the global free energy minimum corresponds to A-A pairs stacked inside the core of the helix with anti-anti conformations in RNA and (high-anti)-(high-anti) conformations in DNA. The next minimum corresponds to anti-syn conformations, whereas syn-syn conformations are higher in energy. Transition rates of the A-A conformations are higher for RNA than DNA. Mechanisms for these various transitions are identified. Additional structural and dynamical aspects of the helical conformations are explored, with a focus on contrasting CAG and GAC duplexes. The neutralizing ion distribution around the noncanonical pairs is described.

Infrared Multiple-Photon Dissociation Action Spectroscopy of the b2+ Ion from PPG: Evidence of Third Residue Affecting b2+ Fragment Structure

Infrared multiple-photon dissociation (IRMPD) action spectroscopy was performed on the b₂⁺ fragment ion from the protonated PPG tripeptide. Comparison of the experimental infrared spectrum with computed spectra for both oxazolone and diketopiperazine structures indicates that the majority of the fragment ion population has an oxazolone structure with the remainder having a diketopiperazine structure. This result is in contrast with a recent study of the IRMPD action spectrum of the PP b₂⁺ fragment ion from PPP, which was found to be nearly 100% diketopiperazine (Martens et al. Int. J. Mass Spectrom. 2015, 377, 179). The diketopiperazine b₂⁺ ion is thermodynamically more stable than the oxazolone but normally requires a trans/cis peptide bond isomerization in the dissociating peptide. Martens et al. showed through IRMPD action spectroscopy that the PPP precursor ion was in a conformation in which the first peptide bond is already in the cis conformation and thus it was energetically favorable to form the thermodynamically-favored diketopiperazine b₂⁺ ion. In the present case, solution-phase NMR spectroscopy and gas-phase IRMPD action spectroscopy show that the PPG precursor ion has its first amide bond in a trans configuration suggesting that the third residue is playing an important role in both the structure of the peptide and the associated ring-closure barriers for oxazolone and diketopiperazine formation. Graphical Abstract ᅟ.

Fast exploration of an optimal path on the multidimensional free energy surface

In a reaction, determination of an optimal path with a high reaction rate (or a low free energy barrier) is important for the study of the reaction mechanism. This is a complicated problem that involves lots of degrees of freedom. For simple models, one can build an initial path in the collective variable space by the interpolation method first and then update the whole path constantly in the optimization. However, such interpolation method could be risky in the high dimensional space for large molecules. On the path, steric clashes between neighboring atoms could cause extremely high energy barriers and thus fail the optimization. Moreover, performing simulations for all the snapshots on the path is also time-consuming. In this paper, we build and optimize the path by a growing method on the free energy surface. The method grows a path from the reactant and extends its length in the collective variable space step by step. The growing direction is determined by both the free energy gradient at the end of the path and the direction vector pointing at the product. With fewer snapshots on the path, this strategy can let the path avoid the high energy states in the growing process and save the precious simulation time at each iteration step. Applications show that the presented method is efficient enough to produce optimal paths on either the two-dimensional or the twelve-dimensional free energy surfaces of different small molecules.

Molecular Dynamics Study of Nitrogen-Pyramidalized Bicyclic β-Proline Oligomers: Length-Dependent Convergence to Organized Structures

In this study, the solution structures of the homooligomers of a conformationally constrained bicyclic proline-type β-amino acid were studied by means of molecular dynamics (MD) calculations in explicit methanol and water using the umbrella sampling method. The ratio of trans-amide and cis-amide was estimated by NMR and the rotational barrier of the amide of acetylated bicyclic amino acid monomer was estimated by two-dimensional (2D) exchange spectroscopy (EXSY) or line-shape analysis. A bias potential was introduced with respect to the amide torsion angle ω to enhance conformational exchange including isomerization of amide bonds by lowering the rotation energy barrier. After determination of reweighting parameters to best reproduce the experimental results of the monomer amide, the free energy profile around the amide torsion angle ω was obtained from the MD trajectory by reweighting of the biased probability density. The MD simulation results support the existence of invertomers of nitrogen-pyramidalized amide. Furthermore, extended structures with a high fraction of trans-amide conformation appear to be increasingly stabilized as the oligomer is elongated, both in methanol and in water. Our conformational analysis of natural and non-natural tertiary-amide-based peptide oligomers indicates that these oligomers preferentially adopt a limited number of conformations.

Effective Riemannian Diffusion Model for Conformational Dynamics of Biomolecular Systems

We present a Riemannian formalism for effective diffusion of biomolecules in collective variable spaces that provides a robust framework for conformational free energy calculation methods. Unlike their Euclidean counterparts, the Riemannian potential of mean force (PMF) and minimum free energy path (MFEP) are invariant under coordinate transformations. The presented formalism can be readily employed to modify the collective variable based enhanced sampling techniques, such as umbrella sampling (US) commonly used in biomolecular simulations, to take into account the role of intrinsic geometry of collective variable space. Although our model is mathematically equivalent to a Euclidean diffusion with a position-dependent diffusion tensor, the Riemannian formulation provides a more convenient framework for free energy calculation methods and path-finding algorithms aimed at characterizing the effective conformational dynamics of biomolecules. A simple three-dimensional toy model and a pentapeptide (met-enkephalin) simulated in an explicit solvent environment are used to illustrate the workings of the formalism and its implementation.

Structure activity relationship study on the peptide hormone preptin, a novel bone-anabolic agent for the treatment of osteoporosis

Preptin is a 34-residue pancreatic hormone shown to be anabolic to bone in vitro and in vivo. The bone activity of preptin resides within the (1-16) N-terminal fragment. Due to its peptidic nature, the truncated fragment of preptin is enzymatically unstable; however it provides an attractive framework for the creation of stable analogues using various peptidomimetic techniques. An alanine scan of preptin (1-16) was undertaken which showed that substitution of Ser at position 3 or Pro at position 14 did not inhibit the proliferative activity of preptin in primary rat osteoblasts (bone-forming cells). Importantly, Ser-3 to Ala substitution also showed a significant activity on osteoblast differentiation in vitro and increased the formation of mineralised bone matrix. Additional modifications with non-proteinogenic amino acids at position 3 improved the stability in liver microsomes, but diminished the osteoblast proliferative activity. In addition, to provide greater structural diversity, a series of macrocyclic preptin (1-16) analogues was synthesised using head-to-tail and head-to-side chain macrolactamisation as well as ring-closing metathesis. However, a detrimental effect on osteoblast activity was observed upon macrocyclisation.

Shape Transitions and Chiral Symmetry Breaking in the Energy Landscape of the Mitotic Chromosome

We derive an unbiased information theoretic energy landscape for chromosomes at metaphase using a maximum entropy approach that accurately reproduces the details of the experimentally measured pairwise contact probabilities between genomic loci. Dynamical simulations using this landscape lead to cylindrical, helically twisted structures reflecting liquid crystalline order. These structures are similar to those arising from a generic ideal homogenized chromosome energy landscape. The helical twist can be either right or left handed so chiral symmetry is broken spontaneously. The ideal chromosome landscape when augmented by interactions like those leading to topologically associating domain formation in the interphase chromosome reproduces these behaviors. The phase diagram of this landscape shows that the helical fiber order and the cylindrical shape persist at temperatures above the onset of chiral symmetry breaking, which is limited by the topologically associating domain interaction strength.

Measuring the Elasticity of Poly-l-Proline Helices with Terahertz Spectroscopy

The rigidity of poly-l-proline is an important contributor to the stability of many protein secondary structures, where it has been shown to strongly influence bulk flexibility. The experimental Young's moduli of two known poly-l-proline helical forms, right-handed all-cis (Form I) and left-handed all-trans (Form II), were determined in the crystalline state by using an approach that combines terahertz time-domain spectroscopy, X-ray diffraction, and solid-state density functional theory. Contrary to expectations, the helices were found to be considerably less rigid than many other natural and synthetic polymers, as well as differing greatly from each other, with Young's moduli of 4.9 and 9.6 GPa for Forms I and II, respectively.

"Wet" Versus "Dry" Folding of Polyproline

When the all-cis polyproline-I helix (PPI, favored in 1-propanol) of polyproline-13 is introduced into water, it folds into the all-trans polyproline-II (PPII) helix through at least six intermediates [Shi, L., Holliday, A.E., Shi, H., Zhu, F., Ewing, M.A., Russell, D.H., Clemmer, D.E.: Characterizing intermediates along the transition from PPI to PPII using ion mobility-mass spectrometry. J. Am. Chem. Soc. 136, 12702-12711 (2014)]. Here, we show that the solvent-free intermediates refold into the all-cis PPI helix with high (>90%) efficiency. Moreover, in the absence of solvent, each intermediate appears to utilize the same small set of pathways observed for the solution-phase PPII → PPI transition upon immersion of PPIIaq in 1-propanol. That folding in solution (under conditions where water is displaced by propanol) and folding in vacuo (where energy required for folding is provided by collisional activation) occur along the same pathway is remarkable. Implicit in this statement is that 1-propanol mimics a "dry" environment, similar to the gas phase. We note that intermediates with structures that are similar to PPIIaq can form PPII under the most gentle activation conditions-indicating that some transitions observed in water (i.e., "wet" folding, are accessible (albeit inefficient) in vacuo. Lastly, these "dry" folding experiments show that PPI (all cis) is favored under "dry" conditions, which underscores the role of water as the major factor promoting preference for trans proline. Graphical Abstract ᅟ.

Calculation of the Local Free Energy Landscape in the Restricted Region by the Modified Tomographic Method

The free energy landscape is the most important information in the study of the reaction mechanisms of the molecules. However, it is difficult to calculate. In a large collective variable space, a molecule must take a long time to obtain the sufficient sampling during the simulation. To save the calculation quantity, decreasing the sampling region and constructing the local free energy landscape is required in practice. However, the restricted region in the collective variable space may have an irregular shape. Simply restricting one or more collective variables of the molecule cannot satisfy the requirement. In this paper, we propose a modified tomographic method to perform the simulation. First, it divides the restricted region by some hyperplanes and connects the centers of hyperplanes together by a curve. Second, it forces the molecule to sample on the curve and the hyperplanes in the simulation and calculates the free energy data on them. Finally, all the free energy data are combined together to form the local free energy landscape. Without consideration of the area outside the restricted region, this free energy calculation can be more efficient. By this method, one can further optimize the path quickly in the collective variable space.