Accelerated Protein Folding Using Greedy-Proximal A*

A discard-and-restart MD algorithm for the sampling of protein intermediate states

We introduce a discard-and-restart molecular dynamics (MD) algorithm tailored for the sampling of realistic protein intermediate states. It aids computational structure-based drug discovery by reducing the simulation times to compute a "quick sketch" of folding pathways by up to 2000×. The algorithm iteratively performs short MD simulations and measures their proximity to a target state via a collective variable loss, which can be defined in a flexible fashion, locally or globally. Using the loss, if the trajectory proceeds toward the target, the MD simulation continues. Otherwise, it is discarded, and a new MD simulation is restarted, with new initial velocities randomly drawn from a Maxwell-Boltzmann distribution. The discard-and-restart algorithm demonstrates efficacy and atomistic accuracy in capturing the folding pathways in several contexts: 1) fast-folding small protein domains, 2) the folding intermediate of the prion protein PrP, and 3) the spontaneous partial unfolding of α-tubulin, a crucial event for microtubule severing. During each iteration of the algorithm, we can perform AI-based analysis of the transitory conformations to find potential binding pockets, which could represent druggable sites. Overall, our algorithm enables systematic and computationally efficient exploration of conformational landscapes, enhancing the design of ligands targeting dynamic protein states.

Combining High-Pressure NMR and Geometrical Sampling to Obtain a Full Topological Description of Protein Folding Landscapes: Application to the Folding of Two MAX Effectors from Magnaporthe oryzae

Despite advances in experimental and computational methods, the mechanisms by which an unstructured polypeptide chain regains its unique three-dimensional structure remains one of the main puzzling questions in biology. Single-molecule techniques, ultra-fast perturbation and detection approaches and improvement in all-atom and coarse-grained simulation methods have greatly deepened our understanding of protein folding and the effects of environmental factors on folding landscape. However, a major challenge remains the detailed characterization of the protein folding landscape. Here, we used high hydrostatic pressure 2D NMR spectroscopy to obtain high-resolution experimental structural information in a site-specific manner across the polypeptide sequence and along the folding reaction coordinate. We used this residue-specific information to constrain Cyana3 calculations, in order to obtain a topological description of the entire folding landscape. This approach was used to describe the conformers populating the folding landscape of two small globular proteins, AVR-Pia and AVR-Pib, that belong to the structurally conserved but sequence-unrelated MAX effectors superfamily. Comparing the two folding landscapes, we found that, in spite of their divergent sequences, the folding pathway of these two proteins involves a similar, inescapable, folding intermediate, even if, statistically, the routes used are different.

Deconvoluting virome-wide antibody epitope reactivity profiles

BACKGROUND

Comprehensive characterization of exposures and immune responses to viral infections is critical to a basic understanding of human health and disease. We previously developed the VirScan system, a programmable phage-display technology for profiling antibody binding to a library of peptides designed to span the human virome. Previous VirScan analytical approaches did not carefully account for antibody cross-reactivity among sequences shared by related viruses or for the disproportionate representation of individual viruses in the library.

METHODS

Here we present the AntiViral Antibody Response Deconvolution Algorithm (AVARDA), a multi-module software package for analyzing VirScan datasets. AVARDA provides a probabilistic assessment of infection with species-level resolution by considering sequence alignment of all library peptides to each other and to all human viruses. We employed AVARDA to analyze VirScan data from a cohort of encephalitis patients with either known viral infections or undiagnosed etiologies. We further assessed AVARDA's utility in associating viral infection with type 1 diabetes and lupus.

FINDINGS

By comparing acute and convalescent sera, AVARDA successfully confirmed or detected encephalitis-associated responses to human herpesviruses 1, 3, 4, 5, and 6, improving the rate of diagnosing viral encephalitis in this cohort by 44%. AVARDA analyses of VirScan data from the type 1 diabetes and lupus cohorts implicated enterovirus and herpesvirus infections, respectively.

INTERPRETATION

AVARDA, in combination with VirScan and other pan-pathogen serological techniques, is likely to find broad utility in the epidemiology and diagnosis of infectious diseases.

FUNDING

This work was made possible by support from the National Institutes of Health (NIH), the US Army Research Office, the Singapore Infectious Diseases Initiative (SIDI), the Singapore Ministry of Health's National Medical Research Council (NMRC) and the Singapore National Research Foundation (NRF).

Molecular mechanism of amyloidogenicity and neurotoxicity of a pro-aggregated tau mutant in the presence of histidine tautomerism via replica-exchange simulation

Considering ΔK280 tau mutation, δε isomer with highest sheet content may accelerate aggregation; generating small compounds to inhibit this would help tp prevent tauopathies.