Gamma-Hemolysin Components: Computational Strategies for LukF-Hlg2 Dimer Reconstruction on a Model Membrane.

The gamma-hemolysin protein is one of the most common pore-forming toxins expressed by the pathogenic bacterium Staphylococcus aureus. The toxin is used by the pathogen to escape the immune system of the host organism, by assembling into octameric transmembrane pores on the surface of the target immune cell and leading to its death by leakage or apoptosis. Despite the high potential risks associated with Staphylococcus aureus infections and the urgent need for new treatments, several aspects of the pore-formation process from gamma-hemolysin are still unclear. These include the identification of the interactions between the individual monomers that lead to the formation of a dimer on the cell membrane, which represents the unit for further oligomerization. Here, we employed a combination of all-atom explicit solvent molecular dynamics simulations and protein-protein docking to determine the stabilizing contacts that guide the formation of a functional dimer. The simulations and the molecular modeling reveal the importance of the flexibility of specific protein domains, in particular the N-terminus, to drive the formation of the correct dimerization interface through functional contacts between the monomers. The results obtained are compared with the experimental data available in the literature.

Fast, Accurate, and System-Specific Variable-Resolution Modeling of Proteins.

In recent years, a few multiple-resolution modeling strategies have been proposed, in which functionally relevant parts of a biomolecule are described with atomistic resolution, with the remainder of the system being concurrently treated using a coarse-grained model. In most cases, the parametrization of the latter requires lengthy reference all-atom simulations and/or the usage of off-shelf coarse-grained force fields, whose interactions have to be refined to fit the specific system under examination. Here, we overcome these limitations through a novel multiresolution modeling scheme for proteins, dubbed coarse-grained anisotropic network model for variable resolution simulations, or CANVAS. This scheme enables a user-defined modulation of the resolution level throughout the system structure; a fast parametrization of the potential without the necessity of reference simulations; and the straightforward usage of the model on the most commonly used molecular dynamics platforms. The method is presented and validated with two case studies, the enzyme adenylate kinase and the therapeutic antibody pembrolizumab, by comparing the results obtained with the CANVAS model against fully atomistic simulations. The modeling software, implemented in Python, is made freely available for the community on a collaborative github repository.

Membrane binding of pore-forming γ-hemolysin components studied at different lipid compositions.

Methicillin-resistant Staphylococcus aureus is among those pathogens currently posing the highest threat to public health. Its host immune evasion strategy is mediated by pore-forming toxins (PFTs), among which the bi-component γ-hemolysin is one of the most common. The complexity of the porogenesis mechanism by γ-hemolysin poses difficulties in the development of antivirulence therapies targeting PFTs from S. aureus, and sparse and apparently contrasting experimental data have been produced. Here, through a large set of molecular dynamics simulations at different levels of resolution, we investigate the first step of pore formation, and in particular the effect of membrane composition on the ability of γ-hemolysin components, LukF and Hlg2, to steadily adhere to the lipid bilayer in the absence of proteinaceous receptors. Our simulations are in agreement with experimental data of γ-hemolysin pore formation on model membranes, which are here explained on the basis of the bilayer properties. Our computational investigation suggests a possible rationale to explain experimental data on phospholipid binding to the LukF component, and to hypothesise a mechanism by which, on purely lipidic bilayers, the stable anchoring of LukF to the cell surface facilitates Hlg2 binding, through the exposure of its N-terminal region. We expect that further insights on the mechanism of transition between soluble and membrane bound-forms and on the role played by the lipid molecules will contribute to the design of antivirulence agents with enhanced efficacy against methicillin-resistant S. aureus infections.

Communication pathways bridge local and global conformations in an IgG4 antibody.

The affinity of an antibody for its antigen is primarily determined by the specific sequence and structural arrangement of the complementarity-determining regions (CDRs). Recent evidence, however, points toward a nontrivial relation between the CDR and distal sites: variations in the binding strengths have been observed upon mutating residues separated from the paratope by several nanometers, thus suggesting the existence of a communication network within antibodies, whose extension and relevance might be deeper than insofar expected. In this work, we test this hypothesis by means of molecular dynamics (MD) simulations of the IgG4 monoclonal antibody pembrolizumab, an approved drug that targets the programmed cell death protein 1 (PD-1). The molecule is simulated in both the apo and holo states, totalling 4 µs of MD trajectory. The analysis of these simulations shows that the bound antibody explores a restricted range of conformations with respect to the apo one, and that the global conformation of the molecule correlates with that of the CDR. These results support the hypothesis that pembrolizumab featues a multi-scale hierarchy of intertwined global and local conformational changes. The analysis pipeline developed in this work is general, and it can help shed further light on the mechanistic aspects of antibody function.

Comparative Molecular Dynamics Investigation of the Electromotile Hearing Protein Prestin.

The mammalian protein prestin is expressed in the lateral membrane wall of the cochlear hair outer cells and is responsible for the electromotile response of the basolateral membrane, following hyperpolarisation or depolarisation of the cells. Its impairment marks the onset of severe diseases, like non-syndromic deafness. Several studies have pointed out possible key roles of residues located in the Transmembrane Domain (TMD) that differentiate mammalian prestins as incomplete transporters from the other proteins belonging to the same solute-carrier (SLC) superfamily, which are classified as complete transporters. Here, we exploit the homology of a prototypical incomplete transporter (rat prestin, rPres) and a complete transporter (zebrafish prestin, zPres) with target structures in the outward open and inward open conformations. The resulting models are then embedded in a model membrane and investigated via a rigorous molecular dynamics simulation protocol. The resulting trajectories are analyzed to obtain quantitative descriptors of the equilibration phase and to assess a structural comparison between proteins in different states, and between different proteins in the same state. Our study clearly identifies a network of key residues at the interface between the gate and the core domains of prestin that might be responsible for the conformational change observed in complete transporters and hindered in incomplete transporters. In addition, we study the pathway of Cl- ions in the presence of an applied electric field towards their putative binding site in the gate domain. Based on our simulations, we propose a tilt and shift mechanism of the helices surrounding the ion binding cavity as the working principle of the reported conformational changes in complete transporters.

From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules.

The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.

Open-Boundary Molecular Mechanics/Coarse-Grained Framework for Simulations of Low-Resolution G-Protein-Coupled Receptor-Ligand Complexes.

G-protein-coupled receptors (GPCRs) constitute as much as 30% of the overall proteins targeted by FDA-approved drugs. However, paucity of structural experimental information and low sequence identity between members of the family impair the reliability of traditional docking approaches and atomistic molecular dynamics simulations for in silico pharmacological applications. We present here a dual-resolution approach tailored for such low-resolution models. It couples a hybrid molecular mechanics/coarse-grained (MM/CG) scheme, previously developed by us for GPCR-ligand complexes, with a Hamiltonian-based adaptive resolution scheme (H-AdResS) for the solvent. This dual-resolution approach removes potentially inaccurate atomistic details from the model while building a rigorous statistical ensemble-the grand canonical one-in the high-resolution region. We validate the method on a well-studied GPCR-ligand complex, for which the 3D structure is known, against atomistic simulations. This implementation paves the way for future accurate in silico studies of low-resolution ligand/GPCRs models.

Open Boundary Simulations of Proteins and Their Hydration Shells by Hamiltonian Adaptive Resolution Scheme.

The recently proposed Hamiltonian adaptive resolution scheme (H-AdResS) allows the performance of molecular simulations in an open boundary framework. It allows changing, on the fly, the resolution of specific subsets of molecules (usually the solvent), which are free to diffuse between the atomistic region and the coarse-grained reservoir. So far, the method has been successfully applied to pure liquids. Coupling the H-AdResS methodology to hybrid models of proteins, such as the molecular mechanics/coarse-grained (MM/CG) scheme, is a promising approach for rigorous calculations of ligand binding free energies in low-resolution protein models. Toward this goal, here we apply for the first time H-AdResS to two atomistic proteins in dual-resolution solvent, proving its ability to reproduce structural and dynamic properties of both the proteins and the solvent, as obtained from atomistic simulations.

Just a Flexible Linker? The Structural and Dynamic Properties of CBP-ID4 Revealed by NMR Spectroscopy.

Here, we present a structural and dynamic description of CBP-ID4 at atomic resolution. ID4 is the fourth intrinsically disordered linker of CREB-binding protein (CBP). In spite of the largely disordered nature of CBP-ID4, NMR chemical shifts and relaxation measurements show a significant degree of α-helix sampling in the protein regions encompassing residues 2-25 and 101-128 (1852-1875 and 1951-1978 in full-length CBP). Proline residues are uniformly distributed along the polypeptide, except for the two α-helical regions, indicating that they play an active role in modulating the structural features of this CBP fragment. The two helical regions are lacking known functional motifs, suggesting that they represent thus-far uncharacterized functional modules of CBP. This work provides insights into the functions of this protein linker that may exploit its plasticity to modulate the relative orientations of neighboring folded domains of CBP and fine-tune its interactions with a multitude of partners.