The nature of the turn in omega loops of proteins

Direct Inhibition of the YAP : TEAD Interaction: An Unprecedented Drug Discovery Challenge

The Hippo pathway, which is key in organ morphogenesis, is frequently deregulated in cancer. The TEAD (TEA domain family member) transcription factors are the most distal elements of this pathway, and their activity is regulated by proteins such as YAP (Yes-associated protein). The identification of inhibitors of the YAP : TEAD interaction is one approach to develop novel anticancer drugs: the first clinical candidate (IAG933) preventing the association between these two proteins by direct competition has just been reported. The discovery of this molecule was particularly challenging because the interface between these two proteins is large (~3500 Å2 buried in complex formation) and made up of distinct contact areas. The most critical of these involves an omega-loop (Ω-loop), a secondary structure element rarely found in protein-protein interactions. This review summarizes how the knowledge gained from structure-function studies of the interaction between the Ω-loop of YAP and TEAD was used to devise the strategy to identify potent low-molecular weight compounds that show a pronounced anti-tumor effect.

Comprehensive assessment of protein loop modeling programs on large-scale datasets: prediction accuracy and efficiency

Protein loops play a critical role in the dynamics of proteins and are essential for numerous biological functions, and various computational approaches to loop modeling have been proposed over the past decades. However, a comprehensive understanding of the strengths and weaknesses of each method is lacking. In this work, we constructed two high-quality datasets (i.e. the General dataset and the CASP dataset) and systematically evaluated the accuracy and efficiency of 13 commonly used loop modeling approaches from the perspective of loop lengths, protein classes and residue types. The results indicate that the knowledge-based method FREAD generally outperforms the other tested programs in most cases, but encountered challenges when predicting loops longer than 15 and 30 residues on the CASP and General datasets, respectively. The ab initio method Rosetta NGK demonstrated exceptional modeling accuracy for short loops with four to eight residues and achieved the highest success rate on the CASP dataset. The well-known AlphaFold2 and RoseTTAFold require more resources for better performance, but they exhibit promise for predicting loops longer than 16 and 30 residues in the CASP and General datasets. These observations can provide valuable insights for selecting suitable methods for specific loop modeling tasks and contribute to future advancements in the field.

Olfactory marker protein contains a leucine-rich domain in the Ω-loop important for nuclear export

Olfactory marker protein (OMP) is a cytosolic protein expressed in mature olfactory receptor neurons (ORNs). OMP modulates cAMP signalling and regulates olfactory sensation and axonal targeting. OMP is a small soluble protein, and passive diffusion between nucleus and cytoplasm is expected. However, OMP is mostly situated in the cytosol and is only sparsely detected in the nuclei of a subset of ORNs, hypothalamic neurons and heterologously OMP-expressing cultured cells. OMP can enter the nucleus in association with transcription factors. However, how OMP is retained in the cytosol at rest is unclear. Because OMP is proposed to affect cell differentiation, it is important to understand how OMP is distributed between cytoplasm and nucleus. To elucidate the structural profile of OMP, we applied several bioinformatics methods to a multiple sequence alignment (MSA) of OMP protein sequences and ranked the evolutionarily conserved residues. In addition to the previously reported cAMP-binding domain, we identified a leucine-rich domain in the Ω-loop of OMP. We introduced mutations into the leucine-rich region and heterologously expressed the mutant OMP in HEK293T cells. Mutations into alanine increased the nuclear distribution of OMP quantified by immunocytochemistry and western blotting. Therefore, we concluded that OMP contains a leucine-rich domain important for nuclear transport.

Thyroxine binding to type III iodothyronine deiodinase

Iodothyronine deiodinases (Dios) are important selenoproteins that control the concentration of the active thyroid hormone (TH) triiodothyronine through regioselective deiodination. The X-ray structure of a truncated monomer of Type III Dio (Dio3), which deiodinates TH inner rings through a selenocysteine (Sec) residue, revealed a thioredoxin-fold catalytic domain supplemented with an unstructured Ω-loop. Loop dynamics are driven by interactions of the conserved Trp207 with solvent in multi-microsecond molecular dynamics simulations of the Dio3 thioredoxin(Trx)-fold domain. Hydrogen bonding interactions of Glu200 with residues conserved across the Dio family anchor the loop’s N-terminus to the active site Ser-Cys-Thr-Sec sequence. A key long-lived loop conformation coincides with the opening of a cryptic pocket that accommodates thyroxine (T₄) through an I⋯Se halogen bond to Sec170 and the amino acid group with a polar cleft. The Dio3-T₄ complex is stabilized by an I⋯O halogen bond between an outer ring iodine and Asp211, consistent with Dio3 selectivity for inner ring deiodination. Non-conservation of residues, such as Asp211, in other Dio types in the flexible portion of the loop sequence suggests a mechanism for regioselectivity through Dio type-specific loop conformations. Cys168 is proposed to attack the selenenyl iodide intermediate to regenerate Dio3 based upon structural comparison with related Trx-fold proteins.

ω-Turn: a novel β-turn mimic in globular proteins stabilized by main-chain to side-chain C−H···O interaction

Mimicry of structural motifs is a common feature in proteins. The 10-membered hydrogen-bonded ring involving the main-chain C − O in a β-turn can be formed using a side-chain carbonyl group leading to Asx-turn. We show that the N − H component of hydrogen bond can be replaced by a C(γ) -H group in the side chain, culminating in a nonconventional C − H···O interaction. Because of its shape this β-turn mimic is designated as ω-turn, which is found to occur ∼ three times per 100 residues. Three residues (i to i + 2) constitute the turn with the C − H···O interaction occurring between the terminal residues, constraining the torsion angles ϕi + 1, ψi + 1, ϕi + 2 and χ'1(i + 2) (using the interacting C(γ) atom). Based on these angles there are two types of ω-turns, each of which can be further divided into two groups. C(β) -branched side-chains, and Met and Gln have high propensities to occur at i + 2; for the last two residues the carbonyl oxygen may participate in an additional interaction involving the S and amino group, respectively. With Cys occupying the i + 1 position, such turns are found in the metal-binding sites. N-linked glycosylation occurs at the consensus pattern Asn-Xaa-Ser/Thr; with Thr at i + 2, the sequence can adopt the secondary structure of a ω-turn, which may be the recognition site for protein modification. Location between two β-strands is the most common occurrence in protein tertiary structure, and being generally exposed ω-turn may constitute the antigenic determinant site. It is a stable scaffold and may be used in protein engineering and peptide design.

Structure of the type III secretion effector protein ExoU in complex with its chaperone SpcU

Disease causing bacteria often manipulate host cells in a way that facilitates the infectious process. Many pathogenic gram-negative bacteria accomplish this by using type III secretion systems. In these complex secretion pathways, bacterial chaperones direct effector proteins to a needle-like secretion apparatus, which then delivers the effector protein into the host cell cytosol. The effector protein ExoU and its chaperone SpcU are components of the Pseudomonas aeruginosa type III secretion system. Secretion of ExoU has been associated with more severe infections in both humans and animal models. Here we describe the 1.92 Å X-ray structure of the ExoU–SpcU complex, a full-length type III effector in complex with its full-length cognate chaperone. Our crystallographic data allow a better understanding of the mechanism by which ExoU kills host cells and provides a foundation for future studies aimed at designing inhibitors of this potent toxin.

MUC4 gene polymorphisms associate with endometriosis development and endometriosis-related infertility

BACKGROUND

Mucin 4 (MUC4) plays an important role in protecting and lubricating the epithelial surface of reproductive tracts, but its role in the pathogenesis of endometriosis is largely unknown.

METHODS

To correlate MUC4 polymorphism with the risk of endometriosis and endometriosis-related infertility, we performed a case-control study of 140 patients and 150 healthy women. Six unique single-nucleotide polymorphisms (SNPs) (rs882605, rs1104760, rs2688513, rs2246901, rs2258447 and rs2291652) were selected for this study. DNA fragments containing the target SNP sites were amplified by polymerase chain reaction using the TaqMan SNP Genotyping Assay System to evaluate allele frequency and distribution of genotype in MUC4 polymorphisms.

RESULTS

Both the T/G genotype of rs882605 and the frequency of haplotype T-T (rs882605 and rs1104760) were higher in patients than in controls and were statistically significant. The frequency of the C allele at rs1104760, the C allele at rs2688513, the G allele at rs2246901 and the A allele at rs2258447 were associated with advanced stage of endometriosis. Moreover, the G allele at rs882605 was verified as a key genetic factor for infertility in patients. Protein sequence analysis indicated that amino acid substitutions by genetic variations at rs882605, rs2688513 and rs2246901 occur in the putative functional loops and the type D von Willebrand factor (VWFD) domain in the MUC4 sequence.

CONCLUSIONS

MUC4 polymorphisms are associated with endometriosis development and endometriosis-related infertility in the Taiwanese population.

Evidence of conformational changes in adsorbed lysozyme molecule on silver colloids

In this article, we discuss metal-protein interactions in the Ag-lysozyme complex by spectroscopic measurements. The analysis of the variation in relative intensities of SERS bands reveals the orientation and the change in conformation of the protein molecules on the Ag surface with time. The interaction kinetics of metal-protein complexes has been analyzed over a period of 3h via Raman measurements. Our analysis indicates that the Ag nanoparticles most likely interact with Trp123 which is in close proximity to Phe34 of the lysozyme molecule.

Estimating the "steric clash" at cis peptide bonds

To account for the scarcity of cis peptide bonds in proteins, especially in nonproline (or secondary amide) cases, a steric-clash argument is often put forward, in a scheme where the R lateral chains are facing parallel one another, and the backbone is kept in an "all- trans"-like arrangement. Although such a steric conflict can be partly relieved through proper adjustment of the backbone dihedral angles, one can try to estimate its associated energy cost. To this end, quantum-chemistry approaches using a differential-torsion protocol and bond-separation-energy analyses are applied to N-ethyl propionamide CH3-CH2-CO-NH-CH2-CH3, regarded as a model capable of exhibiting C beta...C beta interaction as in alanine succession. The calculations provide an increment of 9 kcal/mol, quite close to that obtained in the nearly isostere (gsg) rotamer of n-hexane (10 kcal/mol), suggesting the local effects induced by methyl-methyl contact are similar in both cases. Analogous treatments on larger radicals as encountered in leucine or phenylalanine dimers do not change this increment much, which therefore defines the basic reference per-plaque quota to be overcome along all- cis chains. Explicit modeling indicated it can be reduced by up to a factor of 4.

Inhibition of Ribonuclease A by polyphenols present in green tea

We report the effect of the natural polyphenolic compounds from green tea on the catalytic activity of Ribonuclease A (RNase A). The compounds behave as noncompetitive inhibitors of the protein with inhibition constants ranging from 80-1300 microM. The dissociation constants range from 50-150 microM for the RNase A-polyphenol complexes as determined by ultraviolet (UV) and circular dichroism (CD) studies. We have also investigated the changes in the secondary structure of RNase A on complex formation by CD and Fourier transformed infrared (FTIR) spectroscopy. The presence of the gallate moiety has been shown to be important for the inhibition of enzymatic activity. Docking studies for these compounds indicate that the preferred site of binding is the region encompassing residues 34-39 with possible hydrogen bonding with Lys 7 and Arg 10. Finally we have also looked at changes in the accessible surface area of the interacting residues on complex formation for an insight into the residues involved in the interaction.

Active TEM-1 beta-lactamase mutants with random peptides inserted in three contiguous surface loops

Engineering of alternative binding sites on the surface of an enzyme while preserving the enzymatic activity would offer new opportunities for controlling the activity by binding of non-natural ligands. Loops and turns are the natural substructures in which binding sites might be engineered with this purpose. We have genetically inserted random peptide sequences into three relatively rigid and contiguous loops of the TEM-1 beta-lactamase and assessed the tolerance to insertion by the percentage of active mutants. Our results indicate that tolerance to insertion could not be correlated to tolerance to mutagenesis. A turn between two beta-strands bordering the active site was observed to be tolerant to random mutagenesis but not to insertions. Two rigid loops comprising rather well-conserved amino acid residues tolerated insertions, although with some constraints. Insertions between the N-terminal helix and the first beta-strand generated active libraries if cysteine residues were included at both ends of the insert, suggesting the requirement for a stabilizing disulfide bridge. Random sequences were relatively well accommodated within the loop connecting the final beta-strand to the C-terminal helix, particularly if the wild-type residue was retained at one of the loops' end. This suggests two strategies for increasing the percentage of active mutants in insertion libraries. The amino acid distribution in the engineered loops was analyzed and found to be less biased against hydrophobic residues than in natural medium-sized loops. The combination of these activity-selected libraries generated a huge library containing active hybrid enzymes with all three loops modified.

Ab Initio and Empirical Model MD Simulation Studies of Solvent Effects on the Properties of N-Methylacetamide along a cis−trans Isomerization Pathway

The properties of N-methylacetamide along a cis-trans isomerization pathway described by twisting about the C(O)-N bond are examined at finite temperature both in vacuo and in explicit water solvent. Two distinctly different theoretical descriptions, an ab initio (DFT-BLYP) and an empirical (CHARMM22) model, are studied in order to permit an assessment of the dominant forces active in the system. An analysis of the solvent structure at equilibrium and changes in solvation structure accompanying isomerization is, therefore, given for each model. Many-body polarization effects absent under CHARMM22 but present in the ab initio model are found to have a profound influence on the system. The electronic structure of the NMA molecule predicted by the ab initio method along the reaction coordinate is examined in order to shed further light on changes in peptide "partial-double" bond character [C(O)-N] as isomerization takes place. A new statistical-mechanical interpretation of the entropy change during a chemical reaction is presented to help interpret the thermochemistry of the simple reaction.

Optimization of the GB/SA solvation model for predicting the structure of surface loops in proteins

Implicit solvation models are commonly optimized with respect to experimental data or Poisson-Boltzmann (PB) results obtained for small molecules, where the force field is sometimes not considered. In previous studies, we have developed an optimization procedure for cyclic peptides and surface loops in proteins based on the entire system studied and the specific force field used. Thus, the loop has been modeled by the simplified solvation function E(tot) = E(FF) (epsilon = 2r) + Sigma(i) sigma(i)A(i), where E(FF) (epsilon = nr) is the AMBER force field energy with a distance-dependent dielectric function, epsilon = nr, A(i) is the solvent accessible surface area of atom i, and sigma(i) is its atomic solvation parameter. During the optimization process, the loop is free to move while the protein template is held fixed in its X-ray structure. To improve on the results of this model, in the present work we apply our optimization procedure to the physically more rigorous solvation model, the generalized Born with surface area (GB/SA) (together with the all-atom AMBER force field) as suggested by Still and co-workers (J. Phys. Chem. A 1997, 101, 3005). The six parameters of the GB/SA model, namely, P(1)-P(5) and the surface area parameter, sigma (programmed in the TINKER package) are reoptimized for a "training" group of nine loops, and a best-fit set is defined from the individual sets of optimized parameters. The best-fit set and Still's original set of parameters (where Lys, Arg, His, Glu, and Asp are charged or neutralized) were applied to the training group as well as to a "test" group of seven loops, and the energy gaps and the corresponding RMSD values were calculated. These GB/SA results based on the three sets of parameters have been found to be comparable; surprisingly, however, they are somewhat inferior (e.g, of larger energy gaps) to those obtained previously from the simplified model described above. We discuss recent results for loops obtained by other solvation models and potential directions for future studies.

Characterization of the nonregular regions of proteins by a contortion index

Nonstructured regions in proteins that provide the link between two regular structured regions play a significant role in maintaining the scaffold of the protein. Not only do they act as connectors between two regular secondary structural elements of proteins but they also provide the necessary turn or reversal in the polypeptide chain. This incorporates flexibility in the structure. Thus an understanding of the structural aspects of the nonregular regions is necessary to have a better insight into these features. We can assume the nonregular region to be a contorted polypeptide segment tethered by regular secondary structured regions at both ends. To describe the undulating nature of the nonregular regions, we introduce a parameter called the "contortion index." This index describes how tortuously the region is organized. Our analysis shows that the contortion index is related to other physicochemical parameters and can be used to characterize the nonregular regions of proteins.

V3: HIV's switch-hitter

The third variable region, V3, of the gp120 surface envelope glycoprotein is an approximately 35-residue-long, frequently glycosylated, highly variable, disulfide-bonded structure that has a major influence on HIV-1 tropism. Thus the sequence of V3, directly or indirectly, can determine which coreceptor (CCR5 or CXCR4) is used to trigger the fusion potential of the Env complex, and hence which cells the virus can infect. V3 also influences HIV-1's sensitivity to, and ability to escape from, entry inhibitors that are being developed as antiviral drugs. For some strains, V3 is a prominent target for HIV-1 neutralizing antibodies (NAbs); indeed, for many years it was considered to be the "principal neutralization determinant" (PND). Some efforts to use V3 as a vaccine target continue to this day, despite disappointing progress over more than a decade. Recent findings on the structure, function, antigenicity, and immunogenicity of V3 cast new doubts on the value of this vaccine approach. Here, we review recent advances in the understanding of V3 as a determinant of viral tropism, and discuss how this new knowledge may inform the development of HIV-1 drugs and vaccines.

Expanded turn conformations: characterization and sequence-structure correspondence in alpha-turns with implications in helix folding

Like the beta-turns, which are characterized by a limiting distance between residues two positions apart (i, i+3), a distance criterion (involving residues at positions i and i+4) is used here to identify alpha-turns from a database of known protein structures. At least 15 classes of alpha-turns have been enumerated based on the location in the phi,psi space of the three central residues (i+1 to i+3)-one of the major being the class AAA, where the residues occupy the conventional helical backbone torsion angles. However, moving towards the C-terminal end of the turn, there is a shift in the phi,psi angles towards more negative phi, such that the electrostatic repulsion between two consecutive carbonyl oxygen atoms is reduced. Except for the last position (i+4), there is not much similarity in residue composition at different positions of hydrogen and non-hydrogen bonded AAA turns. The presence or absence of Pro at i+1 position of alpha- and beta-turns has a bearing on whether the turn is hydrogen-bonded or without a hydrogen bond. In the tertiary structure, alpha-turns are more likely to be found in beta-hairpin loops. The residue composition at the beginning of the hydrogen bonded AAA alpha-turn has similarity with type I beta-turn and N-terminal positions of helices, but the last position matches with the C-terminal capping position of helices, suggesting that the existence of a "helix cap signal" at i+4 position prevents alpha-turns from growing into helices. Our results also provide new insights into alpha-helix nucleation and folding.