Experimentally Derived and Computationally Optimized Backbone Conformational Statistics for Blocked Amino Acids

Refined Bonded Terms in Coarse-Grained Models for Intrinsically Disordered Proteins Improve Backbone Conformations

Coarse-grained models designed for intrinsically disordered proteins and regions (IDP/Rs) usually omit some bonded potentials (e.g., angular and dihedral potentials) as a conventional strategy to enhance backbone flexibility. However, a notable drawback of this approach is the generation of inaccurate backbone conformations. Here, we addressed this problem by introducing residue-specific angular, refined dihedral, and correction map (CMAP) potentials, derived based on the statistics from a customized coil database. These bonded potentials were integrated into the existing Mpipi model, resulting in a new model, denoted as the "Mpipi+" model. Results show that the Mpipi+ model can improve backbone conformations. More importantly, it can markedly improve the secondary structure propensity (SSP) based on the experimental chemical shift and, consequently, succeed in capturing transient secondary structures. Moreover, the Mpipi+ model preserves the liquid-liquid phase separation (LLPS) propensities of IDPs.

Chitosan Polysaccharides from a Polarizable Multiscale Approach

We report simulations of chitosan polysaccharides in the aqueous phase, at infinite dilute conditions and zero ionic strength. Those simulations are performed by means of a polarizable multiscale modeling scheme that relies on a polarizable all atom force field to model solutes and on a polarizable solvent coarse grained approach. Force field parameters are assigned only from quantum chemistry ab initio data. We simulate chitosan monomer units, dimers and 50-long chains. Regarding the 50-long chains we simulate three sets of ten randomly built chain replica at three different pH conditions (corresponding to different chain protonation states, the chain degree of deacetylation is 85%). Our simulations show the persistence length of 50-long chitosan chains at strong acidic conditions (pH <5) to be 24 ± 2 nm (at weak/negligible ionic strength conditions), and to be 1 order of magnitude shorter at usual pH conditions. Our simulation data support the most recent simulation and experimental studies devoted to chitosan polysaccharides in the aqueous phase.

Conformational Manifold Sampled by Two Short Linear Motif Segments Probed by Circular Dichroism, Vibrational, and Nuclear Magnetic Resonance Spectroscopy

Disordered protein segments called short linear motifs (SLiM) serve as recognition sites for a variety of biological processes and act as targeting signals, modification, and ligand binding sites. While SLiMs do not adopt one of the known regular secondary structures, the conformational distribution might still reflect the structural propensities of their amino acid residues and possible interactions between them. In the past, conformational analyses of short peptides provided compelling evidence for the notion that individual residues are less conformationally flexible than locally expected for a random coil. Here, we combined various spectroscopies (NMR, IR, vibrational, and UV circular dichroism) to determine the Ramachandran plots of two SLiM motifs, i.e., GRRDSG and GRRTSG. They are two representatives of RxxS motifs that are capable of being phosphorylated by protein kinase A, an enzyme that plays a fundamental role in a variety of biological processes. Our results reveal that the nearest and non-nearest interactions between residues cause redistributions between polyproline II and β-strand basins while concomitantly stabilizing extended relative to turn-forming and helical structures. They also cause shifts in basin positions. With increasing temperature, β-strand populations become more populated at the expense of polyproline II. While molecular dynamics simulations with Amber ff14SB and CHARMM 36m force fields indicate residue-residue interactions, they do not account for the observed structural changes.

Phase Transitions of Associative Biomacromolecules

Multivalent proteins and nucleic acids, collectively referred to as multivalent associative biomacromolecules, provide the driving forces for the formation and compositional regulation of biomolecular condensates. Here, we review the key concepts of phase transitions of aqueous solutions of associative biomacromolecules, specifically proteins that include folded domains and intrinsically disordered regions. The phase transitions of these systems come under the rubric of coupled associative and segregative transitions. The concepts underlying these processes are presented, and their relevance to biomolecular condensates is discussed.

Conformational landscapes of artificial peptides predicted by various force fields: are we ready to simulate β-amino acids?

With the introduction of artificial peptides as antimicrobial agents and organic catalysts, numerous efforts have been made to design foldamers with desirable structures and functions. Computational tools are a helpful proxy for revealing the dynamic structures at atomic resolution and understanding foldamer's complex structure-function relationships. However, the performance of conventional force fields in predicting the structures of artificial peptides has not been systematically evaluated. In this study, we critically assessed three popular force fields, AMBER ff14SB, CHARMM36m, and OPLS-AA/L, in predicting conformational propensities of a β-peptide foldamer at monomer and hexamer levels. Simulation results were compared to those obtained from quantum chemistry calculations and experimental data. We also utilised replica exchange molecular dynamics simulations to investigate the energy landscape of each force field and assess the similarities and differences between force fields. We compared different solvent systems in the AMBER ff14SB and CHARMM36m frameworks and confirmed the unanimous role of hydrogen bonds in shaping energy landscapes. We anticipate that our data will pave the way for further improvements to force fields and for understanding the role of solvents in peptide folding, crystallisation, and engineering.

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)─an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration─physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

AlphaFold and Implications for Intrinsically Disordered Proteins

Accurate predictions of the three-dimensional structures of proteins from their amino acid sequences have come of age. AlphaFold, a deep learning-based approach to protein structure prediction, shows remarkable success in independent assessments of prediction accuracy. A significant epoch in structural bioinformatics was the structural annotation of over 98% of protein sequences in the human proteome. Interestingly, many predictions feature regions of very low confidence, and these regions largely overlap with intrinsically disordered regions (IDRs). That over 30% of regions within the proteome are disordered is congruent with estimates that have been made over the past two decades, as intense efforts have been undertaken to generalize the structure-function paradigm to include the importance of conformational heterogeneity and dynamics. With structural annotations from AlphaFold in hand, there is the temptation to draw inferences regarding the "structures" of IDRs and their interactomes. Here, we offer a cautionary note regarding the misinterpretations that might ensue and highlight efforts that provide concrete understanding of sequence-ensemble-function relationships of IDRs. This perspective is intended to emphasize the importance of IDRs in sequence-function relationships (SERs) and to highlight how one might go about extracting quantitative SERs to make sense of how IDRs function.

Controlled "off-on" fluorescent probe for the specific detection of hyperhomocysteinemia

Hyperhomocysteinemia is an established risk factor for atherosclerosis and vascular disease. Therefore, designing a hyperhomocysteinemia specific probe is of great significance for the early warning of cardiovascular diseases. However, developing probes that can efficiently and specifically recognize homocysteine (Hcy) remains a tremendous challenge. Therefore, we designed an Hcy-specific fluorescent probe (HSFP) with excellent selectivity and anti-interference capability. Interestingly, this probe can automatically “off–on” in water solution, but the fluorescence of HSFP remains “off” when Hcy is present in the solution. The spectroscopic data demonstrated that the fluorescence of HSFP attenuated 13.8 folds toward Hcy in water without interference from other biothiols and amino acids. Furthermore, HSFP can sensitively reflect the change of Hcy content in cells. Therefore, HSFP was further applied to detect hyperhomocysteinemia in vivo with high efficiency. In summary, we have developed an Hcy-specific fluorescent probe to efficiently detect Hcy in vivo and in vitro, which may contribute to basic or clinical research.

An Hcy-specific fluorescent probe (HSFP) with excellent selectivity and anti-interference capability was developed for the detection of hyperhomocysteinemia.

Artificial intelligence techniques for integrative structural biology of intrinsically disordered proteins

We outline recent developments in artificial intelligence (AI) and machine learning (ML) techniques for integrative structural biology of intrinsically disordered proteins (IDP) ensembles. IDPs challenge the traditional protein structure-function paradigm by adapting their conformations in response to specific binding partners leading them to mediate diverse, and often complex cellular functions such as biological signaling, self-organization and compartmentalization. Obtaining mechanistic insights into their function can therefore be challenging for traditional structural determination techniques. Often, scientists have to rely on piecemeal evidence drawn from diverse experimental techniques to characterize their functional mechanisms. Multiscale simulations can help bridge critical knowledge gaps about IDP structure-function relationships-however, these techniques also face challenges in resolving emergent phenomena within IDP conformational ensembles. We posit that scalable statistical inference techniques can effectively integrate information gleaned from multiple experimental techniques as well as from simulations, thus providing access to atomistic details of these emergent phenomena.

q-Canonical Monte Carlo Sampling for Modeling the Linkage between Charge Regulation and Conformational Equilibria of Peptides

The overall charge content and the patterning of charged residues have a profound impact on the conformational ensembles adopted by intrinsically disordered proteins. These parameters can be altered by charge regulation, which refers to the effects of post-translational modifications, pH-dependent changes to charge, and conformational fluctuations that modify the pKa values of ionizable residues. Although atomistic simulations have played a prominent role in uncovering the major sequence-ensemble relationships of IDPs, most simulations assume fixed charge states for ionizable residues. This may lead to erroneous estimates for conformational equilibria if they are linked to charge regulation. Here, we report the development of a new method we term q-canonical Monte Carlo sampling for modeling the linkage between charge regulation and conformational equilibria. The method, which is designed to be interoperable with the ABSINTH implicit solvation model, operates as follows: For a protein sequence with n ionizable residues, we start with all 2n charge microstates and use a criterion based on model compound pKa values to prune down to a subset of thermodynamically relevant charge microstates. This subset is then grouped into mesostates, where all microstates that belong to a mesostate have the same net charge. Conformational distributions, drawn from a canonical ensemble, are generated for each of the charge microstates that make up a mesostate using a method we designate as proton walk sampling. This method combines Metropolis Monte Carlo sampling in conformational space with an auxiliary Markov process that enables interconversions between charge microstates along a mesostate. Proton walk sampling helps identify the most likely charge microstate per mesostate. We then use thermodynamic integration aided by the multistate Bennett acceptance ratio method to estimate the free energies for converting between mesostates. These free energies are then combined with the per-microstate weights along each mesostate to estimate standard state free energies and pH-dependent free energies for all thermodynamically relevant charge microstates. The results provide quantitative estimates of the probabilities and preferred conformations associated with every thermodynamically accessible charge microstate. We showcase the application of q-canonical sampling using two model systems. The results establish the soundness of the method and the importance of charge regulation in systems characterized by conformational heterogeneity.

Improvements to the ABSINTH Force Field for Proteins Based on Experimentally Derived Amino Acid Specific Backbone Conformational Statistics

We present an improved version of the ABSINTH implicit solvation model and force field paradigm (termed ABSINTH-C) by incorporating a grid-based term that bootstraps against experimentally derived and computationally optimized conformational statistics for blocked amino acids. These statistics provide high-resolution descriptions of the intrinsic backbone dihedral angle preferences for all 20 amino acids. The original ABSINTH model generates Ramachandran plots that are too shallow in terms of the basin structures and too permissive in terms of dihedral angle preferences. We bootstrap against the reference optimized landscapes and incorporate CMAP-like residue-specific terms that help us reproduce the intrinsic dihedral angle preferences of individual amino acids. These corrections that lead to ABSINTH-C are achieved by balancing the incorporation of the new residue-specific terms with the accuracies inherent to the original ABSINTH model. We demonstrate the robustness of ABSINTH-C through a series of examples to highlight the preservation of accuracies as well as examples that demonstrate the improvements. Our efforts show how the recent experimentally derived and computationally optimized coil-library landscapes can be used as a touchstone for quantifying errors and making improvements to molecular mechanics force fields.