Conformational Profile Research Articles

Single-cell DNA methylome and 3D genome atlas of the human subcutaneous adipose tissue.

Human subcutaneous adipose tissue (SAT) contains a diverse array of cell-types; however, the epigenomic landscape among the SAT cell-types has remained elusive. Our integrative analysis of single-cell resolution DNA methylation and chromatin conformation profiles (snm3C-seq), coupled with matching RNA expression (snRNA-seq), systematically cataloged the epigenomic, 3D topology, and transcriptomic dynamics across the SAT cell-types. We discovered that the SAT CG methylation (mCG) landscape is characterized by pronounced hyper-methylation in myeloid cells and hypo-methylation in adipocytes and adipose stem and progenitor cells (ASPCs), driving nearly half of the 705,063 detected differentially methylated regions (DMRs). In addition to the enriched cell-type-specific transcription factor binding motifs, we identified TET1 and DNMT3A as plausible candidates for regulating cell-type level mCG profiles. Furthermore, we observed that global mCG profiles closely correspond to SAT lineage, which is also reflected in cell-type-specific chromosome compartmentalization. Adipocytes, in particular, display significantly more short-range chromosomal interactions, facilitating the formation of complex local 3D genomic structures that regulate downstream transcriptomic activity, including those associated with adipogenesis. Finally, we discovered that variants in cell-type level DMRs and A compartments significantly predict and are enriched for variance explained in abdominal obesity. Together, our multimodal study characterizes human SAT epigenomic landscape at the cell-type resolution and links partitioned polygenic risk of abdominal obesity to SAT epigenome.

Molecular modelling of Pro-Gly glyproline and its complexes

In the present work, Pro-Gly glyproline and its biologically active complexes with heparin and with the antithrombin-heparin protein complex have been studied using molecular modelling methods. This molecule is used as a drug and exhibits antithrombotic, anticoagulant immunomodulatory, antiulcer, and antidiabetic effects. The conformational profiles of the dipeptide were studied within the molecular mechanics framework, the geometry and energy of intra- and interresidual interactions were estimated for the most stable states of the investigated dipeptide. A comparative analysis of the extended and folded structures of this molecule optimized by the DFT/B3LYP/6-31+G(d,p) method was carried out, the frontier molecular orbitals were calculated, the surface of the molecular electrostatic potential was studied, the effective charges on the atoms were determined, and the dipole moment and polarizability values were calculated. Molecular docking revealed that the folded structure of Pro-Gly glyproline showed a higher affinity value for both heparin and the antithrombin-heparin protein complex. A model of the pharmacophore of the Pro-Gly molecule for its interaction with a specific receptor was proposed. The obtained data form the basis for the development of effective analogs of glyprolines in the complex with heparin, which are more specific concerning their putative target — antithrombin.

Bi-directional allosteric pathway in NMDA receptor activation and modulation

N-methyl-D-aspartate (NMDA) receptors are ionotropic glutamate receptors involved in learning and memory. NMDA receptors primarily comprise two GluN1 and two GluN2 subunits. The GluN2 subunit dictates biophysical receptor properties, including the extent of receptor activation and desensitization. GluN2A- and GluN2D-containing receptors represent two functional extremes. To uncover the conformational basis of their functional divergence, we utilize single-molecule fluorescence resonance energy transfer to probe the extracellular domains of these receptor subtypes under resting and ligand-bound conditions. We find that the conformational profile of the GluN2 amino-terminal domain correlates with the disparate functions of GluN2A- and GluN2D-containing receptors. Changes at the pre-transmembrane segments inversely correlate with those observed at the amino-terminal domain, confirming direct allosteric communication between these domains. Additionally, binding of a positive allosteric modulator at the transmembrane domain shifts the conformational profile of the amino-terminal domain towards the active state, revealing a bidirectional allosteric pathway between extracellular and transmembrane domains.

The role of hydrogen bonding in the conformational stability of 2-methoxyresorcinol: Insights from theoretical calculations, SERS spectroscopy, and solvent effect

Resorcinol derivatives exhibit unique structural and electronic features. The presence of two hydroxyl groups, along with other substituents, can greatly impact these molecules’ conformational profile and vibrational characteristics. This work investigates molecular and electronic properties of 2-methoxyresorcinol (2-MR) by combining quantum-chemical and spectroscopic techniques. The molecule of 2-MR can exist in the syn, anti, and anti–syn configurations, depending on the orientation of the hydroxyl groups. The syn and anti–syn forms have been verified from DFT, MP2 and MP4 (SDQ) calculations and from infrared/Raman spectroscopy to be stabilized by intramolecular hydrogen bonding that takes place between the hydrogen atoms of the OH groups and the methoxy moiety’s oxygen atom. Such an interaction results in 6–7 kcal/mol more stable states than the corresponding anti-configuration, where only steric forces control the molecule stability. The effect of different solvents on the Raman features has also been accessed. Unlike other solutions, aqueous 2-MR facilitated a split in the OH stretching peak, which can be attributed to the pronounced solute–solvent intermolecular forces through the OH groups. Surface-enhanced Raman scattering (SERS) spectroscopic analysis of trace 2-MR utilizing as-synthesized silver nanoparticles shows a remarkable boost in the Raman signal intensities of several vibrational modes. The SERS detection exhibits a wide dynamic range and a detection limit of 1 × 10−8 M.

Conformational Space Profiling Enhances Generic Molecular Representation for AI-Powered Ligand-Based Drug Discovery.

The molecular representation model is a neural network that converts molecular representations (SMILES, Graph) into feature vectors, and is an essential module applied across a wide range of artificial intelligence-driven drug discovery scenarios. However, current molecular representation models rarely consider the three-dimensional conformational space of molecules, losing sight of the dynamic nature of small molecules as well as the essence of molecular conformational space that covers the heterogeneity of molecule properties, such as the multi-target mechanism of action, recognition of different biomolecules, dynamics in cytoplasm and membrane. In this study, a new model named GeminiMol is proposed to incorporate conformational space profiles into molecular representation learning, which extracts the feature of capturing the complicated interplay between the molecular structure and the conformational space. Although GeminiMol is pre-trained on a relatively small-scale molecular dataset (39290 molecules), it shows balanced and superior performance not only on 67 molecular properties predictions but also on 73 cellular activity predictions and 171 zero-shot tasks (including virtual screening and target identification). By capturing the molecular conformational space profile, the strategy paves the way for rapid exploration of chemical space and facilitates changing paradigms for drug design.

Mechanistic conformational and substrate selectivity profiles emerging in the evolution of enzymes via parallel trajectories

Laboratory evolution studies have demonstrated that parallel evolutionary trajectories can lead to genetically distinct enzymes with high activity towards a non-preferred substrate. However, it is unknown whether such enzymes have convergent conformational dynamics and mechanistic features. To address this question, we use as a model the wild-type Homo sapiens kynureninase (HsKYNase), which is of great interest for cancer immunotherapy. Earlier, we isolated HsKYNase_66 through an unusual evolutionary trajectory, having a 410-fold increase in the kcat/KM for kynurenine (KYN) and reverse substrate selectivity relative to HsKYNase. Here, by following a different evolutionary trajectory we generate a genetically distinct variant, HsKYNase_93D9, that exhibits KYN catalytic activity comparable to that of HsKYNase_66, but instead it is a “generalist” that accepts 3’-hydroxykynurenine (OH-KYN) with the same proficiency. Pre-steady-state kinetic analysis reveals that while the evolution of HsKYNase_66 is accompanied by a change in the rate-determining step of the reactions, HsKYNase_93D9 retains the same catalytic mechanism as HsKYNase. HDX-MS shows that the conformational dynamics of the two enzymes are markedly different and distinct from ortholog prokaryotic enzymes with high KYN activity. Our work provides a mechanistic framework for understanding the relationship between evolutionary mechanisms and phenotypic traits of evolved generalist and specialist enzyme species.

Marine fungal diversity unlocks potent antivirals against monkeypox through methyltransferase inhibition revealed by molecular dynamics and free energy landscape

The escalating threat posed by the Monkeypox virus (MPXV) to global health necessitates the urgent discovery of effective antiviral agents, as there are currently no specific drugs available for its treatment, and existing inhibitors are hindered by toxicity and poor pharmacokinetic profiles. This study aimed to identify potent MPXV inhibitors by screening a diverse library of small molecule compounds derived from marine fungi, focusing on the viral protein VP39, a key methyltransferase involved in viral replication. An extensive virtual screening process identified four promising compounds—CMNPD15724, CMNPD28811, CMNPD30883, and CMNPD18569—alongside a control molecule. Rigorous evaluations, including re-docking, molecular dynamics (MD) simulations, and hydrogen bond analysis, were conducted to assess their inhibitory potential against MPXV VP39. CMNPD15724 and CMNPD30883, in particular, demonstrated a superior binding affinity and stable interactions within the target protein's active site throughout the MD simulations, suggesting a capacity to overcome the limitations associated with sinefungin. The stability of these VP39-compound complexes, corroborated by MD simulations, provided crucial insights into the dynamic behavior of these interactions. Furthermore, Principal Component Analysis (PCA) based free energy landscape assessments offered a detailed understanding of the dynamic conformational changes and energetic profiles underlying these compounds' functional disruption of VP39. These findings establish CMNPD15724, CMNPD28811, CMNPD30883, and CMNPD18569 as promising MPXV inhibitors and highlight marine fungi as a valuable source of novel antiviral agents. These compounds represent potential candidates for further experimental validation, advancing the development of safer and more effective therapeutic options to combat this emerging viral infection.

Structural Mass Spectrometry Captures Residue-Resolved Comprehensive Conformational Rearrangements of a G Protein-Coupled Receptor.

G protein-coupled receptor (GPCR) structural studies with in-solution spectroscopic approaches have offered distinctive insights into GPCR activation and signaling that highly complement those yielded from structural snapshots by crystallography or cryo-EM. While most current spectroscopic approaches allow for probing structural changes at selected residues or loop regions, they are not suitable for capturing a holistic view of GPCR conformational rearrangements across multiple domains. Herein, we develop an approach based on limited proteolysis mass spectrometry (LiP-MS) to simultaneously monitor conformational alterations of a large number of residues spanning both flexible loops and structured transmembrane domains for a given GPCR. To benchmark LiP-MS for GPCR conformational profiling, we studied the adenosine 2A receptor (A2AR) in response to different ligand binding (agonist/antagonist/allosteric modulators) and G protein coupling. Systematic and residue-resolved profiling of A2AR conformational rearrangements by LiP-MS precisely captures structural mechanisms in multiple domains underlying ligand engagement, receptor activation, and allostery, and may also reflect local conformational flexibility. Furthermore, these residue-resolution structural fingerprints of the A2AR protein allow us to readily classify ligands of different pharmacology and distinguish the G protein-coupled state. Thus, our study provides a new structural MS approach that would be generalizable to characterizing conformational transition and plasticity for challenging integral membrane proteins.

Conformational Analysis of 1,3-Difluorinated Alkanes.

Fluorine substitution can have a profound impact on molecular conformation. Here, we present a detailed conformational analysis of how the 1,3-difluoropropylene motif (-CHF-CH2-CHF-) determines the conformational profiles of 1,3-difluoropropane, anti- and syn-2,4-difluoropentane, and anti- and syn-3,5-difluoroheptane. It is shown that the 1,3-difluoropropylene motif strongly influences alkane chain conformation, with a significant dependence on the polarity of the medium. The conformational effect of 1,3-fluorination is magnified upon chain extension, which contrasts with vicinal difluorination. Experimental evidence was obtained from NMR analysis, where polynomial complexity scaling simulation algorithms were necessary to enable J-coupling extraction from the strong second-order spectra, particularly for the large 16-spin systems of the difluorinated heptanes. These results improve our understanding of the conformational control toolkit for aliphatic chains, yield simple rules for conformation population analysis, and demonstrate quantum mechanical time-domain NMR simulations for liquid state systems with large numbers of strongly coupled spins.

A journey in unraveling the enantiorecognition mechanism of 3,5-dinitrobenzoyl-amino acids with two Cinchona alkaloid-based chiral stationary phases: The power of molecular dynamic simulations

BackgroundInnovations in computer hardware and software capabilities have paved the way for advances in molecular modelling techniques and methods, leading to an unprecedented expansion of their potential applications. In contrast to the docking technique, which usually identifies the most stable selector-selectand (SO-SA) complex for each enantiomer, the molecular dynamics (MD) technique enables the consideration of a distribution of the SO-SA complexes based on their energy profile. This approach provides a more truthful representation of the processes occurring within the column. However, benchmark procedures and focused guidelines for computational treatment of enantioselectivity at the molecular level are still missing. ResultsTwenty-eight molecular dynamics simulations were performed to study the enantiorecognition mechanisms of seven N-3,5-dinitrobenzoylated α- and β-amino acids (DNB-AAs), occurring with the two quinine- and quinidine-based (QN-AX and QD-AX) chiral stationary phases (CSPs), under polar-ionic conditions. The MD protocol was optimized in terms of box size, simulation run time, and frame recording frequency. Subsequently, all the trajectories were analyzed by calculating both the type and amount of the interactions engaged by the selectands (SAs) with the two chiral selectors (SOs), as well as the conformational and interaction energy profiles of the formed SA-SO associates. All the MDs were in strict agreement with the experimental enantiomeric elution order and allowed to establish (i) that salt-bridge and H-bond interactions play a pivotal role in the enantiorecognition mechanisms, and (ii) that the π-cation and π-π interactions are the discriminant chemical features between the two SOs in ruling the chiral recognition mechanism. SignificanceThe results of this work clearly demonstrate the high contribution given by MD simulations in the comprehension of the enantiorecognition mechanism with Cinchona alkaloid-based CSPs. However, from this research endeavor it clearly emerged that the MD protocol optimization is crucial for the quality of the produced results.

From Functional Fatty Acids to Potent and Selective Natural-Product-Inspired Mimetics via Conformational Profiling.

Fatty acids play important signaling roles in biology, albeit typically lacking potency or selectivity, due to their substantial conformational flexibility. While being recognized as having properties of potentially great value as therapeutics, it is often the case that the functionally relevant conformation of the natural fatty acid is not known, thereby complicating efforts to develop natural-product-inspired ligands that have similar functional properties along with enhanced potency and selectivity profiles. In other words, without structural information associated with a particular functional relationship and the hopelessly unbiased conformational preferences of the endogenous ligand, one is molecularly ill-informed regarding the precise ligand-receptor interactions that play a role in driving the biological activity of interest. To address this problem, a molecular strategy to query the relevance of distinct subpopulations of fatty acid conformers has been established through "conformational profiling", a process whereby a unique collection of chiral and conformationally constrained fatty acids is employed to deconvolute beneficial structural features that impart natural-product-inspired function. Using oleic acid as an example because it is known to engage a variety of receptors, including GPR40, GPR120, and TLX, a 24-membered collection of mimetics was designed and synthesized. It was then demonstrated that this collection contained members that have enhanced potency and selectivity profiles, with some being clearly biased for engagement of the GPCRs GPR40 and GPR120 while others were identified as potent and selective modulators of the nuclear receptor TLX. A chemical synthesis strategy that exploited the power of modern technology for stereoselective synthesis was critical to achieving success, establishing a common sequence of bond-forming reactions to access a disparate collection of chiral mimetics, whose conformational preferences are impacted by the nature of stereodefined moieties differentially positioned about the C18 skeleton of the parent fatty acid. Overall, this study establishes a foundation to fuel future programs aimed at developing natural-product-inspired fatty acid mimetics as valuable tools in chemical biology and potential therapeutic leads.

4D-QSAR and MIA-QSAR Studies of Aminobenzimidazole Derivatives as Fourth-generation EGFR Inhibitors.

The epidermal growth factor receptor (EGFR) protein has been intensively studied as a therapeutic target for non-small cell lung cancer (NSCLC). The aminobenzimidazole derivatives as the fourth-generation EGFR inhibitors have achieved promising results and overcame EGFR mutations at C797S, del19 and T790M in NSCLC. In order to understand the quantitative structure-activity relationship (QSAR) of aminobenzimidazole derivatives as EGFRdel19 T790M C797S inhibitors, the four-dimensional QSAR (4D-QSAR) and multivariate image analysis (MIA-QSAR) have been performed on the data of 45 known aminobenzimidazole derivatives. The 4D-QSAR descriptors were acquired by calculating the association energies between probes and aligned conformational ensemble profiles (CEP), and the regression models were established by partial least squares (PLS). In order to further understand and verify the 4D-QSAR model, MIA-QSAR was constructed by using chemical structure pictures to generate descriptors and PLS regression. Furthermore, the molecular docking and averaged noncovalent interactions (aNCI) analysis were also performed to further understand the interactions between ligands and the EGFR targets, which was in good agreement with the 4D-QSAR model. The established 4D-QSAR and MIA-QSAR models have strong stability and good external prediction ability. These results will provide theoretical guidance for the research and development of aminobenzimidazole derivatives as new EGFRdel19 T790M C797S inhibitors.

Conformational fingerprinting with Raman spectroscopy reveals protein structure as a translational biomarker of muscle pathology.

Neuromuscular disorders are a group of conditions that can result in weakness of skeletal muscles. Examples include fatal diseases such as amyotrophic lateral sclerosis and conditions associated with high morbidity such as myopathies (muscle diseases). Many of these disorders are known to have abnormal protein folding and protein aggregates. Thus, easy to apply methods for the detection of such changes may prove useful diagnostic biomarkers. Raman spectroscopy has shown early promise in the detection of muscle pathology in neuromuscular disorders and is well suited to characterising the conformational profiles relating to protein secondary structure. In this work, we assess if Raman spectroscopy can detect differences in protein structure in muscle in the setting of neuromuscular disease. We utilise in vivo Raman spectroscopy measurements from preclinical models of amyotrophic lateral sclerosis and the myopathy Duchenne muscular dystrophy, together with ex vivo measurements of human muscle samples from individuals with and without myopathy. Using quantitative conformation profiling and matrix factorisation we demonstrate that quantitative 'conformational fingerprinting' can be used to identify changes in protein folding in muscle. Notably, myopathic conditions in both preclinical models and human samples manifested a significant reduction in α-helix structures, with concomitant increases in β-sheet and, to a lesser extent, nonregular configurations. Spectral patterns derived through non-negative matrix factorisation were able to identify myopathy with a high accuracy (79% in mouse, 78% in human tissue). This work demonstrates the potential of conformational fingerprinting as an interpretable biomarker for neuromuscular disorders.

New model for aspartic acid species in aqueous calcium carbonate growth environments: challenges and perspectives.

The lack of experimental data on the dynamics of aspartic acid species in water for its range of protonation states and the details of their atomic-level interaction with aqueous calcium carbonate species is a driver for accurate force field development. A classical model that is consistent with the few pieces of experimental data available and with first principles calculations has been developed. The complex dynamics of the aspartate anions relevant to biomineralization and calcium carbonate crystal growth has been explored in water, providing a quantitative description of solvation structure and free energies, including conformational free energy profiles and pairing free energies. The model has been used to probe the structure and dynamics of aqueous calcium aspartate homo- and hetero-chiral clusters, confirming their unlikelihood due to weak and water-mediated interactions. This supports the hypothesis that the formation of such clusters, observed while growing vaterite in the presence of acidic chiral amino acids, is favoured by the presence of the crystal surface.

Conserved and divergent gene regulatory programs of the mammalian neocortex

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Metasurface-based conformal linear-to-circular polarization converter for emerging wireless applications

This paper proposes the analysis and design aspects of a metasurface-based reflective type linear-to-circular polarization (LCP) converter in both planar and conformal profiles. The metasurface unit cell consists of a diagonally placed split circular metallic ring on top of a metal backed dielectric substrate. The structure provides LCP conversion with axial ratio (AR) < 3 dB at two frequency bands, which are 27.38–28.61GHz and 33.1–52GHz corresponding to fractional bandwidths of 4.4% and 44.5%, respectively. It exhibits two different handedness at the aforementioned two bands for a TE (or TM) linearly polarized incident wave. Moreover, it shows dual polarized and angularly stable (up to 25°) characteristics. It is observed from the far-field characteristics that the LCP converter with an array of 29 × 29 unit cells (12λ×12λ) produces peak gain of 26dBi, cross polarization level 20 dB below the maximum main lobe level when illuminated by a linearly polarized source at far-field. Subsequently, the same unit cell is used to design a cylindrical conformal metasurface with different radius of curvatures (r). It provides quite stable AR, gain and radiation patterns till r=100mm. The array has been fabricated to experimentally validate the numerical results of both the planar and conformal geometries.

Multi-omic analysis of human kidney tissue identified medulla-specific gene expression patterns

The kidney medulla is a specialized region with important homeostatic functions. It has been implicated in genetic and developmental disorders along with ischemic and drug-induced injuries. Despite its role in kidney function and disease, the medulla's baseline gene expression and epigenomic signatures have not been well described in the adult human kidney. Here we generated and analyzed gene expression (RNA-seq), chromatin accessibility (ATAC-seq), chromatin conformation (Hi-C) and spatial transcriptomic data from the adult human kidney cortex and medulla. Tissue samples were obtained from macroscopically dissected cortex and medulla of tumor-adjacent normal material in nephrectomy specimens from five male patients. We used these carefully annotated specimens to reassign incorrectly labeled samples in the larger public Genotype-Tissue Expression (GTEx) Project, and to extract meaningful medullary gene expression signatures. Using integrated analysis of gene expression, chromatin accessibility and conformation profiles, we found insights into medulla development and function and then validated this by spatial transcriptomics and immunohistochemistry. Thus, our datasets provide a valuable resource for functional annotation of variants from genome-wide association studies and are freely accessible through an epigenome browser portal.

Synthesis And Conformational Behaviors of Unnatural Peptides Alternating Chiral And Achiral α,α-Disubstituted α-Amino Acid Units.

The development of peptidomimetics to modulate the conformational profile of peptides has been extensively studied in the fields of biological and medicinal chemistry. However, large-scale synthesis of peptidomimetics with both an ordered sequence and a controlled secondary structure is highly challenging. In this paper, the framework of peptidomimetics has been designed to be alternating an achiral α,α-disubstituted α-amino acid unit and a chiral α-methylphenylalanine unit. The polymers were synthesized via our invented Ugi reaction-based polycondensation technique. The chiral higher-order structures of the alternating peptides were evaluated mainly through CD spectroscopy. The UV-vis and CD spectra of the polymers in three solvents were systematically measured at various temperatures. The anisotropic factors of CD (gCD ) values were calculated to know the chiroptical response. The results indicate the characteristic conformational behaviors. In a polar solvent, the hydrogen bonds between the N-H group of MePhe unit and the C = O of α,α-diphenylglycine unit outweigh the intraresidue hydrogen bonds in α,α-diphenylglycine unit, leading to the formation of a prevailing preferred-handed 310 -helical conformation. On the other hand, in a less polar solvent, the intrachain hydrogen bonds switch to intraresidue hydrogen bonds in α,α-diphenylglycine unit, which make the polymer adopting a prevailing extended planar C5 -conformation. This article is protected by copyright. All rights reserved.

Insight into the Structure of Victorin, the Host-Selective Toxin from the Oat Pathogen Cochliobolus victoriae. Studies of the Unique Dehydroamino Acid β-Chlorodehydroalanine.

Victorins, a family of peptide toxins, produced by the fungal pathogen Cochliobolus victoriae and responsible for disease of some oat varieties, contain a β-chlorodehydroalanine residue, ΔAla(βCl). To determine the conformational properties of this unique dehydroamino acid, a series of model compounds was studied using X-ray, NMR, and FT-IR methods, supported by theoretical calculations. The ΔAla(βCl) geometrical isomers differ in conformational profile. The isomer Z prefers the helical conformation α (φ, ψ = -61°, -24°), PPII type conformation β (φ, ψ = -47°, 136°), and semiextended conformation β2 (φ, ψ = -116°, 9°) in weakly and more polar solutions. The isomer E prefers mainly the extended conformation C5 (φ, ψ = -177°, 160°), but with an increase of the environment polarity also conformations β (φ, ψ = -44°, 132°) and α (φ, ψ = -53°, -39°). In the most stable conformations the N-H···Cl hydrogen bond (5γ) occurs, created between the chlorine atom of the side chain and the N-H donor of the flanking amide group. The method of synthesis of the β-chlorodehydroalanine residue is proposed, by chlorination of dehydroalanine and then the photoisomerization from the isomer Z to E. The presented results indicate that the assignment of the geometrical isomer of the ΔAla(βCl) residue in naturally occurring victorins still remains an open question, despite being crucial for biological activity.

Tackling Hysteresis in Conformational Sampling: How to Be Forgetful with MEMENTO.

The structure of proteins has long been recognized to hold the key to understanding and engineering their function, and rapid advances in structural biology and protein structure prediction are now supplying researchers with an ever-increasing wealth of structural information. Most of the time, however, structures can only be determined in free energy minima, one at a time. While conformational flexibility may thus be inferred from static end-state structures, their interconversion mechanisms─a central ambition of structural biology─are often beyond the scope of direct experimentation. Given the dynamical nature of the processes in question, many studies have attempted to explore conformational transitions using molecular dynamics (MD). However, ensuring proper convergence and reversibility in the predicted transitions is extremely challenging. In particular, a commonly used technique to map out a path from a starting to a target conformation called steered MD (SMD) can suffer from starting-state dependence (hysteresis) when combined with techniques such as umbrella sampling (US) to compute the free energy profile of a transition. Here, we study this problem in detail on conformational changes of increasing complexity. We also present a new, history-independent approach that we term "MEMENTO" (Morphing End states by Modelling Ensembles with iNdependent TOpologies) to generate paths that alleviate hysteresis in the construction of conformational free energy profiles. MEMENTO utilizes template-based structure modelling to restore physically reasonable protein conformations based on coordinate interpolation (morphing) as an ensemble of plausible intermediates, from which a smooth path is picked. We compare SMD and MEMENTO on well-characterized test cases (the toy peptide deca-alanine and the enzyme adenylate kinase) before discussing its use in more complicated systems (the kinase P38α and the bacterial leucine transporter LeuT). Our work shows that for all but the simplest systems SMD paths should not in general be used to seed umbrella sampling or related techniques, unless the paths are validated by consistent results from biased runs in opposite directions. MEMENTO, on the other hand, performs well as a flexible tool to generate intermediate structures for umbrella sampling. We also demonstrate that extended end-state sampling combined with MEMENTO can aid the discovery of collective variables on a case-by-case basis.