MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology

Modeling pH-Dependent NMR Chemical Shift Perturbations in Peptides

Modeling the pH dependence of protein and peptide chemical shifts outside the range of physiological values (6.5-7) is key to understanding structure-function relationships of these systems. These capabilities are largely not available in current chemical shift prediction software. In this study, we utilize a combination of molecular dynamics and quantum mechanics to investigate the through-space and through-bond contributions of protonation-dependent chemical shift perturbations (CSPs) in model tripeptides. By altering the protonation state of the titratable group in the tripeptides, we observe a notable difference in the conformational ensembles and attendantly compute significant CSPs for all nuclei near the site of protonation. We thus demonstrate the ability to recapitulate experimental pH-dependent CSPs with good agreement (R = 0.85, 0.99, and 0.98 for ¹³C, ¹⁵N, and ¹H, respectively). Broadly, we provide the groundwork for incorporating pH effects into empirical and semiempirical chemical shift predictors.

β-Methylamino-L-alanine substitution of serine in SOD1 suggests a direct role in ALS etiology

Exposure to the environmental toxin β-methylamino-L-alanine (BMAA) is linked to amyotrophic lateral sclerosis (ALS), but its disease-promoting mechanism remains unknown. We propose that incorporation of BMAA into the ALS-linked protein Cu,Zn superoxide dismutase (SOD1) upon translation promotes protein misfolding and aggregation, which has been linked to ALS onset and progression. Using molecular simulation and predictive energetic computation, we demonstrate that substituting any serine with BMAA in SOD1 results in structural destabilization and aberrant dynamics, promoting neurotoxic SOD1 aggregation. We propose that translational incorporation of BMAA into SOD1 is directly responsible for its toxicity in neurodegeneration, and BMAA modification of SOD1 may serve as a biomarker of ALS.

Sequence Tolerance of a Single-Domain Antibody with a High Thermal Stability: Comparison of Computational and Experimental Fitness Profiles

The sequence fitness of a llama single-domain antibody with an unusually high thermal stability is explored by a combined computational and experimental study. Starting with the X-ray crystallographic structure, RosettaBackrub simulations were applied to model sequence-structure tolerance profiles and identify key substitution sites. From the model calculations, an experimental site-directed mutagenesis was used to produce a panel of mutants, and their melting temperatures were determined by thermal denaturation. The results reveal a sequence fitness of an excess stability of approximately 12 °C, a value taken from a decrease in the melting temperature of an electrostatic charge-reversal substitution in the CRD3 without a deleterious effect on the binding affinity to the antigen. The tolerance for the disruption of antigen recognition without loss in the thermal stability was demonstrated by the introduction of a proline in place of a tyrosine in the CDR2, producing a mutant that eliminated binding. To further assist the sequence design and the selection of engineered single-domain antibodies, an assessment of different computational strategies is provided of their accuracy in the detection of substitution "hot spots" in the sequence tolerance landscape.

Driven to near-experimental accuracy by refinement via molecular dynamics simulations

Protein model refinement has been an essential part of successful protein structure prediction. Molecular dynamics simulation-based refinement methods have shown consistent improvement of protein models. There had been progress in the extent of refinement for a few years since the idea of ensemble averaging of sampled conformations emerged. There was little progress in CASP12 because conformational sampling was not sufficiently diverse due to harmonic restraints. During CASP13, a new refinement method was tested that achieved significant improvements over CASP12. The new method intended to address previous bottlenecks in the refinement problem by introducing new features. Flat-bottom harmonic restraints replaced harmonic restraints, sampling was performed iteratively, and a new scoring function and selection criteria were used. The new protocol expanded conformational sampling at reduced computational costs. In addition to overall improvements, some models were refined significantly to near-experimental accuracy.

Structural Model of the ETR1 Ethylene Receptor Transmembrane Sensor Domain

The structure, mechanism of action and copper stoichiometry of the transmembrane sensor domain of the plant ethylene receptor ETR1 and homologs have remained elusive, hampering the understanding on how the perception of the plant hormone ethylene is transformed into a downstream signal. We generated the first structural model of the transmembrane sensor domain of ETR1 by integrating ab initio structure prediction and coevolutionary information. To refine and independently validate the model, we determined protein-related copper stoichiometries on purified receptor preparations and explored the helix arrangement by tryptophan scanning mutagenesis. All-atom molecular dynamics simulations of the dimeric model reveal how ethylene can bind proximal to the copper ions in the receptor, illustrating the initial stages of the ethylene perception process.

Predicting Protein-Protein Interfaces that Bind Intrinsically Disordered Protein Regions

A long-standing goal in biology is the complete annotation of function and structure on all protein-protein interactions, a large fraction of which is mediated by intrinsically disordered protein regions (IDRs). However, knowledge derived from experimental structures of such protein complexes is disproportionately small due, in part, to challenges in studying interactions of IDRs. Here, we introduce IDRBind, a computational method that by combining gradient boosted trees and conditional random field models predicts binding sites of IDRs with performance approaching state-of-the-art globular interface predictions, making it suitable for proteome-wide applications. Although designed and trained with a focus on molecular recognition features, which are long interaction-mediating-elements in IDRs, IDRBind also predicts the binding sites of short peptides more accurately than existing specialized predictors. Consistent with IDRBind's specificity, a comparison of protein interface categories uncovered uniform trends in multiple physicochemical properties, positioning molecular recognition feature interfaces between peptide and globular interfaces.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Ligand binding effects on the activation of the EGFR extracellular domain

The epidermal growth factor receptor (EGFR) is one of the most common target proteins in anti-cancer therapy. The binding of the EGF ligand to the EGFR extracellular domain (EGFR-ECD) promotes its inactive-to-active conformational transition (activation) but the relevant detailed mechanism remains elusive still. Here, the structural characterization and energetics of the EGFR-ECD conformational transition with and without the binding of the EGF are quantitatively explored using an innovative enhanced sampling MD simulation method. Intriguingly, the EGF offers hydrophobic interactions (e.g., EGF residues of Tyr44 and Leu47) and electrostatic interactions (e.g., the EGF residues of Glu5, Asp11, Asp17, and Arg41) to play a dominant role in dragging domain III to close the ligand binding domain gap. Subsequently, the correlation between domains III and II is enhanced through salt-bridges among Glu376, Arg403, and Arg405 from domain III and Glu293, Glu295, and Arg300 from domain II. Finally, the structural bending of domain II is regulated to facilitate the disengagement of domain II from domain IV. In this regard, the functional conformational transition of EGFR-ECD is a consequence of the cooperative motion of protein domains driven by the EGF ligand binding. The present study shows a detailed scenario of the EGF induced activation of EGFR-ECD and provides valuable information for drug discovery targeting the EGFR.

Development and Testing of the OPLS-AA/M Force Field for RNA

Significant improvements have been made to the OPLS-AA force field for modeling RNA. New torsional potentials were optimized based on density functional theory (DFT) scans at the ωB97X-D/6-311++G(d,p) level for potential energy surfaces of the backbone α and γ dihedral angles. In combination with previously reported improvements for the sugar puckering and glycosidic torsion terms, the new force field was validated through diverse molecular dynamics simulations for RNAs in aqueous solution. Results for dinucleotides and tetranucleotides revealed both accurate reproduction of 3 J couplings from NMR and the avoidance of several unphysical states observed with other force fields. Simulations of larger systems with noncanonical motifs showed significant structural improvements over the previous OPLS-AA parameters. The new force field, OPLS-AA/M, is expected to perform competitively with other recent RNA force fields and to be compatible with OPLS-AA models for proteins and small molecules.

Prediction of Membrane Permeation of Drug Molecules by Combining an Implicit Membrane Model with Machine Learning

Lipid membrane permeation of drug molecules was investigated with Heterogeneous Dielectric Generalized Born (HDGB)-based models using solubility-diffusion theory and machine learning. Free energy profiles were obtained for neutral molecules by the standard HDGB and Dynamic HDGB (DHDGB) to account for the membrane deformation upon insertion of drugs. We also obtained hybrid free energy profiles where the neutralization of charged molecules was taken into account upon membrane insertion. The evaluation of the predictions was done against experimental permeability coefficients from Parallel Artificial Membrane Permeability Assays (PAMPA), and effects of partial charge sets, CGenFF, AM1-BCC, and OPLS, on the performance of the predictions were discussed. (D)HDGB-based models improved the predictions over the two-state implicit membrane models, and partial charge sets seemed to have a strong impact on the predictions. Machine learning increased the accuracy of the predictions, although it could not outperform the physics-based approach in terms of correlations.

Glutamine-451 Confers Sensitivity to Oxidative Inhibition and Heme-Thiolate Sulfenylation of Cytochrome P450 4B1

Human cytochrome P450 (P450) family 4 enzymes are involved in the metabolism of fatty acids and the bioactivation of carcinogenic arylamines and toxic natural products, e.g., 4-ipomeanol. These and other drug-metabolizing P450s are redox sensitive, showing a loss of activity resulting from preincubation with H2O2 and recovery with mild reducing agents [Albertolle, M. W., et al. (2017) J. Biol. Chem. 292, 11230-11242]. The inhibition is due to sulfenylation of the heme-thiolate ligand, as determined by chemopreoteomics and spectroscopy. This phenomenon may have implications for chemical toxicity and observed disease-drug interactions, in which the decreased metabolism of P450 substrates occurs in patients with inflammatory diseases (e.g., influenza and autoimmunity). Human P450 1A2 was determined to be redox insensitive. To determine the mechanism underlying the differential redox sensitivity, molecular dynamics (MD) simulations were employed using the crystal structure of rabbit P450 4B1 (Protein Data Bank entry 5T6Q ). In simulating either the thiolate (Cys-S-) or the sulfenic acid (Cys-SOH) at the heme ligation site, MD revealed Gln-451 in either an "open" or "closed" conformation, respectively, between the cytosol and heme-thiolate cysteine. Mutation to either an isosteric leucine (Q451L) or glutamate (Q451E) abrogated the redox sensitivity, suggesting that this "open" conformation allows for reduction of the sulfenic acid and religation of the thiolate to the heme iron. In summary, MD simulations suggest that Gln-451 in P450 4B1 adopts conformations that may stabilize and protect the heme-thiolate sulfenic acid; mutating this residue destabilizes the interaction, producing a redox insensitive enzyme.

Applications of Ensemble Docking in Potential Inhibitor Screening for Mycobacterium tuberculosis Isocitrate Lyase Using a Local Plant Database

Isocitrate lyase (ICL) is a persistent factor for the survival of dormant stage Mycobacterium tuberculosis (MTB), thus a potential drug target for tuberculosis treatment. In this work, ensemble docking approach was used to screen for potential inhibitors of ICL. The ensemble conformations of ICL active site were obtained from molecular dynamics simulation on three dimer form systems, namely the apo ICL, ICL in complex with metabolites (glyoxylate and succinate), and ICL in complex with substrate (isocitrate). Together with the ensemble conformations and the X-ray crystal structures, 22 structures were used for the screening against Malaysian Natural Compound Database (NADI). The top 10 compounds for each ensemble conformation were selected. The number of compounds was then further narrowed down to 22 compounds that were within the Lipinski's Rule of Five for drug-likeliness and were also docked into more than one ensemble conformation. Theses 22 compounds were furthered evaluate using whole cell assay. Some compounds were not commercially available; therefore, plant crude extracts were used for the whole cell assay. Compared to itaconate (the known inhibitor of ICL), crude extracts from Manilkara zapota, Morinda citrifolia, Vitex negundo, and Momordica charantia showed some inhibition activity. The MIC/MBC value were 12.5/25, 12.5/25, 0.78/1.6, and 0.39/1.6 mg/mL, respectively. This work could serve as a preliminary study in order to narrow the scope for high throughput screening in the future.

Allysine modifications perturb tropoelastin structure and mobility on a local and global scale

Elastin provides elastic tissues with resilience through stretch and recoil cycles, and is primarily made of its extensively cross-linked monomer, tropoelastin. Here, we leverage the recently published full atomistic model of tropoelastin to assess how allysine modifications, which are essential to cross-linking, contribute to the dynamics and structural changes that occur in tropoelastin in the context of elastin assembly. We used replica exchange molecular dynamics to generate structural ensembles of allysine containing tropoelastin. We conducted principal component analysis on these ensembles and found that the molecule departs from the canonical structural ensemble. Furthermore, we showed that, while the canonical scissors-twist movement was retained, new movements emerged that deviated from those of the wild type protein, providing evidence for the involvement of a variety of molecular motions in elastin assembly. Additionally, we highlighted secondary structural changes and linked these perturbations to the longevity of specific salt bridges. We propose a model where allysines in tropoelastin contribute to hierarchical elastin assembly through global and local perturbations to molecular structure and dynamics.

•

converting lysine to allysine by lysyl oxidases is needed to generate crosslinks between tropoelastin molecules in order to make elastin

•

structural changes in the intact tropoelastin molecule ensue where modified tropoelastin molecules structurally depart from the canonical ensemble

•

new molecular motions deviate from those of unmodified tropoelastin

•

persistence times of specific salt bridges contribute to these perturbations

•

allysines in tropoelastin contribute to hierarchical elastin assembly through global and local perturbations to molecular structure and dynamics

Impact of branching on the conformational heterogeneity of the lipopolysaccharide from Klebsiella pneumoniae: Implications for vaccine design

Resistance of Klebsiella pneumoniae (KP) to antibiotics has motivated the development of an efficacious KP human vaccine that would not be subject to antibiotic resistance. Klebsiella lipopolysaccharide (LPS) associated O polysaccharide (OPS) types have provoked broad interest as a vaccine antigen as there are only 4 that predominate worldwide (O1, O2a, O3, O5). Klebsiella O1 and O2 OPS are polygalactans that share a common D-Gal-I structure, for which a variant D-Gal-III was recently discovered. To understand the potential impact of this variability on antigenicity, a detailed molecular picture of the conformational differences associated with the addition of the D-Gal-III (1 → 4)-α-Galp branch is presented using enhanced-sampling molecular dynamics simulations. In D-Gal-I two major conformational states are observed while the presence of the 1 → 4 branch in D-Gal-III resulted in only a single dominant extended state. Stabilization of the more folded states in D-Gal-I is due to a O4-H⋯O2 hydrogen bond in the linear backbone that cannot occur in D-Gal-III as the O4 is in the Galp(1 → 4)Galp glycosidic linkage. The impact of branching in D-Gal-III also significantly decreases the accessibility of the monosaccharides in the linear backbone region of D-Gal-I, while the accessibility of the terminal D-Gal-II region of the OPS is not substantially altered. The present results suggest that a vaccine that targets both the D-Gal-I and D-Gal-III LPS can be developed by using D-Gal-III as the antigen combined with cross-reactivity experiments using the Gal-II polysaccharide to assure that this region of the LPS is the primary epitope of the antigen.

Dynamic Effects of CYP2D6 Genetic Variants in a Set of Poor Metaboliser Patients with Infiltrating Ductal Cancer Under Treatment with Tamoxifen

Breast cancer is a group of multigenic diseases. It is the most common cancer diagnosed among women worldwide and is often treated with tamoxifen. Tamoxifen is catalysed by cytochrome P450 2D6 (CYP2D6), and inter-individual variations in the enzyme due to single nucleotide polymorphisms (SNPs) could alter enzyme activity. We evaluated SNPs in patients from Colombia in South America who were receiving tamoxifen treatment for breast cancer. Allelic diversity in the CYP2D6 gene was found in the studied population, with two patients displaying the poor-metaboliser phenotype. Molecular dynamics and trajectory analyses were performed for CYP2D6 from these two patients, comparing it with the common allelic form (CYP2D6*1). Although we found no significant structural change in the protein, its dynamics differ significantly from those of CYP2D6*1, the effect of such differential dynamics resulting in an inefficient enzyme with serious implications for tamoxifen-treated patients, increasing the risk of disease relapse and ineffective treatment.

Benchmarking of different molecular docking methods for protein-peptide docking

Background

Molecular docking studies on protein-peptide interactions are a challenging and time-consuming task because peptides are generally more flexible than proteins and tend to adopt numerous conformations. There are several benchmarking studies on protein-protein, protein-ligand and nucleic acid-ligand docking interactions. However, a series of docking methods is not rigorously validated for protein-peptide complexes in the literature. Considering the importance and wide application of peptide docking, we describe benchmarking of 6 docking methods on 133 protein-peptide complexes having peptide length between 9 to 15 residues. The performance of docking methods was evaluated using CAPRI parameters like FNAT, I-RMSD, L-RMSD.

Result

Firstly, we performed blind docking and evaluate the performance of the top docking pose of each method. It was observed that FRODOCK performed better than other methods with average L-RMSD of 12.46 Å. The performance of all methods improved significantly for their best docking pose and achieved highest average L-RMSD of 3.72 Å in case of FRODOCK. Similarly, we performed re-docking and evaluated the performance of the top and best docking pose of each method. We achieved the best performance in case of ZDOCK with average L-RMSD 8.60 Å and 2.88 Å for the top and best docking pose respectively. Methods were also evaluated on 40 protein-peptide complexes used in the previous benchmarking study, where peptide have length up to 5 residues. In case of best docking pose, we achieved the highest average L-RMSD of 4.45 Å and 2.09 Å for the blind docking using FRODOCK and re-docking using AutoDock Vina respectively.

Conclusion

The study shows that FRODOCK performed best in case of blind docking and ZDOCK in case of re-docking. There is a need to improve the ranking of docking pose generated by different methods, as the present ranking scheme is not satisfactory. To facilitate the scientific community for calculating CAPRI parameters between native and docked complexes, we developed a web-based service named PPDbench (http://webs.iiitd.edu.in/raghava/ppdbench/).

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2449-y) contains supplementary material, which is available to authorized users.

Characterizing the structural and thermodynamic properties of Aβ42 and Aβ40

The self-assembly of amyloid-beta (Aβ) proteins in aqueous extracellular environments is implicated in Alzheimer's disease. Among several alloforms of Aβ proteins differing in sequence length, the 42- and 40-residue forms (Aβ42 and Aβ40) are the most abundant ones in the human body. Although the only difference is the additional I₄₁A₄₂ residues in the C-terminus, Aβ42 exhibits more aggregation tendency and stronger neurotoxicity than Aβ40. Here, we investigate the molecular factors that confer more aggregation potential to Aβ42 than to Aβ40 based on molecular dynamics simulations combined with solvation thermodynamic analyses. It is observed that the most salient structural feature of Aβ42 relative to Aβ40 is the more enhanced β-sheet forming tendency, in particular in the C-terminal region. While such a structural characteristic of Aβ42 will certainly serve to facilitate the formation of aggregate species rich in β-sheet structure, we also detect its interesting thermodynamic consequence. Indeed, we find from the decomposition analysis that the C-terminal region substantially increases the solvation free energy (i.e., overall "hydrophobicity") of Aβ42, which is caused by the dehydration of the backbone moieties showing the enhanced tendency of forming the β-structure. Together with the two additional hydrophobic residues (I₄₁A₄₂), this leads to the higher solvation free energy of Aβ42, implying the larger water-mediated attraction toward the self-assembly. Thus, our computational results provide structural and thermodynamic grounds on why Aβ42 has more aggregation propensity than Aβ40 in aqueous environments.

Structural Insights into Characterizing Binding Sites in Epidermal Growth Factor Receptor Kinase Mutants

Over the last two decades epidermal growth factor receptor (EGFR) kinase has become an important target to treat nonsmall cell lung cancer (NSCLC). Currently, three generations of EGFR kinase-targeted small molecule drugs have been FDA approved. They nominally produce a response at the start of treatment and lead to a substantial survival benefit for patients. However, long-term treatment results in acquired drug resistance and further vulnerability to NSCLC. Therefore, novel EGFR kinase inhibitors that specially overcome acquired mutations are urgently needed. To this end, we carried out a comprehensive study of different EGFR kinase mutants using a structural systems pharmacology strategy. Our analysis shows that both wild-type and mutated structures exhibit multiple conformational states that have not been observed in solved crystal structures. We show that this conformational flexibility accommodates diverse types of ligands with multiple types of binding modes. These results provide insights for designing a new generation of EGFR kinase inhibitor that combats acquired drug-resistant mutations through a multiconformation-based drug design strategy.

Rapid and accurate structure-based therapeutic peptide design using GPU accelerated thermodynamic integration

Peptide-based therapeutics are an alternative to small molecule drugs as they offer superior specificity, lower toxicity, and easy synthesis. Here we present an approach that leverages the dramatic performance increase afforded by the recent arrival of GPU accelerated thermodynamic integration (TI). GPU TI facilitates very fast, highly accurate binding affinity optimization of peptides against therapeutic targets. We benchmarked TI predictions using published peptide binding optimization studies. Prediction of mutations involving charged side-chains was found to be less accurate than for non-charged, and use of a more complex 3-step TI protocol was found to boost accuracy in these cases. Using the 3-step protocol for non-charged side-chains either had no effect or was detrimental. We use the benchmarked pipeline to optimize a peptide binding to our recently discovered cancer target: EME1. TI calculations predict beneficial mutations using both canonical and non-canonical amino acids. We validate these predictions using fluorescence polarization and confirm that binding affinity is increased. We further demonstrate that this increase translates to a significant reduction in pancreatic cancer cell viability.

Nonnative contact effects in protein folding

The folding simulations of three ββα-motifs and β-barrel structured proteins (NTL9, NuG2b, and CspA) were performed to determine the important roles of native and nonnative contacts in protein folding.

Deciphering evolution of immune recognition in antibodies

Background

Antibody, the primary effector molecule of the immune system, evolves after initial encounter with the antigen from a precursor form to a mature one to effectively deal with the antigen. Antibodies of a lineage diverge through antigen-directed isolated pathways of maturation to exhibit distinct recognition potential. In the context of evolution in immune recognition, diversity of antigen cannot be ignored. While there are reports on antibody lineage, structural perspective with respect to diverse recognition potential in a lineage has never been studied. Hence, it is crucial to evaluate how maturation leads to topological tailoring within a lineage enabling them to interact with significantly distinct antigens.

Results

A data-driven approach was undertaken for the study. Global experimental mouse and human antibody-antigen complex structures from PDB were compiled into a coherent database of germline-linked antibodies bound with distinct antigens. Structural analysis of all lineages showed variations in CDRs of both H and L chains. Observations of conformational adaptation made from analysis of static structures were further evaluated by characterizing dynamics of interaction in two lineages, mouse V_H1–84 and human V_H5–51. Sequence and structure analysis of the lineages explained that somatic mutations altered the geometries of individual antibodies with common structural constraints in some CDRs. Additionally, conformational landscape obtained from molecular dynamics simulations revealed that incoming pathogen led to further conformational divergence in the paratope (as observed across datasets) even while maintaining similar overall backbone topology. MM-GB/SA analysis showed binding energies to be in physiological range. Results of the study are coherent with experimental observations.

Conclusions

The findings of this study highlight basic structural principles shaping the molecular evolution of a lineage for significantly diverse antigens. Antibodies of a lineage follow different developmental pathways while preserving the imprint of the germline. From the study, it can be generalized that structural diversification of the paratope is an outcome of natural selection of a conformation from an available ensemble, which is further optimized for antigen interaction. The study establishes that starting from a common lineage, antibodies can mature to recognize a wide range of antigens. This hypothesis can be further tested and validated experimentally.

Electronic supplementary material

The online version of this article (10.1186/s12900-018-0096-1) contains supplementary material, which is available to authorized users.

UV-induced long-lived decays in solvated pyrimidine nucleosides resolved at the MS-CASPT2/MM level

The most relevant 'dark' electronic excited states in DNA/RNA pyrimidine nucleosides are mapped in water employing hybrid MS-CASPT2/MM optimisations with explicit solvation and including the sugar. Conical intersections (CIs) between initially accessed bright ¹ππ* and the lowest energy dark ¹nπ* excited states, involving the lone pair localised on the oxygen and/or nitrogen atoms are characterised. They are found in the vicinities of the Franck-Condon (FC) region and are shown to facilitate non-adiabatic population transfer. The excited state population of the ¹n_Oπ* state, localised in the carbonyl moiety on all pyrimidine nucleosides, is predicted to rapidly evolve to its minimum, displaying non-negligible potential energy barriers along its non-radiative decay, and accounting for the ps signal registered in pump-probe experiments as well as for an efficient population of the triplet state. Cytidine displays an additional ¹n_Nπ* state localised in the N3 atom and that leads to its excited state minimum displaying large potential energy barriers in the pathway connecting to the CI with the ground state. Sugar-to-base hydrogen/proton transfer processes are assessed in solution for the first time, displaying a sizable barrier along its decay and thus being competitive with other slow decay channels in the ps and ns timescales. A unified deactivation scheme for the long-lived channels of pyrimidine nucleosides is delivered, where the ¹n_Oπ* state is found to mediate the long-lived decay in the singlet manifold and act as the doorway for triplet population and thus accounting for the recorded phosphorescence and, more generally, for the transient/photoelectron spectral signals registered up to the ns timescale.

Assessing AMBER force fields for protein folding in an implicit solvent

Molecular dynamics (MD) simulation implemented with a state-of-the-art protein force field and implicit solvent model is an attractive approach to investigate protein folding, one of the most perplexing problems in molecular biology. But how well can force fields developed independently of implicit solvent models work together in reproducing diverse protein native structures and measuring the corresponding folding thermodynamics is not always clear. In this work, we performed enhanced sampling MD simulations to assess the ability of six AMBER force fields (FF99SBildn, FF99SBnmr, FF12SB, FF14ipq, FF14SB, and FF14SBonlysc) as coupled with a recently improved pair-wise GB-Neck2 model in modeling the folding of two helical and two β-sheet peptides. Whilst most of the tested force fields can yield roughly similar features for equilibrium conformational ensembles and detailed folding free-energy profiles for short α-helical TC10b in an implicit solvent, the measured counterparts are significantly discrepant in the cases of larger or β-structured peptides (HP35, 1E0Q, and GTT). Additionally, the calculated folding/unfolding thermodynamic quantities can only partially match the experimental data. Although a combination of the force fields and GB-Neck2 implicit model able to describe all aspects of the folding transitions towards the native structures of all the considered peptides was not identified, we found that FF14SBonlysc coupled with the GB-Neck2 model seems to be a reasonably balanced combination to predict peptide folding preferences.

Role of protein interactions in stabilizing canonical DNA features in simulations of DNA in crowded environments

Background

Cellular environments are highly crowded with biological macromolecules resulting in frequent non-specific interactions. While the effect of such crowding on protein structure and dynamics has been studied extensively, very little is known how cellular crowding affects the conformational sampling of nucleic acids.

Results

The effect of protein crowding on the conformational preferences of DNA (deoxyribonucleic acid) is described from fully atomistic molecular dynamics simulations of systems containing a DNA dodecamer surrounded by protein crowders. From the simulations, it was found that DNA structures prefer to stay in B-like conformations in the presence of the crowders. The preference for B-like conformations results from non-specific interactions of crowder proteins with the DNA sugar-phosphate backbone. Moreover, the simulations suggest that the crowder interactions narrow the conformational sampling to canonical regions of the conformational space.

Conclusions

The overall conclusion is that crowding effects may stabilize the canonical features of DNA that are most important for biological function. The results are complementary to a previous study of DNA in reduced dielectric environments where reduced dielectric environments alone led to a conformational shift towards A-DNA. Such a shift was not observed here suggested that the reduced dielectric response of cellular environments is counteracted by non-specific interactions with protein crowders under in vivo conditions.

Electronic supplementary material

The online version of this article (10.1186/s13628-018-0048-y) contains supplementary material, which is available to authorized users.

Residual Structure Accelerates Binding of Intrinsically Disordered ACTR by Promoting Efficient Folding upon Encounter

Intrinsically disordered proteins (IDPs) often fold into stable structures upon specific binding. The roles of residual structure of unbound IDPs in coupling binding and folding have been under much debate. While many studies emphasize the importance of conformational flexibility for IDP recognition, it was recently demonstrated that stabilization the N-terminal helix of intrinsically disordered ACTR accelerated its binding to another IDP, NCBD of the CREB-binding protein. To understand how enhancing ACTR helicity accelerates binding, we derived a series of topology-based coarse-grained models that mimicked various ACTR mutants with increasing helical contents and reproduced their NCBD binding affinities. Molecular dynamics simulations were then performed to sample hundreds of reversible coupled binding and folding transitions. The results show that increasing ACTR helicity does not alter the baseline mechanism of synergistic folding, which continues to follow "extended conformational selection" with multiple stages of selection and induced folding. Importantly, these coarse-grained models, while only calibrated based on binding thermodynamics, recapitulate the observed kinetic acceleration with increasing ACTR helicity. However, the residual helices do not enhance the association kinetics via more efficient seeding of productive collisions. Instead, they allow the nonspecific collision complexes to evolve more efficiently into the final bound and folded state, which is the primary source of accelerated association kinetics. Meanwhile, reduced dissociation kinetics with increasing ACTR helicity can be directly attributed to smaller entropic cost of forming the bound state. Altogether, this study provides important mechanistic insights into how residual structure may modulate thermodynamics and kinetics of IDP interactions.

Structural motif, topi and its role in protein function and fibrillation

A protein chain is arranged into regions in which the backbone is organized into regular patterns (of conformation and hydrogen bonding) to form the most common secondary structures, α-helix and β-sheet, which are interspersed by turns and more irregular loop regions. A structural motif, topi, is discussed in which a pair of 2-residue segments, each containing hydrogen-bonded five-membered fused-ring motifs, distant in sequence are linked to each other by a hydrogen bond. Though a small motif, it appears to be important in the context of local folding patterns of proteins and occurs near protein active sites. The motif shows quite significant residue preference, and a Cys (or Ser) occupying the second position may further stabilize the motif by forming an additional hydrogen bond across it. Remarkably, topi is found within disease causing misfolded proteins, such as the fibrilled form of Aβ42, and also across the interface between two protein chains. This motif may be an important component of fibrillation and useful for modeling loop regions.

Conformational Selection of a Polyproline Peptide by Ebola Virus VP30

An adaptive temperature-based replica-exchange simulation of a peptide extracted from the Ebola virus nucleoprotein containing a polyproline sequence motif is reported. The simulation results of applying the CHARMM36m force field with a generalized Born solvent model is presented. Conformational heterogeneity is described by potentials of mean force (PMFs) for a set of reaction coordinates that define the topological fold space. Starting from an extended backbone conformation of the peptide observed in an X-ray crystallographic assembly with the Ebola virus protein VP30, the PMFs report a conformational landscape populated by chain excursions to collapsed states with limited transitions to either an extended fold or a canonical polyproline type II helix. Clustering of the conformations and applying an elastic network interpolation model yield a multistep pathway of conformational selection that minimizes the net transition-state cost from the population hub to the bound state. Related difference between the pathway endpoints taken from the PMFs reveal a significant free-energy penalty in reaching a population shift. To evaluate sequence fitness of the Ebola virus peptide in generating probability distributions, two human sequence variants are modeled and are found to produce profiles that show extensive deviations, thus suggesting either dissimilar binding mechanisms or the lack of recognition by VP30.

Molecular dynamics provides insight into how N251A and N251Y mutations in the active site of Bacillus licheniformis RN-01 levansucrase disrupt production of long-chain levan

Produced by levansucrase, levan and levan oligosaccharides (GF_n) have potential applications in food and pharmaceutical industries such as prebiotics, anti-tumor and anti-inflammatory agents. Previous study reported that Bacillus licheniformis RN-01 levansucrase could produce levan oligosaccharides and long-chain levan. However, its N251A and N251Y mutants could effectively produce short-chain oligosaccharides upto GF_3, but they could not produce long-chain levan. We hypothesized that these mutations probably reduced GF₃ binding affinity in levansucrase active site that contains fructosyl-Asp93 intermediate and caused GF₃ to be in an unfavorable orientation for transfructosylation; therefore, levansucrase could not effectively extend GF₃ by one fructosyl residue to produce GF₄ and subsequently long-chain levan. However, these mutations probably did not significantly reduce binding affinity or drastically change orientation of GF₂; therefore, levansucrase could still extend GF₂ to produce GF₃. Using this hypothesis, we employed molecular dynamics to investigate effects of these mutations on GF₂/GF₃ binding in levansucrase active site. Our results reasonably support this hypothesis as N251A and N251Y mutations did not significantly reduce GF₂ binding affinity, as calculated by MM-GBSA technique and hydrogen bond occupations, or drastically change orientation of GF₂ in levansucrase active site, as measured by distance between atoms necessary for transfructosylation. However, these mutations drastically decreased GF₃ binding affinity and caused GF₃ to be in an unfavorable orientation for transfructosylation. Furthermore, the free energy decomposition and hydrogen bond occupation results suggest the importance of Arg255 in GF₂/GF₃ binding in levansucrase active site. This study provides important and novel insight into the effects of N251A and N251Y mutations on GF₂/GF₃ binding in levansucrase active site and how they may disrupt production of long-chain levan. This knowledge could be beneficial in designing levansucrase to efficiently produce levan oligosaccharides with desired length.

Secondary structures transition of tau protein with intrinsically disordered proteins specific force field

Microtubule-associated Tau protein plays a key role in assembling microtubule and modulating the functional organization of the neuron and developing axonal morphology, growth, and polarity. The pathological Tau can aggregate into cross-beta amyloid as one of the hallmarks for Alzheimer's disease (AD). Therefore, one of the top priorities in AD research is to figure out the structural model of Tau aggregation and to screen the inhibitors. The latest generation intrinsically disordered protein specific force field ff14IDPSFF significantly improved the distributions of heterogeneous conformations for intrinsically disordered proteins (IDPs). Here, the molecular dynamics (MD) simulations with three force fields of ff14SB, ff14IDPs, and ff14IDPSFF were employed to investigate the secondary structures transition of Tau (267-312) fragment. The results indicate that ff14IDPSFF can generate more heterogeneous conformers, and the predicted secondary structural distribution is closer to that of the experimental observation. In addition, predicted secondary chemical shifts from ff14IDPSFF are the most approach to those of experiment. Secondary structures transition kinetics for Tau(267-312) with ff14IDPSFF shows that the secondary structures were gradually transformed from α-helix to β-strand and the β-strand located at the regions of the residues 274-280 and residues 305-311. Besides, the driving force for the secondary structures transition of Tau(267-312) is mainly hydrophobic interactions which located at hexa-peptides ²⁷⁵ VQIINK²⁸⁰ and ³⁰⁶ VQIVYK³¹¹ . Secondary structure transition of Tau protein can give insight into the aggregation mechanism for AD.

Computational Insight Into Vitamin K1 ω-Hydroxylation by Cytochrome P450 4F2

Vitamin K₁ (VK1) plays an important role in the modulation of bleeding disorders. It has been reported that ω-hydroxylation on the VK1 aliphatic chain is catalyzed by cytochrome P450 4F2 (CYP4F2), an enzyme responsible for the metabolism of eicosanoids. However, the mechanism of VK1 ω-hydroxylation by CYP4F2 has not been disclosed. In this study, we employed a combination of quantum mechanism (QM) calculations, homology modeling, molecular docking, molecular dynamics (MD) simulations, and combined quantum mechanism/molecular mechanism (QM/MM) calculations to investigate the metabolism profile of VK1 ω-hydroxylation. QM calculations based on the truncated VK1 model show that the energy barrier for ω-hydroxylation is about 6-25 kJ/mol higher than those at other potential sites of metabolism. However, results from the MD simulations indicate that hydroxylation at the ω-site is more favorable than at the other potential sites, which is in accordance with the experimental observation. The evaluation of MD simulations was further endorsed by the QM/MM calculation results. Our studies thus suggest that the active site residues of CYP4F2 play a determinant role in the ω-hydroxylation. Our results provide structural insights into the mechanism of VK1 ω-hydroxylation by CYP4F2 at the atomistic level and are helpful not only for characterizing the CYP4F2 functions but also for looking into the ω-hydroxylation mediated by other CYP4 enzymes.

Discovery of small molecule inhibitors of adenovirus by disrupting E3-19K/HLA-A2 interactions

The binding of the adenovirus (Ad) protein E3-19K with the human leukocyte antigen (HLA) plays an important role in Ad infections, which is the causative agent of a series of gastrointestinal, respiratory and ocular diseases. The objective of this research is to evaluate the essential interactions between E3-19K and HLA-A2 using the X-ray crystal structure of the E3-19K/HLA-A2 complex, and to identify small molecules that could potentially disrupt their binding. Computational methods, including molecular dynamic simulations, MM/GBSA calculations, and computational solvent mapping, were implemented to determine potential binding site(s) for small molecules. The previous experimentally determined hot spot residues, Q54 and E177 in HLA-A2, were also predicted to be the dominant residues for binding to E3-19K by our theoretical calculations. Several other residues were also found to play pivotal roles for the binding of E3-19K with HLA-A2. Residues adjacent to E177, including Q54 and several other residues theoretically predicted to be crucial in HLA-A2 were selected as a potential binding pocket to perform virtual screening with 1200 compounds from the Prestwick library. Seven hits were validated by surface plasmon resonance (SPR) as binders to HLA-A2 as a first step in identifying molecules that can perturb its association with the Ad E3-19K protein.

Atomistic Peptide Folding Simulations Reveal Interplay of Entropy and Long-Range Interactions in Folding Cooperativity

Understanding how proteins fold has remained a problem of great interest in biophysical research. Atomistic computer simulations using physics-based force fields can provide important insights on the interplay of different interactions and energetics and their roles in governing the folding thermodynamics and mechanism. In particular, generalized Born (GB)-based implicit solvent force fields can be optimized to provide an appropriate balance between solvation and intramolecular interactions and successfully recapitulate experimental conformational equilibria for a set of helical and β-hairpin peptides. Here, we further demonstrate that key thermodynamic properties and their temperature dependence obtained from replica exchange molecular dynamics simulations of these peptides are in quantitative agreement with experimental results. Useful lessons can be learned on how the interplay of entropy and sequentially long-range interactions governs the mechanism and cooperativity of folding. These results highlight the great potential of high-quality implicit solvent force fields for studying protein folding and large-scale conformational transitions.

The first crystal structure of a DNA-free nuclear receptor DNA binding domain sheds light on DNA-driven allostery in the glucocorticoid receptor

The glucocorticoid receptor (GR) is a steroid hormone receptor of the nuclear receptor family that regulates gene expression in response to glucocorticoid hormone signaling. Interaction with specific GR DNA binding sequences causes conformational changes in the GR DNA binding domain (DBD) that result in recruitment of specific sets of co-regulators that determine transcriptional outcomes. We have solved the crystal structure of GR DBD in its DNA-free state, the first such crystal structure from any nuclear receptor. In contrast to previous NMR structures, this crystal structure reveals that free GR DBD adopts a conformation very similar to DNA-bound states. The lever arm region is the most variable element in the free GR DBD. Molecular dynamics of the free GR DBD as well as GR DBD bound to activating and repressive DNA elements confirm lever arm flexibility in all functional states. Cluster analysis of lever arm conformations during simulations shows that DNA binding and dimerization cause a reduction in the number of conformations sampled by the lever arm. These results reveal that DNA binding and dimerization drive conformational selection in the GR DBD lever arm region and show how DNA allosterically controls GR structure and dynamics.

Intramolecular Diffusion in α-Synuclein: It Depends on How You Measure It

Intramolecular protein diffusion, the motion of one part of the polypeptide chain relative to another part, is a fundamental aspect of protein folding and may modulate amyloidogenesis of disease-associated intrinsically disordered proteins. Much work has determined such diffusion coefficients using a variety of probes, but there has been an apparent discrepancy between measurements using long-range probes, such as fluorescence resonance energy transfer, and short-range probes, such as Trp-Cys quenching. In this work, we make both such measurements on the same protein, α-synuclein, and confirm that such discrepancy exists. Molecular dynamics simulations suggest that such differences result from a diffusion coefficient that depends on the spatial distance between probes. Diffusional estimates in good quantitative agreement with experiment are obtained by accounting for the distinct distance ranges probed by fluorescence resonance energy transfer and Trp-Cys quenching.

Universal effects of solvent species on the stabilized structure of a protein

We investigate the effects of solvent specificities on the stability of the native structure (NS) of a protein on the basis of our free-energy function (FEF). We use CPB-bromodomain (CBP-BD) and apoplastocyanin (apoPC) as representatives of the protein universe and water, methanol, ethanol, and cyclohexane as solvents. The NSs of CBP-BD and apoPC consist of 66% α-helices and of 35% β-sheets and 4% α-helices, respectively. In order to assess the structural stability of a given protein immersed in each solvent, we contrast the FEF of its NS against that of a number of artificially created, misfolded decoys possessing the same amino-acid sequence but significantly different topology and α-helix and β-sheet contents. In the FEF, we compute the solvation entropy using the morphometric approach combined with the integral equation theories, and the change in electrostatic (ES) energy upon the folding is obtained by an explicit atomistic but simplified calculation. The ES energy change is represented by the break of protein-solvent hydrogen bonds (HBs), formation of protein intramolecular HBs, and recovery of solvent-solvent HBs. Protein-solvent and solvent-solvent HBs are absent in cyclohexane. We are thus able to separately evaluate the contributions to the structural stability from the entropic and energetic components. We find that for both CBP-BD and apoPC, the energetic component dominates in methanol, ethanol, and cyclohexane, with the most stable structures in these solvents sharing the same characteristics described as an association of α-helices. In particular, those in the two alcohols are identical. In water, the entropic component is as strong as or even stronger than the energetic one, with a large gain of translational, configurational entropy of water becoming crucially important so that the relative contents of α-helix and β-sheet and the content of total secondary structures are carefully selected to achieve sufficiently close packing of side chains. If the energetic component is excluded for a protein in water, the priority is given to closest side-chain packing, giving rise to the formation of a structure with very low α-helix and β-sheet contents. Our analysis, which requires minimal computational effort, can be applied to any protein immersed in any solvent and provides robust predictions that are quite consistent with the experimental observations for proteins in different solvent environments, thus paving the way toward a more detailed understanding of the folding process.

Nonproductive Binding Modes as a Prominent Feature of Aβ40 Fiber Elongation: Insights from Molecular Dynamics Simulation

The formation of amyloid fibers has been implicated in a number of neurodegenerative diseases. The growth of amyloid fibers is strongly thermodynamically favorable, but kinetic traps exist where the incoming monomer binds in an incompatible conformation that blocks further elongation. Unfortunately, this process is difficult to follow experimentally at the atomic level. It is also too complex to simulate in full detail and to date has been explored either through coarse-grained simulations, which may miss many important interactions, or full atomic simulations, in which the incoming peptide is constrained to be near the ideal fiber geometry. Here we use an alternate approach starting from a docked complex in which the monomer is from an experimental NMR structure of one of the major conformations in the unbound ensemble, a largely unstructured peptide with the central hydrophobic region in a 3₁₀ helix. A 1000 ns full atomic simulation in explicit solvent shows the formation of a metastable intermediate by sequential, concerted movements of both the fiber and the monomer. A Markov state model shows that the unfolded monomer is trapped at the end of the fiber in a set of interconverting antiparallel β-hairpin conformations. The simulation here may serve as a model for the binding of other non-β-sheet conformations to amyloid fibers.

Comprehensive Peptide Ion Structure Studies Using Ion Mobility Techniques: Part 3. Relating Solution-Phase to Gas-Phase Structures

Molecular dynamics (MD) simulations have been utilized to study peptide ion conformer establishment during the electrospray process. An explicit water model is used for nanodroplets containing a model peptide and hydronium ions. Simulations are conducted at 300 K for two different peptide ion charge configurations and for droplets containing varying numbers of hydronium ions. For all conditions, modeling has been performed until production of the gas-phase ions and the resultant conformers have been compared to proposed gas-phase structures. The latter species were obtained from previous studies in which in silico candidate structures were filtered according to ion mobility and hydrogen-deuterium exchange (HDX) reactivity matches. Results from the present study present three key findings namely (1) the evidence from ion production modeling supports previous structure refinement studies based on mobility and HDX reactivity matching, (2) the modeling of the electrospray process is significantly improved by utilizing initial droplets existing below but close to the calculated Rayleigh limit, and (3) peptide ions in the nanodroplets sample significantly different conformers than those in the bulk solution due to altered physicochemical properties of the solvent. Graphical Abstract ᅟ.

Theoretical Model of the Protochlorophyllide Oxidoreductase from a Hierarchy of Protocols

The enzyme protochlorophyllide oxidoreductase (LPOR) catalyzes the light-driven reduction of protochlorophyllide (Pchlide), a crucial step in chlorophyll biosynthesis. Molecular understanding of the photocatalytic mechanism of LPOR is essential for harnessing light energy to mediate enzymatic reactions. The absence of X-ray crystal structure has promoted the development of LPOR homology models that lack a catalytically competent active site and could not explain the variously reported spectroscopic evidence, including time-resolved optical spectroscopy data. We have refined previous structural models to account for the catalytic active site and the characteristic experimental spectral features of Pchlide binding, including the 26 cm^-1 red shift of the C₁₃₍₁₎ carbonyl stretch vibration in the mid-infrared (IR) and the 12 nm red shift of the Q _x electronic band. A hierarchy of theoretical methods, including homology modeling, molecular dynamics simulations, hybrid quantum mechanics [(TD-)DFT]/molecular mechanics [AMBER] calculations, and computational vibrational and electronic spectroscopies, have been combined in an iterative protocol to reproduce experimental evidence and to predict ultrafast transient IR spectroscopic fingerprints associated with the catalytic process. The successful application to the LPOR enzyme indicates that the presented hierarchical protocol provides a general workflow to protein structure refinement.

The effects of implicit modeling of nonpolar solvation on protein folding simulations

Implicit solvent models, in which the polar and nonpolar solvation free-energies of solute molecules are treated separately, have been widely adopted for molecular dynamics simulation of protein folding. While the development of the implicit models is mainly focused on the methodological improvement and key parameter optimization for polar solvation, nonpolar solvation has been either ignored or described by a simplistic surface area (SA) model. In this work, we performed the folding simulations of multiple β-hairpin and α-helical proteins with varied surface tension coefficients embedded in the SA model to clearly demonstrate the effects of nonpolar solvation treated by a popular SA model on protein folding. The results indicate that the change in the surface tension coefficient does not alter the ability of implicit solvent simulations to reproduce a protein native structure but indeed controls the components of the equilibrium conformational ensemble and modifies the energetic characterization of the folding transition pathway. The suitably set surface tension coefficient can yield explicit solvent simulations and/or experimentally suggested folding mechanism of protein. In addition, the implicit treatment of both polar and nonpolar components of solvation free-energy contributes to the overestimation of the secondary structure in implicit solvent simulations.

PD-L1 Nanobody Competitively Inhibits the Formation of the PD-1/PD-L1 Complex: Comparative Molecular Dynamics Simulations

The anti-PD-L1 monoclonal antibody (mAb) targeting PD-1/PD-L1 immune checkpoint has achieved outstanding results in clinical application and has become one of the most popular anti-cancer drugs. The mechanism of molecular recognition and inhibition of PD-L1 mAbs is not yet clear, which hinders the subsequent antibody design and modification. In this work, the trajectories of PD-1/PD-L1 and nanobody/PD-L1 complexes were obtained via comparative molecular dynamics simulations. Then, a series of physicochemical parameters including hydrogen bond, dihedral angle distribution, pKa value and binding free energy, and so forth, were all comparatively analyzed to investigate the recognition difference between PD-L1 and PD-1 and nanobody. Both LR113 (the amino acid residues in PD-L1 are represented by the lower left sign of L) and LR125 residues of PD-L1 undergo significant conformational change after association with mAbs, which dominates a strong electrostatic interaction. Solvation effect analysis revealed that solvent-water enhanced molecular recognition between PD-L1 and nanobody. By combining the analyses of the time-dependent root mean squared fluctuation (RMSF), free energy landscape, clustering and energy decomposition, the potential inhibition mechanism was proposed that the nanobody competitively and specifically bound to the β-sheet groups of PD-L1, reduced the PD-L1’s flexibility and finally blocked the formation of PD-1/PD-L1 complex. Based on the simulation results, site-directed mutagenesis of ND99 (the amino acid residues in Nano are displayed by the lower left sign of N) and NQ116 in the nanobody may be beneficial for improving antibody activity. This work offers some structural guidance for the design and modification of anticancer mAbs based on the structure of the PD-1/PD-L1 complex.

Order-disorder transition of intrinsically disordered kinase inducible transactivation domain of CREB

Transcription factor cyclic Adenosine monophosphate response-element binding protein plays a critical role in the cyclic AMP response pathway via its intrinsically disordered kinase inducible transactivation domain (KID). KID is one of the most studied intrinsically disordered proteins (IDPs), although most previous studies focus on characterizing its disordered state structures. An interesting question that remains to be answered is how the order-disorder transition occurs at experimental conditions. Thanks to the newly developed IDP-specific force field ff14IDPSFF, the quality of conformer sampling for IDPs has been dramatically improved. In this study, molecular dynamics (MD) simulations were used to study the order-to-disorder transition kinetics of KID based on the good agreement with the experiment on its disordered-state properties. Specifically, we tested four force fields, ff99SBildn, ff99IDPs, ff14IDPSFF, and ff14IDPs in the simulations of KID and found that ff14IDPSFF can generate more diversified disordered conformers and also reproduce more accurate experimental secondary chemical shifts. Kinetics analysis of MD simulations demonstrates that the order-disorder transition of KID obeys the first-order kinetics, and the transition nucleus is I127/L128/L141. The possible transition pathways from the nucleus to the last folded residues were identified as I127-R125-L138-L141-S143-A145 and L128-R125-L138-L141-S143-A145 based on a residue-level dynamical network analysis. These computational studies not only provide testable prediction/hypothesis on the order-disorder transition of KID but also confirm that the ff14IDPSFF force field can be used to explore the correlation between the structure and function of IDPs.

UNRES server for physics-based coarse-grained simulations and prediction of protein structure, dynamics and thermodynamics

A server implementation of the UNRES package (http://www.unres.pl) for coarse-grained simulations of protein structures with the physics-based UNRES model, coined a name UNRES server, is presented. In contrast to most of the protein coarse-grained models, owing to its physics-based origin, the UNRES force field can be used in simulations, including those aimed at protein-structure prediction, without ancillary information from structural databases; however, the implementation includes the possibility of using restraints. Local energy minimization, canonical molecular dynamics simulations, replica exchange and multiplexed replica exchange molecular dynamics simulations can be run with the current UNRES server; the latter are suitable for protein-structure prediction. The user-supplied input includes protein sequence and, optionally, restraints from secondary-structure prediction or small x-ray scattering data, and simulation type and parameters which are selected or typed in. Oligomeric proteins, as well as those containing D-amino-acid residues and disulfide links can be treated. The output is displayed graphically (minimized structures, trajectories, final models, analysis of trajectory/ensembles); however, all output files can be downloaded by the user. The UNRES server can be freely accessed at http://unres-server.chem.ug.edu.pl.

Structure refinement of membrane proteins via molecular dynamics simulations

A refinement protocol based on physics-based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge-based or implicit membrane-based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid-facing residues. Scoring with knowledge-based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane-based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.