Mechanical and thermodynamic properties of Aβ42, Aβ40, and α-synuclein fibrils: a coarse-grained method to complement experimental studies

Durotaxis and Antidurotaxis Droplet Motion onto Gradient Gel-Substrates

The self-sustained motion of fluids on gradient substrates is a spectacular phenomenon, which can be employed and controlled in applications by carefully engineering the substrate properties. Here, we report on a design of a gel substrate with stiffness gradient, which can cause the spontaneous motion of a droplet along (durotaxis) or to the opposite (antidurotaxis) direction of the gradient, depending on the droplet affinity to the substrate. By using extensive molecular dynamics simulations of a coarse-grained model, we find that the mechanisms of the durotaxis and antidurotaxis droplet motion are distinct, require the minimization of the interfacial energy between the droplet and the substrate, and share similarities with those mechanisms previously observed for brush substrates with stiffness gradient. Moreover, durotaxis motion takes place over a wider range of affinities and is generally more efficient (faster motion) than antidurotaxis. Thus, our study points to further possibilities and guidelines for realizing both antidurotaxis and durotaxis motion on the same gradient substrate for applications in microfluidics, energy conservation, and biology.

Wild-Type α-Synuclein Structure and Aggregation: A Comprehensive Coarse-Grained and All-Atom Molecular Dynamics Study

α-Synuclein (α-syn) is a 140 amino acid intrinsically disordered protein (IDP) and the primary component of cytotoxic oligomers implicated in the etiology of Parkinson's disease (PD). While IDPs lack a stable three-dimensional structure, they sample a heterogeneous ensemble of conformations that can, in principle, be assessed through molecular dynamics simulations. However, describing the structure and aggregation of large IDPs is challenging due to force field (FF) accuracy and sampling limitations. To cope with the latter, coarse-grained (CG) FFs emerge as a potential alternative at the expense of atomic detail loss. Whereas CG models can accurately describe the structure of the monomer, less is known about aggregation. The latter is key for assessing aggregation pathways and designing aggregation inhibitor drugs. Herein, we investigate the structure and dynamics of α-syn using different resolution CG (Martini3 and Sirah2) and all-atom (Amber99sb and Charmm36m) FFs to gain insight into the differences and resemblances between these models. The dependence of the magnitude of protein-water interactions and the putative need for enhanced sampling (replica exchange) methods in CG simulations are analyzed to distinguish between force field accuracy and sampling limitations. The stability of the CG models of an α-syn fibril was also investigated. Additionally, α-syn aggregation was studied through umbrella sampling for the CG models and CG/all-atom models for an 11-mer peptide (NACore) from an amyloidogenic domain of α-syn. Our results show that despite the α-syn structures of Martini3 and Sirah2 with enhanced protein-water interactions being similar, major differences exist concerning aggregation. The Martini3 fibril is not stable, and the binding free energy of α-syn and NACore is positive, opposite to Sirah2. Sirah2 peptides in a zwitterionic form, in turn, display termini interactions that are too strong, resulting in end-to-end orientation. Sirah2, with enhanced protein-water interactions and neutral termini, provides, however, a peptide aggregation free energy profile similar to that found with all-atom models. Overall, we find that Sirah2 with enhanced protein-water interactions is suitable for studying protein-protein and protein-drug aggregation.

Exploring Structural Insights of Aβ42 and α-Synuclein Monomers and Heterodimer: A Comparative Study Using Implicit and Explicit Solvent Simulations

Protein misfolding, aggregation, and fibril formation play a central role in the development of severe neurological disorders, including Alzheimer's and Parkinson's diseases. The structural stability of mature fibrils in these diseases is of great importance, as organisms struggle to effectively eliminate amyloid plaques. To address this issue, it is crucial to investigate the early stages of fibril formation when monomers aggregate into small, toxic, and soluble oligomers. However, these structures are inherently disordered, making them challenging to study through experimental approaches. Recently, it has been shown experimentally that amyloid-β 42 (Aβ42) and α-synuclein (α-Syn) can coassemble. This has motivated us to investigate the interaction between their monomers as a first step toward exploring the possibility of forming heterodimeric complexes. In particular, our study involves the utilization of various Amber and CHARMM force-fields, employing both implicit and explicit solvent models in replica exchange and conventional simulation modes. This comprehensive approach allowed us to assess the strengths and weaknesses of these solvent models and force fields in comparison to experimental and theoretical findings, ensuring the highest level of robustness. Our investigations revealed that Aβ42 and α-Syn monomers can indeed form stable heterodimers, and the resulting heterodimeric model exhibits stronger interactions compared to the Aβ42 dimer. The binding of α-Syn to Aβ42 reduces the propensity of Aβ42 to adopt fibril-prone conformations and induces significant changes in its conformational properties. Notably, in AMBER-FB15 and CHARMM36m force fields with the use of explicit solvent, the presence of Aβ42 significantly increases the β-content of α-Syn, consistent with the experiments showing that Aβ42 triggers α-Syn aggregation. Our analysis clearly shows that although the use of implicit solvent resulted in too large compactness of monomeric α-Syn, structural properties of monomeric Aβ42 and the heterodimer were preserved in explicit-solvent simulations. We anticipate that our study sheds light on the interaction between α-Syn and Aβ42 proteins, thus providing the atom-level model required to assess the initial stage of aggregation mechanisms related to Alzheimer's and Parkinson's diseases.

Antidurotaxis Droplet Motion onto Gradient Brush Substrates

Durotaxis motion is a spectacular phenomenon manifesting itself by the autonomous motion of a nano-object between parts of a substrate with different stiffness. This motion usually takes place along a stiffness gradient from softer to stiffer parts of the substrate. Here, we propose a new design of a polymer brush substrate that demonstrates antidurotaxis droplet motion, that is, droplet motion from stiffer to softer parts of the substrate. By carrying out extensive molecular dynamics simulation of a coarse-grained model, we find that antidurotaxis is solely controlled by the gradient in the grafting density of the brush and is favorable for fluids with a strong attraction to the substrate (low surface energy). The driving force of the antidurotaxial motion is the minimization of the droplet-substrate interfacial energy, which is attributed to the penetration of the droplet into the brush. Thus, we anticipate that the proposed substrate design offers a new understanding and possibilities in the area of autonomous motion of droplets for applications in microfluidics, energy conservation, and biology.

The α-Synuclein Monomer May Have Different Misfolding Mechanisms in the Induction of α-Synuclein Fibrils with Different Polymorphs

The aggregation of alpha-synuclein (α-Syn) is closely related to the occurrence of some neurodegenerative diseases such as Parkinson's disease. The misfolding of α-Syn monomer plays a key role in the formation of aggregates and extension of fibril. However, the misfolding mechanism of α-Syn remains elusive. Here, three different α-Syn fibrils (isolated from a diseased human brain, generated by in vitro cofactor-tau induction, and obtained by in vitro cofactor-free induction) were selected for the study. The misfolding mechanisms of α-Syn were uncovered by studying the dissociation of the boundary chains based on the conventional molecular dynamics (MD) and Steered MD simulations. The results showed that the dissociation paths of the boundary chains in the three systems were different. According to the reverse process of dissociation, we concluded that in the human brain system, the binding of the monomer and template starts from the C-terminal and gradually misfolds toward the N-terminal. In the cofactor-tau system, the monomer binding starts from residues 58-66 (contain β3), followed by the C-terminal coil (residues 67-79). Then, the N-terminal coil (residues 36-41) and residues 50-57 (contain β2) bind to the template, followed by residues 42-49 (contain β1). In the cofactor-free system, two misfolding paths were found. One is that the monomer binds to the N/C-terminal (β1/β6) and then binds to the remaining residues. The other one is that the monomer binds sequentially from the C- to N-terminal, similar to the human brain system. Furthermore, in the human brain and cofactor-tau systems, electrostatic interactions (especially from residues 58-66) are the main driving force during the misfolding process, whereas in the cofactor-free system, the contributions of electrostatic and van der Waals interactions are comparable. These results may provide a deeper understanding for the misfolding and aggregation mechanism of α-Syn.

May the force be with you: The role of hyper-mechanostability of the bone sialoprotein binding protein during early stages of Staphylococci infections

The bone sialoprotein-binding protein (Bbp) is a mechanoactive MSCRAMM protein expressed on the surface of Staphylococcus aureus that mediates adherence of the bacterium to fibrinogen-α (Fgα), a component of the bone and dentine extracellular matrix of the host cell. Mechanoactive proteins like Bbp have key roles in several physiological and pathological processes. Particularly, the Bbp: Fgα interaction is important in the formation of biofilms, an important virulence factor of pathogenic bacteria. Here, we investigated the mechanostability of the Bbp: Fgα complex using in silico single-molecule force spectroscopy (SMFS), in an approach that combines results from all-atom and coarse-grained steered molecular dynamics (SMD) simulations. Our results show that Bbp is the most mechanostable MSCRAMM investigated thus far, reaching rupture forces beyond the 2 nN range in typical experimental SMFS pulling rates. Our results show that high force-loads, which are common during initial stages of bacterial infection, stabilize the interconnection between the protein's amino acids, making the protein more "rigid". Our data offer new insights that are crucial on the development of novel anti-adhesion strategies.

Resolution exchange with tunneling for enhanced sampling of protein landscapes

Simulations of protein folding and protein association happen on timescales that are orders of magnitude larger than what can typically be covered in all-atom molecular dynamics simulations. Use of low-resolution models alleviates this problem but may reduce the accuracy of the simulations. We introduce a replica-exchange-based multiscale sampling technique that combines the faster sampling in coarse-grained simulations with the potentially higher accuracy of all-atom simulations. After testing the efficiency of our Resolution Exchange with Tunneling (ResET) in simulations of the Trp-cage protein, an often used model to evaluate sampling techniques in protein simulations, we use our approach to compare the landscape of wild-type and A2T mutant Aβ_{1-42} peptides. Our results suggest a mechanism by that the mutation of a small hydrophobic alanine (A) into a bulky polar threonine (T) may interfere with the self-assembly of Aβ fibrils.

Mapping Mechanostable Pulling Geometries of a Therapeutic Anticalin/CTLA-4 Protein Complex

We used single-molecule AFM force spectroscopy (AFM-SMFS) in combination with click chemistry to mechanically dissociate anticalin, a non-antibody protein binding scaffold, from its target (CTLA-4), by pulling from eight different anchor residues. We found that pulling on the anticalin from residue 60 or 87 resulted in significantly higher rupture forces and a decrease in koff by 2-3 orders of magnitude over a force range of 50-200 pN. Five of the six internal anchor points gave rise to complexes significantly more stable than N- or C-terminal anchor points, rupturing at up to 250 pN at loading rates of 0.1-10 nN s-1. Anisotropic network modeling and molecular dynamics simulations helped to explain the geometric dependency of mechanostability. These results demonstrate that optimization of attachment residue position on therapeutic binding scaffolds can provide large improvements in binding strength, allowing for mechanical affinity maturation under shear stress without mutation of binding interface residues.

Assessing the Stability of Biological Fibrils by Molecular-Scale Simulations

The nanomechanical characterization of several biological fibrils that are the result of protein aggregation via molecular dynamics simulation is nowadays feasible, and together with atomic force microscopy experiments has widened our understanding of the forces in the regime of pN-nN and system sizes of about hundreds of nanometers. Several methodologies have been developed to achieve this target, and they range from the atomistic representation via molecular force fields to coarse-grained strategies that provide comparable results with experiments in a systematic way. In this chapter, we discuss several methodologies for the calculation of mechanical parameters, such as the elastic constants of relevant biological systems. They are presented together with details about parameterization and current limitations. Then, we discuss some of the applications of such methodologies for the description of bacterial filament and β-amyloid systems. Finally, the latest lines of development are discussed.

Explore the role of CR1 genetic variants in late-onset Alzheimer's disease susceptibility

BACKGROUND

Complement component (3b/4b) receptor 1 (CR1) is an interesting candidate gene which has a close connection with Alzheimer's disease, and its polymorphisms have been reported to link to the late-onset Alzheimer's disease (LOAD) susceptibility. However, the findings of these related studies are inconsistent. Objective To explore the effect of CR1 genetic variants in LOAD susceptibility. MethodsWe searched relevant studies for the period up to 1 November 2020. And odds ratios (ORs) and their 95% confidence intervals (CIs) were utilized to assess the strength of the association. In addition, we carried out a case-control association study to assess their genetic association.

RESULTS

Finally, a total of 30 articles with 30108 LOAD cases and 37895 controls were included. Significant allele frequency between LOAD patients and controls was observed in rs3818361 and rs6656401 (rs3818361, T vs. C: OR,1.18; 95% CI, 1.13-1.23; rs6656401, A vs. G: OR, 1.23; 95% CI, 1.10-1.36). Moreover, these results remain significant in subgroup of rs3818361 in Asia or America (OR,1.26; 95% CI,1.06-1.45; OR, 1.18; 95% CI, 1.13-1.24, respectively) and rs6656401 in Europe (OR = 1.26; 95% CI, 1.09-1.42). In addition, the two single nucleotide polymorphisms were proved to significantly increase LOAD risk in the overall population under the dominant model (OR = 1.12; 95% CI, 1.02-1.21; OR = 1.18, 95% CI, 1.15-1.22, respectively). Our case-control study showed that the distribution of rs6656401 genotype was significant (P = 0.000; OR, 6.889; 95% CI, 2.709-17.520), suggesting the A allele of rs6656401 is the risk allele.

CONCLUSION

These available data indicate that rs6656401 in CR1 is significant to increase LOAD risk.

Systematic Search for a Predictor for the Clinical Observables of Alzheimer's Disease

One of the prevailing life-threatening incurable neurodegenerative diseases that are presently endangering human society as a whole, and hence, baffling the entire spectrum of the scientific and pharmaceutical world, is Alzheimer's disease (AD). AD is a manifestation of self-assembly of both wild-type (sporadic) and mutated (familial) forms of the amyloid-β peptide, a proteolytic product of the amyloid precursor protein, where the self-assembly results in the genesis of pathogenic fibrillar aggregates. Currently prevailing diagnostic and hence therapeutic challenges originate from the unavailability of a specific predictor for clinical observables. The continuous emergence of novel pathogenic mutants with unpredictable phenotypes adds immensely to the nonspecific nature of the problem. The current research reports a simple physical parameter, the binding affinity of a protofilament to its protofibril, which predicts the clinical observables of familial AD with astounding accuracy and more importantly, without any adjustable parameters.

An Overview of Several Inhibitors for Alzheimer's Disease: Characterization and Failure

Amyloid beta (Aβ) oligomers are the most neurotoxic aggregates causing neuronal death and cognitive damage. A detailed elucidation of the aggregation pathways from oligomers to fibril formation is crucial to develop therapeutic strategies for Alzheimer's disease (AD). Although experimental techniques rely on the measure of time- and space-average properties, they face severe difficulties in the investigation of Aβ peptide aggregation due to their intrinsically disorder character. Computer simulation is a tool that allows tracing the molecular motion of molecules; hence it complements Aβ experiments, as it allows to explore the binding mechanism between metal ions and Aβ oligomers close to the cellular membrane at the atomic resolution. In this context, integrated studies of experiments and computer simulations can assist in mapping the complete pathways of aggregation and toxicity of Aβ peptides. Aβ oligomers are disordered proteins, and due to a rapid exploration of their intrinsic conformational space in real-time, they are challenging therapeutic targets. Therefore, no good drug candidate could have been identified for clinical use. Our previous investigations identified two small molecules, M30 (2-Octahydroisoquinolin-2(1H)-ylethanamine) and Gabapentin, capable of Aβ binding and inhibiting molecular aggregation, synaptotoxicity, intracellular calcium signaling, cellular toxicity and memory losses induced by Aβ. Thus, we recommend these molecules as novel candidates to assist anti-AD drug discovery in the near future. This review discusses the most recent research investigations about the Aβ dynamics in water, close contact with cell membranes, and several therapeutic strategies to remove plaque formation.

Nanomechanical Stability of Aβ Tetramers and Fibril-like Structures: Molecular Dynamics Simulations

Alzheimer's disease (AD) is a neurodegenerative disorder and one of the main causes of dementia. The disease is associated with amyloid beta (Aβ) peptide aggregation forming initial clusters and then fibril structure and plaques. Other neurodegenerative diseases such as type 2 diabetes, amyotrophic lateral sclerosis, and Parkinson's disease follow a similar mechanism. Therefore, inhibition of Aβ aggregation is considered an effective way to prevent AD. Recent experiments have provided evidence that oligomers are more toxic agents than mature fibrils, prompting researchers to investigate various factors that may influence their properties. One of these factors is nanomechanical stability, which plays an important role in the self-assembly of Aβ and possibly other proteins. This stability is also likely to be related to cell toxicity. In this work, we compare the mechanical stability of Aβ-tetramers and fibrillar structures using a structure-based coarse-grained (CG) approach and all-atom molecular dynamics simulation. Our results support the evidence for an increase in mechanical stability during the Aβ fibrillization process, which is consistent with in vitro AFM characterization of Aβ42 oligomers. Namely, using a CG model, we showed that the Young modulus of tetramers is lower than that of fibrils and, as follows from the experiment, is about 1 GPa. Hydrogen bonds are the dominant contribution to the detachment of one chain from the Aβ fibril fragment. They tend to be more organized along the pulling direction, whereas in the Aβ tetramers no preference is observed.

Rapid Conversion of Amyloid-Beta 1-40 Oligomers to Mature Fibrils through a Self-Catalytic Bimolecular Process

The formation of fibrillar aggregates of the amyloid beta peptide (Aβ) in the brain is one of the hallmarks of Alzheimer's disease (AD). A clear understanding of the different aggregation steps leading to fibrils formation is a keystone in therapeutics discovery. In a recent study, we showed that Aβ40 and Aβ42 form dynamic micellar aggregates above certain critical concentrations, which mediate a fast formation of more stable oligomers, which in the case of Aβ40 are able to evolve towards amyloid fibrils. Here, using different biophysical techniques we investigated the role of different fractions of the Aβ aggregation mixture in the nucleation and fibrillation steps. We show that both processes occur through bimolecular interplay between low molecular weight species (monomer and/or dimer) and larger oligomers. Moreover, we report here a novel self-catalytic mechanism of fibrillation of Aβ40, in which early oligomers generate and deliver low molecular weight amyloid nuclei, which then catalyze the rapid conversion of the oligomers to mature amyloid fibrils. This fibrillation catalytic activity is not present in freshly disaggregated low-molecular weight Aβ40 and is, therefore, a property acquired during the aggregation process. In contrast to Aβ40, we did not observe the same self-catalytic fibrillation in Aβ42 spheroidal oligomers, which could neither be induced to fibrillate by the Aβ40 nuclei. Our results reveal clearly that amyloid fibrillation is a multi-component process, in which dynamic collisions between different interacting species favor the kinetics of amyloid nucleation and growth.

Structural Studies Providing Insights into Production and Conformational Behavior of Amyloid-β Peptide Associated with Alzheimer's Disease Development

Alzheimer's disease is the most common type of neurodegenerative disease in the world. Genetic evidence strongly suggests that aberrant generation, aggregation, and/or clearance of neurotoxic amyloid-β peptides (Aβ) triggers the disease. Aβ accumulates at the points of contact of neurons in ordered cords and fibrils, forming the so-called senile plaques. Aβ isoforms of different lengths are found in healthy human brains regardless of age and appear to play a role in signaling pathways in the brain and to have neuroprotective properties at low concentrations. In recent years, different substances have been developed targeting Aβ production, aggregation, interaction with other molecules, and clearance, including peptide-based drugs. Aβ is a product of sequential cleavage of the membrane glycoprotein APP (amyloid precursor protein) by β- and γ-secretases. A number of familial mutations causing an early onset of the disease have been identified in the APP, especially in its transmembrane domain. The mutations are reported to influence the production, oligomerization, and conformational behavior of Aβ peptides. This review highlights the results of structural studies of the main proteins involved in Alzheimer's disease pathogenesis and the molecular mechanisms by which perspective therapeutic substances can affect Aβ production and nucleation.

Temperature-Dependent Structural Variability of Prion Protein Amyloid Fibrils

Prion protein aggregation into amyloid fibrils is associated with the onset and progression of prion diseases—a group of neurodegenerative amyloidoses. The process of such aggregate formation is still not fully understood, especially regarding their polymorphism, an event where the same type of protein forms multiple, conformationally and morphologically distinct structures. Considering that such structural variations can greatly complicate the search for potential antiamyloid compounds, either by having specific propagation properties or stability, it is important to better understand this aggregation event. We have recently reported the ability of prion protein fibrils to obtain at least two distinct conformations under identical conditions, which raised the question if this occurrence is tied to only certain environmental conditions. In this work, we examined a large sample size of prion protein aggregation reactions under a range of temperatures and analyzed the resulting fibril dye-binding, secondary structure and morphological properties. We show that all temperature conditions lead to the formation of more than one fibril type and that this variability may depend on the state of the initial prion protein molecules.

Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field

We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.

Current Understanding of the Structure, Stability and Dynamic Properties of Amyloid Fibrils

Amyloid fibrils are supramolecular protein assemblies represented by a cross-β structure and fibrous morphology, whose structural architecture has been previously investigated. While amyloid fibrils are basically a main-chain-dominated structure consisting of a backbone of hydrogen bonds, side-chain interactions also play an important role in determining their detailed structures and physicochemical properties. In amyloid fibrils comprising short peptide segments, a steric zipper where a pair of β-sheets with side chains interdigitate tightly is found as a fundamental motif. In amyloid fibrils comprising longer polypeptides, each polypeptide chain folds into a planar structure composed of several β-strands linked by turns or loops, and the steric zippers are formed locally to stabilize the structure. Multiple segments capable of forming steric zippers are contained within a single protein molecule in many cases, and polymorphism appears as a result of the diverse regions and counterparts of the steric zippers. Furthermore, the β-solenoid structure, where the polypeptide chain folds in a solenoid shape with side chains packed inside, is recognized as another important amyloid motif. While side-chain interactions are primarily achieved by non-polar residues in disease-related amyloid fibrils, the participation of hydrophilic and charged residues is prominent in functional amyloids, which often leads to spatiotemporally controlled fibrillation, high reversibility, and the formation of labile amyloids with kinked backbone topology. Achieving precise control of the side-chain interactions within amyloid structures will open up a new horizon for designing useful amyloid-based nanomaterials.

CT-guided radiofrequency treatment of trigeminal neuralgia at different temperatures through foramen rotundus

OBJECTIVE

To compare the efficacy and safety of CT-guided radiofrequency therapy in the treatment of trigeminal neuralgia (TN) through foramen rotundus at different temperatures.

METHODS

A total of 60 patients with TN in our hospital were selected and randomly divided into the control group and the observation group, with 30 cases in each group. CT or MRI was routinely performed before the operation to exclude intracranial space-occupying lesions. The round hole of both groups of patients was targeted under CT guidance. The puncture needle position and depth were determined according to CT images and radiofrequency electrical stimulation. The control group gradually increased the radiofrequency temperature to 75°C, and the observation group gradually increased to 90°C for damage treatment. Both groups of patients were treated with radiofrequency for 360 s. The clinical efficacy was evaluated by BNI classification criteria. At the same time, the visual analog scale (VAS) score, the total score of Pittsburgh sleep quality index (PSQI) and the incidence of complications were compared between the two groups. The Kaplan-Meier method in survival analysis was used to calculate the pain recurrence rate at 1 and 2 years after the operation.

RESULTS

One week after the operation, there was no significant difference in the total effective rate between the two groups (χ²=0.089, P=0.766). After 2 months, the total effective rate of the observation group was higher than that of the control group (χ²=4.043, P=0.044). One day and one week after the operation, the VAS score of the observation group was higher than that of the control group (t=7.365, 6.269; P=0.007, 0.012), and the total score of PSQI was significantly higher than that of the control group (t=8.026, 5.447; P=0.002, 0.015). There was no significant difference in VAS score and PSQI total score between the two groups 1 month and 2 months after operation (P>0.05). There was no significant difference in the incidence of postoperative complications between the two groups (P>0.05). Kaplan-Meier survival analysis was used to calculate the recurrence rate of pain at 1 and 2 years after operation. The recurrence rate of the observation group was significantly lower than that of the control group (χ²=4.219, 4.021; P=0.039, 0.044).

CONCLUSION

CT-guided radiofrequency at 90°C through foramen rotundus is effective in the treatment of trigeminal neuralgia without increasing the incidence of complications.

Molecular Dynamics Simulations of a Catalytic Multivalent Peptide-Nanoparticle Complex

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide-nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure-function relationship, for technological and nanomedical applications.

Enhancing the binding of the β-sheet breaker peptide LPFFD to the amyloid-β fibrils by aromatic modifications: A molecular dynamics simulation study

Alzheimer's is a fatal neurodegenerative disease for which there is no cure at present. The disease is characterized by the presence of plaques in the brains of a patient, which are composed mainly of aggregates of the amyloid-β peptide in the form of β-sheet fibrils. Here, we investigated the possibility of exploiting the superior binding ability of aromatic amino acids to a particular model of the amyloid-β fibrils. which is a difficult target for drug design. The β-sheet breaker peptide LPFFD was modified with aromatic amino acids and its binding to these fibrils was studied. We found that the orientation and the electrostatic complementarity of the modified peptide with respect to the fibrils played a crucial role in determining whether its binding was improved by the aromatic amino acids. The modified LPFFD peptides were able to bind to those fibril residues. which are important in the aggregation of amyloid-β peptides and thus can potentially inhibit the further aggregation of the amyloid-beta peptides by blocking their interactions. We found that the tryptophan modified LPFFD peptides had the best binding affinities. In most cases, the aromatic amino acids in the N-terminus of the modified peptides made more contacts with the fibrils than those in the C-terminus. We also found that increasing the aromatic content did not significantly improve the binding of the LPFFD peptide to the fibrils. Our study can serve as a basis for the design of novel peptide-based drugs for Alzheimer's disease in which aromatic interactions play an important role.

Optimizing Gō-MARTINI Coarse-Grained Model for F-BAR Protein on Lipid Membrane

Coarse-grained (CG) molecular dynamics (MD) simulations allow us to access much larger length and time scales than atomistic MD simulations, providing an attractive alternative to the conventional simulations. Based on the well-known MARTINI CG force field, the recently developed Gō-MARTINI model for proteins describes large-amplitude structural dynamics, which has not been possible with the commonly used elastic network model. Using the Gō-MARTINI model, we conduct MD simulations of the F-BAR Pacsin1 protein on lipid membrane. We observe that structural changes of the non-globular protein are largely dependent on the definition of the native contacts in the Gō model. To address this issue, we introduced a simple cutoff scheme and tuned the cutoff distance of the native contacts and the interaction strength of the Lennard-Jones potentials in the Gō-MARTINI model. With the optimized Gō-MARTINI model, we show that it reproduces structural fluctuations of the Pacsin1 dimer from atomistic simulations. We also show that two Pacsin1 dimers properly assemble through lateral interaction on the lipid membrane. Our work presents a first step towards describing membrane remodeling processes in the Gō-MARTINI CG framework by simulating a crucial step of protein assembly on the membrane.

Characterization of Structural and Energetic Differences between Conformations of the SARS-CoV-2 Spike Protein

The novel coronavirus disease 2019 (COVID-19) pandemic has disrupted modern societies and their economies. The resurgence in COVID-19 cases as part of the second wave is observed across Europe and the Americas. The scientific response has enabled a complete structural characterization of the Severe Acute Respiratory Syndrome-novel Coronavirus 2 (SARS-CoV-2). Among the most relevant proteins required by the novel coronavirus to facilitate the cell entry mechanism is the spike protein. This protein possesses a receptor-binding domain (RBD) that binds the cellular angiotensin-converting enzyme 2 (ACE2) and then triggers the fusion of viral and host cell membranes. In this regard, a comprehensive characterization of the structural stability of the spike protein is a crucial step to find new therapeutics to interrupt the process of recognition. On the other hand, it has been suggested that the participation of more than one RBD is a possible mechanism to enhance cell entry. Here, we discuss the protein structural stability based on the computational determination of the dynamic contact map and the energetic difference of the spike protein conformations via the mapping of the hydration free energy by the Poisson-Boltzmann method. We expect our result to foster the discussion of the number of RBD involved during recognition and the repurposing of new drugs to disable the recognition by discovering new hotspots for drug targets apart from the flexible loop in the RBD that binds the ACE2.

Quantitative determination of mechanical stability in the novel coronavirus spike protein

We report on the novel observation about the gain in nanomechanical stability of the SARS-CoV-2 (CoV2) spike (S) protein in comparison with SARS-CoV from 2002 (CoV1). Our findings have several biological implications in the subfamily of coronaviruses, as they suggest that the receptor binding domain (RBD) (∼200 amino acids) plays a fundamental role as a damping element of the massive viral particle's motion prior to cell-recognition, while also facilitating viral attachment, fusion and entry. The mechanical stability via pulling of the RBD is 250 pN and 200 pN for CoV2 and CoV1 respectively, and the additional stability observed for CoV2 (∼50 pN) might play a role in the increasing spread of COVID-19.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

Differentiating between Inactive and Active States of Rhodopsin by Atomic Force Microscopy in Native Membranes

Membrane proteins, including G protein-coupled receptors (GPCRs), present a challenge in studying their structural properties under physiological conditions. Moreover, to better understand the activity of proteins requires examination of single molecule behaviors rather than ensemble averaged behaviors. Force-distance curve-based AFM (FD-AFM) was utilized to directly probe and localize the conformational states of a GPCR within the membrane at nanoscale resolution based on the mechanical properties of the receptor. FD-AFM was applied to rhodopsin, the light receptor and a prototypical GPCR, embedded in native rod outer segment disc membranes from photoreceptor cells of the retina in mice. Both FD-AFM and computational studies on coarse-grained models of rhodopsin revealed that the active state of the receptor has a higher Young's modulus compared to the inactive state of the receptor. Thus, the inactive and active states of rhodopsin could be differentiated based on the stiffness of the receptor. Differentiating the states based on the Young's modulus allowed for the mapping of the different states within the membrane. Quantifying the active states present in the membrane containing the constitutively active G90D rhodopsin mutant or apoprotein opsin revealed that most receptors adopt an active state. Traditionally, constitutive activity of GPCRs has been described in terms of two-state models where the receptor can achieve only a single active state. FD-AFM data are inconsistent with a two-state model but instead require models that incorporate multiple active states.

The Role of Structural Polymorphism in Driving the Mechanical Performance of the Alzheimer's Beta Amyloid Fibrils

Alzheimer's Disease (AD) is related with the abnormal aggregation of amyloid β-peptides Aβ₁₋₄₀ and Aβ₁₋₄₂, the latter having a polymorphic character which gives rise to U- or S-shaped fibrils. Elucidating the role played by the nanoscale-material architecture on the amyloid fibril stability is a crucial breakthrough to better understand the pathological nature of amyloid structures and to support the rational design of bio-inspired materials. The computational study here presented highlights the superior mechanical behavior of the S-architecture, characterized by a Young's modulus markedly higher than the U-shaped architecture. The S-architecture showed a higher mechanical resistance to the enforced deformation along the fibril axis, consequence of a better interchain hydrogen bonds' distribution. In conclusion, this study, focusing the attention on the pivotal multiscale relationship between molecular phenomena and material properties, suggests the S-shaped Aβ₁₋₄₂ species as a target of election in computational screen/design/optimization of effective aggregation modulators.