Benchmarking of copper(II) LFMM parameters for studying amyloid-β peptides

The Pursuit of the "Inside" of the Amyloid Hypothesis-Is C99 a Promising Therapeutic Target for Alzheimer's Disease?

Aducanumab, co-developed by Eisai (Japan) and Biogen (U.S.), has received Food and Drug Administration approval for treating Alzheimer's disease (AD). In addition, its successor antibody, lecanemab, has been approved. These antibodies target the aggregated form of the small peptide, amyloid-β (Aβ), which accumulates in the patient brain. The "amyloid hypothesis" based therapy that places the aggregation and toxicity of Aβ at the center of the etiology is about to be realized. However, the effects of immunotherapy are still limited, suggesting the need to reconsider this hypothesis. Aβ is produced from a type-I transmembrane protein, Aβ precursor protein (APP). One of the APP metabolites, the 99-amino acids C-terminal fragment (C99, also called βCTF), is a direct precursor of Aβ and accumulates in the AD patient's brain to demonstrate toxicity independent of Aβ. Conventional drug discovery strategies have focused on Aβ toxicity on the "outside" of the neuron, but C99 accumulation might explain the toxicity on the "inside" of the neuron, which was overlooked in the hypothesis. Furthermore, the common region of C99 and Aβ is a promising target for multifunctional AD drugs. This review aimed to outline the nature, metabolism, and impact of C99 on AD pathogenesis and discuss whether it could be a therapeutic target complementing the amyloid hypothesis.

Forcefield evaluation and accelerated molecular dynamics simulation of Zn(II) binding to N-terminus of amyloid-β

We report conventional and accelerated molecular dynamics simulation of Zn(II) bound to the N-terminus of amyloid-β. By comparison against NMR data for the experimentally determined binding mode, we find that certain combinations of forcefield and solvent model perform acceptably in describing the size, shape and secondary structure, and that there is no appreciable difference between implicit and explicit solvent models. We therefore used the combination of ff14SB forcefield and GBSA solvent model to compare the result of different binding modes of Zn(II) to the same peptide, using accelerated MD to enhance sampling and comparing the free peptide simulated in the same way. We show that Zn(II) imparts significant rigidity to the peptide, disrupts the secondary structure and pattern of salt bridges seen in the free peptide, and induces closer contact between residues. Free energy surfaces in 1 or 2 dimensions further highlight the effect of metal coordination on peptide's spatial extent. We also provide evidence that accelerated MD provides improved sampling over conventional MD by visiting as many or more configurations in much shorter simulation times.

The hazardous effects of the environmental toxic gases on amyloid beta-peptide aggregation: A theoretical perspective

Alzheimer's disease (AD) is one of the leading causes of dementia in elderly people. It has been well documented that the exposure to environmental toxins such as CO, CO₂, SO₂ and NO₂ that are present in the air is considered as a hallmark for the progression of Alzheimer's disease. However, their actual mechanism by which environmental toxin triggers the aggregation of Aβ42 peptide at the molecular and atomic levels remain unknown. In this study, molecular dynamics simulation was carried out to study the aggregation mechanism of the Aβ₄₂ peptide due to its interaction of toxic gas (CO, CO₂, SO₂ and NO₂). During the 400 ns simulation, all the Aβ₄₂ interacted toxic gas (CO, CO₂, SO_2, and NO₂) complexes have smaller Root Mean Square Deviation values when compared to the Aβ₄₂ peptide, which shows that the interaction of toxic gases (CO, CO₂, SO_2, and NO₂) would increase the Aβ₄₂ peptide structural stability. The radius of gyration analysis also supports that Aβ₄₂ interacted CO₂ and SO₂ complexes have the minimum value in the range of 0.95 nm and 1.5 nm. It is accounted that the Aβ₄₂ interacted CO₂ and SO₂ complexes have a greater compact structure in comparison to Aβ₄₂ interacted CO and NO₂ complexes. Furthermore, all the Aβ₄₂ interacted toxic gas (CO, CO₂, SO_2, and NO₂) complexes exhibited an enhanced secondary structural probability for coil and turn regions with a reduced α-helix probability, which indicates that the interaction of toxic gases may enhance the toxicity and aggregation of Aβ₄₂.

Theoretical study of copper binding to GHK peptide

We report ligand field molecular mechanics, density functional theory and semi-empirical studies on the binding of Cu(II) to GlyHisLys (GHK) peptide. Following exhaustive conformational searching using molecular mechanics, we show that relative energy and geometry of conformations are in good agreement between GFN2-xTB semi-empirical and B3LYP-D DFT levels. Conventional molecular dynamics simulation of Cu-GHK shows the stability of the copper-peptide binding over 100 ps trajectory. Four equatorial bonds in 3N1O coordination remain stable throughout simulation, while a fifth in apical position from C-terminal carboxylate is more fluxional. We also show that the automated conformer and rotamer search algorithm CREST is able to correctly predict the metal binding position from a starting point consisting of separated peptide, copper and water.

In silico studies of the inhibition mechanism of dengue with papain

Dengue virus is becoming a major global disease; the envelope protein is the major target for vaccine development against Dengue. Nowadays, the attention has focused on developing inhibitors based on Papain is a promising target for treating Dengue. In the present work, the theoretical studies of E-protein(Cys74-Glu79;Lys110)…Papain(Cys25, Asn175 and His159) complexes are analysed by Density Functional Theory (M06-2X/cc-pVDZ) method. Among the E-protein(Cys74-Glu79;Lys110)…Papain(Cys25, Asn175 and Hys159) complexes, E-protein(Glu76)…Papain(Cys25) complex has the highest interaction value of -352.22 kcal/mol. Moreover, the natural bond orbital analysis also supports the above results. The 100 ns Molecular Dynamics simulation reveals that, E-protein(Ala54-Ile129)…Papain(Cys25) complex had the lowest root mean square deviation value of 1 Å compared to the E-protein(Ala54-Ile129)… Papain(Asn175 & His159) complexes. The salt bridge formation between the Asp103 and Lys110 residues are the important stabilizing factor in E-protein(Ala54-Ile129)…Papain(Cys25) complex. This result can extend our knowledge of the functional behaviour of Papain and provides structural insight to target Envelope protein as forthcoming drug targets in Dengue.

Molecular dynamics simulations of copper binding to amyloid-β Glu22 mutants

We report microsecond timescale ligand field molecular dynamics simulations of the copper complexes of three known mutants of the amyloid-β peptide, E22G, E22Q and E22K, alongside the naturally occurring sequence. We find that all three mutants lead to formation of less compact structures than the wild-type: E22Q is the most similar to the native peptide, while E22G and especially E22K are markedly different in size, shape and stability. Turn and coil structures dominate all structures studied but subtle differences in helical and β-sheet distribution are noted, especially in the C-terminal region. The origin of these changes is traced to disruption of key salt bridges: in particular, the Asp23-Lys28 bridge that is prevalent in the wild-type is absent in E22G and E22K, while Lys22 in the latter mutant forms a strong association with Asp23. We surmise that the drastically different pattern of salt bridges in the mutants lead to adoption of a different structural ensemble of the peptide backbone, and speculate that this might affect the ability of the mutant peptides to aggregate in the same manner as known for the wild-type.

Metal Binding to Amyloid-β1-42: A Ligand Field Molecular Dynamics Study

Ligand field molecular mechanics simulation has been used to model the interactions of copper(II) and platinum(II) with the amyloid-β_1-42 peptide monomer. Molecular dynamics over several microseconds for both metalated systems are compared to analogous results for the free peptide. Significant differences in structural parameters are observed, both between Cu and Pt bound systems as well as between free and metal-bound peptide. Both metals stabilize the formation of helices in the peptide as well as reducing the content of β secondary structural elements compared to the unbound monomer. This is in agreement with experimental reports of metals reducing β-sheet structures, leading to formation of amorphous aggregates over amyloid fibrils. The shape and size of the peptide structures also undergo noteworthy change, with the free peptide exhibiting globular-like structure, platinum(II) system adopting extended structures, and copper(II) system resulting in a mixture of conformations similar to both of these. Salt bridge networks exhibit major differences: the Asp23-Lys28 salt bridge, known to be important in fibril formation, has a differing distance profile within all three systems studied. Salt bridges in the metal binding region of the peptide are strongly altered; in particular, the Arg5-Asp7 salt bridge, which has an occurrence of 71% in the free peptide, is reduced to zero in the presence of both metals.