Molecular mechanism of anion permeation through aquaporin 6.

Aquaporins (AQPs) are recognized as transmembrane water channels that facilitate selective water permeation through their monomeric pores. Among the AQP family, AQP6 has an intriguing characteristic as an anion channel, which is allosterically controlled by pH conditions and is eliminated by a single amino acid mutation. However, the molecular mechanism of anion permeation through AQP6 remains unclear. Using molecular dynamics simulations in the presence of a transmembrane voltage utilizing an ion concentration gradient, we show that chloride ions permeate through the pore corresponding to the central axis of the AQP6 homotetramer. Under low pH conditions, a subtle opening of the hydrophobic selectivity filter (SF), located near the extracellular part of the central pore, becomes wetted and enables anion permeation. Our simulations also indicate that a single mutation (N63G) in human AQP6, located at the central pore, significantly reduces anion conduction, consistent with experimental data. Moreover, we demonstrate that the pH-sensing mechanism in which the protonation of H184 and H189 under low pH conditions allosterically triggers the gating of the SF region. These results suggest a unique pH-dependent allosteric anion permeation mechanism in AQP6 and could clarify the role of the central pore in some of the AQP tetramers.

DeepFold: enhancing protein structure prediction through optimized loss functions, improved template features, and re-optimized energy function.

MOTIVATION: Predicting protein structures with high accuracy is a critical challenge for the broad community of life sciences and industry. Despite progress made by deep neural networks like AlphaFold2, there is a need for further improvements in the quality of detailed structures, such as side-chains, along with protein backbone structures. RESULTS: Building upon the successes of AlphaFold2, the modifications we made include changing the losses of side-chain torsion angles and frame aligned point error, adding loss functions for side chain confidence and secondary structure prediction, and replacing template feature generation with a new alignment method based on conditional random fields. We also performed re-optimization by conformational space annealing using a molecular mechanics energy function which integrates the potential energies obtained from distogram and side-chain prediction. In the CASP15 blind test for single protein and domain modeling (109 domains), DeepFold ranked fourth among 132 groups with improvements in the details of the structure in terms of backbone, side-chain, and Molprobity. In terms of protein backbone accuracy, DeepFold achieved a median GDT-TS score of 88.64 compared with 85.88 of AlphaFold2. For TBM-easy/hard targets, DeepFold ranked at the top based on Z-scores for GDT-TS. This shows its practical value to the structural biology community, which demands highly accurate structures. In addition, a thorough analysis of 55 domains from 39 targets with publicly available structures indicates that DeepFold shows superior side-chain accuracy and Molprobity scores among the top-performing groups. AVAILABILITY AND IMPLEMENTATION: DeepFold tools are open-source software available at https://github.com/newtonjoo/deepfold.

Application of conformational space annealing to the protein structure modeling using cryo-EM maps.

Conformational space annealing (CSA), a global optimization method, has been applied to various protein structure modeling tasks. In this paper, we applied CSA to the cryo-EM structure modeling task by combining the python subroutine of CSA (PyCSA) and the fast relax (FastRelax) protocol of PyRosetta. Refinement of initial structures generated from two methods, rigid fitting of predicted structures to the Cryo-EM map and de novo protein modeling by tracing the Cryo-EM map, was performed by CSA. In the refinement of the rigid-fitted structures, the final models showed that CSA can generate reliable atomic structures of proteins, even when large movements of protein domains were required. In the de novo modeling case, although the overall structural qualities of the final models were rather dependent on the initial models, the final models generated by CSA showed improved MolProbity scores and cross-correlation coefficients to the maps. These results suggest that CSA can accomplish flexible fitting and refinement together by sampling diverse conformations effectively and thus can be utilized for cryo-EM structure modeling.

Direct experimental observation of blue-light-induced conformational change and intermolecular interactions of cryptochrome.

Cryptochromes are blue light receptors that mediate circadian rhythm and magnetic sensing in various organisms. A typical cryptochrome consists of a conserved photolyase homology region domain and a varying carboxyl-terminal extension across species. The structure of the flexible carboxyl-terminal extension and how carboxyl-terminal extension participates in cryptochrome's signaling function remain mostly unknown. In this study, we uncover the potential missing link between carboxyl-terminal extension conformational changes and downstream signaling functions. Specifically, we discover that the blue-light induced opening of carboxyl-terminal extension in C. reinhardtii animal-like cryptochrome can structurally facilitate its interaction with Rhythm Of Chloroplast 15, a circadian-clock-related protein. Our finding is made possible by two technical advances. Using single-molecule Förster resonance energy transfer technique, we directly observe the displacement of carboxyl-terminal extension by about 15 Å upon blue light excitation. Combining structure prediction and solution X-ray scattering methods, we propose plausible structures of full-length cryptochrome under dark and lit conditions. The structures provide molecular basis for light active conformational changes of cryptochrome and downstream regulatory functions.

CRFalign: A Sequence-Structure Alignment of Proteins Based on a Combination of HMM-HMM Comparison and Conditional Random Fields.

Sequence-structure alignment for protein sequences is an important task for the template-based modeling of 3D structures of proteins. Building a reliable sequence-structure alignment is a challenging problem, especially for remote homologue target proteins. We built a method of sequence-structure alignment called CRFalign, which improves upon a base alignment model based on HMM-HMM comparison by employing pairwise conditional random fields in combination with nonlinear scoring functions of structural and sequence features. Nonlinear scoring part is implemented by a set of gradient boosted regression trees. In addition to sequence profile features, various position-dependent structural features are employed including secondary structures and solvent accessibilities. Training is performed on reference alignments at superfamily levels or twilight zone chosen from the SABmark benchmark set. We found that CRFalign method produces relative improvement in terms of average alignment accuracies for validation sets of SABmark benchmark. We also tested CRFalign on 51 sequence-structure pairs involving 15 FM target domains of CASP14, where we could see that CRFalign leads to an improvement in average modeling accuracies in these hard targets (TM-CRFalign ≃42.94%) compared with that of HHalign (TM-HHalign ≃39.05%) and also that of MRFalign (TM-MRFalign ≃36.93%). CRFalign was incorporated to our template search framework called CRFpred and was tested for a random target set of 300 target proteins consisting of Easy, Medium and Hard sets which showed a reasonable template search performance.

A Coil-to-Helix Transition Serves as a Binding Motif for hSNF5 and BAF155 Interaction.

Human SNF5 and BAF155 constitute the core subunit of multi-protein SWI/SNF chromatin-remodeling complexes that are required for ATP-dependent nucleosome mobility and transcriptional control. Human SNF5 (hSNF5) utilizes its repeat 1 (RPT1) domain to associate with the SWIRM domain of BAF155. Here, we employed X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, and various biophysical methods in order to investigate the detailed binding mechanism between hSNF5 and BAF155. Multi-angle light scattering data clearly indicate that hSNF5171-258 and BAF155SWIRM are both monomeric in solution and they form a heterodimer. NMR data and crystal structure of the hSNF5171-258/BAF155SWIRM complex further reveal a unique binding interface, which involves a coil-to-helix transition upon protein binding. The newly formed αN helix of hSNF5171-258 interacts with the ß2-α1 loop of hSNF5 via hydrogen bonds and it also displays a hydrophobic interaction with BAF155SWIRM. Therefore, the N-terminal region of hSNF5171-258 plays an important role in tumorigenesis and our data will provide a structural clue for the pathogenesis of Rhabdoid tumors and malignant melanomas that originate from mutations in the N-terminal loop region of hSNF5.

Supramolecular tholos-like architecture constituted by archaeal proteins without functional annotation.

Euryarchaeal genomes encode proteasome-assembling chaperone homologs, PbaA and PbaB, although archaeal proteasome formation is a chaperone-independent process. Homotetrameric PbaB functions as a proteasome activator, while PbaA forms a homopentamer that does not interact with the proteasome. Notably, PbaA forms a complex with PF0014, an archaeal protein without functional annotation. In this study, based on our previous research on PbaA crystal structure, we performed an integrative analysis of the supramolecular structure of the PbaA/PF0014 complex using native mass spectrometry, solution scattering, high-speed atomic force microscopy, and electron microscopy. The results indicated that this highly thermostable complex constitutes ten PbaA and ten PF0014 molecules, which are assembled into a dumbbell-shaped structure. Two PbaA homopentameric rings correspond to the dumbbell plates, with their N-termini located outside of the plates and C-terminal segments left mobile. Furthermore, mutant PbaA lacking the mobile C-terminal segment retained the ability to form a complex with PF0014, allowing 3D modeling of the complex. The complex shows a five-column tholos-like architecture, in which each column comprises homodimeric PF0014, harboring a central cavity, which can potentially accommodate biomacromolecules including proteins. Our findings provide insight into the functional roles of Pba family proteins, offering a novel framework for designing functional protein cages.

Improved Consensus-Fragment Selection in Template-Assisted Prediction of Protein Structures with the UNRES Force Field in CASP13.

The method for protein-structure prediction, which combines the physics-based coarse-grained UNRES force field with knowledge-based modeling, has been developed further and tested in the 13th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP13). The method implements restraints from the consensus fragments common to server models. In this work, the server models to derive fragments have been chosen on the basis of quality assessment; a fully automatic fragment-selection procedure has been introduced, and Dynamic Fragment Assembly pseudopotentials have been fully implemented. The Global Distance Test Score (GDT_TS), averaged over our "Model 1" predictions, increased by over 10 units with respect to CASP12 for the free-modeling category to reach 40.82. Our "Model 1" predictions ranked 20 and 14 for all and free-modeling targets, respectively (upper 20.2% and 14.3% of all models submitted to CASP13 in these categories, respectively), compared to 27 (upper 21.1%) and 24 (upper 18.9%) in CASP12, respectively. For oligomeric targets, the Interface Patch Similarity (IPS) and Interface Contact Similarity (ICS) averaged over our best oligomer models increased from 0.28 to 0.36 and from 12.4 to 17.8, respectively, from CASP12 to CASP13, and top-ranking models of 2 targets (H0968 and T0997o) were obtained (none in CASP12). The improvement of our method in CASP13 over CASP12 was ascribed to the combined effect of the overall enhancement of server-model quality, our success in selecting server models and fragments to derive restraints, and improvements of the restraint and potential-energy functions.

A Mutation in ZNF143 as a Novel Candidate Gene for Endothelial Corneal Dystrophy.

Corneal dystrophies (CDs) are a diverse group of inherited disorders with a heterogeneous genetic background. Here, we report the identification of a novel ZNF143 heterozygous missense mutation in three individuals of the same family with clinical and pathological features that are consistent with endothelial CD. Ophthalmologic examination revealed diffuse corneal clouding and edema with decreased endothelial cell density. Pathological findings showed increased corneal thickness due to edema of basal epithelial cells and stroma, and abnormal metaplastic endothelium with stratified epithelium-like changes. Patients' metaplastic corneal endothelial cells expressed predominantly cytokerain 7, cytokeratin 19, and E-cadherin. Although Sanger sequencing did not detect any mutation associated with endothelial CDs, whole exome sequencing identified the ZNF143 c.937G>C p.(Asp313His) mutation as a candidate gene for our patients' endothelial CD. In-vitro functional studies demonstrated that mutant ZNF143 promoted the mesenchymal-to-epithelial transition; it upregulated the expression of genes associated with epithelialization in human corneal endothelial cells. Additionally, proinflammatory cytokine responsive genes were significantly enriched after mutant ZNF143 transfection, which may contribute to the severe phenotype of the three patients. These findings link a mutation in ZNF143 with endothelial CD for the first time.

Non-sequential protein structure alignment by conformational space annealing and local refinement.

Protein structure alignment is an important tool for studying evolutionary biology and protein modeling. A tool which intensively searches for the globally optimal non-sequential alignments is rarely found. We propose ALIGN-CSA which shows improvement in scores, such as DALI-score, SP-score, SO-score and TM-score over the benchmark set including 286 cases. We performed benchmarking of existing popular alignment scoring functions, where the dependence of the search algorithm was effectively eliminated by using ALIGN-CSA. For the benchmarking, we set the minimum block size to 4 to prevent much fragmented alignments where the biological relevance of small alignment blocks is hard to interpret. With this condition, globally optimal alignments were searched by ALIGN-CSA using the four scoring functions listed above, and TM-score is found to be the most effective in generating alignments with longer match lengths and smaller RMSD values. However, DALI-score is the most effective in generating alignments similar to the manually curated reference alignments, which implies that DALI-score is more biologically relevant score. Due to the high demand on computational resources of ALIGN-CSA, we also propose a relatively fast local refinement method, which can control the minimum block size and whether to allow the reverse alignment. ALIGN-CSA can be used to obtain much improved alignment at the cost of relatively more extensive computation. For faster alignment, we propose a refinement protocol that improves the score of a given alignment obtained by various external tools. All programs are available from http://lee.kias.re.kr.

ConDo: protein domain boundary prediction using coevolutionary information.

MOTIVATION: Domain boundary prediction is one of the most important problems in the study of protein structure and function. Many sequence-based domain boundary prediction methods are either template-based or machine learning (ML) based. ML-based methods often perform poorly due to their use of only local (i.e. short-range) features. These conventional features such as sequence profiles, secondary structures and solvent accessibilities are typically restricted to be within 20 residues of the domain boundary candidate. RESULTS: To address the performance of ML-based methods, we developed a new protein domain boundary prediction method (ConDo) that utilizes novel long-range features such as coevolutionary information in addition to the aforementioned local window features as inputs for ML. Toward this purpose, two types of coevolutionary information were extracted from multiple sequence alignment using direct coupling analysis: (i) partially aligned sequences, and (ii) correlated mutation information. Both the partially aligned sequence information and the modularity of residue-residue couplings possess long-range correlation information. AVAILABILITY AND IMPLEMENTATION: https://github.com/gicsaw/ConDo.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Morphological Evaluation of Meta-stable Oligomers of α-Synuclein with Small-Angle Neutron Scattering.

Amyloidogenesis of α-synuclein (αS) is considered to be a pathological phenomenon related to Parkinson's disease (PD). As a key component to reveal the fibrillation mechanism and toxicity, we have investigated an oligomeric species of αS capable of exhibiting the unit-assembly process leading to accelerated amyloid fibril formation. These oligomers previously shown to exist in a meta-stable state with mostly disordered structure and unable to seed the fibrillation were converted to either temperature-sensitive self-associative oligomers or NaCl-induced non-fibrillating oligomeric species. Despite their transient and disordered nature, the structural information of meta-stable αS oligomers (Meta-αS-Os) was successfully evaluated with small-angle neutron scattering (SANS) technique. By fitting the neutron scattering data with polydisperse Gaussian Coil (pGC) model, Meta-αS-O was analyzed as a sphere with approximate diameter of 100 Å. Its overall shape altered drastically with subtle changes in temperature between 37 °C and 43 °C, which would be responsible for fibrillar polymorphism. Based on their bifurcating property of Meta-αS-Os leading to either on-pathway or off-pathway species, the oligomers could be suggested as a crucial intermediate responsible for the oligomeric diversification and multiple fibrillation processes. Therefore, Meta-αS-Os could be considered as a principal target to control the amyloidogenesis and its pathogenesis.

Use of the UNRES force field in template-assisted prediction of protein structures and the refinement of server models: Test with CASP12 targets.

Knowledge-based methods are, at present, the most effective ones for the prediction of protein structures; however, their results heavily depend on the similarity of a target sequence to those of proteins with known structures. On the other hand, the physics-based methods, although still less accurate and more expensive to execute, are independent of databases and give reasonable results where the knowledge-based methods fail because of weak sequence similarity. Therefore, a plausible approach seems to be the use of knowledge-based methods to determine the sections of the structures that correspond to sufficient sequence similarity and physics-based methods to determine the remaining structure. By participating in the 12th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP12) as the KIAS-Gdansk group, we tested our recently developed hybrid approach, in which protein-structure prediction is carried out by using the physics-based UNRES coarse-grained energy function, with restraints derived from the server models. Best predictions among all groups were obtained for 2 targets and 80% of our models were in the upper 50% of the models submitted to CASP. Our method was also able to exclude, with about 70% confidence, the information from the servers that performed poorly on a given target. Moreover, the method resulted in the best models of 2 refinement targets and performed remarkably well on oligomeric targets.

Data-assisted protein structure modeling by global optimization in CASP12.

In CASP12, 2 types of data-assisted protein structure modeling were experimented. Either SAXS experimental data or cross-linking experimental data was provided for a selected number of CASP12 targets that the CASP12 predictor could utilize for better protein structure modeling. We devised 2 separate energy terms for SAXS data and cross-linking data to drive the model structures into more native-like structures that satisfied the given experimental data as much as possible. In CASP11, we successfully performed protein structure modeling using simulated sparse and ambiguously assigned NOE data and/or correct residue-residue contact information, where the only energy term that folded the protein into its native structure was the term which was originated from the given experimental data. However, the 2 types of experimental data provided in CASP12 were far from being sufficient enough to fold the target protein into its native structure because SAXS data provides only the overall shape of the molecule and the cross-linking contact information provides only very low-resolution distance information. For this reason, we combined the SAXS or cross-linking energy term with our regular modeling energy function that includes both the template energy term and the de novo energy terms. By optimizing the newly formulated energy function, we obtained protein models that fit better with provided SAXS data than the X-ray structure of the target. However, the improvement of the model relative to the 1 modeled without the SAXS data, was not significant. Consistent structural improvement was achieved by incorporating cross-linking data into the protein structure modeling.

Protein structure modeling and refinement by global optimization in CASP12.

For protein structure modeling in the CASP12 experiment, we have developed a new protocol based on our previous CASP11 approach. The global optimization method of conformational space annealing (CSA) was applied to 3 stages of modeling: multiple sequence-structure alignment, three-dimensional (3D) chain building, and side-chain re-modeling. For better template selection and model selection, we updated our model quality assessment (QA) method with the newly developed SVMQA (support vector machine for quality assessment). For 3D chain building, we updated our energy function by including restraints generated from predicted residue-residue contacts. New energy terms for the predicted secondary structure and predicted solvent accessible surface area were also introduced. For difficult targets, we proposed a new method, LEEab, where the template term played a less significant role than it did in LEE, complemented by increased contributions from other terms such as the predicted contact term. For TBM (template-based modeling) targets, LEE performed better than LEEab, but for FM targets, LEEab was better. For model refinement, we modified our CASP11 molecular dynamics (MD) based protocol by using explicit solvents and tuning down restraint weights. Refinement results from MD simulations that used a new augmented statistical energy term in the force field were quite promising. Finally, when using inaccurate information (such as the predicted contacts), it was important to use the Lorentzian function for which the maximal penalty arising from wrong information is always bounded.

Methods for estimation of model accuracy in CASP12.

Methods to reliably estimate the quality of 3D models of proteins are essential drivers for the wide adoption and serious acceptance of protein structure predictions by life scientists. In this article, the most successful groups in CASP12 describe their latest methods for estimates of model accuracy (EMA). We show that pure single model accuracy estimation methods have shown clear progress since CASP11; the 3 top methods (MESHI, ProQ3, SVMQA) all perform better than the top method of CASP11 (ProQ2). Although the pure single model accuracy estimation methods outperform quasi-single (ModFOLD6 variations) and consensus methods (Pcons, ModFOLDclust2, Pcomb-domain, and Wallner) in model selection, they are still not as good as those methods in absolute model quality estimation and predictions of local quality. Finally, we show that when using contact-based model quality measures (CAD, lDDT) the single model quality methods perform relatively better.

An artificially constructed dimer through deformation of a short zinc-binding loop.

We have analyzed the crystal structure of the dimeric form of d-glycero-d-manno-heptose-1,7-bisphosphate phosphatase from Burkholderia thailandensis (BtGmhB), catalyzing the removal of the phosphate at the 7 position of d-glycero-d-manno-heptose-1,7-bisphosphate. The crystal structure of BtGmhB revealed a dimeric form caused by a disruption of a short zinc-binding loop. The dimeric BtGmhB structure was induced by triggering the loss of Zn2+via the protonation of cysteine residues at pH 4.8 of the crystallization condition. Similarly, the addition of EDTA also causes the dimerization of BtGmhB. It appears there are two dimeric forms in solution with and without the disulfide bridge mediated by Cys95. The disulfide-free dimer produced by the loss of Zn2+ in the short zinc-binding loop is further converted to a stable disulfide-bonded dimer in vitro. Though the two dimeric forms are reversible, both of them are inactive due to a deformation of the active site. Single and triple mutant experiments confirmed the presence of two dimeric forms in vitro. Phosphatase assay results showed that only a zinc-bound monomeric form contains catalytic activity in contrast to the inactive zinc-free dimeric forms. The monomer-to-dimer transition caused by the loss of Zn2+ observed in this study is an example of reversal phenomenon caused by artificial proteins containing protein engineered zinc-finger motifs where the monomer-to-dimer transitions occurred in the presence of Zn2+. Therefore, this unusual dimerization process may be applicable to designing proteins possessing a short zinc-binding loop with a novel regulatory role.

Ergodicity and model quality in template-restrained canonical and temperature/Hamiltonian replica exchange coarse-grained molecular dynamics simulations of proteins.

Molecular simulations restrained to single or multiple templates are commonly used in protein-structure modeling. However, the restraints introduce additional barriers, thus impairing the ergodicity of simulations, which can affect the quality of the resulting models. In this work, the effect of restraint types and simulation schemes on ergodicity and model quality was investigated by performing template-restrained canonical molecular dynamics (MD), multiplexed replica-exchange molecular dynamics, and Hamiltonian replica exchange molecular dynamics (HREMD) simulations with the coarse-grained UNRES force field on nine selected proteins, with pseudo-harmonic log-Gaussian (unbounded) or Lorentzian (bounded) restraint functions. The best ergodicity was exhibited by HREMD. It has been found that non-ergodicity does not affect model quality if good templates are used to generate restraints. However, when poor-quality restraints not covering the entire protein are used, the improved ergodicity of HREMD can lead to significantly improved protein models. © 2017 Wiley Periodicals, Inc.

Architecture of the type IV coupling protein complex of Legionella pneumophila.

Many bacteria, including Legionella pneumophila, rely on the type IV secretion system to translocate a repertoire of effector proteins into the hosts for their survival and growth. Type IV coupling protein (T4CP) is a hexameric ATPase that links translocating substrates to the transenvelope secretion conduit. Yet, how a large number of effector proteins are selectively recruited and processed by T4CPs remains enigmatic. DotL, the T4CP of L. pneumophila, contains an ATPase domain and a C-terminal extension whose function is unknown. Unlike T4CPs involved in plasmid DNA translocation, DotL appeared to function by forming a multiprotein complex with four other proteins. Here, we show that the C-terminal extension of DotL interacts with DotN, IcmS, IcmW and an additionally identified subunit LvgA, and that this pentameric assembly binds Legionella effector proteins. We determined the crystal structure of this assembly and built an architecture of the T4CP holocomplex by combining a homology model of the ATPase domain of DotL. The holocomplex is a hexamer of a bipartite structure composed of a membrane-proximal ATPase domain and a membrane-distal substrate-recognition assembly. The presented information demonstrates the architecture and functional dissection of the multiprotein T4CP complexes and provides important insights into their substrate recruitment and processing.

Protein Loop Structure Prediction Using Conformational Space Annealing.

We have developed a protein loop structure prediction method by combining a new energy function, which we call EPLM (energy for protein loop modeling), with the conformational space annealing (CSA) global optimization algorithm. The energy function includes stereochemistry, dynamic fragment assembly, distance-scaled finite ideal gas reference (DFIRE), and generalized orientation- and distance-dependent terms. For the conformational search of loop structures, we used the CSA algorithm, which has been quite successful in dealing with various hard global optimization problems. We assessed the performance of EPLM with two widely used loop-decoy sets, Jacobson and RAPPER, and compared the results against the DFIRE potential. The accuracy of model selection from a pool of loop decoys as well as de novo loop modeling starting from randomly generated structures was examined separately. For the selection of a nativelike structure from a decoy set, EPLM was more accurate than DFIRE in the case of the Jacobson set and had similar accuracy in the case of the RAPPER set. In terms of sampling more nativelike loop structures, EPLM outperformed EDFIRE for both decoy sets. This new approach equipped with EPLM and CSA can serve as the state-of-the-art de novo loop modeling method.