Predicting the quaternary structure of protein complex is an important problem. Inter-chain residue-residue contact prediction can provide useful information to guide the ab initio reconstruction of quaternary structures. However, few methods have been developed to build quaternary structures from predicted inter-chain contacts. Here, we introduce a gradient descent optimization algorithm (GD) to build quaternary structures of protein dimers utilizing inter-chain contacts as distance restraints. We evaluate GD on several datasets of homodimers and heterodimers using true or predicted contacts. GD consistently performs better than a simulated annealing method and a Markov Chain Monte Carlo simulation method. Using true inter-chain contacts as input, GD can reconstruct high-quality structural models for homodimers and heterodimers with average TM-score ranging from 0.92 to 0.99 and average interface root mean square distance (I-RMSD) from 0.72 Å to 1.64 Å. On a dataset of 115 homodimers, using predicted inter-chain contacts as input, the average TM-score of the structural models built by GD is 0.76. For 46% of the homodimers, high-quality structural models with TM-score >= 0.9 are reconstructed from predicted contacts. There is a strong correlation between the quality of the reconstructed models and the precision and recall of predicted contacts. If the precision or recall of predicted contacts is >20%, GD can reconstruct good models for most homodimers, indicating only a moderate precision or recall of inter-chain contact prediction is needed to build good structural models for most homodimers. Moreover, the accuracy of reconstructed models positively correlates with the contact density in dimers.

Lignin is one of the world’s most abundant organic polymers, and 2-pyrone-4,6-dicarboxylate lactonase (LigI) catalyzes the hydrolysis of 2-pyrone-4,6-dicarboxylate (PDC) in the degradation of lignin. The pH has profound effects on enzyme catalysis and therefore we studied this in the context of LigI. We found that changes of the pH mostly affects surface residues, while the residues at the active site are more subject to changes of the surrounding microenvironment. In accordance with this, a high pH facilitates the deprotonation of the substrate. Detailed free energy calculations by the empirical valence bond (EVB) approach revealed that the overall hydrolysis reaction is more likely when the three active site histidines (His31, His33 and His180) are protonated at the ɛ site, however, protonation at the δ site may be favored during specific steps of reaction. Our studies have uncovered the determinant role of the protonation state of the active site residues His31, His33 and His180 in the hydrolysis of PDC.

We have investigated the pressure- and temperature-induced conformational changes associated with the low complexity domain of hnRNP A1, an RNA-binding protein able to phase separate in response to cellular stress. Solution NMR spectra of the hnRNP A1 low-complexity domain fused with protein-G B1 domain were collected from 1 to 2,500 bar and from 268 K to 290 K. While the GB1 domain shows the typical pressure-induced and cold temperature-induced unfolding expected for small globular domains, the low-complexity domain of hnRNP A1 exhibits unusual pressure and temperature dependences. We observed that the low-complexity domain is pressure sensitive, undergoing a major conformational transition within the prescribed pressure range. Remarkably, this transition has the inverse temperature dependence of a typical folding-unfolding transition. Our results suggest the presence of a low-lying extended, and fully solvated state(s) of the low-complexity domain that may play a role in phase separation. This study highlights the exquisite sensitivity of solution NMR spectroscopy to observe subtle conformational changes and illustrates how pressure perturbation can be used to determine the properties of metastable conformational ensembles.

Protein structure networks (PSNs) have long been used to provide a coarse yet meaningful representation of protein structure, dynamics, and internal communication pathways. An important question is what criteria should be applied to construct the network so that to include relevant interresidue contacts while avoiding unnecessary connections. To address this issue we systematically considered varying residue distance cutoff length and the probability threshold for contact formation to construct PSNs based on atomistic molecular dynamics in order to assess the amount of mutual information within the resulting representations. We found that the minimum in mutual information is universally achieved at the cutoff length of 5 Å, irrespective of the applied contact formation probability threshold in all considered, distinct proteins. Assuming that the optimal PSNs should be characterised by the least amount of redundancy, which corresponds to the minimum in mutual information, this finding suggests an objective criterion for cutoff distance and supports the existing preference towards its customary selection around 5 Å length, typically based to date on heuristic criteria.

The denatured state of several proteins has been shown to display transient structures that are relevant for folding, stability and aggregation. To detect them by nuclear magnetic resonance (NMR) spectroscopy, the denatured state must be stabilized by chemical agents or changes in temperature. This makes the environment different from that experienced in biologically relevant processes. Using high-resolution heteronuclear NMR spectroscopy, we have characterized several denatured states of a monomeric variant of HIV-1 protease induced by different concentrations of urea, guanidinium chloride and acetic acid. We have extrapolated the chemical shifts and the relaxation parameters to the denaturant-free denatured state at native conditions, showing that they converge to the same values. Subsequently, we characterized the conformational properties of this biologically relevant denatured state under native conditions by advanced molecular dynamics simulations and validated the results by comparison to experimental data. We show that the denatured state of HIV-1 protease under native conditions displays rich patterns of transient native and non-native structures, which could be of relevance to its guidance through a complex folding process.

Multi-domain proteins are not only formed through natural evolution but can also be generated by recombinant DNA technology. Because many fusion proteins can enhance the selectivity of cell targeting, these artificially produced molecules, called multi-specific biologics, are promising drug candidates, especially for immunotherapy. Moreover, the rational design of domain linkers in fusion proteins is becoming an essential step toward a quantitative understanding of the dynamics in these biopharmaceutics. We developed a computational framework to characterize the impacts of peptide linkers on the dynamics of multi-specific biologics. We constructed a benchmark containing six types of linkers that represent various lengths and degrees of flexibility and used them to connect two natural proteins as a test system. The microsecond dynamics of these proteins generated from Anton were projected onto a coarse-grained conformational space. The similarity of dynamics among different proteins in this low-dimensional space was further analyzed by a neural network model. Finally, hierarchical clustering was applied to place linkers into different subgroups based on the neural network classification results. The clustering results suggest that the length of linkers used to spatially separate different functional modules plays the most important role in regulating the dynamics of this fusion protein. Given the same number of amino acids, linker flexibility functions as a regulator of protein dynamics. In summary, we illustrated that a new computational strategy can be used to study the dynamics of multi-domain fusion proteins by a combination of long timescale molecular dynamics simulation, coarse-grained modeling, and artificial intelligence.

One way in which trichocyte keratin intermediate filament proteins (keratins) and keratin associated proteins (KAPs) differ from their epithelial equivalents is in their higher levels of cysteine residues. Interactions between these cysteine residues within a mammalian fiber, and the putative regular organization of interactions (i.e., types of disulfide bond) are likely important for defining fiber mechanical properties, and thus biological functionality of hairs. Here we extend a previous study of cysteine accessibility under different levels of exposure to reducing compounds to explore a finer set of levels associated with interactions between keratins and KAPs. We found that most of the cysteines in the KAPs were close to either the N- or C- terminal domains of these proteins. The most accessible cysteines in keratins were present in the head or tail domains indicating their function in readily forming intermolecular bonds with KAPs. Some of the more buried cysteines in keratins were discovered either close to or within the rod region in positions previously identified in human epithelial keratins as being involved in crosslinking between the heterodimers of the tetramer. Our present study therefore provides a deeper understanding of the accessibility of disulfides especially in keratins and thus proves that there is some specificity to the disulfide bond interactions leading to these intermolecular bonds stabilizing the fiber structure.

Acetylcholinesterase (AChE) is the crucial enzyme in the central nervous system. It is the target of various organophosphorus nerve agents and pesticides, and the inhibition of AChE is a therapeutic strategy for the treatment of various neurological-related diseases. The Glu202 is a key residue adjacent to the catalytic His447 and plays important role in catalysis. Although the Glu202 has long been considered as negatively charged in many studies, more and more evidences support a protonated Glu202. However, Glu202 is freely accessible by solvent, and thus it seems more reasonable for Glu202 to majorly take the deprotonated state. In the present work, we carried out a series of molecular dynamics simulations with the Glu202 adopting different protonation states. Our results show that the protonated Glu202 is important in maintaining the key hydrogen bond network that supports the catalytic triad, whereas the deprotonated Glu202 results in the collapse of the key hydrogen bond network which consequently destabilizes the catalytic His447. We also notice that different protonation states of Glu202 merely alters the binding mode of ACh. However, since the catalytic His447 is disrupted if Glu202 is deprotonated, His447 can not facilitate the nucleophilic attack performed by Ser203. Therefore, the catalytic efficiency of ACh hydrolysis should be remarkably decreased if Glu202 is deprotonated. Our findings suggest that, when designing and developing highly active AChE inhibitors or proposing mechanistic hypotheses for AChE-catalyzed reactions, the protonated state of Glu202 should be considered.

Multimeric protein complexes are molecular apparatuses to regulate biological systems and often determine their fate. Among proteins forming such molecular assemblies, amyloid proteins have drawn attention over a half-century since amyloid fibril formation of these proteins is supposed to be a common pathogenic cause for neurodegenerative diseases. This process is triggered by the accumulation of fibril-like aggregates, while the microscopic mechanisms are mostly elusive due to technical limitation of experimental methodologies in individually observing each of diverse aggregate species in the aqueous solution. We then addressed this problem by employing atomistic molecular dynamics simulations for the paradigmatic amyloid protein, amyloid-β (1-42) (Aβ 42). Seven different dimeric forms of oligomeric Aβ 42 fibril-like aggregate in aqueous solution, ranging from tetramer to decamer, were considered. We found additive effects of the size of these fibril-like aggregates on their thermodynamic stability and have clarified kinetic suppression of protomer-protomer dissociation reactions at and beyond the point of pentamer dimer formation. This observation was obtained from the specific combination of the Aβ 42 protomer structure and the physicochemical condition that we here examined, while it is worthwhile to recall that several amyloid fibrils take dimeric forms of their protomers. We could thus conclude that the stable formation of fibril-like protomer dimer should be involved in a turning point where rapid growth of amyloid fibrils is triggered.

Scoring docking solutions is a difficult task, and many methods have been developed for this purpose. In docking, only a handful of the hundreds of thousands of models generated by docking algorithms are acceptable, causing difficulties when developing scoring functions. Today’s best scoring functions can significantly increase the number of top-ranked models but still fails for most targets. Here, we examine the possibility of utilising predicted residues on a protein-protein interface to score docking models generated during the scan stage of a docking algorithm. Many methods have been developed to infer the portions of a protein surface that interact with another protein, but most have not been benchmarked using docking algorithms. Different interface prediction methods are systematically tested for scoring >300.000 low-resolution rigid-body template free docking decoys. Overall we find that BIPSPI is the best method to identify interface amino acids and score docking solutions. Further, using BIPSPI provides better docking results than state of the art scoring functions, with >12% of first ranked docking models being acceptable. Additional experiments indicated precision as a high-importance metric when estimating interface prediction quality, focusing on docking constraints production. We also discussed several limitations for the adoption of interface predictions as constraints in a docking protocol.

Chaperonin Hsp60, as a protein found in all organisms, is of great interest in medicine, since it is present in many tissues and can be used both as a drug and as an object of targeted therapy. Hence, Hsp60 deserves a fundamental comparative analysis to assess its evolutionary characteristics. It was found that the percent identity of Hsp60 amino acid sequences both within and between phyla was not high enough to identify Hsp60s as highly conserved proteins. In turn, their amino acid composition remained relatively constant. At the same time, the analysis of the nucleotide sequences showed that GC content in the Hsp60 genes was comparable to or greater than the genomic values, which may indicate a high resistance to mutations due to tight control of the nucleotide composition by DNA repair systems. Natural selection plays a dominant role in the evolution of Hsp60 genes. The degree of mutational pressure affecting the Hsp60 genes is quite low, and its direction does not depend on taxonomy. Interestingly, for the Hsp60 genes from Chordata, Arthropoda, and Proteobacteria the exact direction of mutational pressure could not be determined. However, upon further division into classes, it was found that the direction of the mutational pressure for Hsp60 genes from Fish differs from that for other chordates. The direction of the mutational pressure affects the synonymous codon usage bias. The number of high and low represented codons increases with increasing GC content, which can improve codon usage.

SR/RS dipeptide repeats vary in both length and position and are phosphorylated by SR protein kinases (SRPKs). PIM-1L, the long isoform of PIM-1 kinase, the splicing of which has been implicated in acute myeloid leukemia, contains a domain which consists largely of repeating SR/RS and SH/HS dipeptides (SR/SH-rich). In order to extend our knowledge on the specificity and cellular functions of SRPK1, here we investigate whether PIM-1L could act as substrate of SRPK1 by a combination of biochemical and computational approaches. Our biochemical data showed that the SR/SH-rich domain of PIM-1L was able to associate with SRPK1, yet it could not act as a substrate but, instead, inactivated the kinase. In line with our biochemical data, molecular modeling followed by a microsecond-scale all-atom molecular dynamics (MD) simulation suggests that the SR/SH-rich domain acts as a pseudo-docking peptide that binds to the same acidic docking-groove used in other SRPK1 interactions and induces inactive SRPK1 conformations. Comparative community network analysis of the MD trajectories, unraveled the dynamic architecture of apo SRPK1 and notable alterations of allosteric communications upon PIM-1L peptide binding. This analysis also allowed us to identify key SRPK1 residues, including unique ones, with a pivotal role in mediating allosteric signal propagation within the kinase core. Interestingly, most of the identified amino acids correspond to cancer-associated amino acid changes, validating our results. In total, this work provides insights not only on the details of SRPK1 inhibition by the PIM-1L SR/SH-domain, but also contributes to an in-depth understanding of SRPK1 regulation.

NADPH:protochlorophyllide (Pchlide) oxidoreductase (POR) is a key enzyme of chlorophyll biosynthesis in angiosperms. It is one of few known photoenzymes, which catalyzes the light-activated trans-reduction of the C17-C18 double bond of Pchlide’s porphyrin ring. Due to the light requirement, dark-grown angiosperms cannot synthesize chlorophyll. No crystal structure of POR is available, so to improve understanding of the protein’s three-dimensional structure, its dimerization, and binding of ligands (both the cofactor NADPH and substrate Pchlide), we computationally investigated the sequence and structural relationships among homologous proteins identified through database searches. The results indicate that α4 and α7 helices of monomers form the interface of POR dimers. On the basis of conserved residues, we predicted 11 functionally important amino acids that play important roles in POR binding to NADPH. Structural comparison of available crystal structures revealed that they participate in formation of binding pockets that accommodate the Pchlide ligand, and that five atoms of the closed tetrapyrrole are involved in non-bonding interactions. However, we detected no clear pattern in the physico-chemical characteristics of the amino acids they interact with. Thus, we hypothesize that interactions of these atoms in the Pchlide porphyrin ring are important to hold the ligand within the POR binding site. Analysis of Pchlide binding in POR by molecular docking and PELE simulations revealed that the orientation of the nicotinamide group is important for Pchlide binding. These findings highlight the complexity of interactions of porphyrin-containing ligands with proteins, and we suggest that fit-inducing processes play important roles in POR-Pchlide interactions.

Amyloid beta (Aβ) peptides, a major contributor to Alzheimers disease, occur in differing lengths, each of which forms a multitude of assembly types. The most toxic soluble oligomers are formed by Aβ42; some of which have antiparallel β-sheets. Previously, our group proposed molecular models of Aβ42 hexamers in which the C-terminus third of the peptide (S3) forms an antiparallel 6-stranded β-barrel that is surrounded by an antiparallel barrel formed by the more polar N-terminus (S1) and middle (S2) portions. These hexamers were proposed to act as seeds from which dodecamers, octadecamers, both smooth and beaded annular protofibrils, and transmembrane channels form. Since then, numerous aspects of our models have been supported by experimental findings. Recently, NMR-based structures have been proposed for Aβ42 tetramers and octamers, and NMR studies have been reported for oligomers composed of ~ 32 monomers. Here we propose a range of concentric β-barrel models and compare their dimensions to image-averaged electron micrographs of both beaded annular protofibrils (bAPFs) and smooth annular protofibrils (sAPFs) of Aβ42. The smaller oligomers have 6, 8, 12, 16, and 18 monomers. These beads string together to form necklace-like bAPFs. These gradually morph into sAPFs in which a S3 β-barrel is shielded on one or both sides by β-barrels formed from S1 and S2 segments.

Amyloid beta (Aβ of Alzheimer’s disease) and α-synuclein (α-Syn of Parkinson’s disease) form large fibrils. Evidence is increasing however that much smaller oligomers are more toxic and that these oligomers can form transmembrane ion channels. We have proposed previously that Aβ42 oligomers, annular protofibrils, and ion channels adopt concentric β-barrel molecular structures. Here we extend that hypothesis to the superfamily of α, β, and γ-synucleins. Our models of numerous Synuclein oligomers, annular protofibrils, tubular protofibrils, lipoproteins, and ion channels were developed to be consistent with sizes, shapes, molecular weights, and secondary structures of assemblies as determined by EM and other studies. The models have the following features: 1) all subunits have identical structures and interactions; 2) they are consistent with conventional β-barrel theory; 3) the distance between walls of adjacent β-barrels is between 0.6 and 1.2 nm; 4) hydrogen bonds, salt bridges, interactions among aromatic side-chains, burial and tight packing of hydrophobic side-chains, and aqueous solvent exposure of hydrophilic side-chains are relatively optimal; and 5) residues that are identical among distantly related homologous proteins cluster in the interior of most oligomers whereas residues that are hypervariable are exposed on protein surfaces. Atomic scale models of some assemblies were developed.

σ factors are essential parts of bacterial RNA polymerase (RNAP) as they allow to recognize promotor sequences and initiate transcription. Domain 1.1 of vegetative σ factors occupies the primary channel of RNAP and also prevents binding of the σ factor to promoter DNA alone. Here, we show that domain 1.1 of Bacillus subtilis σ A exists in two structurally distinct variants in dynamic equilibrium. The major conformation at room temperature is represented by a previously reported well-folded structure solved by nuclear magnetic resonance (NMR), but 4% of the protein molecules are present in a less thermodynamically favorable state. We show that this population increases with temperature and may represent as much as 20% at 43.5 ◦ C. We characterized the minor state of the domain 1.1 using specialized methods of NMR. We found that, in contrast to the major state, the detected minor state is partially unfolded. Its propensity to form secondary structure elements is especially decreased for the first and third α helices, while the second α helix and β strand close to the C-terminus are more stable. In summary, this study reveals conformational dynamics of domain 1.1 and provides a basis for studies of its interaction with RNAP and effects on transcription regulation.

Histidine kinases (HKs) are a central part of bacterial environmental-sensing two component systems, providing their hosts with the ability to respond to a wide range of physical and chemical signals, being a paradigmatic example of protein diversity and molecular evolution. HK are multidomain proteins consisting of at least a sensor domain, dimerization and phosphorylation domain (DHp) and a catalytic domain. They work as homodimers, and the existence of two different autophosphorylation mechanisms (cis and trans), has been proposed as relevant for pathway specificity. Although several HKs have been intensively studied, a precise sequence-to-structure explanation of why and how either cis or trans phosphorylation occurs, is still unavailable, nor is there any evolutionary analysis on the subject. In this work we show that AlphaFold can accurately determine whether a HK dimerizes in a cis or trans structure. We then use it to explore the molecular determinants of the phosphorylation mechanism. We conclude that it is the difference in lengths of the helices surrounding the DHp loop that determines the mechanism. We also show that very small changes in these helices can cause a mechanism switch, providing these systems with a way to diverge rapidly and adapt to new stimuli.

The bacterial adhesin FimH is a model for the study of protein allostery because its structure has been resolved in multiple configurations, including the active and the inactive state. FimH consists of a pilin domain (PD) that anchors it to the rest of the fimbria and an allosterically regulated lectin domain (LD) that binds mannose on the surface of infected cells. Under normal conditions, the two domains are docked to each other and LD binds mannose weakly. However, in the presence of tensile force generated by shear the domains separate and conformational changes propagate across LD resulting in a stronger bond to mannose. Recently, the crystallographic structure of a variant of FimH has been resolved, called FimH FocH, where PD contains 10 mutations near the inter-domain interface. Although the X-ray structures of FimH and FimH FocH are almost identical, experimental evidence shows that FimH FocH is activated even in the absence of shear. Here, molecular dynamics simulations combined with the Jarzinski equality were used to investigate the discrepancy between the crystallographic structures and the functional assays. The results indicate that the free energy barrier of the unbinding process between LD and PD is drastically reduced in FimH FocH. Rupture of an inter-domain hydrogen bond involving R166 constitutes a rate limiting step of the domains separation process and occurs more readily in FimH FocH than FimH. In conclusion, the mutations in FimH FocH shift the equilibrium towards an equal occupancy of bound and unbound states for LD and PD by reducing a rate limiting step.