Predicting the Affinity of Peptides to Major Histocompatibility Complex Class II by Scoring Molecular Dynamics Simulations.

Predicting the binding affinity of peptides able to interact with major histocompatibility complex (MHC) molecules is a priority for researchers working in the identification of novel vaccines candidates. Most available approaches are based on the analysis of the sequence of peptides of known experimental affinity. However, for MHC class II receptors, these approaches are not very accurate, due to the intrinsic flexibility of the complex. To overcome these limitations, we propose to estimate the binding affinity of peptides bound to an MHC class II by averaging the score of the configurations from finite-temperature molecular dynamics simulations. The score is estimated for 18 different scoring functions, and we explored the optimal manner for combining them. To test the predictions, we considered eight peptides of known binding affinity. We found that six scoring functions correlate with the experimental ranking of the peptides significantly better than the others. We then assessed a set of techniques for combining the scoring functions by linear regression and logistic regression. We obtained a maximum accuracy of 82% for the predicted sign of the binding affinity using a logistic regression with optimized weights. These results are potentially useful to improve the reliability of in silico protocols to design high-affinity binding peptides for MHC class II receptors

Protein physics by advanced computational techniques: conformational sampling and folded state discrimination

Proteins are essential parts of organisms and participate in virtually every process within cells. Many proteins are enzymes that catalyze biochemical reactions and are vital to metabolism. Proteins also have structural or mechanical functions, such as actin and myosin in muscle that are in charge of motion and locomotion of cells and organisms. Others proteins are important for transporting materials, cell signaling, immune response, and several other functions. Proteins are the main building blocks of life. A protein is a polymer chain of amino acids whose sequence is defined in a gene: three nucleo type basis specify one out of the 20 natural amino acids. All amino acids possess common structural features. They have an \u3b1-carbon to which an amino group, a carboxyl group, a hydrogen atom and a variable side chain are attached. In a protein, the amino acids are linked together by peptide bonds between the carboxyl and amino groups of adjacent residues..

BioEM: GPU-accelerated computing of Bayesian inference of electron microscopy images

In cryo-electron microscopy (EM), molecular structures are determined from large numbers of projection images of individual particles. To harness the full power of this single-molecule information, we use the Bayesian inference of EM (BioEM) formalism. By ranking structural models using posterior probabilities calculated for individual images, BioEM in principle addresses the challenge of working with highly dynamic or heterogeneous systems not easily handled in traditional EM reconstruction. However, the calculation of these posteriors for large numbers of particles and models is computationally demanding. Here we present highly parallelized, GPU-accelerated computer software that performs this task efficiently. Our flexible formulation employs CUDA, OpenMP, and MPI parallelization combined with both CPU and GPU computing. The resulting BioEM software scales nearly ideally both on pure CPU and on CPU+GPU architectures, thus enabling Bayesian analysis of tens of thousands of images in a reasonable time. The general mathematical framework and robust algorithms are not limited to cryo-electron microscopy but can be generalized for electron tomography and other imaging experiments

Microscopic Theory, Analysis, and Interpretation of Conductance Histograms in Molecular Junctions

Molecular electronics break-junction experiments are widely used to investigate fundamental physics and chemistry at the nanoscale. Reproducibility in these experiments relies on measuring conductance on thousands of freshly formed molecular junctions, yielding a broad histogram of conductance events. Experiments typically focus on the most probable conductance, while the information content of the conductance histogram has remained unclear. Here, we develop a theory for the conductance histogram by merging the theory of force-spectroscopy with molecular conductance. To do so, we propose a microscopic model of the junction evolution under the modulation of external mechanical forces and combine it with the statistics of junction rupture and formation. Our formulation focuses on contributions to the conductance dispersion that emerge due to changes in the conductance during mechanical manipulation. The final shape of the histogram is determined by the statistics of junction rupture and formation. The procedure yields analytical equations for the conductance histogram in terms of parameters that describe the free-energy profile of the junction, its mechanical manipulation, and the ability of the molecule to transport charge. All physical parameters that define our microscopic model can be extracted from separate conductance and rupture force measurements on molecular junctions. Our theory accurately fits experimental conductance histograms and augments the information content that can be extracted from experiments. Further, the predicted behavior with respect to physical parameters can be used to design experiments with narrower conductance distribution and to test the range of validity of the model

Pensamiento computacional en la resolución de problemas de ecuaciones cuadráticas en estudiantes de secundaria de una institución pública, Lima – 2022

La presente investigación tuvo como objetivo determinar la influencia del pensamiento computacional en la resolución de problemas de ecuaciones cuadráticas en estudiantes de secundaria de una institución pública, Lima – 2022. El tipo de investigación básica, diseño no experimental, de corte transversal, subtipo, correlacional causal. La muestra estuvo conformada por 110 estudiantes de educación básica obtenida a través del muestreo no probabilístico por conveniencia. Se utilizó como técnica la prueba de rendimiento y como instrumento una prueba. La validez del instrumento se obtuvo a través de un juicio de expertos y para determinar la confiabilidad se utilizó los coeficientes de KR (20) y el Alpha de Cronbach, obteniéndose para la variable pensamiento computacional 0,80; y para la resolución de problemas de ecuaciones cuadráticas en estudiantes de una institución pública, 0,90. Para el procesamiento de datos se utilizó el SPSS_AMOS. Para determinar el grado de influencia entre las variables se utilizó el modelo de ecuaciones estructurales a través del método de estimación de la distribución libre asintótica. Los resultados indican que existe una fuerte influencia positiva y significativa del pensamiento computacional sobre la resolución de problemas de ecuaciones cuadráticas. Asimismo, se evidencia una fuerte influencia positiva y significativa del pensamiento computacional sobre las dimensiones comprender, el problema, configurar un plan, ejecutar el plan y mirar hacia atrás

PARCE: Protocol for Amino acid Refinement through Computational Evolution

The in silico design of peptides and proteins as binders is useful for diagnosis and therapeutics due to their low adverse effects and major specificity. To select the most promising candidates, a key matter is to understand their interactions with protein targets. In this work, we present PARCE, an open source Protocol for Amino acid Refinement through Computational Evolution that implements an advanced and promising method for the design of peptides and proteins. The protocol performs a random mutation in the binder sequence, then samples the bound conformations using molecular dynamics simulations, and evaluates the protein-protein interactions from multiple scoring. Finally, it accepts or rejects the mutation by applying a consensus criterion based on binding scores. The procedure is iterated with the aim to explore efficiently novel sequences with potential better affinities toward their targets. We also provide a tutorial for running and reproducing the methodology

Omicron Mutations Increase Interdomain Interactions and Reduce Epitope Exposure in the SARS-CoV-2 Spike

Omicron BA.1 is a highly infectious variant of SARS-CoV-2 that carries more than thirty mutations on the spike protein in comparison to the Wuhan wild type (WT). Some of the Omicron mutations, located on the receptor binding domain (RBD), are exposed to the surrounding solvent and are known to help evade immunity. However, the impact of buried mutations on the RBD conformations and on the mechanics of the spike opening is less evident. Here, we use all-atom molecular dynamics (MD) simulations with metadynamics to characterize the thermodynamic RBD-opening ensemble, identifying significant differences between WT and Omicron. Specifically, the Omicron mutations S371L, S373P, and S375F make more RBD interdomain contacts during the spike's opening. Moreover, Omicron takes longer to reach the transition state than WT. It stabilizes up-state conformations with fewer RBD epitopes exposed to the solvent, potentially favoring immune or antibody evasion.© 2023 The Author(s)

Impact of Structural Observables From Simulations to Predict the Effect of Single-Point Mutations in MHC Class II Peptide Binders

The prediction of peptide binders to Major Histocompatibility Complex (MHC) class II receptors is of great interest to study autoimmune diseases and for vaccine development. Most approaches predict the affinities using sequence-based models trained on experimental data and multiple alignments from known peptide substrates. However, detecting activity differences caused by single-point mutations is a challenging task. In this work, we used interactions calculated from simulations to build scoring matrices for quickly estimating binding differences by single-point mutations. We modelled a set of 837 peptides bound to an MHC class II allele, and optimized the sampling of the conformations using the Rosetta backrub method by comparing the results to molecular dynamics simulations. From the dynamic trajectories of each complex, we averaged and compared structural observables for each amino acid at each position of the 9°mer peptide core region. With this information, we generated the scoring-matrices to predict the sign of the binding differences. We then compared the performance of the best scoring-matrix to different computational methodologies that range in computational costs. Overall, the prediction of the activity differences caused by single mutated peptides was lower than 60% for all the methods. However, the developed scoring-matrix in combination with existing methods reports an increase in the performance, up to 86% with a scoring method that uses molecular dynamics