Backbone Fluctuations Research Articles

Heterogeneous Slowdown of Dynamics in the Condensate of an Intrinsically Disordered Protein.

The high concentration of proteins and other biological macromolecules inside biomolecular condensates leads to dense and confined environments, which can affect the dynamic ensembles and the time scales of the conformational transitions. Here, we use atomistic molecular dynamics (MD) simulations of the intrinsically disordered low complexity domain (LCD) of the human fused in sarcoma (FUS) RNA-binding protein to study how self-crowding inside a condensate affects the dynamic motions of the protein. We found a heterogeneous retardation of the protein dynamics in the condensate with respect to the dilute phase, with large-amplitude motions being strongly slowed by up to 2 orders of magnitude, whereas small-scale motions, such as local backbone fluctuations and side-chain rotations, are less affected. The results support the notion of a liquid-like character of the condensates and show that different protein motions respond differently to the environment.

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations.

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

Antiangiogenic potential of phytochemicals from Clerodendrum inerme (L.) Gaertn investigated through in silico and quantum computational methods.

Suppressing vascular endothelial growth factor (VEGF), its receptor (VEGFR2), and the VEGF/VEGFR2 signaling cascade system to inhibit angiogenesis has emerged as a possible cancer therapeutic target. The present work was designed to discover and evaluate bioactive phytochemicals from the Clerodendrum inerme (L.) Gaertn plant for their anti-angiogenic potential. Molecular docking of twenty-one phytochemicals against the VEGFR-2 (PDB ID: 3VHE) protein was performed, followed by ADMET profiling and molecular docking simulations. These investigations unveiled two hit compounds, cirsimaritin (- 12.29kcal/mol) and salvigenin (- 12.14kcal/mol), with the highest binding energy values when compared to the reference drug, Sorafenib (- 15.14kcal/mol). Furthermore, only nine phytochemicals (cirsimaritin and salvigenin included) obeyed Lipinski's rule of five and passed ADMET filters. Molecular dynamics simulations run over 100ns revealed that the protein-ligand complexes remained stable with minimal backbone fluctuations. The binding free energy values of cirsimaritin (- 52.35kcal/mol) and salvigenin (- 55.89kcal/mol), deciphered by MM-GBSA analyses, further corroborated the docking interactions. The HOMO-LUMO band energy gap (ΔE) was calculated using density-functional theory (DFT) and substantiated using density of state (DOS) spectra. The chemical reactivity analyses revealed that salvigenin exhibited the highest chemical softness value (6.384eV), the lowest hardness value (0.07831eV), and the lowest ΔE value (0.1566eV), which implies salvigenin was less stable and chemically more reactive than cirsimaritin and sorafenib. These findings provide further evidence that cirsimaritin and salvigenin have the ability to prevent angiogenesis and the development of cancer. Nevertheless, more in vitro and in vivo confirmation is necessary.

Influence of counterion type on the scattering of a semiflexible polyelectrolyte.

Understanding the influence of counterion and backbone solvation on the conformational and thermodynamic properties of polyelectrolytes in solution is one of the main open challenges in polyelectrolyte science. To address this problem, we study the scattering from semidilute solutions of a semiflexible polyelectrolyte, carboxymethyl cellulose (CMC) with alkaline and tetra-alkyl-ammonium (TAA) counterions in aqueous media using small-angle neutron scattering (SANS), and small-angle X-ray scattering (SAXS), which allow us to probe concentration fluctuations of the polymer backbone and counterions. In SAXS, the calculated contrast arises primarily from the polymer backbone for both alkaline and TAA salts of CMC. In SANS, however, the contrast is dominated by the counterions for the TAA salts and the polymer backbone for the alkaline salts. Solutions are found to display a correlation peak in their scattering function, which at low concentrations is independent of counterion type. At moderate salt concentrations (c ≳ 0.1 M), the peak positions obtained from SANS and SAXS for the CMC salts with the TAA counterions differ. This divergence suggests a decoupling in the lengthscale over which the couterions and the polymer fluctuate. Upturns in the scattering intensity in the low-q region signal the presence of long-ranged compositional inhomogeneities in the solutions. The strength of these decreases with increasing counterion-solvent interaction strength, as measured by the viscosity B coefficient, and are strongest for the corresponding sodium salt of CMC.

Queuine as a potential multi-target drug for alzheimer’s disease: insights from protein dynamics

Alzheimer’s disease (AD) is a neurodegenerative disorder with a complex pathogenesis. One promising approach to treating AD is simultaneously targeting multiple aspects of the disease using a multi-target drug (MTD). In this study, multi-target drug (MTD) potential of the nutraceutical molecule Queuine was explored using molecular docking and molecular dynamics (MD) simulations with five different protein targets engaged in AD: AChE, beta-site amyloid precursor protein cleaving enzyme-1 (BACE-1), N-methyl-D-aspartate receptor (NMDAR), monoamine oxidase A (MAO-A), and Synapsin III. Queuine revealed significant binding affinities, the docking scores being −10.1, −5.97, −5.63, −8.40, and −10.56 kcal/mol for AChE, BACE-1, NMDAR, MAO-A, and Synapsin III, respectively. MD simulations showed that Queuine formed stable complexes and preserved its stability throughout the simulation, the backbone fluctuations remaining within 2.5 Å specifically in the case of the BACE-1. Elastic network model simulations and principal component analysis were carried out to illustrate the dynamics of the protein systems. Significant hinge-bending and twisting-type motions that may be relevant to function were observed around the dimerization interfaces or binding sites. Structural clustering based on PCA analysis and cross-correlation maps demonstrated that Queuine binding altered the protein dynamics more drastically in the case of highly mobile proteins NMDAR and MAO-A. We propose that the neuroprotective effect of Queuine may stem from its prominent inhibitory action on enzymes BACE-1 and AChE. Our results suggest that Queuine may serve as a promising MTD candidate for the treatment of AD. Communicated by Ramaswamy H. Sarma

Nanodisc-embedded cytochrome P450 P3A4 binds diverse ligands by distributing conformational dynamics to its flexible elements

Cytochrome P450 3A4 (CYP3A4) metabolizes a wide range of drugs and toxins. Interactions of CYP3A4 with ligands are difficult to predict due to promiscuity and conformational flexibility. To better understand CYP3A4 conformational responses to ligands we use hydrogen deuterium exchange mass spectrometry (HDX-MS) to investigate the effect of ligands on nanodisc-embedded CYP3A4. For a subset of CYP3A4-ligand complexes, differences in the low-frequency modes derived by principal component analyses of molecular dynamics trajectories mirrored the HDX-MS results. The effects of ligands are distributed to flexible elements of CYP3A4 between stretches of secondary structure. The largest effects occur in the F- and G-helices, where most ligands increase the flexibility of the F-helix and connecting loops and decrease the flexibility of the C-term of the G-helix. Most ligands affect the E-F-G, CD and HI regions of the protein. Ligand-dependent differences are observed in the A"-A' loop, BC region, E-helix, K-β1 region, proximal loop, and C-term loop. Correlated HDX responses were observed in the CD region and the C-term of the G-helix that were most pronounced for Type II ligands. Collectively, the HDX and molecular dynamics results suggest that CYP3A4 accommodates diverse binding partners by propagating local backbone fluctuations from the binding site onto the flexible regions of the enzyme via long-range interactions that are differentially modulated by ligands. In contrast to the paradigm wherein ligands decrease protein dynamics at their binding site, a wide range of ligands modestly increase CYP3A4 dynamics throughout the protein including effects remote from the active site.

Structure and Dynamics of Three Escherichia coli NfsB Nitro-Reductase Mutants Selected for Enhanced Activity with the Cancer Prodrug CB1954

Escherichia coli NfsB has been studied extensively for its potential for cancer gene therapy by reducing the prodrug CB1954 to a cytotoxic derivative. We have previously made several mutants with enhanced activity for the prodrug and characterised their activity in vitro and in vivo. Here, we determine the X-ray structure of our most active triple and double mutants to date, T41Q/N71S/F124T and T41L/N71S. The two mutant proteins have lower redox potentials than wild-type NfsB, and the mutations have lowered activity with NADH so that, in contrast to the wild-type enzyme, the reduction of the enzyme by NADH, rather than the reaction with CB1954, has a slower maximum rate. The structure of the triple mutant shows the interaction between Q41 and T124, explaining the synergy between these two mutations. Based on these structures, we selected mutants with even higher activity. The most active one contains T41Q/N71S/F124T/M127V, in which the additional M127V mutation enlarges a small channel to the active site. Molecular dynamics simulations show that the mutations or reduction of the FMN cofactors of the protein has little effect on its dynamics and that the largest backbone fluctuations occur at residues that flank the active site, contributing towards its broad substrate range.

Studies of the local assembly dynamics of the T4 bacteriophage helicase loading protein (gp59) and single-stranded (ss) DNA binding protein (gp32) at and near (iCy3)2-labeled ss-dsDNA junctions by ensemble and single-molecule polarization sweep microscopy.

The T4 bacteriophage DNA replication system (replisome) provides an excellent model for studying DNA replication, as it can self-assemble in vitro and carries out the same basic functions as those present in higher organisms. During the elongation phase of DNA replication, eight different T4 proteins, including the single-stranded DNA binding protein (gp32) and the helicase loading protein (gp59), recognize and bind to non-base-sequence specific sites to form the replisome complex. Replisome assembly is facilitated by DNA ‘breathing’, which is defined as thermally-induced DNA backbone and base fluctuations that transiently expose its template bases in suitable conformations for protein binding. Understanding the local conformations of sugar-phosphate backbones near ss-dsDNA junctions is important for elucidating the assembly mechanisms for proteins that bind and function specifically at replication junctions. To site-specifically label the DNA junctions, an exciton-coupled dimer probe, (iCy3)2, was inserted into positions at and near ss-dsDNA junction to report on sugar-phosphate backbone structure and dynamics. The optical properties of the (iCy3)2 dimer-labeled ss-dsDNA constructs were characterized using absorbance and circular dichroism (CD) spectroscopy, and the local dynamics were studied using polarization-sweep single-molecule fluorescence (PS-SMF) microscopy. We applied a relatively simple exciton-coupling model to calculate the (Cy3)2 absorbance and CD spectra. Optimization of this model to obtain a best-fit to the experimental data permits us to determine the position-specific structural parameters associated with the (iCy3)2 dimer-labeled ss-dsDNA constructs, both in the presence and absence of proteins. To interpret our PS-SMF data, we applied a kinetic network model to determine the minimum number of distinct conformational macrostates present, in addition to elucidating the time scales of their interconversion.

Using single molecule polarization-sweep spectroscopy to monitor site-specific fluctuations of internal (iCy3)2 dimer-labeled DNA constructs at ss-ds DNA junctions.

Local fluctuations of the sugar-phosphate backbones of DNA (a form of DNA ‘breathing’) play key roles in protein-DNA assembly and enzymatic function. In this work, we monitored spectroscopic signals from single-molecules of DNA fork constructs labeled with two cyanine optical probes [(iCy3)2], which were rigidly inserted into the sugar-phosphate backbones at opposite positions of complementary single-strands. We thus studied the local conformational fluctuations of the DNA sugar-phosphate backbones at specific (iCy3)2 dimer-labeled positions relative to a single-stranded (ss) double-stranded (ds) DNA junction. These experiments are based on a newly-developed single-molecule spectroscopic method that uses a linearly polarized continuous wave laser in which the polarization direction is continuously rotated at radio frequencies. The emitted fluorescence from the single-molecule sample contains information about conformational changes of the exciton-coupled (iCy3)2 dimer probe, which monitor conformational fluctuations of the vicinal DNA sugar-phosphate backbones. Our results indicate that the local conformation of the DNA sugar-phosphate backbones at positions near ss-dsDNA junctions adopts four topologically-relevant macrostates. The probability that a particular macrostate is occupied depends on backbone position. We apply a kinetic network approach to interpret our observations of DNA backbone fluctuations at and near ss-dsDNA junctions, and we propose conformational assignments based on exciton-coupling models for the macrostates that we observe. This polarization-sweep single-molecule fluorescence method and analysis can be applied to study protein binding and macromolecular function at and near ss-dsDNA junctions, thus providing site-specific information about the free energy landscapes at various functionally relevant positions within such macromolecular complexes.

Effects of PCSK9 missense variants on molecular conformation and biological activity in transfected HEK293FT cells

PCSK9 gain-of-function (GOF) variants increase degradation of low-density lipoprotein receptor (LDLR) and are potentially associated with Familial Hypercholesterolemia (FH). This study aimed to explore the effects of PCSK9 missense variants on protein structure and interactions with LDLR using molecular modeling analyses and in vitro functional studies. Variants in FH-related genes were identified in a Brazilian FH cohort using an exon-target gene sequencing strategy. Eight PCSK9 missense variants in pro- [p.(E32K) and p.(E57K)], catalytic [p.(R237W), p.(P279T) and p.(A443T)], and C-terminal histidine-cysteine rich (CHR) [p.(R469W), p.(Q619P) and p.(R680Q)] domains were identified. Molecular dynamics analyses revealed that GOF variants p.(E32K) and p.(R469W) increased extreme motions in PCSK9 amino acid backbone fluctuations and affected Hbond and water bridge interactions between the pro-domain and CM1 region of the CHR domain. HEK293FT cells transfected with plasmids carrying p.(E32K) and p.(R469W) variants reduced LDLR expression (8.7 % and 14.8 %, respectively) compared to wild type (p < 0.05) but these GOF variants did not affect PCSK9 expression and secretion. The missense variants p.(P279T) and p.(Q619P) also reduced protein stability and altered Hbond interactions. In conclusion, PCSK9 p.(E32K), p.(R469W), p.(P279T) and p.(Q619P) variants disrupt intramolecular interactions that are essential for PCSK9 structural conformation and biological activity and may have a potential role in FH pathogenesis.

Measurement of Protein Dynamics from Site Directed Cu(II) Labeling

AbstractThis review describes the use of Electron Paramagnetic Resonance (EPR) to measure residue specific dynamics in proteins with a specific focus on Cu(II)‐based spin labels. First, we outline approaches used to measure protein motion by nitroxide‐based spin labels. Here, we describe conceptual details and outline challenges that limit the use of nitroxide spin labels to solvent‐exposed α‐helical sites. The bulk of this review showcases the use of newly developed Cu(II)‐based protein labels. In this approach, the strategic mutation of native residues on a protein to generate two neighboring Histidine residues (i.e., the dHis motif) is exploited to enable a rigid site‐selective binding of a Cu(II) complex. The chelation of the Cu(II) complex to dHis directly anchors the Cu(II) spin label to the protein backbone. The improvement in rigidity expands both the spin‐labeling toolkit as well as the resolution of many EPR measurements. We describe how EPR measurements of the Cu(II) label directly reflect backbone motion and fluctuations. The EPR are complemented by Molecular Dynamics simulations. Finally, the dHis motif provides access to the measurement of site‐specific dynamics at both α‐helices and β‐sheets. The review outlines the limitations of the dHis method and provides an outlook for future developments.

Interdisciplinary in silico studies to understand in-depth molecular level mechanism of drug resistance involving NS3-4A protease of HCV

Hepatitis C virus (HCV) causes hepatitis, a life-threatening disease responsible for liver cirrhosis. Urgent measures have been taken to develop therapeutics against this deadly pathogen. NS3/4A protease is an extremely important target. A series of inhibitors have been developed against this viral protease including Faldaprevir. Unfortunately, the error-prone viral RNA polymerase causes the emergence of resistance, thereby causing reduced effectiveness of those peptidomimetic inhibitors. Among the drug resistant variants, three single amino acid residues (R155, A156 and D168) are notable for their presence in clinical isolates and also their effectivity against most of the known inhibitors in clinical development. Therefore, it is crucial to understand the mechanistic role of those drug resistant variants while designing potent novel inhibitors. In this communication, we have deeply analyzed through using in silico studies to understand the molecular mechanism of alteration of inhibitor binding between wild type and its R155K, A156V and D168V variants. Principal component analysis was carried to identify the backbone fluctuations of important residues in HCV NS3/4A responsible for the inhibitor binding and maintaining drug resistance. Free energy landscape as a function of the principal components has been used to identify the stability and conformation of the key residues that regulate inhibitor binding and their impact in developing drug resistance. Our findings are consistent with the trend of experimental results. The observations are also true in case of other Faldaprevir-like peptidomimetic inhibitors. Understanding this binding mechanism would be significant for the development of novel inhibitors with less susceptibility towards drug resistance. Communicated by Ramaswamy H. Sarma

Tuberculosis and HIV responses threatened by nCOVID-19: A situation prompting an in silico investigation of reported MbtA inhibitors for combined inhibition of SARS-CoV-2 and HIV-TB co-infection

The menace of infectious diseases has constantly been a reason of concern for humankind since time immemorial. As evident by the name, infectious diseases can infect a huge population within a short period, leading to an eruption of pandemics and epidemics. The present human era is fortunate enough to have a wide array of readily available drugs that help cure and prevent various diseases. Moreover, the scientific community has always responded to the needs of society through its drug discovery and development programs. The co-existence of multiple diseases calls forth the scientific community to design and develop drugs that could have a broad spectrum of activity. In this perspective, our goal was to investigate the potential of reported MbtA inhibitors (antitubercular molecules) in inhibiting HIV-1 RT and nCovid-19-RdRp and eventually leading to the identification of a multi-targeted ligand (triple co-infection inhibitor). In this study, the primary success was attained by capitalizing on the structure-based virtual screening drug discovery approach. Results were quite promising. Molecular docking results showed that GV17 interacted strongly with the active site residues of both the target proteins (HIV-1 RT and nCOVID-19-RdRp). Moreover, the docking score of GV17 was more than that of the internal ligands of both the target proteins, which indicates a firm binding. Molecular dynamics further validated these results as identical amino acid residues were observed in the protein’s docked pose with the ligand. The detailed atomic interactions of ligand GV17 with the protein residues have been discussed. Overall, the protein–ligand complexes remained stable throughout the simulation, and the system’s backbone fluctuations were modest. MM-GBSA analysis revealed free binding energy of − 72.30 ± 7.85 kcal/mol and − 65.40 ± 7.25 kcal/mol for 1RT2 and 7BV2, respectively. The more negative binding energy indicates a stronger affinity of GV17 with both the receptors. GV17 also gave satisfactory predictive in silico ADMET results. Overall, this computational study identified GV17 as a potential HIT molecule and findings can open up a new avenue to explore and develop inhibitors against nCOVID-19-HIV-TB triple-infections.

Evolutionary conservation of protein dynamics: insights from all-atom molecular dynamics simulations of ‘peptidase’ domain of Spt16

Protein function is encoded in its sequence, manifested in its three-dimensional structure, and facilitated by its dynamics. Studies have suggested that protein structures with higher sequence similarity could have more similar patterns of dynamics. However, such studies of protein dynamics within and across protein families typically rely on coarse-grained models, or approximate metrics like crystallographic B-factors. This study uses µs scale molecular dynamics (MD) simulations to explore the conservation of dynamics among homologs of ∼50 kDa N-terminal module of Spt16 (Spt16N). Spt16N from Saccharomyces cerevisiae (Sc-Spt16N) and three of its homologs with 30–40% sequence identities were available in the PDB. To make our data-set more comprehensive, the crystal structure of an additional homolog (62% sequence identity with Sc-Spt16N) was solved at 1.7 Å resolution. Cumulative MD simulations of 6 µs were carried out on these Spt16N structures and on two additional protein structures with varying degrees of similarity to it. The simulations revealed that correlation in patterns of backbone fluctuations vary linearly with sequence identity. This trend could not be inferred using crystallographic B-factors. Further, normal mode analysis suggested a similar pattern of inter-domain (inter-lobe) motions not only among Spt16N homologs, but also in the M24 peptidase structure. On the other hand, MD simulation results highlighted conserved motions that were found unique for Spt16N protein, this along with electrostatics trends shed light on functional aspects of Spt16N. Communicated by Ramaswamy H. Sarma.

Low-frequency collective motion of DNA-binding domain defines p53 function.

Most mutations in the DNA-binding domain (DBD) of p53 inactivate or rescue the protein function interacting with the minor groove of DNA. However, how the conformation changes propagating from the mutation sites result in distinct molecular recognition is still not well understood. As the protein mobility is an intrinsic property encrypted in its primary structure, we examined if different structures of wild-type and mutant p53 core domains display any unique patterns of intrinsic mobility. Normal mode calculation was employed to characterize the collective dynamics of DBD in p53 monomer and tetramer as well as their mutants. Intriguingly, the low-frequency collective motions of DBD show similar patterns between the wild-type protein and the rescued mutants. The analysis on atomic backbone fluctuations and low-frequency vibration mode statistics does further support the correlation between the intrinsic collective motion of DBD and the p53 protein function. The mutations in the DBD influence the low-frequency vibration of the p53 tetramer via the change of the collective motions among its four monomers. These findings thus provide new insights for understanding the physical mechanism of p53 protein structure-function relationship and help find the small molecule drug to modulate protein dynamic for disease therapy.

Pseudokinases repurpose flexibility signatures associated with the protein kinase fold for noncatalytic roles.

The bilobal protein kinase-like fold in pseudokinases lack one or more catalytic residues, conserved in canonical protein kinases, and are considered enzymatically deficient. Tertiary structures of pseudokinases reveal that their loops topologically equivalent to activation segments of kinases adopt contracted configurations, which is typically extended in active conformation of kinases. Herein, anisotropic network model based normal mode analysis (NMA) was conducted on 51 active conformation structures of protein kinases and 26 crystal structures of pseudokinases. Our observations indicate that although backbone fluctuation profiles are similar for individual kinase-pseudokinase families, low intensity mean square fluctuations in pseudo-activation segment and other sub-structures impart rigidity to pseudokinases. Analyses of collective motions from functional modes reveal that pseudokinases, compared to active kinases, undergo distinct conformational transitions using the same structural fold. All-atom NMA of protein kinase-pseudokinase pairs from each family, sharing high amino acid sequence identities, yielded distinct community clusters, partitioned by residues exhibiting highly correlated fluctuations. It appears that atomic fluctuations from equivalent activation segments guide community membership and network topologies for respective kinase and pseudokinase. Our findings indicate that such adaptations in backbone and side-chain fluctuations render pseudokinases competent for catalysis-independent roles.

Tomatidine and Patchouli Alcohol as Inhibitors of SARS-CoV-2 Enzymes (3CLpro, PLpro and NSP15) by Molecular Docking and Molecular Dynamics Simulations.

Considering the current dramatic and fatal situation due to the high spreading of SARS-CoV-2 infection, there is an urgent unmet medical need to identify novel and effective approaches for prevention and treatment of Coronavirus disease (COVID 19) by re-evaluating and repurposing of known drugs. For this, tomatidine and patchouli alcohol have been selected as potential drugs for combating the virus. The hit compounds were subsequently docked into the active site and molecular docking analyses revealed that both drugs can bind the active site of SARS-CoV-2 3CLpro, PLpro, NSP15, COX-2 and PLA2 targets with a number of important binding interactions. To further validate the interactions of promising compound tomatidine, Molecular dynamics study of 100 ns was carried out towards 3CLpro, NSP15 and COX-2. This indicated that the protein-ligand complex was stable throughout the simulation period, and minimal backbone fluctuations have ensued in the system. Post dynamic MM-GBSA analysis of molecular dynamics data showed promising mean binding free energy 47.4633 ± 9.28, 51.8064 ± 8.91 and 54.8918 ± 7.55 kcal/mol, respectively. Likewise, in silico ADMET studies of the selected ligands showed excellent pharmacokinetic properties with good absorption, bioavailability and devoid of toxicity. Therefore, patchouli alcohol and especially, tomatidine may provide prospect treatment options against SARS-CoV-2 infection by potentially inhibiting virus duplication though more research is guaranteed and secured.

Bridging the Bridging Imidazolate in the Bimetallic Center of the Cu/Zn SOD1 and ALS.

Metallation status of human Cu/Zn superoxide dismutase 1 (SOD1) plays a pivotal role in the pathogenesis of amyotrophic lateral sclerosis (ALS). All of the amino acids found in the bimetallic center have been associated with ALS except for two positions. H63 which forms the bridging imidazolate ion in the bimetallic center and K136 which is not directly involved in coordination but located in the bimetallic center were not reported to be mutated in any of the identified ALS cases. In this study, we investigated the structure and flexibility of five SOD1 variants by using classical molecular dynamics simulations. These variants include three substitutions on the non-ALS-linked positions; H63A, H63R, K136A and ALS-linked positions; G37R, H46R/H48D. We have generated four systems for each variant differing in metallation and presence of the intramolecular disulfide bond. Overall, a total of 24 different dimers including the wild-type were generated and simulated at two temperatures, 298 and 400 K. We have monitored backbone mobility, fluctuations and compactness of the dimer structures to assess whether the hypothetical mutations would behave similar to the ALS-linked variants. Results showed that particularly two mutants, H63R and K136A, drastically affected the dimer dynamics by increasing the fluctuations of the metal binding loops compared with the control mutations. Further, these variants resulted in demetallation of the dimers, highlighting probable ALS toxicity that could be elicited by the SOD1 variants of H63R and K136A. Overall, this study bridges two putative SOD1 positions in the metallic center and ALS, underlining the potential use of atomistic simulations for studying disease variants.

Met35 Oxidation Hinders Aβ25-35 Peptide Aggregation within the Dimyristoylphosphatidylcholine Bilayer.

Using all-atom explicit solvent replica exchange molecular dynamics simulations, we studied the aggregation of oxidized (ox) Aβ25-35 peptides into dimers mediated by the zwitterionic dimyristoylphosphatidylcholine (DMPC) lipid bilayer. By comparing oxAβ25-35 aggregation with that observed for reduced and phosphorylated Aβ25-35 peptides, we elucidated plausible impact of post-translational modifications on cytotoxicity of Aβ peptides involved in Alzheimer's disease. We found that Met35 oxidation reduces helical propensity in oxAβ25-35 peptides bound to the lipid bilayer and enhances backbone fluctuations. These factors destabilize the wild-type head-to-tail dimer interface and lower the aggregation propensity. Met35 oxidation diversifies aggregation pathways by adding monomeric species to the bound conformational ensemble. The oxAβ25-35 dimer becomes partially expelled from the DMPC bilayer and as a result inflicts limited disruption to the bilayer structure compared to wild-type Aβ25-35. Interestingly, the effect of Ser26 phosphorylation is largely opposite, as it preserves the wild-type head-to-tail aggregation interface and strengthens, not weakens, aggregation propensity. The differing effects can be attributed to the sequence locations of these post-translational modifications, since in contrast to Ser26 phosphorylation, Met35 oxidation directly affects the wild-type C-terminal aggregation interface. A comparison with experimental data is provided.

Submillisecond Conformational Transitions of Short Single-Stranded DNA Lattices by Photon Correlation Single-Molecule Förster Resonance Energy Transfer.

Thermally driven conformational fluctuations (or "breathing") of DNA play important roles in the function and regulation of the "macromolecular machinery of genome expression." Fluctuations in double-stranded (ds) DNA are involved in the transient exposure of pathways to protein binding sites within the DNA framework, leading to the binding of regulatory proteins to single-stranded (ss) DNA templates. These interactions often require that the ssDNA sequences, as well as the proteins involved, assume transient conformations critical for successful binding. Here, we use microsecond-resolved single-molecule Förster resonance energy transfer (smFRET) experiments to investigate the backbone fluctuations of short [oligo(dT)n] templates within DNA constructs that also serve as models for ss-dsDNA junctions. Such junctions, together with the attached ssDNA sequences, are involved in interactions with the ssDNA binding (ssb) proteins that control and integrate the functions of DNA replication complexes. We analyze these data using a chemical network model based on multiorder time-correlation functions and probability distribution functions that characterize the kinetic and thermodynamic behavior of the system. We find that the oligo(dT)n tails of ss-dsDNA constructs interconvert, on submillisecond time scales, between three macrostates with distinctly different end-to-end distances. These are (i) a "compact" macrostate that represents the dominant species at equilibrium; (ii) a "partially extended" macrostate that exists as minority species; and (iii) a "highly extended" macrostate that is present in trace amounts. We propose a model for ssDNA secondary structure that advances our understanding of how spontaneously formed nucleic acid conformations may facilitate the activities of ssDNA-associating proteins.