Systematic validation of protein force fields against experimental data

Differentially heterogeneous hydration environment of the familial mutants of α-synuclein

The behavior of hydration water around familial Parkinson's disease linked mutants of α-synuclein may be linked to the early-onset of Parkinson's disease. For the first time, this study compares the local structure and dynamics of hydration water around different segments of some of the natural mutants of α-synuclein, i.e., E46K, G51D, A30P, and A53E, with that of the wild-type protein through explicit water MD simulations. The results show that the C-terminal segments of the fast aggregating mutants such as E46K and A30P are less exposed to water, while those of the slow aggregating ones such as A53E and G51D are more exposed to water relative to that of the wild-type protein. In addition, the water molecules are found to be more ordered around the C-terminal segment of the A53E and G51D mutants as compared to the wild-type protein. This is due to an increase in the overall charge of α-syn upon A53E and G51D mutations. The translational and rotational motions of water molecules in the hydration shell of the C-terminal segment of A53E and G51D mutants are found to be faster relative to that of the wild-type protein. This study validates the differential hydration environment around the C-terminal segment for the causative and protective mutants of α-synuclein.

Toward Gaussian Process Regression Modeling of a Urea Force Field

FFLUX is a next-generation, machine-learnt force field built on three cornerstones: quantum chemical topology, Gaussian process regression, and (high-rank) multipolar electrostatics. It is capable of performing molecular dynamics with near-quantum accuracy at a lower computational cost than standard ab initio molecular dynamics. Previous work with FFLUX was concerned with water and formamide. In this study, we go one step further and challenge FFLUX to model urea, a larger and more flexible system. In result, we have trained urea models at the B3LYP/aug-cc-pVTZ level of theory, with a mean absolute error of 0.4 kJ mol-1 and a maximum prediction error below 7.0 kJ mol-1. To test their performance in molecular dynamics simulations, two sets of FFLUX geometry optimizations were carried out: 5 dimers corresponding to energy minima and 75 random dimers. The 5 dimers were recovered with a root-mean-square deviation below 0.1 Å with respect to their ab initio references. Out of the 75 random dimers, 68% converged to the qualitatively same dimer as those obtained at the ab initio level. Furthermore, we have ranked the 5 FFLUX-optimized dimers in the order of their relative FFLUX single-point energies and compared them with the ab initio method. The energy ranking fully agreed but for one crossover between two successive minima. Finally, we have demonstrated the importance of geometry-dependent (i.e., flexible) multipole moments, showing that the lack of multipole moment flexibility can lead to average errors in the total intermolecular electrostatic energy of more than 2 orders of magnitude.

Identifying the minimal sets of distance restraints for FRET-assisted protein structural modeling

Proteins naturally occur in crowded cellular environments and interact with other proteins, nucleic acids, and organelles. Since most previous experimental protein structure determination techniques require that proteins occur in idealized, non-physiological environments, the effects of realistic cellular environments on protein structure are largely unexplored. Recently, Förster resonance energy transfer (FRET) has been shown to be an effective experimental method for investigating protein structure in vivo. Inter-residue distances measured in vivo can be incorporated as restraints in molecular dynamics (MD) simulations to model protein structural dynamics in vivo. Since most FRET studies only obtain inter-residue separations for a small number of amino acid pairs, it is important to determine the minimum number of restraints in the MD simulations that are required to achieve a given root-mean-square deviation (RMSD) from the experimental structural ensemble. Further, what is the optimal method for selecting these inter-residue restraints? Here, we implement several methods for selecting the most important FRET pairs and determine the number of pairsN r that are needed to induce conformational changes in proteins between two experimentally determined structures. We find that enforcing only a small fraction of restraints,N r / N ≲ 0.08 , where N is the number of amino acids, can induce the conformational changes. These results establish the efficacy of FRET-assisted MD simulations for atomic scale structural modeling of proteins in vivo.

In silico study on graphene quantum dots modified with various functional groups inhibiting α‑synuclein dimerization

HYPOTHESIS

Graphene quantum dots (GQDs) with various functional groups are hypothesized to inhibit the α-synuclein (αS) dimerization, a crucial step in Parkinson's disease pathogenesis. The potential of differently functionalized GQDs is systematically explored.

EXPERIMENTS

All-atom replica-exchange molecular dynamics simulations (accumulating to 75.6 μs) in explicit water were performed to study the dimerization of the αS non-amyloid component region and the influence of GQDs modified with various functional groups. Conformation ensemble, binding behavior, and free energy analysis were conducted.

FINDINGS

All studied GQDs inhibit β-sheet and backbone hydrogen bond formation in αS dimers, leading to looser oligomeric conformations. Charged GQDs severely impede the growth of extended β-sheets by providing extra contact surface. GQD binding primarily disrupts αS inter-peptide interactions through π-π stacking, CH-π interactions, and for charged GQDs, additionally through salt-bridge and hydrogen bonding interactions. GQD-COO- showed the most optimal inhibitory effect, binding mode, and intensity, which holds promise for the development of nanomedicines targeting amyloid aggregation in neurodegenerative diseases.

Tips and Tricks in the Modeling of Supramolecular Peptide Assemblies

Supramolecular peptide assemblies (SPAs) hold promise as materials for nanotechnology and biomedicine. Although their investigation often entails adapting experimental techniques from their protein counterparts, SPAs are fundamentally distinct from proteins, posing unique challenges for their study. Computational methods have emerged as indispensable tools for gaining deeper insights into SPA structures at the molecular level, surpassing the limitations of experimental techniques, and as screening tools to reduce the experimental search space. However, computational studies have grappled with issues stemming from the absence of standardized procedures and relevant crystal structures. Fundamental disparities between SPAs and protein simulations, such as the absence of experimentally validated initial structures and the importance of the simulation size, number of molecules, and concentration, have compounded these challenges. Understanding the roles of various parameters and the capabilities of different models and simulation setups remains an ongoing endeavor. In this review, we aim to provide readers with guidance on the parameters to consider when conducting SPA simulations, elucidating their potential impact on outcomes and validity.

The Legionella collagen-like protein employs a distinct binding mechanism for the recognition of host glycosaminoglycans

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual trimer arrangement with a positively charged external surface and negatively charged solvent exposed internal cavity. Through molecular dynamics simulations, we show how the glycosaminoglycan chondroitin-4-sulphate associates with the Lcl-CTD surface via distinct binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate-binding mechanism.

Advanced computational approaches to understand protein aggregation

Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.

On the Validation of Protein Force Fields Based on Structural Criteria

Molecular dynamics simulations depend critically on the quality of the force field used to describe the interatomic interactions and the extent to which it has been validated for use in a specific application. Using a curated test set of 52 high-resolution structures, 39 derived from X-ray diffraction and 13 solved using NMR, we consider the extent to which different parameter sets of the GROMOS protein force field can be distinguished based on comparing a range of structural criteria, including the number of backbone hydrogen bonds, the number of native hydrogen bonds, polar and nonpolar solvent-accessible surface area, radius of gyration, the prevalence of secondary structure elements, J-coupling constants, nuclear Overhauser effect (NOE) intensities, positional root-mean-square deviations (RMSD), and the distribution of backbone ϕ and ψ dihedral angles. It is shown that while statistically significant differences between the average values of individual metrics could be detected, these were in general small. Furthermore, improvements in agreement in one metric were often offset by loss of agreement in another. The work establishes a framework and test set against which protein force fields can be validated. It also highlights the danger of inferring the relative quality of a given force field based on a small range of structural properties or small number of proteins.

Transient interdomain interactions modulate the monomeric structural ensemble and oligomerization landscape of Huntingtin Exon 1

Polyglutamine expansion (≥ 36 residues) within the N-terminal exon-1 of Huntingtin (Httex1) leads to Huntington's disease, a neurogenerative condition marked by the presence of intranuclear Htt inclusions. Notably, the polyglutamine tract in Httex1 is flanked by an N-terminal coiled-coil domain - N17 (17 amino acids), which undergoes self-association to promote the formation of soluble Httex1 oligomers and brings the aggregation-prone polyQ tracts in close spatial proximity. However, the mechanisms underlying the subsequent conversion of soluble oligomers into insoluble β-rich aggregates with increasing polyQ length, remain unclear. Current knowledge suggests that expansion of the polyQ tract increases its helicity, and this favors its oligomerization and aggregation. In addition, studies utilizing conformation-specific antibodies and a stable coiled-coil heterotetrametric system fused to polyQ indicate that domain "cross-talk" (i.e., interdomain interactions) may be necessary to efficiently promote the emergence of toxic conformations (in monomers and oligomers) and fibrillar aggregation. Here, we performed extensive atomistic molecular dynamics (MD) simulations (aggregate time ∼ 0.7 ms) of N17-polyQ fragments to uncover the interplay between structural transformation and domain "cross-talk" on the monomeric structural ensemble and oligomerization landscape of Httex1. Our simulation ensembles of N17-polyQ monomers validated against 13 C NMR chemical shifts indicated that in addition to elevated α-helicity, polyQ expansion also favors transient, interdomain (N17-polyQ) interactions which result in the emergence of β-conformations. Further, interdomain interactions decreased the overall stability of N17-mediated dimers by counteracting the stabilizing effect of increased α-helicity and promoted a heterogenous oligomerization landscape on the sub-microsecond timescale. Overall, our study uncovers the significance of domain "cross-talk" in modulating the monomeric conformational ensemble and oligomerization landscape of Httex1 to favor the formation of amyloid aggregates.

N4BP1 functions as a dimerization-dependent linear ubiquitin reader which regulates TNF signalling

Signalling through TNFR1 modulates proinflammatory gene transcription and programmed cell death, and its impairment causes autoimmune diseases and cancer. NEDD4-binding protein 1 (N4BP1) is a critical suppressor of proinflammatory cytokine production that acts as a regulator of innate immune signalling and inflammation. However, our current understanding about the molecular properties that enable N4BP1 to exert its suppressive potential remain limited. Here, we show that N4BP1 is a novel linear ubiquitin reader that negatively regulates NFκB signalling by its unique dimerization-dependent ubiquitin-binding module that we named LUBIN. Dimeric N4BP1 strategically positions two non-selective ubiquitin-binding domains to ensure preferential recognition of linear ubiquitin. Under proinflammatory conditions, N4BP1 is recruited to the nascent TNFR1 signalling complex, where it regulates duration of proinflammatory signalling in LUBIN-dependent manner. N4BP1 deficiency accelerates TNFα-induced cell death by increasing complex II assembly. Under proapoptotic conditions, caspase-8 mediates proteolytic processing of N4BP1, resulting in rapid degradation of N4BP1 by the 26 S proteasome, and acceleration of apoptosis. In summary, our findings demonstrate that N4BP1 dimerization creates a novel type of ubiquitin reader that selectively recognises linear ubiquitin which enables the timely and coordinated regulation of TNFR1-mediated inflammation and cell death.

Hierarchical Framework for Predicting Entropies in Bottom-Up Coarse-Grained Models

The thermodynamic entropy of coarse-grained (CG) models stands as one of the most important properties for quantifying the missing information during the CG process and for establishing transferable (or extendible) CG interactions. However, performing additional CG simulations on top of model construction often leads to significant additional computational overhead. In this work, we propose a simple hierarchical framework for predicting the thermodynamic entropies of various molecular CG systems. Our approach employs a decomposition of the CG interactions, enabling the estimation of the CG partition function and thermodynamic properties a priori. Starting from the ideal gas description, we leverage classical perturbation theory to systematically incorporate simple yet essential interactions, ranging from the hard sphere model to the generalized van der Waals model. Additionally, we propose an alternative approach based on multiparticle correlation functions, allowing for systematic improvements through higher-order correlations. Numerical applications to molecular liquids validate the high fidelity of our approach, and our computational protocols demonstrate that a reduced model with simple energetics can reasonably estimate the thermodynamic entropy of CG models without performing any CG simulations. Overall, our findings present a systematic framework for estimating not only the entropy but also other thermodynamic properties of CG models, relying solely on information from the reference system.

Effect of caffeine on the aggregation of amyloid-β-A 3D RISM study

Alzheimer's disease is a detrimental neurological disorder caused by the formation of amyloid fibrils due to the aggregation of amyloid-β peptide. The primary therapeutic approaches for treating Alzheimer's disease are targeted to prevent this amyloid fibril formation using potential inhibitor molecules. The discovery of such inhibitor molecules poses a formidable challenge to the design of anti-amyloid drugs. This study investigates the effect of caffeine on dimer formation of the full-length amyloid-β using a combined approach of all-atom, explicit water molecular dynamics simulations and the three-dimensional reference interaction site model theory. The change in the hydration free energy of amyloid-β dimer, with and without the inhibitor molecules, is calculated with respect to the monomeric amyloid-β, where the hydration free energy is decomposed into energetic and entropic components, respectively. Dimerization is accompanied by a positive change in the partial molar volume. Dimer formation is spontaneous, which implies a decrease in the hydration free energy. However, a reverse trend is observed for the dimer with inhibitor molecules. It is observed that the negatively charged residues primarily contribute for the formation of the amyloid-β dimer. A residue-wise decomposition reveals that hydration/dehydration of the side-chain atoms of the charged amino acid residues primarily contribute to dimerization.

A comprehensive computational study to explore promising natural bioactive compounds targeting glycosyltransferase MurG in Escherichia coli for potential drug development

Peptidoglycan is a carbohydrate with a cross-linked structure that protects the cytoplasmic membrane of bacterial cells from damage. The mechanism of peptidoglycan biosynthesis involves the main synthesizing enzyme glycosyltransferase MurG, which is known as a potential target for antibiotic therapy. Many MurG inhibitors have been recognized as MurG targets, but high toxicity and drug-resistant Escherichia coli strains remain the most important problems for further development. In addition, the discovery of selective MurG inhibitors has been limited to the synthesis of peptidoglycan-mimicking compounds. The present study employed drug discovery, such as virtual screening using molecular docking, drug likeness ADMET proprieties predictions, and molecular dynamics (MD) simulation, to identify potential natural products (NPs) for Escherichia coli. We conducted a screening of 30,926 NPs from the NPASS database. Subsequently, 20 of these compounds successfully passed the potency, pharmacokinetic, ADMET screening assays, and their validation was further confirmed through molecular docking. The best three hits and the standard were chosen for further MD simulations up to 400 ns and energy calculations to investigate the stability of the NPs-MurG complexes. The analyses of MD simulations and total binding energies suggested the higher stability of NPC272174. The potential compounds can be further explored in vivo and in vitro for promising novel antibacterial drug discovery.

Unraveling motion in proteins by combining NMR relaxometry and molecular dynamics simulations: A case study on ubiquitin

Nuclear magnetic resonance (NMR) relaxation experiments shine light onto the dynamics of molecular systems in the picosecond to millisecond timescales. As these methods cannot provide an atomically resolved view of the motion of atoms, functional groups, or domains giving rise to such signals, relaxation techniques have been combined with molecular dynamics (MD) simulations to obtain mechanistic descriptions and gain insights into the functional role of side chain or domain motion. In this work, we present a comparison of five computational methods that permit the joint analysis of MD simulations and NMR relaxation experiments. We discuss their relative strengths and areas of applicability and demonstrate how they may be utilized to interpret the dynamics in MD simulations with the small protein ubiquitin as a test system. We focus on the aliphatic side chains given the rigidity of the backbone of this protein. We find encouraging agreement between experiment, Markov state models built in the χ1/χ2 rotamer space of isoleucine residues, explicit rotamer jump models, and a decomposition of the motion using ROMANCE. These methods allow us to ascribe the dynamics to specific rotamer jumps. Simulations with eight different combinations of force field and water model highlight how the different metrics may be employed to pinpoint force field deficiencies. Furthermore, the presented comparison offers a perspective on the utility of NMR relaxation to serve as validation data for the prediction of kinetics by state-of-the-art biomolecular force fields.

Approaching the catalytic mechanism of protein lysine methyltransferases by biochemical and simulation techniques

Protein lysine methyltransferases (PKMTs) transfer up to three methyl groups to the side chains of lysine residues in proteins and fulfill important regulatory functions by controlling protein stability, localization and protein/protein interactions. The methylation reactions are highly regulated, and aberrant methylation of proteins is associated with several types of diseases including neurologic disorders, cardiovascular diseases, and various types of cancer. This review describes novel insights into the catalytic machinery of various PKMTs achieved by the combined application of biochemical experiments and simulation approaches during the last years, focusing on clinically relevant and well-studied enzymes of this group like DOT1L, SMYD1-3, SET7/9, G9a/GLP, SETD2, SUV420H2, NSD1/2, different MLLs and EZH2. Biochemical experiments have unraveled many mechanistic features of PKMTs concerning their substrate and product specificity, processivity and the effects of somatic mutations observed in PKMTs in cancer cells. Structural data additionally provided information about the substrate recognition, enzyme-substrate complex formation, and allowed for simulations of the substrate peptide interaction and mechanism of PKMTs with atomistic resolution by molecular dynamics and hybrid quantum mechanics/molecular mechanics methods. These simulation technologies uncovered important mechanistic details of the PKMT reaction mechanism including the processes responsible for the deprotonation of the target lysine residue, essential conformational changes of the PKMT upon substrate binding, but also rationalized regulatory principles like PKMT autoinhibition. Further developments are discussed that could bring us closer to a mechanistic understanding of catalysis of this important class of enzymes in the near future. The results described here illustrate the power of the investigation of enzyme mechanisms by the combined application of biochemical experiments and simulation technologies.

Understanding the cell: Future views of structural biology

Determining the structure and mechanisms of all individual functional modules of cells at high molecular detail has often been seen as equal to understanding how cells work. Recent technical advances have led to a flush of high-resolution structures of various macromolecular machines, but despite this wealth of detailed information, our understanding of cellular function remains incomplete. Here, we discuss present-day limitations of structural biology and highlight novel technologies that may enable us to analyze molecular functions directly inside cells. We predict that the progression toward structural cell biology will involve a shift toward conceptualizing a 4D virtual reality of cells using digital twins. These will capture cellular segments in a highly enriched molecular detail, include dynamic changes, and facilitate simulations of molecular processes, leading to novel and experimentally testable predictions. Transferring biological questions into algorithms that learn from the existing wealth of data and explore novel solutions may ultimately unveil how cells work.

Applications of Molecular Dynamics Simulations in Drug Discovery

In the current drug development process, molecular dynamics (MD) simulations have proven to be very useful. This chapter provides an overview of the current applications of MD simulations in drug discovery, from detecting protein druggable sites and validating drug docking outcomes to exploring protein conformations and investigating the influence of mutations on its structure and functions. In addition, this chapter emphasizes various strategies to improve the conformational sampling efficiency in molecular dynamics simulations. With a growing computer power and developments in the production of force fields and MD techniques, the importance of MD simulations in helping the drug development process is projected to rise significantly in the future.

Advances in Drug Discovery and Design using Computer-aided Molecular Modeling

Computer-aided molecular modeling is a rapidly emerging technology that is being used to accelerate the discovery and design of new drug therapies. It involves the use of computer algorithms and 3D structures of molecules to predict interactions between molecules and their behavior in the body. This has drastically improved the speed and accuracy of drug discovery and design. Additionally, computer-aided molecular modeling has the potential to reduce costs, increase the quality of data, and identify promising targets for drug development. Through the use of sophisticated methods, such as virtual screening, molecular docking, pharmacophore modeling, and quantitative structure-activity relationships, scientists can achieve higher levels of efficacy and safety for new drugs. Moreover, it can be used to understand the activity of known drugs and simplify the process of formulating, optimizing, and predicting the pharmacokinetics of new and existing drugs. In conclusion, computer-aided molecular modeling is an effective tool to rapidly progress drug discovery and design by predicting the interactions between molecules and anticipating the behavior of new drugs in the body.

Atomistic molecular simulations of Aβ-Zn conformational ensembles

The amyloid-forming Aβ peptide is able to interact with metal cations to form very stable complexes that influence fibril formation and contribute to the onset of Alzheimer's disease. Multiple structures of peptides derived from Aβ in complex with different metals have been resolved experimentally to provide an atomic-level description of the metal-protein interactions. However, Aβ is intrinsically disordered, and hence more amenable to an ensemble description. Molecular dynamics simulations can now reach the timescales needed to generate ensembles for these type of complexes. However, this requires accurate force fields both for the protein and the protein-metal interactions. Here we use state-of-the-art methods to generate force field parameters for the Zn(II) cations in a set of complexes formed by different Aβ variants and combine them with the Amber99SB*-ILDN optimized force field. Upon comparison of NMR experiments with the simulation results, further optimized with a Bayesian/Maximum entropy approach, we provide an accurate description of the molecular ensembles for most Aβ-metal complexes. We find that the resulting conformational ensembles are more heterogeneous than the NMR models deposited in the Protein Data Bank.

Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Despite advances in artificial intelligence methods, protein folding remains in many ways an enigma to be solved. Accurate computation of protein folding energetics could help drive fields such as protein and drug design and genetic interpretation. However, the challenge of calculating the state functions governing protein folding from first-principles remains unaddressed. We present here a simple approach that allows us to accurately calculate the energetics of protein folding. It is based on computing the energy of the folded and unfolded states at different temperatures using molecular dynamics simulations. From this, two essential quantities (ΔH and ΔCp) are obtained and used to calculate the conformational stability of the protein (ΔG). With this approach, we have successfully calculated the energetics of two- and three-state proteins, representatives of the major structural classes, as well as small stability differences (ΔΔG) due to changes in solution conditions or variations in an amino acid residue.

The Legionella collagen-like protein employs a unique binding mechanism for the recognition of host glycosaminoglycans

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans (GAGs) on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual dynamic trimer arrangement with a positively charged external surface and a negatively charged solvent exposed internal cavity. Through Molecular Dynamics (MD) simulations, we show how the GAG chondroitin-4-sulphate associates with the Lcl-CTD surface via unique binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate binding mechanism.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) make up the largest superfamily of human membrane proteins and represent primary targets of ∼1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon the binding of positive and negative allosteric modulators (PAMs and NAMs). The mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy prOfiling Workflow (GLOW). GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence and absence of the modulator. DL and free energy calculations revealed significantly reduced dynamic fluctuations and conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G-protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes. Therefore, GPCR allostery exhibits a dynamic "conformational selection" mechanism. In the absence of available modulator-bound structures as for most current GPCRs, it is critical to use a structural ensemble of representative GPCR conformations rather than a single structure for compound docking ("ensemble docking"), which will potentially improve structure-based design of novel allosteric drugs of GPCRs.

Abnormal myosin post-translational modifications and ATP turnover time associated with human congenital myopathy-related RYR1 mutations

AIM

Conditions related to mutations in the gene encoding the skeletal muscle ryanodine receptor 1 (RYR1) are genetic muscle disorders and include congenital myopathies with permanent weakness, as well as episodic phenotypes such as rhabdomyolysis/myalgia. Although RYR1 dysfunction is the primary mechanism in RYR1-related disorders, other downstream pathogenic events are less well understood and may include a secondary remodeling of major contractile proteins. Hence, in the present study, we aimed to investigate whether congenital myopathy-related RYR1 mutations alter the regulation of the most abundant contractile protein, myosin.

METHODS

We used skeletal muscle tissues from five patients with RYR1-related congenital myopathy and compared those with five controls and five patients with RYR1-related rhabdomyolysis/myalgia. We then defined post-translational modifications on myosin heavy chains (MyHCs) using LC/MS. In parallel, we determined myosin relaxed states using Mant-ATP chase experiments and performed molecular dynamics (MD) simulations.

RESULTS

LC/MS revealed two additional phosphorylations (Thr1309-P and Ser1362-P) and one acetylation (Lys1410-Ac) on the β/slow MyHC of patients with congenital myopathy. This method also identified six acetylations that were lacking on MyHC type IIa of these patients (Lys35-Ac, Lys663-Ac, Lys763-Ac, Lys1171-Ac, Lys1360-Ac, and Lys1733-Ac). MD simulations suggest that modifying myosin Ser1362 impacts the protein structure and dynamics. Finally, Mant-ATP chase experiments showed a faster ATP turnover time of myosin heads in the disordered-relaxed conformation.

CONCLUSIONS

Altogether, our results suggest that RYR1 mutations have secondary negative consequences on myosin structure and function, likely contributing to the congenital myopathic phenotype.

In silico design and identification of new peptides for mitigating hIAPP aggregation in type 2 diabetes

The aberrant misfolding and self-aggregation of human islet amyloid polypeptide (hIAPP or amylin) into cytotoxic aggregates are implicated in the pathogenesis of type 2 diabetes (T2D). Among various inhibitors, short peptides derived from the amyloidogenic regions of hIAPP have been employed as hIAPP aggregation inhibitors due to their low immunogenicity, biocompatibility, and high chemical diversity. Recently, hIAPP fragment HSSNN18-22 was identified as an amyloidogenic sequence and displayed higher antiproliferative activity to RIN-5F cells. Various hIAPP aggregation inhibitors have been designed by chemical modifications of the highly amyloidogenic sequence (NFGAIL) of hIAPP. In this work, a library of pentapeptides based on fragment HSSNN18-22 was designed and assessed for their efficacy in blocking hIAPP aggregation using an integrated computational screening approach. The binding free energy calculations by molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method identified HSSQN and HSSNQ that bind to hIAPP monomer with a binding affinity of -21.25 ± 4.90 and -19.73 ± 3.10 kcal/mol, respectively, which is notably higher as compared to HSSNN (-11.90 ± 4.12 kcal/mol). The sampling of the non aggregation-prone helical conformation was notably increased from 23.5 ± 3.0 in the hIAPP monomer to 38.1 ± 3.6, and 33.8 ± 3.0% on the incorporation of HSSQN, and HSSNQ, respectively, which indicate reduced aggregation propensity of hIAPP monomer. The pentapeptides, HSSQN and HSSNQ, identified as hIAPP aggregation inhibitors in this work can be further conjugated with various metal chelating peptides to yield more efficacious and clinically relevant multifunctional modulators for targeting various pathological hallmarks of T2D.

Optimizing molecular potential models by imposing kinetic constraints with path reweighting

Empirical force fields employed in molecular dynamics simulations of complex systems are often optimized to reproduce experimentally determined structural and thermodynamic properties. In contrast, experimental knowledge about the interconversion rates between metastable states in such systems is hardly ever incorporated in a force field due to a lack of an efficient approach. Here, we introduce such a framework based on the relationship between dynamical observables, such as rate constants, and the underlying molecular model parameters using the statistical mechanics of trajectories. Given a prior ensemble of molecular dynamics trajectories produced with imperfect force field parameters, the approach allows for the optimal adaption of these parameters such that the imposed constraint of equally predicted and experimental rate constant is obeyed. To do so, the method combines the continuum path ensemble maximum caliber approach with path reweighting methods for stochastic dynamics. When multiple solutions are found, the method selects automatically the combination that corresponds to the smallest perturbation of the entire path ensemble, as required by the maximum entropy principle. To show the validity of the approach, we illustrate the method on simple test systems undergoing rare event dynamics. Next to simple 2D potentials, we explore particle models representing molecular isomerization reactions and protein-ligand unbinding. Besides optimal interaction parameters, the methodology gives physical insights into what parts of the model are most sensitive to the kinetics. We discuss the generality and broad implications of the methodology.

Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein

Molecular dynamics (MD) is an extremely powerful, highly effective, and widely used approach to understanding the nature of chemical processes in atomic details for proteins. The accuracy of results from MD simulations is highly dependent on force fields. Currently, molecular mechanical (MM) force fields are mainly utilized in MD simulations because of their low computational cost. Quantum mechanical (QM) calculation has high accuracy, but it is exceedingly time consuming for protein simulations. Machine learning (ML) provides the capability for generating accurate potential at the QM level without increasing much computational effort for specific systems that can be studied at the QM level. However, the construction of general machine learned force fields, needed for broad applications and large and complex systems, is still challenging. Here, general and transferable neural network (NN) force fields based on CHARMM force fields, named CHARMM-NN, are constructed for proteins by training NN models on 27 fragments partitioned from the residue-based systematic molecular fragmentation (rSMF) method. The NN for each fragment is based on atom types and uses new input features that are similar to MM inputs, including bonds, angles, dihedrals, and non-bonded terms, which enhance the compatibility of CHARMM-NN to MM MD and enable the implementation of CHARMM-NN force fields in different MD programs. While the main part of the energy of the protein is based on rSMF and NN, the nonbonded interactions between the fragments and with water are taken from the CHARMM force field through mechanical embedding. The validations of the method for dipeptides on geometric data, relative potential energies, and structural reorganization energies demonstrate that the CHARMM-NN local minima on the potential energy surface are very accurate approximations to QM, showing the success of CHARMM-NN for bonded interactions. However, the MD simulations on peptides and proteins indicate that more accurate methods to represent protein-water interactions in fragments and non-bonded interactions between fragments should be considered in the future improvement of CHARMM-NN, which can increase the accuracy of approximation beyond the current mechanical embedding QM/MM level.

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.

Development and Benchmarking of Open Force Field 2.0.0: The Sage Small Molecule Force Field

We introduce the Open Force Field (OpenFF) 2.0.0 small molecule force field for drug-like molecules, code-named Sage, which builds upon our previous iteration, Parsley. OpenFF force fields are based on direct chemical perception, which generalizes easily to highly diverse sets of chemistries based on substructure queries. Like the previous OpenFF iterations, the Sage generation of OpenFF force fields was validated in protein-ligand simulations to be compatible with AMBER biopolymer force fields. In this work, we detail the methodology used to develop this force field, as well as the innovations and improvements introduced since the release of Parsley 1.0.0. One particularly significant feature of Sage is a set of improved Lennard-Jones (LJ) parameters retrained against condensed phase mixture data, the first refit of LJ parameters in the OpenFF small molecule force field line. Sage also includes valence parameters refit to a larger database of quantum chemical calculations than previous versions, as well as improvements in how this fitting is performed. Force field benchmarks show improvements in general metrics of performance against quantum chemistry reference data such as root-mean-square deviations (RMSD) of optimized conformer geometries, torsion fingerprint deviations (TFD), and improved relative conformer energetics (ΔΔE). We present a variety of benchmarks for these metrics against our previous force fields as well as in some cases other small molecule force fields. Sage also demonstrates improved performance in estimating physical properties, including comparison against experimental data from various thermodynamic databases for small molecule properties such as ΔHmix, ρ(x), ΔGsolv, and ΔGtrans. Additionally, we benchmarked against protein-ligand binding free energies (ΔGbind), where Sage yields results statistically similar to previous force fields. All the data is made publicly available along with complete details on how to reproduce the training results at https://github.com/openforcefield/openff-sage.

Fitting Force Field Parameters to NMR Relaxation Data

We present an approach to optimize force field parameters using time-dependent data from NMR relaxation experiments. To do so, we scan parameters in the dihedral angle potential energy terms describing the rotation of the methyl groups in proteins and compare NMR relaxation rates calculated from molecular dynamics simulations with the modified force fields to deuterium relaxation measurements of T4 lysozyme. We find that a small modification of C^γ methyl groups improves the agreement with experiments both for the protein used to optimize the force field and when validating using simulations of CI2 and ubiquitin. We also show that these improvements enable a more effective a posteriori reweighting of the MD trajectories. The resulting force field thus enables more direct comparison between simulations and side-chain NMR relaxation data and makes it possible to construct ensembles that better represent the dynamics of proteins in solution.

Structural insights into SARS-CoV-2 main protease conformational plasticity

The spread of different severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants underscores the need for insights into the structural properties of its structural and non-structural proteins. The highly conserved homo-dimeric chymotrypsin-like protease (3CL M^PRO ), belonging to the class of cysteine hydrolases, plays an indispensable role in processing viral polyproteins that are involved in viral replication and transcription. Studies have successfully demonstrated the role of M^PRO as an attractive drug target for designing antiviral treatments because of its importance in the viral life cycle. Herein, we report the structural dynamics of six experimentally solved structures of M^PRO (i.e., 6LU7, 6M03, 6WQF, 6Y2E, 6Y84, and 7BUY including both ligand-free and ligand-bound states) at different resolutions. We have employed a structure-based balanced forcefield, CHARMM36m through state-of-the-art all-atoms molecular dynamics simulations at µ-seconds scale at room temperature (303K) and pH 7.0 to explore their structure-function relationship. The helical domain-III responsible for dimerization mostly contributes to the altered conformational states and destabilization of M^PRO . A keen observation of the high degree of flexibility in the P5 binding pocket adjoining domain II-III highlights the reason for observation of conformational heterogeneity among the structural ensembles of M^PRO . We also observe a differential dynamics of the catalytic pocket residues His41, Cys145, and Asp187, which may lead to catalytic impairment of the monomeric proteases. Among the highly populated conformational states of the six systems, 6LU7 and 7M03 forms the most stable and compact M^PRO conformation with intact catalytic site and structural integrity. Altogether, our findings from this extensive study provides a benchmark to identify physiologically relevant structures of such promising drug targets for structure-based drug design and discovery of potent drug-like compounds having clinical potential.

Aggregation propensities of proteins with varying degrees of disorder

The hydration thermodynamics of a globular protein (AcP), three intrinsically disordered protein regions (1CD3, 1MVF, 1F0R) and a fully disordered protein (α-synuclein) is studied by an approach that combines an all-atom explicit water molecular dynamics simulations and three-dimensional reference interaction site model (3D-RISM) theory. The variation in hydration free energy with percentage disorder of the selected proteins is investigated through its nonelectrostatic and electrostatic components. A decrease in hydration free energy is observed with an increase in percentage disorder, indicating favorable interactions of the disordered proteins with the solvent. This confirms the role of percentage disorder in determining the aggregation propensity of proteins which is measured in terms of the hydration free energy in addition to their respective mean net charge and mean hydrophobicity. The hydration free energy is decoupled into energetic and entropic terms. A residue-wise decomposition analysis of the hydration free energy for the selected proteins is evaluated. The decomposition shows that the disordered regions contribute more than the ordered ones for the intrinsically disordered protein regions. The dominant role of electrostatic interactions is confirmed from the residue-wise decomposition of the hydration free energy. The results depict that the negatively charged residues contribute more to the total hydration free energy for the proteins with negative mean net charge, while the positively charged residues contribute more for proteins with positive mean net charge.

Relaxation dynamics measure the aggregation propensity of amyloid-β and its mutants

Atomistic molecular dynamics simulations are employed to investigate the global and segmental relaxation dynamics of the amyloid-β protein and its causative and protective mutants. Amyloid-β exhibits significant global/local dynamics that span a broad range of length and time scales due to its intrinsically disordered nature. The relaxation dynamics of the amyloid-β protein and its mutants is quantitatively correlated with its experimentally measured aggregation propensity. The protective mutant has slower relaxation dynamics, whereas the causative mutants exhibit faster global dynamics compared with that of the wild-type amyloid-β. The local dynamics of the amyloid-β protein or its mutants is governed by a complex interplay of the charge, hydrophobicity, and change in the molecular mass of the mutated residue.

A Comparison of Bonded and Nonbonded Zinc(II) Force Fields with NMR Data

Classical molecular dynamics (MD) simulations are widely used to inspect the behavior of zinc(II)-proteins at the atomic level, hence the need to properly model the zinc(II) ion and the interaction with its ligands. Different approaches have been developed to represent zinc(II) sites, with the bonded and nonbonded models being the most used. In the present work, we tested the well-known zinc AMBER force field (ZAFF) and a recently developed nonbonded force field (NBFF) to assess how accurately they reproduce the dynamic behavior of zinc(II)-proteins. For this, we selected as benchmark six zinc-fingers. This superfamily is extremely heterogenous in terms of architecture, binding mode, function, and reactivity. From repeated MD simulations, we computed the order parameter (S2) of all backbone N-H bond vectors in each system. These data were superimposed to heteronuclear Overhauser effect measurements taken by NMR spectroscopy. This provides a quantitative estimate of the accuracy of the FFs in reproducing protein dynamics, leveraging the information about the protein backbone mobility contained in the NMR data. The correlation between the MD-computed S2 and the experimental data indicated that both tested FFs reproduce well the dynamic behavior of zinc(II)-proteins, with comparable accuracy. Thus, along with ZAFF, NBFF represents a useful tool to simulate metalloproteins with the advantage of being extensible to diverse systems such as those bearing dinuclear metal sites.

Influence of force field choice on the conformational landscape of rat and human islet amyloid polypeptide

Human islet amyloid polypeptide (hIAPP) is a naturally occurring, intrinsically disordered protein (IDP) whose abnormal aggregation into toxic soluble oligomers and insoluble amyloid fibrils is a pathological feature in type-2 diabetes. Rat IAPP (rIAPP) differs from hIAPP by only six amino acids yet has a reduced tendency to aggregate or form fibrils. The structures of the monomeric forms of IAPP are difficult to characterize due to their intrinsically disordered nature. Molecular dynamics simulations can provide a detailed characterization of the monomeric forms of rIAPP and hIAPP in near-physiological conditions. In this work, the conformational landscapes of rIAPP and hIAPP as a function of secondary structure content were predicted using well-tempered bias exchange metadynamics simulations. Several combinations of commonly used biomolecular force fields and water models were tested. The predicted conformational preferences of both rIAPP and hIAPP are typical of IDPs, exhibiting dominant random coil structures but showing a low propensity for transient α-helical conformations. Predicted nuclear magnetic resonance Cα chemical shifts reveal different preferences with each force field towards certain conformations, with AMBERff99SBnmr2/TIP4Pd showing the best agreement with the experiment. Comparisons of secondary structure content demonstrate residue-specific differences between hIAPP and rIAPP that may reflect their different aggregation propensities.

Protein Design Using Physics Informed Neural Networks

The inverse protein folding problem, also known as protein sequence design, seeks to predict an amino acid sequence that folds into a specific structure and performs a specific function. Recent advancements in machine learning techniques have been successful in generating functional sequences, outperforming previous energy function-based methods. However, these machine learning methods are limited in their interoperability and robustness, especially when designing proteins that must function under non-ambient conditions, such as high temperature, extreme pH, or in various ionic solvents. To address this issue, we propose a new Physics-Informed Neural Networks (PINNs)-based protein sequence design approach. Our approach combines all-atom molecular dynamics simulations, a PINNs MD surrogate model, and a relaxation of binary programming to solve the protein design task while optimizing both energy and the structural stability of proteins. We demonstrate the effectiveness of our design framework in designing proteins that can function under non-ambient conditions.

Modern non-polarizable force fields diverge in modeling the enzyme-substrate complex of a canonical serine protease

Classical molecular dynamics simulation is a powerful and established method of modern computational chemistry. Being able to obtain accurate information on molecular behavior is crucial to get valuable insights into structure-function relationships that translate into fundamental findings and practical applications. Active sites of enzymes are known to be particularly intricate, therefore, simpler non-polarizable force fields may provide an inaccurate description. In this work, we addressed this hypothesis in a case of a canonical serine triad protease trypsin in its complex with a substrate-mimicking inhibitor. We tested six modern and popular force fields to find that significantly diverging results may be obtained. Amber FB-15 and OPLS-AA/M turned out to model the active site incorrectly. Amber ff19sb and ff15ipq demonstrated mixed performance. The best performing force fields were CHARMM36m and Amber ff99sb-ildn, therefore, they are recommended for use with this and related systems. We speculate that a similar lack of cross-force field convergence may be characteristic of other enzymatic systems. Therefore, we advocate for careful consideration of different force fields in any study within the field of computational enzymology.

An Imbalance in the Force: The Need for Standardized Benchmarks for Molecular Simulation

Force fields (FFs) for molecular simulation have been under development for more than half a century. As with any predictive model, rigorous testing and comparisons of models critically depends on the availability of standardized data sets and benchmarks. While such benchmarks are rather common in the fields of quantum chemistry, this is not the case for empirical FFs. That is, few benchmarks are reused to evaluate FFs, and development teams rather use their own training and test sets. Here we present an overview of currently available tests and benchmarks for computational chemistry, focusing on organic compounds, including halogens and common ions, as FFs for these are the most common ones. We argue that many of the benchmark data sets from quantum chemistry can in fact be reused for evaluating FFs, but new gas phase data is still needed for compounds containing phosphorus and sulfur in different valence states. In addition, more nonequilibrium interaction energies and forces, as well as molecular properties such as electrostatic potentials around compounds, would be beneficial. For the condensed phases there is a large body of experimental data available, and tools to utilize these data in an automated fashion are under development. If FF developers, as well as researchers in artificial intelligence, would adopt a number of these data sets, it would become easier to compare the relative strengths and weaknesses of different models and to, eventually, restore the balance in the force.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Structural Perturbation of Monomers Determines the Amyloid Aggregation Propensity of Calcitonin Variants

Human calcitonin (hCT) is a polypeptide hormone that participates in calcium-phosphorus metabolism. Irreversible aggregation of 32-amino acid hCT into β-sheet-rich amyloid fibrils impairs physiological activity and increases the risk of medullary carcinoma of the thyroid. Amyloid-resistant hCT derivatives substituting critical amyloidogenic residues are of particular interest for clinical applications as therapeutic drugs against bone-related diseases. Uncovering the aggregation mechanism of hCT at the molecular level, therefore, is important for the design of amyloid-resistant hCT analogues. Here, we investigated the aggregation dynamics of hCT, non-amyloidogenic salmon calcitonin (sCT), and two hCT analogues with reduced aggregation tendency─TL-hCT and phCT─using long timescale discrete molecular dynamics simulations. Our results showed that hCT monomers mainly adopted unstructured conformations with dynamically formed helices around the central region. hCT self-assembled into helix-rich oligomers first, followed by a conformational conversion into β-sheet-rich oligomers with β-sheets formed by residues 10-30 and stabilized by aromatic and hydrophobic interactions. Our simulations confirmed that TL-hCT and phCT oligomers featured more helices and fewer β-sheets than hCT. Substitution of central aromatic residues with leucine in TL-hCT and replacing C-terminal hydrophobic residue with hydrophilic amino acid in phCT only locally suppressed β-sheet propensities in the central region and C-terminus, respectively. Having mutations in both central and C-terminal regions, sCT monomers and dynamically formed oligomers predominantly adopted helices, confirming that both central aromatic and C-terminal hydrophobic residues played important roles in the fibrillization of hCT. We also observed the formation of β-barrel intermediates, postulated as the toxic oligomers in amyloidosis, for hCT but not for sCT. Our computational study depicts a complete picture of the aggregation dynamics of hCT and the effects of mutations. The design of next-generation amyloid-resistant hCT analogues should consider the impact on both amyloidogenic regions and also take into account the amplification of transient β-sheet population in monomers upon aggregation.

The Effect of Arginine on the Phase Stability of Aqueous Hen Egg-White Lysozyme Solutions

The effect of arginine on the phase stability of the hen egg-white lysozyme (HEWL) has been studied via molecular dynamics computer simulations, as well as experimentally via cloud-point temperature determination. The experiments show that the addition of arginine increases the stability of the HEWL solutions. The computer simulation results indicate that arginine molecules tend to self-associate. If arginine residues are located on the protein surface, the free arginine molecules stay in their vicinity and prevent the way protein molecules "connect" through them to form clusters. The results are not sensitive to a particular force field and suggest a possible microscopic mechanism of the stabilizing role of arginine as an excipient.

Comparative study of binding pocket structure and dynamics in cardiac and skeletal myosin

The development of small molecule myosin modulators has seen an increased effort in recent years due to their possible use in the treatment of cardiac and skeletal myopathies. Omecamtiv mecarbil (OM) is the first-in-class cardiac myotrope and the first to enter clinical trials. Its selectivity toward slow/beta-cardiac myosin lies at the heart of its function; however, little is known about the underlying reasons for selectivity to this isoform as opposed to other closely related ones such as fast-type skeletal myosins. In this work, we compared the structure and dynamics of the OM binding site in cardiac and in fasttype IIa skeletal myosin to identify possible reasons for OM selectivity. We found that the different shape, size, and composition of the binding pocket in skeletal myosin directly affects the binding mode and related affinity of OM, which is potentially a result of weaker interactions and less optimal molecular recognition. Moreover, we identified a side pocket adjacent to the OM binding site that shows increased accessibility in skeletal myosin compared with the cardiac isoform. These findings could pave the way to the development of skeletal-selective compounds that can target this region of the protein and potentially be used to treat congenital myopathies where muscle weakness is related to myosin loss of function.

Molecular Dynamics Methods for Antibody Design

Complex and coordinated dynamics are closely connected with protein functions, including the binding of antibodies to antigens. Knowledge of such dynamics could improve the design of antibodies. Molecular dynamics (MD) simulations provide a "computational microscope" that can resolve atomic motions and inform antibody design efforts.

Computer Aided Development of Nucleic Acid Applications in Nanotechnologies

Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.

Molecular Dynamics Simulation of the Thermal Behavior of Hydroxyapatite

Hydroxyapatite (HAP) is the main mineral component of bones and teeth. Due to its biocompatibility, HAP is widely used in medicine as a filler that replaces parts of lost bone and as an implant coating that promotes new bone growth. The modeling and calculations of the structure and properties of HAP showed that various structural defects have a significant effect on the properties of the material. By varying these structural heterogeneities, it is possible to increase the biocompatibility of HAP. An important role here is played by OH group vacancies, which are easily formed when these hydroxyl groups leave OH channels of HAP. In this case, the temperature dependence of the concentration of OH ions, which also determines the thermal behavior of HAP, is important. To study the evaporation of OH ions from HAP structures with increasing temperatures, molecular dynamics simulation (MDS) methods were used in this work. As a program for MDS modeling, we used the PUMA-CUDA software package. The initial structure of HAP, consisting of 4 × 4 × 2 = 32 unit cells of the hexagonal HAP phase, surrounded by a 15-Å layer of water was used in the modelling. Multiple and statistically processed MDS, running calculations in the range of 700-1400 K, showed that active evaporation of OH ions begins at the temperature of 1150 K. The analysis of the obtained results in comparison with those available in the literature data shows that these values are very close to the experiments. Thus, this MDS approach demonstrates its effective applicability and shows good results in the study of the thermal behavior of HAP.

Exploiting Sequence-Dependent Rotamer Information in Global Optimization of Proteins

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

Comparison of the United- and All-Atom Representations of (Halo)alkanes Based on Two Condensed-Phase Force Fields Optimized against the Same Experimental Data Set

The level of accuracy that can be achieved by a force field is influenced by choices made in the interaction-function representation and in the relevant simulation parameters. These choices, referred to here as functional-form variants (FFVs), include for example the model resolution, the charge-derivation procedure, the van der Waals combination rules, the cutoff distance, and the treatment of the long-range interactions. Ideally, assessing the effect of a given FFV on the intrinsic accuracy of the force-field representation requires that only the specific FFV is changed and that this change is performed at an optimal level of parametrization, a requirement that may prove extremely challenging to achieve in practice. Here, we present a first attempt at such a comparison for one specific FFV, namely the choice of a united-atom (UA) versus an all-atom (AA) resolution in a force field for saturated acyclic (halo)alkanes. Two force-field versions (UA vs AA) are optimized in an automated way using the CombiFF approach against 961 experimental values for the pure-liquid densities ρ_liq and vaporization enthalpies ΔH_vap of 591 compounds. For the AA force field, the torsional and third-neighbor Lennard-Jones parameters are also refined based on quantum-mechanical rotational-energy profiles. The comparison between the UA and AA resolutions is also extended to properties that have not been included as parameterization targets, namely the surface-tension coefficient γ, the isothermal compressibility κ_T, the isobaric thermal-expansion coefficient α_P, the isobaric heat capacity c_P, the static relative dielectric permittivity ϵ, the self-diffusion coefficient D, the shear viscosity η, the hydration free energy ΔG_wat, and the free energy of solvation ΔG_che in cyclohexane. For the target properties ρ_liq and ΔH_vap, the UA and AA resolutions reach very similar levels of accuracy after optimization. For the nine other properties, the AA representation leads to more accurate results in terms of η; comparably accurate results in terms of γ, κ_T, α_P, ϵ, D, and ΔG_che; and less accurate results in terms of c_P and ΔG_wat. This work also represents a first step toward the calibration of a GROMOS-compatible force field at the AA resolution.

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.

A comprehensive review on preparation, structure-activities relationship, and calcium bioavailability of casein phosphopeptides

Calcium is one of the important elements for human health. Calcium deficiencies can lead to numerous diseases. Calcium chelating peptides have shown potential application in the management of calcium deficiencies. Casein phosphopeptides (CPP) are phosphoseryl-containing fragments of casein by enzymatic hydrolysis or fermentation during manufacture of milk products as well as during intestinal digestion. An increasing number of CPP with the ability to facilitate and enhance the bioavailability of calcium are being discovered and identified. In this review, 249 reported CPP derived from four types of bovine casein (αs1, αs2, β and κ) were collected, and the amino acid sequence and phosphoserine group information were sorted out. This review outlines the current enzyme hydrolysis, detection methods, purification, structure-activity relationship and mechanism of intestinal calcium absorption in vitro and in vivo as well as application of CPP.