Computational and experimental insights into the mechanism of substrate recognition and feedback inhibition of protoporphyrinogen oxidase

A New V361A Mutation in Amaranthus palmeri PPX2 Associated with PPO-Inhibiting Herbicide Resistance

Weeds resistant to PPO-inhibiting herbicides threaten the profitability of crop producers relying on this chemistry. In Amaranthus palmeri, mutations at G210 (∆G210) and R128 (R128G/M) of the PPX2 gene were reported to confer PPO-inhibitor resistance. Here, A. palmeri samples from nine states in America, having survived a field application of a PPO-inhibitor, were genotyped to determine the prevalence of these mutations. Less than 5% of the 1828 A. palmeri plants screened contained the ∆G210 mutation. Of the plants lacking ∆G210, a R128 substitution was only found in a single plant. An A. palmeri population from Alabama without mutations at G210 or R128 had a resistance ratio of 3.1 to 3.5 for fomesafen. Of the candidate PPX2 mutations identified in this population, only V361A conferred resistance to lactofen and fomesafen in a transformed bacterial strain. This is the first report of the V361A substitution of PPX2 conferred PPO-inhibiting herbicide resistance in any plant species. Future molecular screens of PPO-inhibitor resistance in A. palmeri and other species should encompass the V361A mutation of PPX2 to avoid false-negative results.

Identifying novel selective PPO inhibitors through structure-based virtual screening and bio-evaluation

Protoporphyrinogen oxidase (PPO) is a key enzyme in chlorophyll and heme biosynthesis, and the development of its inhibitors is of great importance both in the pharmaceutical and pesticide industries. However, the currently developed PPO inhibitors have insignificant bio-selectivity and have a serious impact on non-target organisms. In this study, a docking-based virtual screening approach combined with bio-activity testing was used to obtain novel selective inhibitors of PPO. The results of the bio-activity test showed that thirteen compounds showed 10-fold selectivity over human PPO. And the best selective compound, ZINC70338, has a K _i value of 2.21 μM for Nicotiana tabacum PPO and >113-fold selectivity for human PPO. The selectivity mechanism of ZINC70338 in different species of PPO was then analyzed by molecular dynamics simulations to provide a design basis and theoretical guidance for the design of novel selective inhibitors.

Discovery of novel phenoxypyridine as promising protoporphyrinogen IX oxidase inhibitors

Protoporphyrinogen oxidase (PPO, EC 1.3.3.4) is a significant target for the discovery of novel bleaching herbicides. Starting from the active fragments of several known commercial herbicides, a series of PPO inhibitors with diphenyl ether scaffolds were designed and synthesized by substructure splicing and bioisosterism methods. The greenhouse herbicidal activity and the PPO inhibitory activity in vitro were measured. The results showed that the novel synthesized compounds have good PPO inhibitory activity, and the IC₅₀ value against corn PPO ranges from 0.032 ± 0.008 mg/L to 3.245 ± 0.247 mg/L. Among all target compounds, compound P2 showed the best herbicidal activity, with a half inhibitory concentration (IC₅₀) of 0.032 ± 0.008 mg/L. In addition, the molecular docking results showed that the benzene ring part of compound P2 can form a π-π stacking with PHE-392, and the trifluoromethyl group and ARG-98 form two hydrogen bonds. Crop safety experiments and cumulative concentration analysis experiments indicated that compound P2 can be used for weed control in rice, wheat, soybean and corn. Therefore, compound P2 can be selected to develop potential lead compounds for novel PPO inhibitors.

Can double PPO mutations exist in the same allele and are such mutants functional?

BACKGROUND

Resistance to protoporphyrinogen oxidase (PPO)-inhibiting herbicides is endowed primarily by target-site mutations at the PPX2 gene that compromise binding of the herbicide to the catalytic domain. In Amaranthus spp. PPX2, the most prevalent target mutations are deletion of the G210 codon, and the R128G and G339A substitutions. These mutations strongly affect the dynamic of the PPO2 binding pocket, resulting in reduced affinity with the ligand. Here we investigated the likelihood of co-occurrence of the most widespread target site mutations in the same PPX2 allele.

RESULTS

Plants carrying R128G+/+ ΔG210+/-, where + indicates presence of the mutation, were crossed with each other. The PPX2 of the offspring was subjected to pyrosequencing and E. coli-based Sanger sequencing to determine mutation frequencies and allele co-occurrence. The data show that R128G ΔG210 can occur in one allele only; the second allele carries only one mutation. Double mutation in both alleles is less likely because of significant loss of enzyme activity. The segregation of offspring populations derived from a cross between heterozygous plants carrying ΔG210 G399A also showed no co-occurrence in the same allele. The offspring exhibited the expected mutation distribution patterns with few exceptions.

CONCLUSIONS

Homozygous double-mutants are not physiologically viable. Double-mutant plants can only exist in a heterozygous state. Alternatively, if two mutations are detected in one plant, each mutation would occur in a separate allele. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Design, Synthesis, and Herbicidal Activity of Diphenyl Ether Derivatives Containing a Five-Membered Heterocycle

Protoporphyrinogen oxidase (PPO, EC 1.3.3.4) is an important target for discovering novel herbicides, and it causes bleaching symptoms by inhibiting the synthesis of chlorophyll and heme. In this study, the active fragments of several commercial herbicides were joined by substructure splicing and bioisosterism, and a series of novel diphenyl ether derivatives containing five-membered heterocycles were synthesized. The greenhouse herbicidal activity and the PPO inhibitory activity in vitro were discussed in detail. The results showed that most compounds had good PPO inhibitory activity, and target compounds containing trifluoromethyl groups tended to have higher activity. Among them, compound G4 showed the best inhibitory activity, with a half-maximal inhibitory concentration (IC₅₀) of 0.0468 μmol/L, which was approximately 3 times better than that of oxyfluorfen (IC₅₀ = 0.150 μmol/L). In addition, molecular docking indicated that compound G4 formed obvious π-π stacking interactions and hydrogen bond interactions with PHE-392 and ARG-98, respectively. Remarkably, compound G4 had good safety for corn, wheat, rice, and soybean, and the cumulative concentration in crops was lower than that of oxyfluorfen. Therefore, compound G4 can be used to develop potential lead compounds for novel PPO inhibitors.

Functional PPO2 mutations: co-occurrence in one plant or the same ppo2 allele of herbicide-resistant Amaranthus palmeri in the US mid-south

BACKGROUND

Protoporphyrinogen IX oxidase 2 (PPO2) inhibitors are important for the management of glyphosate- and acetolactate synthase-resistant Palmer amaranth [Amaranthus palmeri (S.) Wats.]. The evolving resistance to PPO inhibitors is of great concern. We surveyed the evolution of resistance to fomesafen in the US Mid-south and determined its correlation with the known functional PPO2 target-site mutations (TSM).

RESULTS

The 167 accessions analyzed were grouped into five categories, four resistant (147) and one susceptible (20). Arkansas accessions constituted 100% of the susceptible group while the Missouri accessions comprised 60% of the most resistant category. The majority of Mississippi accessions (88%) clustered in the high-survival-high-injury category, manifesting an early-stage resistance evolution. One hundred and fifteen accessions were genotyped for four known TSMs; 74% of accessions carried at least one TSM. The most common single TSM was ΔG210 (18% of accessions) and the predominant double mutation was ΔG210 + G399A (17%). Other mutations are likely less favorable, hence are rare. All TSMs were detected in three accessions. Further examination revealed that 9 and two individuals carried G399A + G210 and G399A + R128G TSM in the same allele, respectively. The existence of these combinations is supported by molecular modeling.

CONCLUSIONS

Resistance to PPO inhibitors is widespread across the Mid-southern USA. Highly resistant field populations have plants with multiple mutations. G399A is the most prone to co-occur with other ppo2 mutations in the same allele. Mutation at R128 in the configuration of the PPO2 catalytic domain restrains the co-occurrence of R128G with ΔG210, making ΔG210 + G399A the most plausible, tolerable functional mutation combination to co-occur in the same ppo2 allele.

Conformation of the Intermediates in the Reaction Catalyzed by Protoporphyrinogen Oxidase: An In Silico Analysis

Protoporphyrinogen oxidase (PPO) is a critical enzyme across life as the last common step in the synthesis of many metalloporphyrins. The reaction mechanism of PPO was assessed in silico and the unstructured loop near the binding pocket was investigated. The substrate, intermediates, and product were docked in the catalytic domain of PPO using a modified Autodock method, introducing flexibility in the macrocycles. Sixteen PPO protein sequences across phyla were aligned and analyzed with Phyre2 and ProteinPredict to study the unstructured loop from residue 204–210 in the H. sapiens structure. Docking of the substrate, intermediates, and product all resulted in negative binding energies, though the substrate had a lower energy than the others by 40%. The α-H of C10 was found to be 1.4 angstroms closer to FAD than the β-H, explaining previous reports of the reaction occurring on the meso face of the substrate. A lack of homology in sequence or length in the unstructured loop indicates a lack of function for the protein reaction. This docking study supports a reaction mechanism proposed previously whereby all hydride abstractions occur on the C10 of the tetrapyrrole followed by tautomeric rearrangement to prepare the intermediate for the next reaction.

Design, Synthesis, and Molecular Mechanism Studies of N-Phenylisoxazoline-thiadiazolo[3,4-a]pyridazine Hybrids as Protoporphyrinogen IX Oxidase Inhibitors

Protoporphyrinogen oxidase (PPO, EC 1.3.3.4) is an important target for green agrochemical discovery. Herein, a novel N-phenylisoxazoline-thiadiazolo[3,4-a]pyridazine herbicidal active scaffold was designed by the scaffold hybridization strategy. Systematic structural optimization enabled the discovery of a series of derivatives with excellent weed control at 9.375-150 g ai/ha by the post-emergent application. Some derivatives exhibited improved Nicotiana tabacum PPO (NtPPO)-inhibitory activity than fluthiacet-methyl. Of these, 2b, with K_i = 21.8 nM, displayed higher weed control than fluthiacet-methyl at the rate of 12-75 g ai/ha, and selective to maize at 75 g ai/ha. In planta, 2b was converted into a bioactive metabolite 5 (K_i = 4.6 nM), which exhibited 4.6-fold more potency than 2b in inhibiting the activity of NtPPO. Molecular dynamics simulation explained that 5 formed stronger π-π interaction with Phe392 than that of 2b. This work not only provides a promising lead compound for weed control in maize fields but is also helpful to understand the molecular mechanism and basis of the designed hybrids.

Review on Structures of Pesticide Targets

Molecular targets play important roles in agrochemical discovery. Numerous pesticides target the key proteins in pathogens, insect, or plants. Investigating ligand-binding pockets and/or active sites in the proteins' structures is usually the first step in designing new green pesticides. Thus, molecular target structures are extremely important for the discovery and development of such pesticides. In this manuscript, we present a review of the molecular target structures, including those of antiviral, fungicidal, bactericidal, insecticidal, herbicidal, and plant growth-regulator targets, currently used in agrochemical research. The data will be helpful in pesticide design and the discovery of new green pesticides.

A novel mutation A212T in chloroplast Protoporphyrinogen oxidase (PPO1) confers resistance to PPO inhibitor Oxadiazon in Eleusine indica

BACKGROUND

Protoporphyrinogen oxidase (PPO) with two isoforms, chloroplast-targeted (PPO1) and mitochondrial-targeted (PPO2), catalyzes a step in the biosynthesis of chlorophyll and heme. PPO1 and PPO2 are herbicide target sites of PPO-inhibiting herbicides. Target-site mutations conferring resistance to PPO inhibitors have all thus far been in PPO2. Oxadiazon is a unique PPO inhibitor utilized for preemergence Eleusine indica control. In this research, we evaluated the response of two previously confirmed oxadiazon-resistant and susceptible E. indica biotypes to other PPO inhibitors and identified the resistance mechanism in two oxadiazon-resistant E. indica biotypes.

RESULTS

Two E. indica biotypes were resistant to oxadiazon, but not to other structurally unrelated PPO inhibitors, such as lactofen, flumioxazin and sulfentrazone. A novel mutation A212T was identified in the chloroplast-targeted PPO1, conferring resistance to oxadiazon in a heterologous expression system. Computational structural modeling provided a mechanistic explanation for reduced herbicide binding to the variant protein: the presence of a methyl group of threonine 212 changes the PPO1 active site and produces repulsive electrostatic interactions that repel oxadiazon from the binding pocket.

CONCLUSION

The novel A212T mutation in PPO1 conferring resistance specifically to PPO inhibitor oxadiazon was characterized. This is the first evidence of the direct role of PPO1 in the PPO mode of action, and the first evidence of evolved resistance in PPO1. © 2019 Society of Chemical Industry.

Investigating target-site resistance mechanism to the PPO-inhibiting herbicide fomesafen in waterhemp and interspecific hybridization of Amaranthus species using next generation sequencing

BACKGROUND

Waterhemp (Amaranthus tuberculatus (Moq.) J. D. Sauer) is one of the most pernicious weeds in cropping systems of the USA due to its evolved resistance against several herbicide sites-of-action, including protoporphyrinogen oxidase inhibitors (PPO-R). Currently, the only source of PPO-R documented in waterhemp is ΔG210 of PPX2. Gene flow may not only lead to a transfer of herbicide-resistant alleles, but also produce a hybrid genotype more competitively fit than one or both parents. However, investigating gene flow of Amaranthus species has been of interest in the past two decades with limited evidence.

RESULTS

Here, a high-throughput MiSeq amplicon sequencing method was used to investigate alterations of the PPX2 gene in 146 PPO-R waterhemp populations across five Midwest states of the USA. Five R128 codons of PPX2, novel to waterhemp, were found including AGG (R), GGA (G), GGG (G), AAA (K) and ATA (I). R128G, R128I, and R128K were found in 11, 3, and 2 populations, respectively. R128G and R128I, but not R128K, conferred fomesafen resistance in a bacterial system. Sequence alignment of the R128 region of PPX2 identified a tumble pigweed (Amaranthus albus)-type and Palmer amaranth (Amaranthus palmeri)-type PPX2 allele to be present and widespread in the surveyed waterhemp populations, thus providing strong evidence of gene flow between Amaranthus species.

CONCLUSION

Using a next-generation sequencing method, we identified two PPO target-site mutations R128G/I novel to waterhemp and provided evidence of gene flow of Amaranthus species in a large group of screened waterhemp populations from five Midwest states of the USA. © 2019 Society of Chemical Industry.

Characterization of HemY-type protoporphyrinogen IX oxidase genes from cyanobacteria and their functioning in transgenic Arabidopsis

We investigated the functions of two cyanobacterial HemY protoporphyrinogen IX oxidase (PPO) genes with in vitro and in vivo assays and evaluated their applicability as resistance traits to PPO-inhibiting herbicides. We isolated HemY-type protoporphyrinogen IX oxidase (PPO) genes from cyanobacteria, OnPPO gene from Oscillatoria nigro-viridis PCC7112 and HaPPO gene from Halothece sp. PCC7418. The alignment of amino acid sequences as well as phylogenetic analyses conducted showed that OnPPO and HaPPO are classified as HemY-type PPO and are more closely related to plastidic PPOs than to mitochondrial PPOs. The PPO-deficient Escherichia coli BT3 strain, which requires heme supplementation, could obtain normal growth in the absence of heme supplementation when complemented with OnPPO and HaPPO. The enzyme assays of OnPPO, HaPPO, and Arabidopsis thaliana PPO1 (AtPPO1) proteins each revealed different kinetic properties in terms of catalytic efficiency, substrate affinity, and the degree of inhibition by PPO inhibitors. In particular, the catalytic efficiencies (k_cat/Km) of OnPPO and HaPPO were approximately twofold higher than that of AtPPO1. The elution profiles of all three PPOs, acquired by size-exclusion chromatography, showed only a single peak with a molecular weight of approximately 52-54 kDa, which corresponds to a monomeric form. Moreover, functional complementation with OnPPO and HaPPO in AtPPO1-silenced Arabidopsis resulted in restored growth, whereas AtPPO1-silenced wild type Arabidopsis suffered necrotic death. In addition, we observed that overexpression of OnPPO and HaPPO in Arabidopsis conferred resistance to the PPO-inhibiting herbicides tiafenacil and saflufenacil. These results suggest that two HemY-type PPOs of cyanobacteria can functionally substitute for plastidic PPO activity in Arabidopsis and can enhance resistance to tiafenacil and saflufenacil.

Effects of herbicide on non-target microorganisms: Towards a new class of biomarkers?

Conventional agriculture still relies on the general use of agrochemicals (herbicides, fungicides and insecticides) to control various pests (weeds, fungal pathogens and insects), to ensure the yield of crop and to feed a constantly growing population. The generalized use of pesticides in agriculture leads to the contamination of soil and other connected environmental resources. The persistence of pesticide residues in soil is identified as a major threat for in-soil living organisms that are supporting an important number of ecosystem services. Although authorities released pesticides on the market only after their careful and thorough evaluation, the risk assessment for in-soil living organisms is unsatisfactory, particularly for microorganisms for which pesticide toxicity is solely considered by one global test measuring N mineralization. Recently, European Food Safety Authority (EFSA) underlined the lack of standardized methods to assess pesticide ecotoxicological effects on soil microorganisms. Within this context, there is an obvious need to develop innovative microbial markers sensitive to pesticide exposure. Biomarkers that reveal direct effects of pesticides on microorganisms are often viewed as the panacea. Such biomarkers can only be developed for pesticides having a mode of action inhibiting a specific enzyme not only found in the targeted organisms but also in microorganisms which are considered as "non-target organisms" by current regulations. This review explores possible ways of innovation to develop such biomarkers for herbicides. We scanned the herbicide classification by considering the mode of action, the targeted enzyme and the ecotoxicological effects of each class of active substance in order to identify those that can be tracked using sensitive microbial markers.

A Novel Single-Site Mutation in the Catalytic Domain of Protoporphyrinogen Oxidase IX (PPO) Confers Resistance to PPO-Inhibiting Herbicides

Protoporphyrinogen oxidase (PPO)-inhibiting herbicides are used to control weeds in a variety of crops. These herbicides inhibit heme and photosynthesis in plants. PPO-inhibiting herbicides are used to control Amaranthus palmeri (Palmer amaranth) especially those with resistance to glyphosate and acetolactate synthase (ALS) inhibiting herbicides. While investigating the basis of high fomesafen-resistance in A. palmeri, we identified a new amino acid substitution of glycine to alanine in the catalytic domain of PPO2 at position 399 (G399A) (numbered according to the protein sequence of A. palmeri). G399 is highly conserved in the PPO protein family across eukaryotic species. Through combined molecular, computational, and biochemical approaches, we established that PPO2 with G399A mutation has reduced affinity for several PPO-inhibiting herbicides, possibly due to steric hindrance induced by the mutation. This is the first report of a PPO2 amino acid substitution at G399 position in a field-selected weed population of A. palmeri. The mutant A. palmeri PPO2 showed high-level in vitro resistance to different PPO inhibitors relative to the wild type. The G399A mutation is very likely to confer resistance to other weed species under selection imposed by the extensive agricultural use of PPO-inhibiting herbicides.

Two new PPX2 mutations associated with resistance to PPO-inhibiting herbicides in Amaranthus palmeri

BACKGROUND

Resistance to herbicides that inhibit protoporphyrinogen oxidase (PPO) is a widespread and growing problem for weed managers across the midwestern and midsouthern United States. In Amaranthus spp., this resistance is known to be conferred by a glycine deletion at the 210^th amino acid (ΔG210) in PPO2. Preliminary analysis indicated that the ΔG210 mutation did not fully account for observed resistance to PPO inhibitors in two Amaranthus palmeri populations from Tennessee and one from Arkansas.

RESULTS

Sequencing PPX2 cDNA from six resistant plants uncovered two new mutations at the R98 site (R98G and R98M), a site previously found to endow PPO-inhibitor resistance in Ambrosia artemisiifolia. Sequencing of this region from additional plants sprayed with 264 g fomesafen ha^-1 showed the presence of one or both R98 mutations in a subset of the resistant plants from all three populations. No plants sensitive to fomesafen contained either mutation. A derived cleaved amplified polymorphic sequence (dCAPS) assay to test for the presence of these mutations in A. palmeri was developed.

CONCLUSION

Two new mutations of PPX2 (R98G, R98M) likely confer resistance to PPO-inhibitors in A. palmeri, and can be rapidly identified using a dCAPS assay. © 2017 Society of Chemical Industry.

Mechanism of Developmental Effects in Rats Caused by an N-Phenylimide Herbicide: Transient Fetal Anemia and Sequelae during Mid-to-Late Gestation

BACKGROUND

Rat developmental toxicity including embryolethality and teratogenicity (mainly ventricular septal defects [VSDs] and wavy ribs) was produced by an N-phenylimide herbicide that inhibits protoporphyrinogen oxidase (PPO) common to chlorophyll and heme biosynthesis. Major characteristics of the developmental toxicity included species difference between rats and rabbits, compound-specific difference among structurally similar herbicides, and sensitive period. Protoporphyrin accumulation in treated fetuses closely correlated with the major characteristics. Iron deposits in erythroblastic mitochondria and degeneration of erythroblasts were observed in treated rat fetuses. In this study we investigated fetal anemia and subsequent developmental effects in rats, and inhibition of PPO in rats, rabbits, and humans by the herbicides in vitro.

METHODS

Fetuses were treated on gestational day (GD) 12 and removed on GDs 13 through 20. All litters were examined externally. One half of litters were examined for blood and skeletal development, and the other half for interventricular foramen closure. Effects on PPO were determined in mitochondria from embryos and adult livers.

RESULTS

Fetal anemia in rats was evident on GDs 13 through 16. Subsequently, enlarged heart, delayed closure of the foramen, reduced serum protein, and retarded rib ossification were observed. In vitro PPO inhibition exhibited species- and compound-specific differences corresponding to the developmental toxicity.

CONCLUSION

We propose that developmental toxicity results from PPO inhibition in primitive erythroblasts, causing transient fetal anemia followed by death. Compensatory enlargement of the fetal heart results in failure of interventricular foramen closure and VSD. Reduced serum protein leads to delayed ossification and wavy ribs.

Molecular phylogeny and intricate evolutionary history of the three isofunctional enzymes involved in the oxidation of protoporphyrinogen IX

Tetrapyrroles such as heme and chlorophyll are essential for biological processes, including oxygenation, respiration, and photosynthesis. In the tetrapyrrole biosynthesis pathway, protoporphyrinogen IX oxidase (Protox) catalyzes the formation of protoporphyrin IX, the last common intermediate for the biosynthesis of heme and chlorophyll. Three nonhomologous isofunctional enzymes, HemG, HemJ, and HemY, for Protox have been identified. To reveal the distribution and evolution of the three Protox enzymes, we identified homologs of each along with other heme biosynthetic enzymes by whole-genome clustering across three domains of life. Most organisms possess only one of the three Protox types, with some exceptions. Detailed phylogenetic analysis revealed that HemG is mostly limited to γ-Proteobacteria whereas HemJ may have originated within α-Proteobacteria and transferred to other Proteobacteria and Cyanobacteria. In contrast, HemY is ubiquitous in prokaryotes and is the only Protox in eukaryotes, so this type may be the ancestral Protox. Land plants have a unique HemY homolog that is also shared by Chloroflexus species, in addition to the main HemY homolog originating from Cyanobacteria. Meanwhile, organisms missing any Protox can be classified into two groups; those lacking most heme synthetic genes, which necessarily depend on external heme supply, and those lacking only genes involved in the conversion of uroporphyrinogen III into heme, which would use a precorrin2-dependent alternative pathway. However, hemN encoding coproporphyrinogen IX oxidase was frequently found in organisms lacking Protox enzyme, which suggests a unique role of this gene other than in heme biosynthesis.

Understanding resistance mechanism of protoporphyrinogen oxidase-inhibiting herbicides: insights from computational mutation scanning and site-directed mutagenesis

The potential of protoporphyrinogen oxidase (PPO) to develop resistance against five PPO-inhibiting herbicides has been studied using computational mutation scanning (CMS) protocol, leading to valuable insights into the resistance mechanisms and structure-resistance relationship of the PPO inhibitors. The calculated shifts in the binding free energies caused by the mutations correlated very well with those derived from the corresponding experimental data obtained from site-directed mutagenesis of PPO, leading to valuable insights into the resistance mechanisms of PPO inhibitors. The calculated entropy change was related to the conformational flexibility of the inhibitor, which demonstrated that inhibitors with appropriate conformational flexibility may inhibit both the wild type and mutants simultaneously. The reasonable correlation between the computational and experimental data further validate that CMS protocol is valuable for predicting resistance associated with amino acid mutations on target proteins.

Binding structures and energies of the human neonatal Fc receptor with human Fc and its mutants by molecular modeling and dynamics simulations

Homology modeling and molecular dynamics simulations have been carried out to model the detailed structures of the human neonatal Fc receptor (FcRn) binding with the wild-type Fc of human immunoglobulin G1 (IgG1) and its various mutants. Based on the modeled human FcRn-Fc binding structures, it has been proposed that the protein-protein binding interface is composed of three subsites. The first subsite is a hydrophobic core where residue I39 of human Fc can be accommodated very well, and the other two subsites are all composed of critical salt bridges between human FcRn and human Fc. All of the modeled structures and the calculated binding energies are qualitatively consistent with the available experimental data, suggesting that the modeled human FcRn-Fc binding structures are reasonable. The modeled human FcRn-Fc binding structure may be valuable for future rational design of novel mutants of human Fc and Fc-fused therapeutic proteins with a potentially higher binding affinity for human FcRn and, thus, a longer in vivo half-life in humans.