Modeling Protein Structure from its Sequence

In silico study on the effects of disulfide bonds in ORF8 of SARS-CoV-2

The COVID-19 epidemic, caused by virus SARS-CoV-2, has been a pandemic and threatening everyone's health in the past two years. In SARS-CoV-2, ORF8 is one of the most important accessory...

RASSF effectors couple diverse RAS subfamily GTPases to the Hippo pathway

Small guanosine triphosphatases (GTPases) of the RAS superfamily signal by directly binding to multiple downstream effector proteins. Effectors are defined by a folded RAS-association (RA) domain that binds exclusively to GTP-loaded (activated) RAS, but the binding specificities of most RA domains toward more than 160 RAS superfamily GTPases have not been characterized. Ten RA domain family (RASSF) proteins comprise the largest group of related effectors and are proposed to couple RAS to the proapoptotic Hippo pathway. Here, we showed that RASSF1-6 formed complexes with the Hippo kinase ortholog MST1, whereas RASSF7-10 formed oligomers with the p53-regulating effectors ASPP1 and ASPP2. Moreover, only RASSF5 bound directly to activated HRAS and KRAS, and RASSFs did not augment apoptotic induction downstream of RAS oncoproteins. Structural modeling revealed that expansion of the RASSF effector family in vertebrates included amino acid substitutions to key residues that direct GTPase-binding specificity. We demonstrated that the tumor suppressor RASSF1A formed complexes with the RAS-related GTPases GEM, REM1, REM2, and the enigmatic RASL12. Furthermore, interactions between RASSFs and RAS GTPases blocked YAP1 nuclear localization. Thus, these simple scaffolds link the activation of diverse RAS family small G proteins to Hippo or p53 regulation.

Actin Cytoskeletal Reorganization Function of JRAB/MICAL-L2 Is Fine-tuned by Intramolecular Interaction between First LIM Zinc Finger and C-terminal Coiled-coil Domains

JRAB/MICAL-L2 is an effector protein of Rab13, a member of the Rab family of small GTPase. JRAB/MICAL-L2 consists of a calponin homology domain, a LIM domain, and a coiled-coil domain. JRAB/MICAL-L2 engages in intramolecular interaction between the N-terminal LIM domain and the C-terminal coiled-coil domain, and changes its conformation from closed to open under the effect of Rab13. Open-form JRAB/MICAL-L2 induces the formation of peripheral ruffles via an interaction between its calponin homology domain and filamin. Here, we report that the LIM domain, independent of the C-terminus, is also necessary for the function of open-form JRAB/MICAL-L2. In mechanistic terms, two zinc finger domains within the LIM domain bind the first and second molecules of actin at the minus end, potentially inhibiting the depolymerization of actin filaments (F-actin). The first zinc finger domain also contributes to the intramolecular interaction of JRAB/MICAL-L2. Moreover, the residues of the first zinc finger domain that are responsible for the intramolecular interaction are also involved in the association with F-actin. Together, our findings show that the function of open-form JRAB/MICAL-L2 mediated by the LIM domain is fine-tuned by the intramolecular interaction between the first zinc finger domain and the C-terminal domain.

Outer membrane phospholipase A's roles in Helicobacter pylori acid adaptation

Background

The pH of the human gastric mucosa varies around 2.5 so that only bacteria with strong acidic stress tolerance can colonize it. The ulcer causing Helicobacter pylori thrives in the gastric mucosa. We analyse the roles of the key outer membrane protein OMPLA in its roles in acid tolerance.

Results

The homology model of Helicobacter pylori outer membrane phospholipase A (OMPLA) reveals a twelve stranded β-barrel with a pore that allows molecules to pass with a diameter up to 4 Å. Structure based multiple sequence alignments revealed the functional roles of many amino acids, and led to the suggestion that OMPLA has multiple functions. Besides its role as phospholipase it lets urea enter and ammonium exit the periplasm. Combined with an extensive literature study, our work leads to a comprehensive model for H. pylori’s acid tolerance. This model is based on the conversion of urea into ammonium, and it includes multiple roles for OMPLA and involves two hitherto little studied membrane channels in the OMPLA operon.

Conclusion

The three-dimensional model of OMPLA predicts a transmembrane pore that can aid H. pylori’s acid tolerance through urea influx and ammonium efflux. After urea passes through OMPLA into the periplasm, it passes through the pH-gated inner membrane channel UreI into the cytoplasm where urease hydrolyses it into NH₃ and CO₂. Most of the NH₃ becomes NH₄⁺ that is likely to need an inner membrane channel to reach the periplasm. Two genes that are co-regulated with OMPLA in gastric Helicobacter operons could aid this transport. The NH₄⁺ that might leave the cell through the OMPLA pore has been implicated in H. pylor’s pathogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13099-017-0184-y) contains supplementary material, which is available to authorized users.

Local alignment refinement using structural assessment

Homology modeling is the most commonly used technique to build a three-dimensional model for a protein sequence. It heavily relies on the quality of the sequence alignment between the protein to model and related proteins with a known three dimensional structure. Alignment quality can be assessed according to the physico-chemical properties of the three dimensional models it produces.

In this work, we introduce fifteen predictors designed to evaluate the properties of the models obtained for various alignments. They consist of an energy value obtained from different force fields (CHARMM, ProsaII or ANOLEA) computed on residue selected around misaligned regions. These predictors were evaluated on ten challenging test cases. For each target, all possible ungapped alignments are generated and their corresponding models are computed and evaluated.

The best predictor, retrieving the structural alignment for 9 out of 10 test cases, is based on the ANOLEA atomistic mean force potential and takes into account residues around misaligned secondary structure elements. The performance of the other predictors is significantly lower. This work shows that substantial improvement in local alignments can be obtained by careful assessment of the local structure of the resulting models.