Searching for Structure: Characterizing the Protein Conformational Landscape with Clustering-Based Algorithms.

The identification and characterization of the main conformations from a protein population are a challenging and inherently high-dimensional problem. Here, we evaluate the performance of the Secondary sTructural Ensembles with machine LeArning (StELa) double-clustering method, which clusters protein structures based on the relationship between the φ and ψ dihedral angles in a protein backbone and the secondary structure of the protein, thus focusing on the local properties of protein structures. The classification of states as vectors composed of the clusters' indices arising naturally from the Ramachandran plot is followed by the hierarchical clustering of the vectors to allow for the identification of the main features of the corresponding free energy landscape (FEL). We compare the performance of StELa with the established root-mean-squared-deviation (RMSD)-based clustering algorithm, which focuses on global properties of protein structures and with Combinatorial Averaged Transient Structure (CATS), the combinatorial averaged transient structure clustering method based on distributions of the φ and ψ dihedral angle coordinates. Using ensembles of conformations from molecular dynamics simulations of intrinsically disordered proteins (IDPs) of various lengths (tau protein fragments) or short fragments from a globular protein, we show that StELa is the clustering method that identifies many of the minima and relevant energy states around the minima from the corresponding FELs. In contrast, the RMSD-based algorithm yields a large number of clusters that usually cover most of the FEL, thus being unable to distinguish between states, while CATS does not sample well the FELs for long IDPs and fragments from globular proteins.

The quaternary question: Determining allostery in spastin through dynamics classification learning and bioinformatics.

The nanomachine from the ATPases associated with various cellular activities superfamily, called spastin, severs microtubules during cellular processes. To characterize the functionally important allostery in spastin, we employed methods from evolutionary information, to graph-based networks, to machine learning applied to atomistic molecular dynamics simulations of spastin in its monomeric and the functional hexameric forms, in the presence or absence of ligands. Feature selection, using machine learning approaches, for transitions between spastin states recognizes all the regions that have been proposed as allosteric or functional in the literature. The analysis of the composition of the Markov State Model macrostates in the spastin monomer, and the analysis of the direction of change in the top machine learning features for the transitions, indicate that the monomer favors the binding of ATP, which primes the regions involved in the formation of the inter-protomer interfaces for binding to other protomer(s). Allosteric path analysis of graph networks, built based on the cross-correlations between residues in simulations, shows that perturbations to a hub specific for the pre-hydrolysis hexamer propagate throughout the structure by passing through two obligatory regions: the ATP binding pocket, and pore loop 3, which connects the substrate binding site to the ATP binding site. Our findings support a model where the changes in the terminal protomers due to the binding of ligands play an active role in the force generation in spastin. The secondary structures in spastin, which are found to be highly degenerative within the network paths, are also critical for feature transitions of the classification models, which can guide the design of allosteric effectors to enhance or block allosteric signaling.

Design, Rationalization, and Automation of a Catalytic Sensing Mechanism for Homogeneous SERS Biosensors.

The current pandemic has shown that we need sensitive and deployable diagnostic technologies. Surface-enhanced Raman scattering (SERS) sensors can be an ideal solution for developing such advanced point-of-need (PON) diagnostic tests. Homogeneous (reagentless) SERS sensors work by directly responding to the target without any processing step, making them capable for simple one-pot assays, but their limitation is the achievable sensitivity, insufficient compared to what is needed for sensing of viral biomarkers. Noncovalent DNA catalysis mechanisms have been recently exploited for catalytic amplification in SERS assays. These advances used catalytic hairpin assembly (CHA) and other DNA self-assembly processes to develop sensing mechanisms with improved sensitivities. However, these mechanisms have not been used in OFF-to-ON homogeneous sensors, and they often target the same biomarker, likely due to the complexity of the mechanism design. There is still a strong need for a catalytic SERS sensor with a homogeneous mechanism and a rationalization of the catalytic sensing mechanism to translate this sensing strategy to different targets and applications. We developed and investigated a homogeneous SERS sensing mechanism that uses catalytic amplification based on DNA self-assembly. We systematically investigated the role of three domains in the fuel strand (internal loop, stem, and toehold), which drives the catalytic mechanism. The thermodynamic parameters determined in our studies were used to build an algorithm for automated design of catalytic sensors that we validated on target sequences associated with malaria and SARS-CoV-2 strains. With our mechanism, we were able to achieve an amplification level of 20-fold for conventional DNA and of 36-fold using locked nucleic acids (LNAs), with corresponding improvements observed in the sensor limit of detection (LOD). We also show a single-base sequence specificity for a sensor targeting a sequence associated with the omicron variant, tested against a delta variant target. This work on catalytic amplification of homogeneous SERS sensors has the potential to enable the use of this sensing modality in new applications, such as infectious disease surveillance, by improving the LOD while conserving the sensor's homogeneous character.

Microtubule Severing Enzymes Oligomerization and Allostery: A Tale of Two Domains.

Severing proteins are nanomachines from the AAA+ (ATPases associated with various cellular activities) superfamily whose function is to remodel the largest cellular filaments, microtubules. The standard AAA+ machines adopt hexameric ring structures for functional reasons, while being primarily monomeric in the absence of the nucleotide. Both major severing proteins, katanin and spastin, are believed to follow this trend. However, studies proposed that they populate lower-order oligomers in the presence of cofactors, which are functionally relevant. Our simulations show that the preferred oligomeric assembly is dependent on the binding partners and on the type of severing protein. Essential dynamics analysis predicts that the stability of an oligomer is dependent on the strength of the interface between the helical bundle domain (HBD) of a monomer and the convex face of the nucleotide binding domain (NBD) of a neighboring monomer. Hot spots analysis found that the region consisting of the HBD tip and the C-terminal (CT) helix is the only common element between the allosteric networks responding to nucleotide, substrate, and intermonomer binding. Clustering analysis indicates the existence of multiple pathways for the transition between the secondary structure of the HBD tip in monomers and the structure(s) it adopts in oligomers.

Exploring the Effect of Mechanical Anisotropy of Protein Structures in the Unfoldase Mechanism of AAA+ Molecular Machines.

Essential cellular processes of microtubule disassembly and protein degradation, which span lengths from tens of µm to nm, are mediated by specialized molecular machines with similar hexameric structure and function. Our molecular simulations at atomistic and coarse-grained scales show that both the microtubule-severing protein spastin and the caseinolytic protease ClpY, accomplish spectacular unfolding of their diverse substrates, a microtubule lattice and dihydrofolate reductase (DHFR), by taking advantage of mechanical anisotropy in these proteins. Unfolding of wild-type DHFR requires disruption of mechanically strong ß-sheet interfaces near each terminal, which yields branched pathways associated with unzipping along soft directions and shearing along strong directions. By contrast, unfolding of circular permutant DHFR variants involves single pathways due to softer mechanical interfaces near terminals, but translocation hindrance can arise from mechanical resistance of partially unfolded intermediates stabilized by ß-sheets. For spastin, optimal severing action initiated by pulling on a tubulin subunit is achieved through specific orientation of the machine versus the substrate (microtubule lattice). Moreover, changes in the strength of the interactions between spastin and a microtubule filament, which can be driven by the tubulin code, lead to drastically different outcomes for the integrity of the hexameric structure of the machine.

Microtubule assembly and disassembly dynamics model: Exploring dynamic instability and identifying features of Microtubules' Growth, Catastrophe, Shortening, and Rescue.

Microtubules (MTs), a cellular structure element, exhibit dynamic instability and can switch stochastically from growth to shortening; but the factors that trigger these processes at the molecular level are not understood. We developed a 3D Microtubule Assembly and Disassembly DYnamics (MADDY) model, based upon a bead-per-monomer representation of the αß-tubulin dimers forming an MT lattice, stabilized by the lateral and longitudinal interactions between tubulin subunits. The model was parameterized against the experimental rates of MT growth and shortening, and pushing forces on the Dam1 protein complex due to protofilaments splaying out. Using the MADDY model, we carried out GPU-accelerated Langevin simulations to access dynamic instability behavior. By applying Machine Learning techniques, we identified the MT characteristics that distinguish simultaneously all four kinetic states: growth, catastrophe, shortening, and rescue. At the cellular 25 µM tubulin concentration, the most important quantities are the MT length L , average longitudinal curvature κ long , MT tip width w , total energy of longitudinal interactions in MT lattice U long , and the energies of longitudinal and lateral interactions required to complete MT to full cylinder U long add and U lat add . At high 250 µM tubulin concentration, the most important characteristics are L , κ long , number of hydrolyzed αß-tubulin dimers n hyd and number of lateral interactions per helical pitch n lat in MT lattice, energy of lateral interactions in MT lattice U lat , and energy of longitudinal interactions in MT tip u long . These results allow greater insights into what brings about kinetic state stability and the transitions between states involved in MT dynamic instability behavior.

Prevalence of early loss of primary molar and its impact in schoolchildren's quality of life.

BACKGROUND: Evidence regarding the impact of premature loss of primary molars on oral health-related quality of life (OHRQoL) is lacking. AIM: To evaluate the prevalence of premature loss of primary molars in schoolchildren and its impact on OHRQoL. DESIGN: We randomly selected 667 children aged 8-9 years enrolled in public and private schools from Aracaju (Brazil). The presence of untreated dental caries and early loss of primary molars were recorded, and its impacts on OHRQoL evaluated using Child Perceptions Questionnaire (CPQ8-10 ). Influence of some variables on tooth loss and OHRQoL was assessed by multilevel Poisson's regression analysis (P < 0.05). RESULTS: The prevalence for early tooth loss was 65.4% (95% confidence interval 51.1%-77.3%). Children with untreated dental caries and attending public schools presented higher prevalence of early tooth loss. Moreover, early loss of primary molars provoked negative impact on OHRQoL considering total score of CPQ8-10 and domains "oral symptoms," "functional limitations" and "emotional well-being." Children presenting early teeth loss presented significantly higher scores of CPQ8-10 , as well as females and children with dental caries in other teeth. CONCLUSIONS: The prevalence of early loss of primary molars is high, and it impacts the OHRQoL of children aged 8-9 years.