Homodimeric Protein Research Articles

The Integrity of the Intradimer Interface of the Hepatitis B Virus Capsid Protein Dimer Regulates Capsid Self-Assembly.

During the hepatitis B virus lifecycle, 120 copies of homodimeric capsid protein assemble around a copy of reverse transcriptase and viral RNA and go on to produce an infectious virion. Assembly needs to be tightly regulated by protein conformational change to ensure symmetry, fidelity, and reproducibility. Here, we show that structures at the intradimer interface regulate conformational changes at the distal interdimer interface and so regulate assembly. A pair of interacting charged residues, D78 from each monomer, conspicuously located at the top of a four-helix bundle that forms the intradimer interface, were mutated to serine to disrupt communication between the two monomers. The mutation slowed assembly and destabilized the dimer to thermal and chemical denaturation. Mutant dimers showed evidence of transient partial unfolding based on the appearance of new proteolytically sensitive sites. Though the mutant dimer was less stable, the resulting capsids were as stable as the wildtype, based on assembly and thermal denaturation studies. Cryo-EM image reconstructions of capsid indicated that the subunits adopted an "open" state more usually associated with a free dimer and that the spike tips were either disordered or highly flexible. Molecular dynamics simulations provide mechanistic explanations for these results, suggesting that D78 stabilizes helix 4a, which forms part of the intradimer interface, by capping its N-terminus and hydrogen-bonding to nearby residues, whereas the D78S mutation disrupts these interactions, leading to partial unwinding of helix 4a. This in turn weakens the connection from helix 4 and the intradimer interface to helix 5, which forms the interdimer interface.

Influence of protein (human galectin-3) design on aspects of lectin activity

The concept of biomedical significance of the functional pairing between tissue lectins and their glycoconjugate counterreceptors has reached the mainstream of research on the flow of biological information. A major challenge now is to identify the principles of structure–activity relationships that underlie specificity of recognition and the ensuing post-binding processes. Toward this end, we focus on a distinct feature on the side of the lectin, i.e. its architecture to present the carbohydrate recognition domain (CRD). Working with a multifunctional human lectin, i.e. galectin-3, as model, its CRD is used in protein engineering to build variants with different modular assembly. Hereby, it becomes possible to compare activity features of the natural design, i.e. CRD attached to an N-terminal tail, with those of homo- and heterodimers and the tail-free protein. Thermodynamics of binding disaccharides proved full activity of all proteins at very similar affinity. The following glycan array testing revealed maintained preferential contact formation with N-acetyllactosamine oligomers and histo-blood group ABH epitopes irrespective of variant design. The study of carbohydrate-inhibitable binding of the test panel disclosed up to qualitative cell-type-dependent differences in sections of fixed murine epididymis and especially jejunum. By probing topological aspects of binding, the susceptibility to inhibition by a tetravalent glycocluster was markedly different for the wild-type vs the homodimeric variant proteins. The results teach the salient lesson that protein design matters: the type of CRD presentation can have a profound bearing on whether basically suited oligosaccharides, which for example tested positively in an array, will become binding partners in situ. When lectin-glycoconjugate aggregates (lattices) are formed, their structural organization will depend on this parameter. Further testing (ga)lectin variants will thus be instrumental (i) to define the full range of impact of altering protein assembly and (ii) to explain why certain types of design have been favored during the course of evolution, besides opening biomedical perspectives for potential applications of the novel galectin forms.

In Vitro Selection of New DNA Aptamers for Human Vascular Endothelial Growth Factor 165

Two DNA aptamers that bind the heparin-binding domain (HBD) of the human vascular endothelial growth factor 165 (VEGF-165) have been previously reported. Although VEGF-165 is a homodimeric protein and the two aptamers have different sequences and secondary structures, the aptamers appear to occupy the same binding site and cannot form a 2 : 1 aptamer/protein complex, thus making them unsuitable for creating a higher-affinity dimeric DNA aptamer. This has motivated us to conduct a new in vitro selection experiment to search for new VEGF-165-binding DNA aptamers with different properties. We undertook a multistream selection strategy in which the concentration of VEGF-165 was varied significantly. We carried out 11 rounds of selection, and next-generation sequencing was conducted for every round in each stream. From comprehensive sequence analysis, we identified four classes of DNA aptamers, of which two were reported before, but two are new DNA aptamers. One of the new aptamers exhibits a unique property that has never been observed before: it is capable of forming the 2 : 1 aptamer/protein complex with VEGF-165. This work has expanded the repertoire of VEGF-165-binding DNA aptamers and creates a possibility to engineer a higher affinity homodimeric aptamer for VEGF-165.

ERH proteins: connecting RNA processing to tumorigenesis?

With the development of -omics approaches, the scientific community is now submerged by a wealth of information that can be used to analyze various parameters: the degree of protein sequence conservation, protein 3D structures as well as RNA and protein expression levels in various benign and tumor tissues, during organism development or upon exposure to chemicals such as endocrine disrupters. However, if such information can be used to identify genes with potentially important biological function, additional studies are needed to deeply characterize their cellular function in model organisms. Here, we discuss the case of such a gene: ERH, encoding a highly conserved homodimeric protein found in unicellular eukaryotes, plants and metazoan, of yet unknown biological function, which might be linked to mRNA metabolism and that is emerging as important for cell migration and metastasis.

Identification of potential anti-leishmanial agents using computational investigation and biological evaluation against trypanothione reductase

Trypanothione reductase of Leishmania donovani is a flavin adenine dinucleotide containing homodimeric protein essential for parasite survival. The flavoenzyme utilizes nicotinamide adenine dinucleotide phosphate in the reaction to convert oxidized trypanothione to reduced trypanothione which is further used up by tryparedoxin/tryparedoxin peroxidase system to neutralize the reactive oxygen species generated by the macrophages. Some of the drugs previously reported against the disease include sodium stibogluconate, miltefosine and amphotericin B. However, due to the resistance and toxicity problem associated with these molecules, there is an urgent need to develop new drugs against L. donovani. Trypanothione reductase of L. donovani is one such essential target whose inhibition could lead to a decline in parasite growth. In this work, we have performed a computational studies using Maybridge library of chemical compounds to identify potential inhibitors of Trypanothione reductase of L. donovani. Structure-based virtual screening method in combination with molecular docking was employed to identify and prioritize 30 compounds which were further subjected to molecular dynamics simulation. Ten compounds which showed stable ligand root-mean-square deviation plot, c-alpha backbone and root-mean-square fluctuation were considered for trypanothione reductase inhibition assay and subsequent inhibition studies of parasite growth. Enzyme inhibition assay resulted in shortlisting of four compounds that were found to inhibit Trypanothione reductase of L. donovani. Subsequently, the anti-leishmanial screening highlighted one compound as the potential anti-leishmanial agent, with IC50 value of 15.2 µM, that can be further optimised with medicinal chemistry efforts to improve its activity. Communicated by Ramaswamy H. Sarma

A DNA Vaccine That Encodes an Antigen-Presenting Cell-Specific Heterodimeric Protein Protects against Cancer and Influenza

Immunogenicity of DNA vaccines can be increased by constructing the DNA in such a way that it encodes secreted homodimeric fusion proteins that target antigen-presenting cells (APCs). In this study, we have developed novel APC-targeting vaccine molecules with an increased flexibility due to introduction of a heterodimerization motif. The heterodimeric proteins permit four different fusions within a single molecule, thus allowing expression of two different APC-targeting moieties and two different antigens. Two types of heterodimeric fusion proteins were developed that employed either the ACID/BASE or the Barnase/Barstar motifs, respectively. The ACID/BASE heterodimeric vaccines conferred protection against challenges with either influenza virus or tumor cells in separate preclinical models. The ACID/BASE motif was flexible since a large number of different targeting moieties and antigens could be introduced with maintenance of specificity, antigenicity, and secretion. APC-targeting ACID/BASE vaccines expressing two different antigens induced antibody and T cell responses against either of the two antigens. Heterodimeric ACID/BASE DNA vaccines were of approximately the same potency as previously reported homodimeric DNA vaccines. The flexibility and potency of the ACID/BASE format suggest that it could be a useful platform for DNA vaccines that encode APC-targeting fusion proteins.

A Novel Ribonuclease with HIV-1 Reverse Transcriptase Inhibitory Activity Purified from the Field Blewit Mushroom Lepista personata (Agaricomycetes).

A ribonuclease was purified from dry fruiting bodies of the wild edible mushroom Lepista personata (LPR) to 259-fold with a specific activity of 280 U/mg. The purification protocol involved ion-exchange chromatography on DEAE-cellulose, CM-cellulose and SP-sepharose, followed by size exclusion chromatography on Superdex 75. LPR is a homodimeric protein with a molecular weight of 27.8 kDa as determined by SDS-PAGE and by gel filtration. Three inner peptide sequences for LPR were obtained by LC-MS-MS analysis. It demonstrated the optimum pH of 4.0 and temperature optimum of 60°C. The specificity ribonuclease potencies order toward polyhomoribonucleotides was poly C > poly A > poly G > poly U. The ribonuclease inhibited HIV-1 reverse transcriptase with an IC50 of 0.53 μM.

Comprehensive in silico screening and molecular dynamics studies of missense mutations in Sjogren-Larsson syndrome associated with the ALDH3A2 gene.

Sjögren-Larsson syndrome (SLS) is an autoimmune disorder inherited in an autosomal recessive pattern. To date, 80 missense mutations have been identified in association with the Aldehyde Dehydrogenase 3 Family Member A2 (ALDH3A2) gene causing SLS. Disruption of the function of ALDH3A2 leads to excessive accumulation of fat in the cells, which interferes with the normal function of protective membranes or materials that are necessary for the body to function normally. We retrieved 54 missense mutations in the ALDH3A2 from the OMIM, UniProt, dbSNP, and HGMD databases that are known to cause SLS. These mutations were examined with various in silico stability tools, which predicted that the mutations p.S308N and p.R423H that are located at the protein-protein interaction domains are the most destabilizing. Furthermore, to determine the atomistic-level differences within the protein-protein interactions owing to mutations, we performed macromolecular simulation (MMS) using GROMACS to validate the motion patterns and dynamic behavior of the biological system. We found that both mutations (p.S380N and p.R423H) had significant effects on the protein-protein interaction and disrupted the dimeric interactions. The computational pipeline provided in this study helps to elucidate the potential structural and functional differences between the ALDH3A2 native and mutant homodimeric proteins, and will pave the way for drug discovery against specific targets in the SLS patients.

Catalytic Cycle of Glycoside Hydrolase BglX from Pseudomonas aeruginosa and Its Implications for Biofilm Formation.

BglX is a heretofore uncharacterized periplasmic glycoside hydrolase (GH) of the human pathogen Pseudomonas aeruginosa. X-ray analysis identifies it as a protein homodimer. The two active sites of the homodimer comprise catalytic residues provided by each monomer. This arrangement is seen in <2% of the hydrolases of known structure. In vitro substrate profiling shows BglX is a catalyst for β-(1→2) and β-(1→3) saccharide hydrolysis. Saccharides with β-(1→4) or β-(1→6) bonds, and the β-(1→4) muropeptides from the cell-wall peptidoglycan, are not substrates. Additional structural insights from X-ray analysis (including structures of a mutant enzyme-derived Michaelis complex, two transition-state mimetics, and two enzyme-product complexes) enabled the comprehensive description of BglX catalysis. The half-chair (4H3) conformation of the transition-state oxocarbenium species, the approach of the hydrolytic water molecule to the oxocarbenium species, and the stepwise release of the two reaction products were also visualized. The substrate pattern for BglX aligns with the [β-(1→2)-Glc]x and [β-(1→3)-Glc]x periplasmic osmoregulated periplasmic glucans, and possibly with the Psl exopolysaccharides, of P. aeruginosa. Both polysaccharides are implicated in biofilm formation. Accordingly, we show that inactivation of the bglX gene of P. aeruginosa PAO1 attenuates biofilm formation.

Anti Mullerian hormone levels in women of central Madhya Pradesh- A retrospective analysis

Introduction: Anti Mullerian Hormone(AMH) is homodimeric protein that is classified under the superfamily of transforming growth factor beta(TGF-b ). The granulose cells of the primary, preantral and small antral follicles of the ovary express this hormone. AMH levels remain fairly constant throughout the cycle which makes it a favourable and suitable marker for determining the ovarian reserve in women of reproductive age. The present retrospective study was carried out with the aim of studying the AMH levels in women upto 50 years of age coming to our diagnostic centre for AMH screening for different reasons Materials and Methods: This was a retrospective study conducted between January to August 2019 in women upto 50 years of age coming to our diagnostic centre for AMH screening test. They were divided into further classified as 10.0 ng /ml .The mean and median values were calculated for each group. Results: AMH levels below 1 ng/ml were observed in 102/400(25.5%) patients. AMH levels between 1-4 ng/ml was observed in 171/400 (42.75%) women. High AMH values between 4-8 ng/ml were observed in 90/400(22.5%) women out of which66/90(73.33%)were in 21-30 years age range, followed by 18.89% in 31-40 years, 6.67% below 20 years and 1.11?ove 40 years of age Conclusion: Since AMH levels do not vary with the days of the menstrual cycle in women, it is one of the best markers available in the current scenario for assessing the ovarian reserve in women of reproductive age group. Our study emphasizes the observation that PCOS is a growing health problem in women of reproductive age due to several factors like sedentary life style, lack of exercise and unhealthy eating habits which can be easily managed by lifestyle changes. Keywords: Anti Mullerian hormone, LH, FSH, PCOS.

Unraveling Determinants of Affinity Enhancement in Dimeric Aptamers for a Dimeric Protein

High-affinity aptamers can be derived de novo by using stringent conditions in SELEX (Systematic Evolution of Ligands by EXponential enrichment) experiments or can be engineered post SELEX via dimerization of selected aptamers. Using electrophoretic mobility shift assays, we studied a series of heterodimeric and homodimeric aptamers, constructed from two DNA aptamers with distinct primary sequences and secondary structures, previously isolated for VEGF-165, a homodimeric protein. We investigated four factors envisaged to impact the affinity of a dimeric aptamer to a dimeric protein: (1) length of the linker between two aptamer domains, (2) linking orientation, (3) binding-site compatibility of two component aptamers in a heterodimeric aptamer, and (4) steric acceptability of the two identical aptamers in a homodimeric aptamer. All heterodimeric aptamers for VEGF-165 were found to exhibit monomeric aptamer-like affinity and the lack of affinity enhancement was attributed to binding-site overlap by the constituent aptamers. The best homodimeric aptamer showed 2.8-fold better affinity than its monomeric unit (Kd = 13.6 ± 2.7 nM compared to 37.9 ± 14 nM), however the barrier to further affinity enhancement was ascribed to steric interference of the constituent aptamers. Our findings point to the need to consider the issues of binding-site compatibility and spatial requirement of aptamers for the development of dimeric aptamers capable of bivalent recognition. Thus, determinants highlighted herein should be assessed in future multimerization efforts.

Analysis of Heterodimeric "Mutual Synergistic Folding"-Complexes.

Several intrinsically disordered proteins (IDPs) are capable to adopt stable structures without interacting with a folded partner. When the folding of all interacting partners happens at the same time, coupled with the interaction in a synergistic manner, the process is called Mutual Synergistic Folding (MSF). These complexes represent a discrete subset of IDPs. Recently, we collected information on their complexes and created the MFIB (Mutual Folding Induced by Binding) database. In a previous study, we compared homodimeric MSF complexes with homodimeric and monomeric globular proteins with similar amino acid sequence lengths. We concluded that MSF homodimers, compared to globular homodimeric proteins, have a greater solvent accessible main-chain surface area on the contact surface of the subunits, which becomes buried during dimerization. The main driving force of the folding is the mutual shielding of the water-accessible backbones, but the formation of further intermolecular interactions can also be relevant. In this paper, we will report analyses of heterodimeric MSF complexes. Our results indicate that the amino acid composition of the heterodimeric MSF monomer subunits slightly diverges from globular monomer proteins, while after dimerization, the amino acid composition of the overall MSF complexes becomes more similar to overall amino acid compositions of globular complexes. We found that inter-subunit interactions are strengthened, and additionally to the shielding of the solvent accessible backbone, other factors might play an important role in the stabilization of the heterodimeric structures, likewise energy gain resulting from the interaction of the two subunits with different amino acid compositions. We suggest that the shielding of the β-sheet backbones and the formation of a buried structural core along with the general strengthening of inter-subunit interactions together could be the driving forces of MSF protein structural ordering upon dimerization.

Oligomeric state of the ZIKV E protein defines protective immune responses

The current leading Zika vaccine candidates in clinical testing are based on live or killed virus platforms, which have safety issues, especially in pregnant women. Zika subunit vaccines, however, have shown poor performance in preclinical studies, most likely because the antigens tested do not display critical quaternary structure epitopes present on Zika E protein homodimers that cover the surface of the virus. Here, we produce stable recombinant E protein homodimers that are recognized by strongly neutralizing Zika specific monoclonal antibodies. In mice, the dimeric antigen stimulate strongly neutralizing antibodies that target epitopes that are similar to epitopes recognized by human antibodies following natural Zika virus infection. The monomer antigen stimulates low levels of E-domain III targeting neutralizing antibodies. In a Zika challenge model, only E dimer antigen stimulates protective antibodies, not the monomer. These results highlight the importance of mimicking the highly structured flavivirus surface when designing subunit vaccines.

Differential modulation of energy landscapes of cyclic AMP receptor protein (CRP) as a regulatory mechanism for class II CRP-dependent promoters

The Escherichia coli cAMP receptor protein, CRP, is a homodimeric global transcription activator that employs multiple mechanisms to modulate the expression of hundreds of genes. These mechanisms require different interfacial interactions among CRP, RNA, and DNA of varying sequences. The involvement of such a multiplicity of interfaces requires a tight control to ensure the desired phenotype. CRP-dependent promoters can be grouped into three classes. For decades scientists in the field have been puzzled over the differences in mechanisms between class I and II promoters. Using a new crystal structure, IR spectroscopy, and computational analysis, we defined the energy landscapes of WT and 14 mutated CRPs to determine how a homodimeric protein can distinguish nonpalindromic DNA sequences and facilitate communication between residues located in three different activation regions (AR) in CRP that are ∼30 Å apart. We showed that each mutation imparts differential effects on stability among the subunits and domains in CRP. Consequently, the energetic landscapes of subunits and domains are different, and CRP is asymmetric. Hence, the same mutation can exert different effects on ARs in class I or II promoters. The effect of a mutation is transmitted through a network by long-distance communication not necessarily relying on physical contacts between adjacent residues. The mechanism is simply the sum of the consequences of modulating the synchrony of dynamic motions of residues at a distance, leading to differential effects on ARs in different subunits. The computational analysis is applicable to any system and potentially with predictive capability.

Negative Cooperativity in NAD(P)H Quinone Oxidoreductase 1 (NQO1).

NAD(P)H quinone oxidoreductase-1 (NQO1) is a homodimeric protein that acts as a detoxifying enzyme or as a chaperone protein. Dicourmarol interacts with NQO1 at the NAD(P)H binding site and can both inhibit enzyme activity and modulate the interaction of NQO1 with other proteins. We show that the binding of dicoumarol and related compounds to NQO1 generates negative cooperativity between the monomers. This does not occur in the presence of the reducing cofactor, NAD(P)H, alone. Alteration of Gly150 (but not Gly149 or Gly174) abolished the dicoumarol-induced negative cooperativity. Analysis of the dynamics of NQO1 with the Gaussian network model indicates a high degree of collective motion by monomers and domains within NQO1. Ligand binding is predicted to alter NQO1 dynamics both proximal to the ligand binding site and remotely, close to the second binding site. Thus, drug-induced modulation of protein motion might contribute to the biological effects of putative inhibitors of NQO1.

Identification and Molecular Characterization of Superoxide Dismutases Isolated From A Scuticociliate Parasite: Physiological Role in Oxidative Stress

Philasterides dicentrarchi is a free-living microaerophilic scuticociliate that can become a facultative parasite and cause a serious parasitic disease in farmed fish. Both the free-living and parasitic forms of this scuticociliate are exposed to oxidative stress associated with environmental factors and the host immune system. The reactive oxygen species (ROS) generated by the host are neutralized by the ciliate by means of antioxidant defences. In this study we aimed to identify metalloenzymes with superoxide dismutase (SOD) activity capable of inactivating the superoxide anion (•O2−) generated during induction of oxidative stress. P. dicentrarchi possesses the three characteristic types of SOD isoenzymes in eukaryotes: copper/zinc-SOD, manganese-SOD and iron-SOD. The Cu/Zn-SOD isoenzymes comprise three types of homodimeric proteins (CSD1-3) of molecular weight (MW) 34–44 kDa and with very different AA sequences. All Cu/Zn-SODs are sensitive to NaCN, located in the cytosol and in the alveolar sacs, and one of them (CSD2) is extracellular. Mn- and Fe-SOD transcripts encode homodimeric proteins (MSD and FSD, respectively) in their native state: a) MSD (MW 50 kDa) is insensitive to H2O2 and NaN3 and is located in the mitochondria; and b) FSD (MW 60 kDa) is sensitive to H2O2, NaN3 and the polyphenol trans-resveratrol and is located extracellularly. Expression of SOD isoenzymes increases when •O2− is induced by ultraviolet (UV) irradiation, and the increase is proportional to the dose of energy applied, indicating that these enzymes are actively involved in cellular protection against oxidative stress.

Crystal structure and catalytic mechanism of the essential m1G37 tRNA methyltransferase TrmD from Pseudomonas aeruginosa

The tRNA (m1G37) methyltransferase TrmD catalyzes m1G formation at position 37 in many tRNA isoacceptors and is essential in most bacteria, which positions it as a target for antibiotic development. In spite of its crucial role, little is known about TrmD in Pseudomonas aeruginosa (PaTrmD), an important human pathogen. Here we present detailed structural, substrate, and kinetic properties of PaTrmD. The mass spectrometric analysis confirmed the G36G37-containing tRNAs Leu(GAG), Leu(CAG), Leu(UAG), Pro(GGG), Pro(UGG), Pro(CGG), and His(GUG) as PaTrmD substrates. Analysis of steady-state kinetics with S-adenosyl-l-methionine (SAM) and tRNALeu(GAG) showed that PaTrmD catalyzes the two-substrate reaction by way of a ternary complex, while isothermal titration calorimetry revealed that SAM and tRNALeu(GAG) bind to PaTrmD independently, each with a dissociation constant of 14 ± 3 µM. Inhibition by the SAM analog sinefungin was competitive with respect to SAM (Ki = 0.41 ± 0.07 µM) and uncompetitive for tRNA (Ki = 6.4 ± 0.8 µM). A set of crystal structures of the homodimeric PaTrmD protein bound to SAM and sinefungin provide the molecular basis for enzyme competitive inhibition and identify the location of the bound divalent ion. These results provide insights into PaTrmD as a potential target for the development of antibiotics.

Structural and Functional Properties of the Capsid Protein of Dengue and Related Flavivirus.

Dengue, West Nile and Zika, closely related viruses of the Flaviviridae family, are an increasing global threat, due to the expansion of their mosquito vectors. They present a very similar viral particle with an outer lipid bilayer containing two viral proteins and, within it, the nucleocapsid core. This core is composed by the viral RNA complexed with multiple copies of the capsid protein, a crucial structural protein that mediates not only viral assembly, but also encapsidation, by interacting with host lipid systems. The capsid is a homodimeric protein that contains a disordered N-terminal region, an intermediate flexible fold section and a very stable conserved fold region. Since a better understanding of its structure can give light into its biological activity, here, first, we compared and analyzed relevant mosquito-borne Flavivirus capsid protein sequences and their predicted structures. Then, we studied the alternative conformations enabled by the N-terminal region. Finally, using dengue virus capsid protein as main model, we correlated the protein size, thermal stability and function with its structure/dynamics features. The findings suggest that the capsid protein interaction with host lipid systems leads to minor allosteric changes that may modulate the specific binding of the protein to the viral RNA. Such mechanism can be targeted in future drug development strategies, namely by using improved versions of pep14-23, a dengue virus capsid protein peptide inhibitor, previously developed by us. Such knowledge can yield promising advances against Zika, dengue and closely related Flavivirus.

Positional Distributions of the Tethered Modules in Nitric Oxide Synthase: Monte Carlo Calculations and Pulsed EPR Measurements.

The nitric oxide synthase (NOS) enzyme consists of multiple domains connected by flexible random coil tethers. In a catalytic cycle, the NOS domains move within the limits determined by the length and flexibility of the interdomain tethers and form docking complexes with each other. This process represents a key component of the electron transport from the flavin adenine dinucleotide/reduced nicotinamide adenine dinucleotide phosphate binding domain to the catalytic heme centers located in the oxygenase domain. Studying the conformational behavior of NOS is therefore imperative for a full understanding of the overall catalytic mechanism. In this work, we have investigated the equilibrium positional distributions of the NOS domains and the bound calmodulin (CaM) by using Monte Carlo calculations of the NOS conformations. As a main experimental reference, we have used the magnetic dipole interaction between a bifunctional spin label attached to T34C/S38C mutant CaM and the NOS heme centers, which was measured by pulsed electron paramagnetic resonance. In general, the calculations of the conformational distributions allow one to determine the range and statistics of positions occupied by the tethered protein domains, assess the crowding effect of the multiple domains on each other, evaluate the accessibility of various potential domain docking sites, and estimate the interaction energies required to achieve target populations of the docked states. In the particular application described here, we have established the specific mechanisms by which the bound CaM facilitates the flavin mononucleotide (FMN)/heme interdomain docking in NOS. We have also shown that the intersubunit FMN/heme domain docking and electron transfer in the homodimeric NOS protein are dictated by the existing structural makeup of the protein. Finally, from comparison of the calculated and experimental docking probabilities, the characteristic stabilization energies for the CaM/heme domain and the FMN domain/heme domain docking complexes have been estimated as -4.5kT and -10.5kT, respectively.

A Single Reactive Noncanonical Amino Acid Is Able to Dramatically Stabilize Protein Structure.

A p-isothiocyanate phenylalanine mutant of the homodimeric protein homoserine o-succinyltransferase (MetA) was isolated in a temperature dependent selection from a library of metA mutants containing noncanonical amino acids. This mutant protein has a dramatic increase of 24 °C in thermal stability compared to the wild type protein. Peptide mapping experiments revealed that the isothiocyanate group forms a thiourea cross-link to the N terminal proline of the other monomer, despite the two positions being >30 Å apart in the X-ray crystal structure of the wild type protein. These results show that an expanded set of building blocks beyond the canonical 20 amino acids can lead to significant changes in the properties of proteins.