RESUMEN
The characterization of the statistical ensemble of conformations of intrinsically disordered regions (IDRs) is a great challenge both from experimental and computational points of view. In this respect, a number of protocols have been developed using molecular dynamics (MD) simulations to sample the huge conformational space of the molecule. In this work, we consider one of the best methods available, replica exchange solute tempering (REST), as a reference to compare the results obtained using this method with the results obtained using other methods, in terms of experimentally measurable quantities. Along with the methods assessed, we propose here a novel protocol called probabilistic MD chain growth (PMD-CG), which combines the flexible-meccano and hierarchical chain growth methods with the statistical data obtained from tripeptide MD trajectories as the starting point. The system chosen for testing is a 20-residue region from the C-terminal domain of the p53 tumor suppressor protein (p53-CTD). Our results show that PMD-CG provides an ensemble of conformations extremely quickly, after suitable computation of the conformational pool for all peptide triplets of the IDR sequence. The measurable quantities computed on the ensemble of conformations agree well with those based on the REST conformational ensemble.
RESUMEN
The study of proteins with intrinsically disordered regions (IDRs) has emerged as an active field of research due to their intriguing nature. Although IDRs lack a well-defined folded structure, they play important functional roles in cells, following biological mechanisms different from those of the traditional structured proteins. Consequently, it has been necessary to re-design experimental and theoretical methods in order to face the challenges introduced by the dynamic nature of IDRs. In this work, we present an accurate and cost-effective method to study the conformational dynamics of IDRs based on the use of residue-local probabilistic expressions that characterize the conformational ensembles obtained from finite-temperature molecular dynamics (MD) simulations. It is shown that the good performance and the high convergence rates achieved with our method are independent of the IDR lengths, since the method takes advantage of the major influence of the identity and conformation of the nearest residue neighbors on the amino-acid conformational preferences to evaluate the IDR conformational ensembles. This allows us to characterize the conformational space of IDRs using a reduced number of probabilities which can be obtained from comparatively short MD simulations or experimental databases. To exemplify the usefulness of our approach, we present an application to directly detect Molecular Recognition Features (MoRFs) in an IDR domain of the protein p53, and to follow the time evolution of the thermodynamic magnitudes of this system during its exploration of the conformational space.
Asunto(s)
Proteínas Intrínsecamente Desordenadas , Conformación Proteica , Proteínas Intrínsecamente Desordenadas/química , Simulación de Dinámica Molecular , Probabilidad , AminoácidosRESUMEN
Herein, we compared the ability of linear and cyclic peptides generated in silico to target different protein sites: internal pockets and solvent-exposed sites. We selected human lysozyme (HuL) as a model target protein combined with the computational evolution of linear and cyclic peptides. The sequence evolution of these peptides was based on the PARCE algorithm. The generated peptides were screened based on their aqueous solubility and HuL binding affinity. The latter was evaluated by means of scoring functions and atomistic molecular dynamics (MD) trajectories in water, which allowed prediction of the structural features of the protein-peptide complexes. The computational results demonstrated that cyclic peptides constitute the optimal choice for solvent exposed sites, while both linear and cyclic peptides are capable of targeting the HuL pocket effectively. The most promising binders found in silico were investigated experimentally by surface plasmon resonance (SPR), nuclear magnetic resonance (NMR), and electrospray ionization mass spectrometry (ESI-MS) techniques. All tested peptides displayed dissociation constants in the micromolar range, as assessed by SPR; however, both NMR and ESI-MS suggested multiple binding modes, at least for the pocket binding peptides. A detailed NMR analysis confirmed that both linear and cyclic pocket peptides correctly target the binding site they were designed for.
Asunto(s)
Ligandos , Simulación de Dinámica Molecular , Muramidasa/química , Péptidos/química , Algoritmos , Secuencia de Aminoácidos , Sitios de Unión , Muramidasa/metabolismo , Resonancia Magnética Nuclear Biomolecular , Péptidos/metabolismo , Péptidos Cíclicos/química , Péptidos Cíclicos/metabolismo , Unión Proteica , Espectrometría de Masa por Ionización de Electrospray , Resonancia por Plasmón de SuperficieRESUMEN
The bottom-up design of smart nanodevices largely depends on the accuracy by which each of the inherent nanometric components can be functionally designed with predictive methods. Here, we present a rationally designed, self-assembled nanochip capable of capturing a target protein by means of pre-selected binding sites. The sensing elements comprise computationally evolved peptides, designed to target an arbitrarily selected binding site on the surface of beta-2-Microglobulin (ß2m), a globular protein that lacks well-defined pockets. The nanopatterned surface was generated by an atomic force microscopy (AFM)-based, tip force-driven nanolithography technique termed nanografting to construct laterally confined self-assembled nanopatches of single stranded (ss)DNA. These were subsequently associated with an ssDNA-peptide conjugate by means of DNA-directed immobilization, therefore allowing control of the peptide's spatial orientation. We characterized the sensitivity of such peptide-containing systems against ß2m in solution by means of AFM-based differential topographic imaging and surface plasmon resonance (SPR) spectroscopy. Our results show that the confined peptides are capable of specifically capturing ß2m from the surface-liquid interface with micromolar affinity, hence providing a viable proof-of-concept for our approach to peptide design.
Asunto(s)
Biología Computacional/métodos , ADN de Cadena Simple/metabolismo , Péptidos/metabolismo , Microglobulina beta-2/metabolismo , Sitios de Unión/genética , Técnicas Biosensibles/métodos , ADN de Cadena Simple/química , ADN de Cadena Simple/genética , Humanos , Cinética , Microscopía de Fuerza Atómica/métodos , Simulación de Dinámica Molecular , Péptidos/química , Péptidos/genética , Unión Proteica , Resonancia por Plasmón de Superficie/métodos , Microglobulina beta-2/química , Microglobulina beta-2/genéticaRESUMEN
The SERS-based detection of protein sequences with single-residue sensitivity suffers from signal dominance of aromatic amino acid residues and backbones, impeding detection of non-aromatic amino acid residues. Herein, we trap a gold nanoparticle in a plasmonic nanohole to generate a single SERS hot spot for single-molecule detection of 2 similar polypeptides (vasopressin and oxytocin) and 10 distinct amino acids that constitute the 2 polypeptides. Significantly, both aromatic and non-aromatic amino acids are detected and discriminated at the single-molecule level either at individual amino acid molecules or within the polypeptide chains. Correlated with molecular dynamics simulations, our results suggest that the signal dominance due to large spatial occupancy of aromatic rings of the polypeptide sidechains on gold surfaces can be overcome by the high localization of the single hot spot. The superior spectral and spatial discriminative power of our approach can be applied to single-protein analysis, fingerprinting, and sequencing.
Asunto(s)
Aminoácidos/química , Péptidos/química , Espectrometría Raman/métodos , Oro/química , Límite de Detección , Nanopartículas del Metal/química , Simulación de Dinámica MolecularRESUMEN
Mutation protocols are a key tool in computational biophysics for modelling unknown side chain conformations. In particular, these protocols are used to generate the starting structures for molecular dynamics simulations. The accuracy of the initial side chain and backbone placement is crucial to obtain a stable and quickly converging simulation. In this work, we assessed the performance of several mutation protocols in predicting the most probable conformer observed in finite temperature molecular dynamics simulations for a set of protein-peptide crystals differing only by single-point mutations in the peptide sequence. Our results show that several programs which predict well the crystal conformations fail to predict the most probable finite temperature configuration. Methods relying on backbone-dependent rotamer libraries have, in general, a better performance, but even the best protocol fails in predicting approximately 30% of the mutations.
Asunto(s)
Aminoácidos/química , Mutación , Temperatura , Secuencia de Aminoácidos , Modelos Moleculares , Simulación de Dinámica MolecularRESUMEN
Nanobodies offer a viable alternative to antibodies for engineering high affinity binders. Their small size has an additional advantage: it allows exploiting computational protocols for optimizing their biophysical features, such as the binding affinity. The efficient prediction of this quantity is still considered a daunting task especially for modelled complexes. We show how molecular dynamics can successfully assist in the binding affinity prediction of modelled nanobody-protein complexes. The approximate initial configurations obtained by in silico design must undergo large rearrangements before achieving a stable conformation, in which the binding affinity can be meaningfully estimated. The scoring functions developed for the affinity evaluation of crystal structures will provide accurate estimates for modelled binding complexes if the scores are averaged over long finite temperature molecular dynamics simulations.
Asunto(s)
Complejo Antígeno-Anticuerpo/química , Simulación de Dinámica Molecular , Proteínas/inmunología , Anticuerpos de Cadena Única/inmunología , Secuencia de Aminoácidos , Afinidad de Anticuerpos , Complejo Antígeno-Anticuerpo/metabolismo , Humanos , Muramidasa/química , Muramidasa/inmunología , Estructura Terciaria de Proteína , Proteínas/química , Receptor ErbB-2/química , Receptor ErbB-2/metabolismo , Alineación de Secuencia , TemperaturaRESUMEN
Despite the increasing evidence that conformational transitions in peptides and proteins are driven by specific vibrational energy pathways along the molecule, the current experimental techniques of analysis do as yet not allow to study these biophysical processes in terms of anisotropic energy flows. Computational methods offer a complementary approach to obtain a more detailed understanding of the vibrational and conformational dynamics of these systems. Accordingly, in this work we investigate jointly the vibrational energy distribution and the conformational dynamics of trialanine peptide in water solution at room temperature by applying the Instantaneous Normal Mode analysis to the results derived from equilibrium molecular dynamics simulations. It is shown that conformational changes in trialanine are triggered by the vibrational energy accumulated in the low-frequency modes of the molecule, and that excitation is caused exclusively by thermal fluctuations of the solute-solvent system, thus excluding the possibility of an intramolecular vibrational energy redistribution process.
Asunto(s)
Péptidos/química , Termodinámica , Modelos Moleculares , Simulación de Dinámica Molecular , Oligopéptidos/química , Conformación Proteica , Solventes , Vibración , Agua/químicaRESUMEN
The oriented immobilization of proteins, key for the development of novel responsive biomaterials, relies on the availability of effective probes. These are generally provided by standard approaches based on in vivo maturation and in vitro selection of antibodies and/or aptamers. These techniques can suffer technical problems when a non-immunogenic epitope needs to be targeted. Here we propose a strategy to circumvent this issue by in silico design. In our method molecular binders, in the form of cyclic peptides, are computationally evolved by stochastically exploring their sequence and structure space to identify high-affinity peptides for a chosen epitope of a target globular protein: here a solvent-exposed site of ß2-microglobulin (ß2m). Designed sequences were screened by explicit solvent molecular dynamics simulations (MD) followed by experimental validation. Five candidates gave dose-response surface plasmon resonance signals with dissociation constants in the micromolar range. One of them was further analyzed by means of isothermal titration calorimetry, nuclear magnetic resonance, and 250 ns of MD. Atomic-force microscopy imaging showed that this peptide is able to immobilize ß2m on a gold surface. In short, we have shown by a variety of experimental techniques that it is possible to capture a protein through an epitope of choice by computational design.
Asunto(s)
Técnicas de Química Analítica/métodos , Simulación por Computador , Péptidos Cíclicos/química , Proteínas/aislamiento & purificación , Epítopos/química , Modelos Químicos , Simulación de Dinámica Molecular , Péptidos Cíclicos/metabolismoRESUMEN
The chaperonin complex GroEL-GroES is able to accelerate the folding process of knotted proteins considerably. However, the folding mechanism inside the chaperonin cage is elusive. Here we use a combination of lattice and off-lattice Monte Carlo simulations of simple Go models to study the effect of physical confinement and local flexibility on the folding process of protein model systems embedding a trefoil knot in their native structure. This study predicts that steric confinement plays a specific role in the folding of knotted proteins by increasing the knotting probability for very high degrees of confinement. This effect is observed for protein MJ0366 even above the melting temperature for confinement sizes compatible with the size of the GroEL/GroES chaperonin cage. An enhanced local flexibility produces the same qualitative effects on the folding process. In particular, we observe that knotting probability increases up to 40% in the transition state of protein MJ0366 when flexibility is enhanced. This is underlined by a structural change in the transition state, which becomes devoid of helical content. No relation between the knotting mechanism and flexibility was found in the context of the off-lattice model adopted in this work.
Asunto(s)
Proteínas Bacterianas/química , Modelos Moleculares , Proteínas Bacterianas/metabolismo , Chaperonina 60/metabolismo , Cinética , Método de Montecarlo , Pliegue de Proteína , Termodinámica , Temperatura de TransiciónRESUMEN
Using analytical excited-state gradients, vibrational normal modes have been calculated at the minimum of the electronic excited-state potential energy surfaces for a set of extended conjugated molecules with different coupling between them. Molecular model systems composed of units of polyphenylene ethynylene (PPE), polyphenylenevinylene (PPV), and naphthacene/pentacene (NP) have been considered. In all cases except the NP model, the influence of the nonadiabatic coupling on the excited-state equilibrium normal modes is revealed as a unique highest frequency adiabatic vibrational mode that overlaps with the coupling vector. This feature is removed by using a locally diabatic representation in which the effect of NA interaction is removed. Comparison of the original adiabatic modes with a set of vibrational modes computed in the locally diabatic representation demonstrates that the effect of nonadiabaticity is confined to only a few modes. This suggests that the nonadiabatic character of a molecular system may be detected spectroscopically by identifying these unique state-specific high frequency vibrational modes.
Asunto(s)
Modelos Moleculares , Análisis Espectral/métodos , Vibración , Algoritmos , Naftacenos/química , Polímeros/químicaRESUMEN
The folding properties of a protein whose native structure contains a 52 knot are investigated by means of extensive Monte Carlo simulations of a simple lattice model and compared with those of a 31 knot. A 52 knot embedded in the native structure enhances the kinetic stability of the carrier lattice protein in a way that is clearly more pronounced than in the case of the 31 knot. However, this happens at the expense of a severe loss in folding efficiency, an observation that is consistent with the relative abundance of 31 and 52 knots in the Protein Data Bank. The folding mechanism of the 52 knot shares with that of the 31 knot the occurrence of a threading movement of the chain terminus that lays closer to the knotted core. However, co-concomitant knotting and folding in the 52 knot occurs with negligible probability, in sharp contrast to what is observed for the 31 knot. The study of several single point mutations highlights the importance in the folding of knotted proteins of the so-called structural mutations (i.e., energetic perturbations of native interactions between residues that are critical for knotting but not for folding). On the other hand, the present study predicts that mutations that perturb the folding transition state may significantly enhance the kinetic stability of knotted proteins provided they involve residues located within the knotted core.
Asunto(s)
Enfermedades por Prión/genética , Conformación Proteica , Proteínas/química , Cristalografía por Rayos X , Bases de Datos de Proteínas , Cinética , Simulación de Dinámica Molecular , Método de Montecarlo , Enfermedades por Prión/metabolismo , Enfermedades por Prión/patología , Pliegue de Proteína , Proteínas/genética , TermodinámicaRESUMEN
In the framework of the rational design of macromolecules capable of binding to a specific target for biosensing applications, we here further develop an evolutionary protocol designed to optimize the binding affinity of protein binders. In particular we focus on the optimization of the binding portion of small antibody fragments known as nanobodies (or VHH) and choose the hen egg white lysozyme (HEWL) as our target. By implementing a replica exchange scheme for this optimization, we show that an initial hit is not needed and similar solutions can be found by either optimizing an already known anti-HEWL VHH or a randomly selected binder (here a VHH selective towards another macromolecule). While we believe that exhaustive searches of the mutation space are most appropriate when only few key residues have to be optimized, in case a lead binder is not available the proposed evolutionary algorithm should be instead the method of choice.
Asunto(s)
Fragmentos de Inmunoglobulinas , Anticuerpos de Dominio Único , Animales , Fragmentos de Inmunoglobulinas/genética , Mutación , Anticuerpos de Dominio Único/química , PollosRESUMEN
Nanobodies (VHHs) are engineered fragments of the camelid single-chain immunoglobulins. The VHH domain contains the highly variable segments responsible for antigen recognition. VHHs can be easily produced as recombinant proteins. Their small size is a good advantage for in silico approaches. Computer methods represent a valuable strategy for the optimization and improvement of their binding affinity. They also allow for epitope selection offering the possibility to design new VHHs for regions of a target protein that are not naturally immunogenic. Here we present an in silico mutagenic protocol developed to improve the binding affinity of nanobodies together with the first step of their in vitro production. The method, already proven successful in improving the low Kd of a nanobody hit obtained by panning, can be employed for the ex novo design of antibody fragments against selected protein target epitopes.
Asunto(s)
Anticuerpos de Dominio Único , Afinidad de Anticuerpos , Anticuerpos de Dominio Único/química , Epítopos , Proteínas Recombinantes/genéticaRESUMEN
BACKGROUND: P-21-activated kinases (PAKs) are protein serine/threonine kinases, part of the RAS/mitogen-activated protein kinase pathway. PAK1 is highly expressed in the central nervous system and crucially involved in neuronal migration and brain developmental processes. Recently, de novo heterozygous missense variants in PAK1 have been identified as an ultrarare cause of pediatric neurodevelopmental disorders. METHODS: We report a series of children affected with postnatal macrocephaly, neurodevelopmental impairment, and drug-resistant epilepsy. Repeated electroencephalographic (EEG) and video-EEG evaluations were performed over a two- to 10-year period during follow-up to delineate electroclinical histories. Genetic sequencing studies and computational evaluation of the identified variants were performed in our patient cohort. RESULTS: We identified by whole-exome sequencing three novel de novo variants in PAK1 (NM_001128620: c.427A>G, p.Met143Val; c.428T>C, p.Met143Thr; c.428T>A, p.Met143Lys) as the underlying cause of the disease in our families. The three variants affected the same highly conserved Met143 residue within the cysteine-rich inhibitor of PAK1 (CRIPaK) domain, which was identified before as a PAK1 inhibitor target. Computational studies suggested a defective autoinhibition presumably due to impaired PAK1 autoregulation as a result of the recurrent substitution. CONCLUSIONS: We delineated the electroclinical phenotypes of PAK1-related neurological disorders and highlight a novel mutational hotspot that may involve defective autoinhibition of the PAK1 protein. The three novel variants affecting the same hotspot residue within the CRIPaK domain highlight potentially impaired PAK1-CRIPaK interaction as a novel disease mechanism. These findings shed light on possible future treatments targeted at the CRIPaK domain, to modulate PAK1 activity and function.
Asunto(s)
Trastornos del Neurodesarrollo , Quinasas p21 Activadas , Niño , Humanos , Quinasas p21 Activadas/genética , Quinasas p21 Activadas/química , Quinasas p21 Activadas/metabolismo , Proteínas Serina-Treonina Quinasas/genética , Mutación/genética , Trastornos del Neurodesarrollo/genética , Mutación MissenseRESUMEN
The nonadiabatic excited-state molecular dynamics (NA-ESMD) method and excited-state instantaneous normal modes (ES-INMs) analyses have been applied to describe the state-specific vibrations that participate in the unidirectional energy transfer between the coupled chromophores in a branched dendrimeric molecule. Our molecule is composed of two-, three-, and four-ring linear poly(phenyleneethynylene) (PPE) units linked through meta-substitutions. After an initial laser excitation, an ultrafast sequential S(3) â S(2) â S(1) electronic energy transfer from the shortest to longest segment takes place. During each S(n) â S(n-1) (n = 3, 2) transition, ES-INM(S(n)) and ES-INM(S(n-1)) analyses have been performed on S(n) and S(n-1) states, respectively. Our results reveal a unique vibrational mode localized on the S(n) state that significantly matches with the corresponding nonadiabatic coupling vector d(n,(n-1)). This mode also corresponds to the highest frequency ES-INM(S(n)) and it is seen mainly during the electronic transitions. Furthermore, its absence as a unique ES-INM(S(n-1)) reveals that state-specific vibrations play the main role in the efficiency of the unidirectional S(n) â S(n-1) electronic and vibrational energy funneling in light-harvesting dendrimers.
Asunto(s)
Dendrímeros/química , Simulación de Dinámica Molecular , Polímeros/química , Vibración , Transferencia de EnergíaRESUMEN
Antibodies have become the Swiss Army tool for molecular biology and nanotechnology. Their outstanding ability to specifically recognise molecular antigens allows their use in many different applications from medicine to the industry. Moreover, the improvement of conventional structural biology techniques (e.g., X-ray, NMR) as well as the emergence of new ones (e.g., Cryo-EM), have permitted in the last years a notable increase of resolved antibody-antigen structures. This offers a unique opportunity to perform an exhaustive structural analysis of antibody-antigen interfaces by employing the large amount of data available nowadays. To leverage this factor, different geometric as well as chemical descriptors were evaluated to perform a comprehensive characterization.
RESUMEN
Computational peptide design is useful for therapeutics, diagnostics, and vaccine development. To select the most promising peptide candidates, the key is describing accurately the peptide-target interactions at the molecular level. We here review a computational peptide design protocol whose key feature is the use of all-atom explicit solvent molecular dynamics for describing the different peptide-target complexes explored during the optimization. We describe the milestones behind the development of this protocol, which is now implemented in an open-source code called PARCE. We provide a basic tutorial to run the code for an antibody fragment design example. Finally, we describe three additional applications of the method to design peptides for different targets, illustrating the broad scope of the proposed approach.
Asunto(s)
Simulación de Dinámica Molecular , Péptidos , Péptidos/química , SolventesRESUMEN
Protein-protein docking typically consists of the generation of putative binding conformations, which are subsequently ranked by fast heuristic scoring functions. The simplicity of these functions allows for computational efficiency but has severe repercussions on their discrimination capabilities. In this work, we show the effectiveness of suitable descriptors calculated along short scaled molecular dynamics runs in recognizing the nearest-native bound conformation among a set of putative structures generated by the HADDOCK tool for eight protein-protein systems.
Asunto(s)
Simulación de Dinámica Molecular , Conformación Molecular , Simulación del Acoplamiento Molecular , Unión Proteica , Conformación Proteica , Proteínas/metabolismoRESUMEN
In vivo clinical applications of nanobodies (VHHs) require molecules that induce minimal immunoresponse and therefore possess sequences as similar as possible to the human VH domain. Although the relative sequence variability in llama nanobodies has been used to identify scaffolds with partially humanized signature, the transformation of the Camelidae hallmarks in the framework2 still represents a major problem. We assessed a set of mutants in silico and experimentally to elucidate what is the contribution of single residues to the VHH stability and how their combinations affect the mutant nanobody stability. We described at molecular level how the interaction among residues belonging to different structural elements enabled a model llama nanobody (C8WT, isolated from a naïve library) to be functional and maintain its stability, despite the analysis of its primary sequence would classify it as aggregation-prone. Five chimeras formed by grafting CDRs isolated from different nanobodies into C8WT scaffold were successfully expressed as soluble proteins and both tested clones preserved their antigen binding specificity. We identified a nanobody with human hallmarks that seems suitable for humanizing selected camelid VHHs by grafting heterologous CDRs in its scaffold and could serve for the preparation of a synthetic library of human-like single domains.