A CRITERION FOR ORBITAL HYBRIDIZATION AND CHARGE DISTRIBUTION IN CHEMICAL BONDS

Backbone Torsion Angle Determination Using Proton Detected Magic-Angle Spinning Nuclear Magnetic Resonance

Protein torsion angles define the backbone secondary structure of proteins. Magic-angle spinning (MAS) NMR methods using carbon detection have been developed to measure torsion angles by determining the relative orientation between two anisotropic interactions─dipolar coupling or chemical shift anisotropy. Here we report a new proton-detection based method to determine the backbone torsion angle by recoupling NH and CH dipolar couplings within the HCANH pulse sequence, for protonated or partly deuterated samples. We demonstrate the efficiency and precision of the method with microcrystalline chicken α spectrin SH3 protein and the influenza A matrix 2 (M2) membrane protein, using 55 or 90 kHz MAS. For M2, pseudo-4D data detect a turn between transmembrane and amphipathic helices.

Serotonin and Melatonin Show Different Modes of Action on Aβ42 Protofibril Destabilization

Alzheimer's disease (AD) is associated with the aberrant self-assembly of amyloid-β (Aβ) protein into fibrillar deposits. The disaggregation of Aβ fibril is believed as one of the major therapeutic strategies for treating AD. Previous experimental studies reported that serotonin (Ser), one of the indoleamine neurotransmitters, and its derivative melatonin (Mel) are able to disassemble preformed Aβ fibrils. However, the fibril-disruption mechanisms are unclear. As the first step to understand the underlying mechanism, we investigated the interactions of Ser and Mel molecules with the LS-shaped Aβ₄₂ protofibril by performing a total of nine individual 500 ns all-atom molecular dynamics (MD) simulations. The simulations demonstrate that both Ser and Mel molecules disrupt the local β-sheet structure, destroy the salt bridges between K28 side chain and A42 COO^-, and consequently destabilize the global structure of Aβ₄₂ protofibril. The Mel molecule exhibits a greater binding capacity than the Ser molecule. Intriguingly, we find that Ser and Mel molecules destabilize Aβ₄₂ protofibril through different modes of action. Ser preferentially binds with the aromatic residues in the N-terminal region through π-π stacking interactions, while Mel binds not only with the N-terminal aromatic residues but also with the C-terminal hydrophobic residues via π-π and hydrophobic interactions. This work reveals the disruptive mechanisms of Aβ₄₂ protofibril by Ser and Mel molecules and provides useful information for designing drug candidates against AD.

Impact of Mutations on the Conformational Transition from α-Helix to β-Sheet Structures in Arctic-Type Aβ40: Insights from Molecular Dynamics Simulations

The amyloid-β (Aβ) protein aggregation into toxic oligomers and fibrils has been recognized as a key player in the pathogenesis of Alzheimer's disease. Recent experiments reported that a double alanine mutation (L17A/F19A) in the central hydrophobic core (CHC) region of [G22]Aβ₄₀ (familial Arctic mutation) diminished the self-assembly propensity of [G22]Aβ₄₀. However, the molecular mechanism behind the decreased aggregation tendency of [A17/A19/G22]Aβ₄₀ is not well understood. Herein, we carried out molecular dynamics simulations to elucidate the structure and dynamics of [G22]Aβ₄₀ and [A17/A19/G22]Aβ₄₀. The results for the secondary structure analysis reveal a significantly increased amount of the helical content in the CHC and C-terminal region of [A17/A19/G22]Aβ₄₀ as compared to [G22]Aβ₄₀. The bending free-energy analysis of D23-K28 salt bridge suggests that the double alanine mutation in the CHC region of [G22]Aβ₄₀ has the potential to reduce the fibril formation rate by 0.57 times of [G22]Aβ₄₀. Unlike [G22]Aβ₄₀, [A17/A19/G22]Aβ₄₀ largely sampled helical conformation, as determined by the minimum energy conformations extracted from the free-energy landscape. The present study provided atomic level details into the experimentally observed diminished aggregation tendency of [A17/A19/G22]Aβ₄₀ as compared to [G22]Aβ₄₀.

How Does the Mono-Triazole Derivative Modulate Aβ42 Aggregation and Disrupt a Protofibril Structure: Insights from Molecular Dynamics Simulations

Clinical studies have identified that abnormal self-assembly of amyloid-β (Aβ) peptide into toxic fibrillar aggregates is associated with the pathology of Alzheimer's disease (AD). The most acceptable therapeutic approach to stop the progression of AD is to inhibit the formation of β-sheet-rich structures. Recently, we designed and evaluated a series of novel mono-triazole derivatives 4(a-x), where compound 4v was identified as the most potent inhibitor of Aβ₄₂ aggregation and disaggregates preformed Aβ₄₂ fibrils significantly. Moreover, 4v strongly averts the Cu²⁺-induced Aβ₄₂ aggregation and disaggregates the preformed Cu²⁺-induced Aβ₄₂ fibrils, halts the generation of reactive oxygen species, and shows neuroprotective effects in SH-SY5Y cells. However, the underlying molecular mechanism of inhibition of Aβ₄₂ aggregation by 4v and disaggregation of preformed Aβ₄₂ fibrils remains obscure. In this work, molecular dynamics (MD) simulations have been performed to explore the conformational ensemble of the Aβ₄₂ monomer and a pentameric protofibril structure of Aβ₄₂ in the presence of 4v. The MD simulations highlighted that 4v binds preferentially at the central hydrophobic core region of the Aβ₄₂ monomer and chains D and E of the Aβ₄₂ protofibril. The dictionary of secondary structure of proteins analysis indicated that 4v retards the conformational conversion of the helix-rich structure of the Aβ₄₂ monomer into the aggregation-prone β-sheet conformation. The binding free energy calculated by the molecular mechanics Poisson-Boltzmann surface area method revealed an energetically favorable process with ΔG _binding = -44.9 ± 3.3 kcal/mol for the Aβ₄₂ monomer-4v complex. The free energy landscape analysis highlighted that the Aβ₄₂ monomer-4v complex sampled conformations with significantly higher helical contents (35 and 49%) as compared to the Aβ₄₂ monomer alone (17%). Compound 4v displayed hydrogen bonding with Gly37 (chain E) and π-π interactions with Phe19 (chain D) of the Aβ₄₂ protofibril. Further, the per-residue binding free energy analysis also highlighted that Phe19 (chain D) and Gly37 (chain E) of the Aβ₄₂ protofibril showed the maximum contribution in the binding free energy. The decreased binding affinity and residue-residue contacts between chains D and E of the Aβ₄₂ protofibril in the presence of 4v indicate destabilization of the Aβ₄₂ protofibril structure. Overall, the structural information obtained through MD simulations indicated that 4v stabilizes the native helical conformation of the Aβ₄₂ monomer and persuades a destabilization in the protofibril structure of Aβ₄₂. The results of the study will be useful in the rational design of potent inhibitors against amyloid aggregation.

Photochemical Transformation of some 3-benzyloxy-2-(benzo[b]thiophen-2-yl)-4Hchromen-4-ones: A Remote Substituent Effect

The effect on photochemical transformations of the substituents present remotely from the reaction site in 3-benzyloxy-2-(benzo[b]thiophen-2-yl)-4H-chromen-4-ones has been determined. The structure(s) of the substrates and photoproducts were established by spectroscopic techniques (UV, IR, and NMR). The substituents had profound effects on product yield and distribution. Electron withdrawing groups (EWGs) on the benzenoid moiety of the chromenone nucleus increased the yield of the photoproducts whereas electron donating groups (EDGs) decreased the yield. These results may be attributed to “state switching” of the substituents during excitation.

Progressive Stereo Locking (PSL): A Residual Dipolar Coupling Based Force Field Method for Determining the Relative Configuration of Natural Products and Other Small Molecules

Establishing the relative configuration of a bioactive natural product represents the most challenging part in determining its structure. Residual dipolar couplings (RDCs) are sensitive probes of the relative spatial orientation of internuclear vectors. We adapted a force field structure calculation methodology to allow free sampling of both R and S configurations of the stereocenters of interest. The algorithm uses a floating alignment tensor in a simulated annealing protocol to identify the conformations and configurations that best fit experimental RDC and distance restraints (from NOE and J-coupling data). A unique configuration (for rigid molecules) or a very small number of configurations (for less rigid molecules) of the structural models having the lowest chiral angle energies and reasonable magnitudes of the alignment tensor are provided as the best predictions of the unknown configuration. For highly flexible molecules, the progressive locking of their stereocenters into their statistically dominant R or S state dramatically reduces the number of possible relative configurations. The result is verified by checking that the same configuration is obtained by initiating the locking from different regions of the molecule. For all molecules tested having known configurations (with conformations ranging from mostly rigid to highly flexible), the method accurately determined the correct configuration.

Impact of modified ribose sugars on nucleic acid conformation and function

The modification of the ribofuranose in nucleic acids is a widespread method of manipulating the activity of nucleic acids. These alterations, however, impact the local conformation and chemical reactivity of the sugar. Changes in the conformation and dynamics of the sugar moiety alter the local and potentially global structure and plasticity of nucleic acids, which in turn contributes to recognition, binding of ligands and enzymatic activity of proteins. This review article introduces the conformational properties of the (deoxy)ribofuranose ring and then explores sugar modifications and how they impact local and global structure and dynamics in nucleic acids.

Finding Our Way in the Dark Proteome

The traditional structure-function paradigm has provided significant insights for well-folded proteins in which structures can be easily and rapidly revealed by X-ray crystallography beamlines. However, approximately one-third of the human proteome is comprised of intrinsically disordered proteins and regions (IDPs/IDRs) that do not adopt a dominant well-folded structure, and therefore remain "unseen" by traditional structural biology methods. This Perspective considers the challenges raised by the "Dark Proteome", in which determining the diverse conformational substates of IDPs in their free states, in encounter complexes of bound states, and in complexes retaining significant disorder requires an unprecedented level of integration of multiple and complementary solution-based experiments that are analyzed with state-of-the art molecular simulation, Bayesian probabilistic models, and high-throughput computation. We envision how these diverse experimental and computational tools can work together through formation of a "computational beamline" that will allow key functional features to be identified in IDP structural ensembles.

Comparison of structure determination methods for intrinsically disordered amyloid-β peptides

Intrinsically disordered proteins (IDPs) represent a new frontier in structural biology since the primary characteristic of IDPs is that structures need to be characterized as diverse ensembles of conformational substates. We compare two general but very different ways of combining NMR spectroscopy with theoretical methods to derive structural ensembles for the disease IDPs amyloid-β 1-40 and amyloid-β 1-42, which are associated with Alzheimer's Disease. We analyze the performance of de novo molecular dynamics and knowledge-based approaches for generating structural ensembles by assessing their ability to reproduce a range of NMR experimental observables. In addition to the comparison of computational methods, we also evaluate the relative value of different types of NMR data for refining or validating the IDP structural ensembles for these important disease peptides.

Coarse-grained simulations of protein aggregation

Protein aggregation is believed to be responsible for a number of human diseases and limits the yields of pharmaceutical proteins during production. Computer simulations can be used to develop novel experimentally testable hypotheses pertaining to aggregation. While all-atom simulations with explicit solvent are too computationally intensive to address the multitude of relevant time scales, coarse-grained models make it possible to observe the transition of monomers to an equilibrium containing aggregates. Here, we provide the reader with background information and a list of steps for setting up, performing, and analyzing computer simulations of aggregating coarse-grained (CG) proteins.

Homogeneous and heterogeneous tertiary structure ensembles of amyloid-β peptides

The interplay of modern molecular simulation and high-quality nuclear magnetic resonance (NMR) experiments has reached a fruitful stage for quantitative characterization of structural ensembles of disordered peptides. Amyloid-β 1-42 (Aβ42), the primary peptide associated with Alzheimer's disease, and fragments such as Aβ21-30 are both classified as intrinsically disordered peptides (IDPs). We use a variety of NMR observables to validate de novo molecular dynamics simulations in explicit water to characterize the tertiary structure ensemble of Aβ42 and Aβ21-30 from the perspective of their classification as IDPs. Unlike the Aβ21-30 fragment that conforms to expectations of an IDP that is primarily extended, we find that Aβ42 samples conformations reflecting all possible secondary structure categories and spans the range of IDP classifications from collapsed structured states to highly extended conformations, making it an IDP with a far more heterogeneous tertiary ensemble.

Structure and dynamics of the Abeta(21-30) peptide from the interplay of NMR experiments and molecular simulations

We combine molecular dynamics simulations and new high-field NMR experiments to describe the solution structure of the Abeta(21-30) peptide fragment that may be relevant for understanding structural mechanisms related to Alzheimer's disease. By using two different empirical force-field combinations, we provide predictions of the three-bond scalar coupling constants ((3)J(H(N)H(alpha))), chemical-shift values, (13)C relaxation parameters, and rotating-frame nuclear Overhauser effect spectroscopy (ROESY) crosspeaks that can then be compared directly to the same observables measured in the corresponding NMR experiment of Abeta(21-30). We find robust prediction of the (13)C relaxation parameters and medium-range ROESY crosspeaks by using new generation TIP4P-Ew water and Amber ff99SB protein force fields, in which the NMR validates that the simulation yields both a structurally and dynamically correct ensemble over the entire Abeta(21-30) peptide. Analysis of the simulated ensemble shows that all medium-range ROE restraints are not satisfied simultaneously and demonstrates the structural diversity of the Abeta(21-30) conformations more completely than when determined from the experimental medium-range ROE restraints alone. We find that the structural ensemble of the Abeta(21-30) peptide involves a majority population (approximately 60%) of unstructured conformers, lacking any secondary structure or persistent hydrogen-bonding networks. However, the remaining minority population contains a substantial percentage of conformers with a beta-turn centered at Val24 and Gly25, as well as evidence of the Asp23 to Lys28 salt bridge important to the fibril structure. This study sets the stage for robust theoretical work on Abeta(1-40) and Abeta(1-42), for which collection of detailed NMR data on the monomer will be more challenging because of aggregation and fibril formation on experimental timescales at physiological conditions. In addition, we believe that the interplay of modern molecular simulation and high-quality NMR experiments has reached a fruitful stage for characterizing structural ensembles of disordered peptides and proteins in general.

Perturbation of Conjugation in Allylic Lithium Compounds Due to Stereochemical Control of Internal Lithium Coordination: Crystallography, NMR, and Calculational Studies

Several allylic lithium compounds have been prepared with ligands tethered at C(2). These are with (CH(3)OCH(2)CH(2))(2)NCH(2)-, 6, 1-TMS 5, 1,3-bis(TMS) 8, and 1,1,3-tris(TMS) 9. Allylic lithiums with (CH(3)OCH(2)CH(2))(2)NCH(2)C(CH(3))(2)-, are 10, 1-TMS 11, and 1,3-bis(TMS), 12 compounds with -C(CH(3))(2)CH(2)N-((S)-(2-methoxymethyl)-pyrrolidino) at C(2) 13, 1-TMS 14, and 1,3-bis(TMS) 15. In the solid state, 8-10 and 12 are monomers, 6 and 13 are Li-bridged dimers, and 5 and 7 are polymers. In solution (NMR data), 5, 7-12, 14, and 15 are monmeric, and 6 is a dimer. All samples show lithium to be closest to one of the terminal allyl carbons in the crystal structures and to exhibit one-bond (13)C-(7)Li or (13)C(1)-(7)Li spin coupling, for the former typically ca. 3 Hz and for the latter 6-8 Hz. In every structure, the C(1)-C(2) allyl bond is longer than the C(2)-C(3) bond, and both lie between those for solvated delocalized and unsolvated localized allylic lithium compounds, respectively, as is also the case for the terminal allyl (13)C NMR shifts. Lithium lies 40-70 degrees off the axis perpendicular to the allyl plane at C(1). These effects are variable, so the trend is that the differences between the C(1)-C(2) and C(2)-C(3) bond lengths, (13)delta(3)-(13)delta(1) values, and the (13)C(1)-(7)Li or (13)C-(6)Li coupling constants all increase with decreasing values of the torsional angle that C(1)-Li makes with respect to the allyl plane.