Water-membrane partition thermodynamics of an amphiphilic lipopeptide: an enthalpy-driven hydrophobic effect

GTP-Bound N-Ras Conformational States and Substates Are Modulated by Membrane and Point Mutation

Oncogenic Ras proteins are known to present multiple conformational states, as reported by the great variety of crystallographic structures. The GTP-bound states are grouped into two main states: the "inactive" state 1 and the "active" state 2. Recent reports on H-Ras have shown that state 2 exhibits two substates, directly related to the orientation of Tyr32: toward the GTP-bound pocket and outwards. In this paper, we show that N-Ras exhibits another substate of state 2, related to a third orientation of Tyr32, toward Ala18 and parallel to the GTP-bound pocket. We also show that this substate is highly sampled in the G12V mutation of N-Ras and barely present in its wild-type form, and that the G12V mutation prohibits the sampling of the GTPase-activating protein (GAP) binding substate, rendering this mutation oncogenic. Furthermore, using molecular dynamics simulations, we explore the importance of the membrane on N-Ras' conformational state dynamics and its strong influence on Ras protein stability. Moreover, the membrane has a significant influence on the conformational (sub)states sampling of Ras. This, in turn, is of crucial importance in the activation/deactivation cycle of Ras, due to the binding of guanine nucleotide exchange factor proteins (GEFs)/GTPase-activating proteins (GAPs).

Conformational ensemble-dependent lipid recognition and segregation by prenylated intrinsically disordered regions in small GTPases

We studied diverse prenylated intrinsically disordered regions (PIDRs) of Ras and Rho family small GTPases using long timescale atomistic molecular dynamics simulations in an asymmetric model membrane of phosphatidylcholine (PC) and phosphatidylserine (PS) lipids. Here we show that conformational plasticity is a key determinant of lipid sorting by polybasic PIDRs and provide evidence for lipid sorting based on both headgroup and acyl chain structures. We further show that conformational ensemble-based lipid recognition is generalizable to all polybasic PIDRs, and that the sequence outside the polybasic domain (PBD) modulates the conformational plasticity, bilayer adsorption, and interactions of PIDRs with membrane lipids. Specifically, we find that palmitoylation, the ratio of basic to acidic residues, and the hydrophobic content of the sequence outside the PBD significantly impact the diversity of conformational substates and hence the extent of conformation-dependent lipid interactions. We thus propose that the PBD is required but not sufficient for the full realization of lipid sorting by prenylated PBD-containing membrane anchors, and that the membrane anchor is not only responsible for high affinity membrane binding but also directs the protein to the right target membrane where it participates in lipid sorting.

Conformational ensemble dependent lipid recognition and segregation by prenylated intrinsically disordered regions in small GTPases

We studied diverse prenylated intrinsically disordered regions (PIDRs) of Ras and Rho family small GTPases using long timescale atomistic molecular dynamics simulations in an asymmetric model membrane of phosphatidylcholine (PC) and phosphatidylserine (PS) lipids. We show that conformational plasticity is a key determinant of lipid sorting by polybasic PIDRs and provide evidence for lipid sorting based on both headgroup and acyl chain structures. We further show that conformational ensemble-based lipid recognition is generalizable to all polybasic PIDRs, and that the sequence outside the polybasic domain (PBD) modulates the conformational plasticity, bilayer adsorption, and interactions of PIDRs with membrane lipids. Specifically, we found that palmitoylation, the ratio of basic to acidic residues, and the hydrophobic content of the sequence outside the PBD significantly impact the diversity of conformational substates and hence the extent of conformation-dependent lipid interactions. We thus propose that the PBD is required but not sufficient for the full realization of lipid sorting by prenylated PBD-containing membrane anchors, and that the membrane anchor is not only responsible for high affinity membrane binding but also directs the protein to the right target membrane where it participates in lipid sorting.

Engineered biomembrane-derived nanoparticles for nanoscale theranostics

Currently, biological membrane-derived nanoparticles (NPs) have shown enormous potential as drug delivery vehicles due to their outstanding biomimetic properties. To make these NPs more adaptive to complex biological systems, some methods have been developed to modify biomembranes and endow them with more functions while preserving their inherent natures. In this review, we introduce five common approaches used for biomembrane decoration: membrane hybridization, the postinsertion method, chemical methods, metabolism engineering and gene engineering. These methods can functionalize a series of biomembranes derived from red blood cells, white blood cells, tumor cells, platelets, exosomes and so on. Biomembrane engineering could markedly facilitate the targeted drug delivery, treatment and diagnosis of cancer, inflammation, immunological diseases, bone diseases and Alzheimer's disease. It is anticipated that these membrane modification techniques will advance biomembrane-derived NPs into broader applications in the future.

Remodeling of the Plasma Membrane by Surface-Bound Protein Monomers and Oligomers: The Critical Role of Intrinsically Disordered Regions

The plasma membrane (PM) of cells is a dynamic structure whose morphology and composition is in constant flux. PM morphologic changes are particularly relevant for the assembly and disassembly of signaling platforms involving surface-bound signaling proteins, as well as for many other mechanochemical processes that occur at the PM surface. Surface-bound membrane proteins (SBMP) require efficient association with the PM for their function, which is often achieved by the coordinated interactions of intrinsically disordered regions (IDRs) and globular domains with membrane lipids. This review focuses on the role of IDR-containing SBMPs in remodeling the composition and curvature of the PM. The ability of IDR-bearing SBMPs to remodel the Gaussian and mean curvature energies of the PM is intimately linked to their ability to sort subsets of phospholipids into nanoclusters. We therefore discuss how IDRs of many SBMPs encode lipid-binding specificity or facilitate cluster formation, both of which increase their membrane remodeling capacity, and how SBMP oligomers alter membrane shape by monolayer surface area expansion and molecular crowding.

Phosphatidylserine and Phosphatidylethanolamine Asymmetry Have a Negligible Effect on the Global Structure, Dynamics, and Interactions of the KRAS Lipid Anchor

The intrinsically disordered C-terminus of the prominent oncogenic protein KRAS-4B (KRAS) selectively interacts and clusters with phosphatidylserine (PS) lipids in the plasma membrane (PM). This 11-residue segment, called tK, contains a polybasic domain (PBD) of six contiguous lysine residues and a farnesylated cysteine. Previous molecular dynamics (MD) simulation studies of tK in phosphatidylcholine (PC)/PS bilayers have suggested that backbone conformational dynamics modulate tK-PS interactions. These simulations have been conducted in symmetric membranes whereas the PM is compositionally asymmetric, with the inner leaflet, where KRAS is localized, being enriched with PS and phosphatidylethanolamine (PE) lipids. To examine if bilayer asymmetry affects tK conformational dynamics and interaction with lipids, we conducted two 10 μs long MD simulations of tK bound to a PC/PS and a PC/PS/PE bilayer in which the PS and PE lipids are distributed in one leaflet. We found that, first, these compositional asymmetries caused differences in acyl chain dynamics between leaflets, but the equilibrium structural and dynamic properties of the two asymmetric bilayers are similar; second, in both systems tK is highly dynamic and samples at least two distinct conformational states; third, PS-tK hydrogen-bonding interactions vary with peptide backbone conformations, and lysine side chains in the PBD predominantly interact with the serine oxygens of PS. These results are in good agreement with previous observations of tK in symmetric membranes. The effects of POPS asymmetry or the presence of POPE on tK are limited to modulating the relative contribution of individual side chains to interactions with lipids and redistributing conformational substates. Additional observations include the larger flexibility of tK in the current simulations, which we attribute to the longer duration of the simulations and the use of the CHARMM36m force field, which more accurately models intrinsically disordered peptides such as tK.

Elucidation of the key role of Pt···Pt interactions in the directional self-assembly of platinum(II) complexes

Molecular self-assembly provides a bottom-up platform to design supramolecular functional materials, attracting numerous interests in material sciences. The utilization of platinum(II) complexes as building blocks of supramolecular assemblies opens up the unique noncovalent Pt···Pt interaction as one of the driving forces, imparting the supramolecular materials with rich spectroscopic features. However, the exact role of Pt···Pt interactions in molecular assembly remains elusive. The current study combines experimental and computational techniques to elucidate the role of Pt···Pt interactions in the self-assembly process of a representative amphiphilic platinum(II) complex. This work demonstrates the directional role of Pt···Pt interactions in assisting the molecular assembly in an anisotropic manner, achieving the formation of ordered self-assembled structures.

Here, we report the use of an amphiphilic Pt(II) complex, K[Pt{(O₃SCH₂CH₂CH₂)₂bzimpy}Cl] (PtB), as a model to elucidate the key role of Pt···Pt interactions in directing self-assembly by combining temperature-dependent ultraviolet-visible (UV-Vis) spectroscopy, stopped-flow kinetic experiments, quantum mechanics (QM) calculations, and molecular dynamics (MD) simulations. Interestingly, we found that the self-assembly mechanism of PtB in aqueous solution follows a nucleation-free isodesmic model, as revealed by the temperature-dependent UV-Vis experiments. In contrast, a cooperative growth is found for the self-assembly of PtB in acetone–water (7:1, vol/vol) solution, which is further verified by the stopped-flow experiments, which clearly indicates the existence of a nucleation phase in the acetone–water (7:1, vol/vol) solution. To reveal the underlying reasons and driving forces for these self-assembly processes, we performed QM calculations and show that the Pt···Pt interactions arising from the interaction between the p_z and d_z² orbitals play a crucial role in determining the formation of ordered self-assembled structures. In subsequent oligomer MD simulations, we demonstrate that this directional Pt···Pt interaction can indeed facilitate the formation of linear structures packed in a helix-like fashion. Our results suggest that the self-assembly of PtB in acetone–water (7:1, vol/vol) solution is predominantly driven by the directional noncovalent Pt···Pt interaction, leading to the cooperative growth and the formation of fibrous nanostructures. On the contrary, the self-assembly in aqueous solution forms spherical nanostructures of PtB, which is primarily due to the predominant contribution from the less directional hydrophobic interactions over the directional Pt···Pt and π−π interactions that result in an isodesmic growth.

Interactions between Pea Protein Isolate and Carboxymethylcellulose in Neutral and Acid Aqueous Systems

Pea protein isolate (PPI), as an emerging plant protein, has gradually aroused the attention of the public, but the PPI, especially high-concentration PPI's low stability in the acidic aqueous system, was still a problem that limited its application. In this research, we investigated the interactions between relatively high concentrations of PPI (3.0%) and carboxymethylcellulose (CMC, 0-0.5%) in neutral and acid aqueous systems to explore the change of the phase behavior and stability of PPI as affected by CMC. It showed that the stability of PPI in the aqueous systems strongly depended on the CMC concentration, especially at the acidic aqueous systems. At neutral aqueous system, a certain amount addition of CMC into the PPI caused serious phase separation. While stable PPI solutions can be obtained at a narrow region around pH 4.5 to 5.5 by adding different amounts of CMC. The enhancement in the electrostatic repulsion and steric hindrance between the newly formed PPI-CMC biopolymers, as well as the increase in bulk viscosity with the adding of CMC at pH 4.5, contributed to the higher stability of PPI in acidic aqueous systems.

Programming Cell-Cell Interactions through Non-genetic Membrane Engineering

The ability to direct targeted intercellular interactions has the potential to enable and expand the use of cell-based therapies for regenerative medicine, tissue engineering, and immunotherapy. While genetic engineering approaches have proven effective, these techniques are not amenable to all cell types and often yield permanent modifications with potentially long-lasting adverse effects, restricting their application. To circumvent these limitations, there is intense interest in developing non-genetic methods to modify cell membranes with functional groups that will enable the recognition of target cells. While many such techniques have been developed, relatively few have been applied to directing specific cell-cell interactions. This review details these non-genetic membrane engineering approaches-namely, hydrophobic membrane insertion, chemical modification, liposome fusion, metabolic engineering, and enzymatic remodeling-and summarizes their major applications. Based on this analysis, perspective is provided on the ideal features of these systems with an emphasis on the potential for clinical translation.

Interaction of cyclic and linear Labaditin peptides with anionic and zwitterionic micelles

Conformational changes of the cyclic (Lo) peptide Labaditin (VWTVWGTIAG) and its linear analogue (L1) promoted by presence of anionic sodium dodecyl sulfate (SDS) and zwitterionic L-α-Lysophosphatidylcholine (LPC) micelles were investigated. Results from λ(max) blue-shift of tryptophan fluorescence emission combined with Stern-Volmer constants values and molecular dynamics (MD) simulations indicated that L1 interacts with SDS micelles to a higher extent than does Lo. Further, the MD simulation demonstrated that both Lo and L1 interact similarly with LPC micelles, being preferentially located at the micelle/water interface. The peptide-micelle interaction elicits conformational changes in the peptides. Lo undergoes limited modifications and presents unordered structure in both LPC and SDS micelles. On the other hand, L1 displays a random-coil structure in aqueous medium, pH 7.0, and it acquires a β-structure upon interaction with SDS and LPC, albeit with structural differences in each medium.

Overview of simulation studies on the enzymatic activity and conformational dynamics of the GTPase Ras

Over the last 40 years, we have learnt a great deal about the Ras onco-proteins. These intracellular molecular switches are essential for the function of a variety of physiological processes, including signal transduction cascades responsible for cell growth and proliferation. Molecular simulations and free energy calculations have played an essential role in elucidating the conformational dynamics and energetics underlying the GTP hydrolysis reaction catalysed by Ras. Here we present an overview of the main lessons from molecular simulations on the GTPase reaction and conformational dynamics of this important anti-cancer drug target. In the first part, we summarise insights from quantum mechanical and combined quantum mechanical/molecular mechanical simulations as well as other free energy methods and highlight consensus viewpoints as well as remaining controversies. The second part provides a very brief overview of new insights emerging from large-scale molecular dynamics simulations. We conclude with a perspective regarding future studies of Ras where computational approaches will likely play an active role.

Deformation of a Two-domain Lipid Bilayer due to Asymmetric Insertion of Lipid-modified Ras Peptides

Ras proteins are attached to the inner leaflet of the plasma membrane via a lipid-modified anchor. Membrane-bound Ras proteins laterally segregate into nanoscale signaling platforms called nanoclusters. It has been shown that the membrane domain preference of Ras nanoclusters varies with the nature of lipidation but their effect on the membrane has not been well understood. To investigate the effect of Ras insertion on membrane structure, we carried out numerous coarse-grained molecular dynamics (CGMD) simulations on a two-domain DPPC/DLiPC/cholesterol lipid bilayer in which different number and type of H-ras peptides were attached on one side. We have shown previously that this lipid mixture forms co-existing liquid-ordered/liquid-disordered (L o /L d ) domains and that different H-ras peptides form clusters that variously accumulate at the L o /L d regions or the boundary between them. Here we show that asymmetric insertion of each of these peptides induces a vertical relative displacement of the domains and deforms the bilayer, with the domain boundary serving as the center of deformation. The extent of the deformation, however, varies with the type and number of lipid modification. This is because the number and type of the Ras lipid tails determines the degree to which the stress caused by asymmetric peptide insertion is relieved by inter-leaflet cholesterol transfer and lipid tilt. In addition, we have characterized the mechanism of bilayer deformation based on the collective effect of the Ras peptides on inter-leaflet surface area, pressure profile and line tension differences. This allowed us to elucidate how Ras lipid modification affects membrane geometry and how a two-domain bilayer adjusts its shape through boundary deformation. The result contributes to a better understanding of Ras signaling platforms and highlights some of the mechanisms by which a multi-domain membrane responds to external perturbation.

Lessons from computer simulations of Ras proteins in solution and in membrane

BACKGROUND

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

SCOPE OF REVIEW

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

MAJOR CONCLUSIONS

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

GENERAL SIGNIFICANCE

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.

Aggregation of lipid-anchored full-length H-Ras in lipid bilayers: simulations with the MARTINI force field

Lipid-anchored Ras oncoproteins assemble into transient, nano-sized substructures on the plasma membrane. These substructures, called nanoclusters, were proposed to be crucial for high-fidelity signal transmission in cells. However, the molecular basis of Ras nanoclustering is poorly understood. In this work, we used coarse-grained (CG) molecular dynamics simulations to investigate the molecular mechanism by which full-length H-ras proteins form nanoclusters in a model membrane. We chose two different conformations of H-ras that were proposed to represent the active and inactive state of the protein, and a domain-forming model bilayer made up of di16:0-PC (DPPC), di18:2-PC (DLiPC) and cholesterol. We found that, irrespective of the initial conformation, Ras molecules assembled into a single large aggregate. However, the two binding modes, which are characterized by the different orientation of the G-domain with respect to the membrane, differ in dynamics and organization during and after aggregation. Some of these differences involve regions of Ras that are important for effector/modulator binding, which may partly explain observed differences in the ability of active and inactive H-ras nanoclusters to recruit effectors. The simulations also revealed some limitations in the CG force field to study protein assembly in solution, which we discuss in the context of proposed potential avenues of improvement.

Staurosporines disrupt phosphatidylserine trafficking and mislocalize Ras proteins

Oncogenic mutant Ras is frequently expressed in human cancers, but no anti-Ras drugs have been developed. Since membrane association is essential for Ras biological activity, we developed a high content assay for inhibitors of Ras plasma membrane localization. We discovered that staurosporine and analogs potently inhibit Ras plasma membrane binding by blocking endosomal recycling of phosphatidylserine, resulting in redistribution of phosphatidylserine from plasma membrane to endomembrane. Staurosporines are more active against K-Ras than H-Ras. K-Ras is displaced to endosomes and undergoes proteasomal-independent degradation, whereas H-Ras redistributes to the Golgi and is not degraded. K-Ras nanoclustering on the plasma membrane is also inhibited. Ras mislocalization does not correlate with protein kinase C inhibition or induction of apoptosis. Staurosporines selectively abrogate K-Ras signaling and proliferation of K-Ras-transformed cells. These results identify staurosporines as novel inhibitors of phosphatidylserine trafficking, yield new insights into the role of phosphatidylserine and electrostatics in Ras plasma membrane targeting, and validate a new target for anti-Ras therapeutics.

Formation and domain partitioning of H-ras peptide nanoclusters: effects of peptide concentration and lipid composition

Experiments have shown that homologous Ras proteins containing different lipid modification, which is required for membrane binding, form nonoverlapping nanoclusters on the plasma membrane. However, the physical basis for clustering and lateral organization remains poorly understood. We have begun to tackle this issue using coarse-grained molecular dynamics simulations of the H-ras lipid anchor (tH), a triply lipid-modified heptapeptide embedded in a domain-forming mixed lipid bilayer [Janosi L. et al. Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 8097]. Here we use the same simulation approach to investigate the effect of peptide concentration and bilayer composition on the clustering and lateral distribution of tH. We found no major difference in the clustering behavior of tH above a certain concentration. However, the simulations predict the existence of a critical concentration below which tH does not form nanoclusters. Moreover, our data demonstrate that cholesterol enhances the stability of tH nanoclusters but is not required for their formation. Finally, analyses of peptide distributions and partition free energies allowed us to quantitatively describe how clustering facilitates the accumulation of tH at the interface between ordered and disordered domains of the simulated bilayer systems. These thermodynamic insights represent some of the key elements for a comprehensive understanding of the molecular basis for the formation and stability of Ras signaling platforms.

What drives the clustering of membrane-bound Ras?

The dynamic assembly and lateral organization of Ras proteins on the plasma membrane has been the focus of much research in recent years. It has been shown that different isoforms of Ras proteins share a nearly identical catalytic domain, yet form distinct and non-overlapping nanoclusters. Though this difference in the clustering behavior of Ras proteins has been attributed largely to their different C terminal lipid modification, its precise physical basis was not determined. Recently, we used computer simulations to study the mechanism by which the triply lipid-modified membrane-anchor of H-ras, and its partially de-lipidated variants, form nanoclusters in a model lipid bilayer. We found that the specific nature of the lipid modification is less important for cluster formation, but plays a key role for the domain-specific distribution of the nanoclusters. Here we provide additional details on the interplay between bilayer structure perturbation and peptide-peptide association that provide the physical driving force for clustering. We present some thoughts about how enthalpic (i.e., interaction) and entropic effects might regulate nanocluster size and stability.

Segregation of negatively charged phospholipids by the polycationic and farnesylated membrane anchor of Kras

The Kras protein, a member of the Ras family of bio-switches that are frequently mutated in cancer and developmental disorders, becomes functional when anchored to the inner surface of the plasma membrane. It is well known that membrane attachment involves the farnesylated and poylcationic C-terminus of the protein. However, little is known about the structure of the complex and the specific protein-lipid interactions that are responsible for the binding. On the basis of data from extensive (>0.55 μs) molecular dynamics simulations of multiple Kras anchors in bilayers of POPC/POPG lipids (4:1 ratio), we show that, as expected, Kras is tethered to the bilayer surface by specific lysine-POPG salt bridges and by nonspecific farnesyl-phospholipid van der Waals interactions. Unexpectedly, however, only the C-terminal five of the eight Kras Lys side chains were found to directly interact with the bilayer, with the N-terminal ones staying in water. Furthermore, the positively charged Kras anchors pull the negatively charged POPG lipids together, leading to the clustering of the POPG lipids around the proteins. This selective Kras-POPG interaction is directly related to the specific geometry of the backbone, which exists in two major conformational states: 1), a stable native-like ensemble of structures characterized by an extended geometry with a pseudohelical turn; and 2), less stable nonnative ensembles of conformers characterized by severely bent geometries. Finally, although the interface-bound anchor has little effect on the overall structure of the bilayer, it induces local thinning within a persistence length of ∼12 Å. Our results thus go beyond documenting how Kras attaches to a mixed bilayer of charged and neutral lipids; they highlight a fascinating process of protein-induced lipid sorting coupled with the (re)shaping of a surface-bound protein by the host lipids.

An antibody as surrogate receptor reveals determinants of activity of an innate immune peptide antibiotic

Drug discovery initiatives often depend critically on knowledge of ligand-receptor interactions. However, the identity or structure of the target receptor may not be known in every instance. The concept of receptor surrogate, a molecular environment mimic of natural receptor, may prove beneficial under such circumstances. Here, we demonstrate the potential of monoclonal antibodies (mAbs) to act as surrogate receptors for a class of innate immune peptide antibiotics, a strategy that can help comprehend their action mechanism and identify chemical entities crucial for activity. A panel of antibody surrogates was raised against indolicidin, a tryptophan-rich cationic broad spectrum antimicrobial peptide of innate immune origin. Employing an elegant combination of thermodynamics, crystallography, and molecular modeling, interactions of the peptide with a high affinity anti-indolicidin monoclonal antibody were analyzed and were used to identify a motif that contained almost the entire antibiotic activity of native indolicidin. The analysis clarified the interaction of the peptide with previously proposed targets such as bacterial cell membrane and DNA and could further be correlated with antimicrobial compounds whose actions involve varied other mechanisms. These features suggest a multipronged assault pathway for indolicidin. Remarkably, the anti-indolicidin mAb surrogate was able to isolate additional independent bactericidal sequences from a random peptide library, providing compelling evidence as to the physiological relevance of surrogate receptor concept and suggesting applications in receptor-based pharmacophore research.

Induced fit, conformational selection and independent dynamic segments: an extended view of binding events

Single molecule and NMR measurements of protein dynamics increasingly uncover the complexity of binding scenarios. Here, we describe an extended conformational selection model that embraces a repertoire of selection and adjustment processes. Induced fit can be viewed as a subset of this repertoire, whose contribution is affected by the bond types stabilizing the interaction and the differences between the interacting partners. We argue that protein segments whose dynamics are distinct from the rest of the protein ('discrete breathers') can govern conformational transitions and allosteric propagation that accompany binding processes and, as such, might be more sensitive to mutational events. Additionally, we highlight the dynamic complexity of binding scenarios as they relate to events such as aggregation and signalling, and the crowded cellular environment.

Absolute Single-Molecule Entropies from Quasi-Harmonic Analysis of Microsecond Molecular Dynamics: Correction Terms and Convergence Properties

The convergence properties of the absolute single-molecule configurational entropy and the correction terms used to estimate it are investigated using microsecond molecular dynamics simulation of a peptide test system and an improved methodology. The results are compared with previous applications for systems of diverse chemical nature. It is shown that (i) the effect of anharmonicity is small, (ii) the effect of pairwise correlation is typically large, and (iii) the latter affects to a larger extent the entropy estimate of thermodynamic states characterized by a higher motional correlation. The causes of such deviations from a quasi-harmonic behavior are explained. This improved approach provides entropies also for molecular systems undergoing conformational transitions and characterized by highly frustrated energy surfaces, thus not limited to systems sampling a single quasi-harmonic basin. Overall, this study emphasizes the need for extensive phase-space sampling in order to obtain a reliable estimation of entropic contributions.