Conformational Dynamics of Matrix Metalloproteinase-1·Triple-Helical Peptide Complexes

Identification and Structural Characterization of a Novel COL3A1 Gene Duplication in a Family With Vascular Ehlers-Danlos Syndrome

BACKGROUND

Vascular Ehlers-Danlos syndrome (vEDS) is caused by alterations in the COL3A1 gene, typically involving missense variants that replace glycine residues. In contrast, short in-frame insertions, deletions, and duplications are rare and pose significant challenges for investigation.

METHODS

The histological examination of vascular tissue from a 26-year-old man, who died from a common iliac artery aneurysm and whose mother died at age 60 from an abdominal aortic dissection, strongly suggested a diagnosis of Ehler-Danlos type IV. Ex vivo collagen phenotype assessment, molecular analysis, and in silico structural studies of type III collagen were subsequently performed.

RESULTS

Ex vivo analysis of the patient's fibroblasts revealed altered collagen synthesis, whereas the molecular testing identified a novel 18-nucleotide in-frame duplication (c.2868_2885dup-GGGTCTTGCAGGACCACC) in the COL3A1 gene, resulting in a six-amino acid insertion, p.(Leu958_Gly963dup). Structural investigation indicated that this duplication led to a local perturbation of the collagen triple helix near a metalloproteinase cleavage site.

CONCLUSION

This study highlights the pathogenic role of a novel in-frame duplication in the COL3A1 gene, demonstrating how this seemingly benign alteration significantly compromises collagen turnover and contributes to the development of vEDS.

Novel Insights into the Catalytic Mechanism of Collagenolysis by Zn(II)-Dependent Matrix Metalloproteinase-1

Collagen hydrolysis, catalyzed by Zn(II)-dependent matrix metalloproteinases (MMPs), is a critical physiological process. Despite previous computational investigations into the catalytic mechanisms of MMP-mediated collagenolysis, a significant knowledge gap in understanding remains regarding the influence of conformational sampling and entropic contributions at physiological temperature on enzymatic collagenolysis. In our comprehensive multilevel computational study, employing quantum mechanics/molecular mechanics (QM/MM) metadynamics (MetD) simulations, we aimed to bridge this gap and provide valuable insights into the catalytic mechanism of MMP-1. Specifically, we compared the full enzyme-substrate complex in solution, clusters in solution, and gas-phase to elucidate insights into MMP-1-catalyzed collagenolysis. Our findings reveal significant differences in the catalytic mechanism when considering thermal effects and the dynamic evolution of the system, contrasting with conventional static potential energy surface QM/MM reaction path studies. Notably, we observed a significant stabilization of the critical tetrahedral intermediate, attributed to contributions from conformational flexibility and entropy. Moreover, we found that protonation of the scissile bond nitrogen occurs via proton transfer from a Zn(II)-coordinated hydroxide rather than from a solvent water molecule. Following C-N bond cleavage, the C-terminus remains coordinated to the catalytic Zn(II), while the N-terminus forms a hydrogen bond with a solvent water molecule. Subsequently, the release of the C-terminus is facilitated by the coordination of a water molecule. Our study underscores the pivotal role of protein conformational dynamics at physiological temperature in stabilizing the transition state of the rate-limiting step and key intermediates, compared to the corresponding reaction in solution. These fundamental insights into the mechanism of collagen degradation provide valuable guidance for the development of MMP-1-specific inhibitors.

Engineering metalloproteinase inhibitors: tissue inhibitors of metalloproteinases or antibodies, that is the question

Targeting metalloproteinases (MPs) has been the center of attention for developing therapeutics due to their contribution to a wide range of diseases, including cancer, cardiovascular, neurodegenerative disease, and preterm labor. Protein-based MP inhibitors offer higher stability and selectivity, which is critical for developing efficient therapeutics with low off-target effects. Tissue inhibitors of metalloproteinases (TIMPs), natural inhibitors of MPs, and antibodies provide excellent protein scaffolds for engineering selective or multispecific MP inhibitors. Advances in protein engineering and design techniques, such as rational design and directed evolution using yeast display to develop potent MP inhibitors, are discussed, including but not limited to loop grafting, swapping, and counterselective selection.

Catalytic Mechanism of Collagen Hydrolysis by Zinc(II)-Dependent Matrix Metalloproteinase-1

Human matrix metalloproteinase-1 (MMP-1) is a zinc(II)-dependent enzyme that catalyzes collagenolysis. Despite the availability of extensive experimental data, the mechanism of MMP-1-catalyzed collagenolysis remains poorly understood due to the lack of experimental structure of a catalytically productive enzyme-substrate complex of MMP-1. In this study, we apply molecular dynamics and combined quantum mechanics/molecular mechanics to reveal the reaction mechanism of MMP-1 based on a computationally modeled structure of the catalytically competent complex of MMP-1 that contains a large triple-helical peptide substrate. Our proposed mechanism involves the participation of an auxiliary (second) water molecule (wat2) in addition to the zinc(II)-coordinated water (wat1). The reaction initiates through a proton transfer to Glu219, followed by a nucleophilic attack by a zinc(II)-coordinated hydroxide anion nucleophile at the carbonyl carbon of the scissile bond, leading to the formation of a tetrahedral intermediate (IM2). The process continues with a hydrogen-bond rearrangement to facilitate proton transfer from wat2 to the amide nitrogen of the scissile bond and, finally, C-N bond cleavage. The calculations indicate that the rate-determining step is the water-mediated nucleophilic attack with an activation energy barrier of 22.3 kcal/mol. Furthermore, the calculations show that the hydrogen-bond rearrangement/proton-transfer step can proceed in a consecutive or concerted manner, depending on the conformation of the tetrahedral intermediate, with the consecutive mechanism being energetically preferable. Overall, the study reveals the crucial role of a second water molecule and the dynamics for effective MMP-1-catalyzed collagenolysis.

Molecular Dynamics Simulations of Matrix Metalloproteinase 13 and the Analysis of the Specificity Loop and the S1'-Site

The specificity loop of Matrix Metalloproteinases (MMPs) is known to regulate recognition of their substrates, and the S1'-site surrounded by the loop is a unique place to address the selectivity of ligands toward each MMP. Molecular dynamics (MD) simulations of apo-MMP-13 and its complex forms with various ligands were conducted to identify the role of the specificity loop for the ligand binding to MMP-13. The MD simulations showed the dual role of T247 as a hydrogen bond donor to the ligand, as well as a contributor to the formation of the van der Waal surface area, with T245 and K249 on the S1'-site. The hydrophobic surface area mediated by T247 blocks the access of water molecules to the S1'-site of MMP-13 and stabilizes the ligand in the site. The F252 residue is flexible in order to search for the optimum location in the S1'-site of the apo-MMP-13, but once a ligand binds to the S1'-site, it can form offset π-π or edge-to-π stacking interactions with the ligand. Lastly, H222 and Y244 provide the offset π-π and π-CH(Cβ) interactions on each side of the phenyl ring of the ligand, and this sandwiched interaction could be critical for the ligand binding to MMP-13.

Mechanism of the Early Catalytic Events in the Collagenolysis by Matrix Metalloproteinase-1

Metalloproteinase-1 (MMP-1) catalyzed collagen degradation is essential for a wide variety of normal physiological processes, while at the same time contributing to several diseases in humans. Therefore, a comprehensive understanding of this process is of great importance. Although crystallographic and spectroscopic studies provided fundamental information about the structure and function of MMP-1, the precise mechanism of collagen degradation especially considering the complex and flexible structure of the substrate, remains poorly understood. In addition, how the protein environment dynamically reorganizes at the atomic scale into a catalytically active state capable of collagen hydrolysis remains unknown. In this study, we applied experimentally-guided multiscale molecular modeling methods including classical molecular dynamics (MD), well-tempered (WT) classical metadynamics (MetD), combined quantum mechanics/molecular mechanics (QM/MM) MD and QM/MM MetD simulations to explore and characterize the early catalytic events of MMP-1 collagenolysis. Importantly the study provided a complete atomic and dynamic description of the transition from the open to the closed form of the MMP-1•THP complex. Notably, the formation of catalytically active Michaelis complex competent for collagen cleavage was characterized. The study identified the changes in the coordination state of the catalytic zinc(II) associated with the conformational transformation and the formation of catalytically productive ES complex. Our results confirm the essential role of the MMP-1 catalytic domain's α-helices (hA, hB and hC) and the linker region in the transition to the catalytically competent ES complex. Overall, the results provide unique mechanistic insight into the conformational transformations and associated changes in the coordination state of the catalytic zinc(II) that would be important for the design of effective MMP-1 inhibitors.

Structure of Vibrio collagenase VhaC provides insight into the mechanism of bacterial collagenolysis

The collagenases of Vibrio species, many of which are pathogens, have been regarded as an important virulence factor. However, there is little information on the structure and collagenolytic mechanism of Vibrio collagenase. Here, we report the crystal structure of the collagenase module (CM) of Vibrio collagenase VhaC and the conformation of VhaC in solution. Structural and biochemical analyses and molecular dynamics studies reveal that triple-helical collagen is initially recognized by the activator domain, followed by subsequent cleavage by the peptidase domain along with the closing movement of CM. This is different from the peptidolytic mode or the proposed collagenolysis of Clostridium collagenase. We propose a model for the integrated collagenolytic mechanism of VhaC, integrating the functions of VhaC accessory domains and its collagen degradation pattern. This study provides insight into the mechanism of bacterial collagenolysis and helps in structure-based drug design targeting of the Vibrio collagenase.

The collagenolytic mechanism of Vibrio collagenase, a virulence factor, remains unclear. Here, the authors report the structure of Vibrio collagenase VhaC and propose the mechanism for collagen recognition and degradation, providing new insight into bacterial collagenolysis.

Effects of the Nature of the Metal Ion, Protein and Substrate on the Catalytic Center in Matrix Metalloproteinase-1: Insights from Multilevel MD, QM/MM and QM Studies

Matrix metalloproteinase-1 (MMP-1) is a Zn(II) dependent endopeptidase involved in the degradation of collagen, the most abundant structural protein in the extracellular matrix of connective tissues and the human body. Herein we performed a multilevel computational analysis including molecular dynamics (MD), combined quantum mechanics/molecular mechanics (QM/MM), and quantum mechanics (QM) calculations to characterize the structure and geometry of the catalytic Zn(II) within the MMP-1 protein environment in comparison to crystallographic and spectroscopic data. The substrate's removal fine-tuned impact on the conformational dynamics and geometry of the catalytic Zn(II) center was also explored. Finally, the study examined the effect of substituting catalytic Zn(II) by Co(II) on the overall structure and dynamics of the MMP-1 THP complex and specifically on the geometry of the catalytic metal center. Overall our QM/MM and QM studies were in good agreement with the MM description of the Zn(II) centers in the MD simulations.

A synergy between the catalytic and structural Zn(II) ions and the enzyme and substrate dynamics underlies the structure-function relationships of matrix metalloproteinase collagenolysis

Matrix metalloproteinases (MMPs) are Zn(II) dependent endopeptidases involved in the degradation of collagen. Unbalanced collagen breakdown results in numerous pathological conditions, including cardiovascular and neurodegenerative diseases and tumor growth and invasion. Matrix metalloproteinase-1 (MMP-1) is a member of the MMPs family. The enzyme contains catalytic and structural Zn(II) ions. Despite many studies on the enzyme, there is little known about the synergy between the two Zn(II) metal ions and the enzyme and substrate dynamics in MMP-1 structure-function relationships. We performed a computational study of the MMP-1•triple-helical peptide (THP) enzyme•substrate complex to provide this missing insight. Our results revealed Zn(II) ions' importance in modulating the long-range correlated motions in the MMP-1•THP complex. Overall, our results reveal the importance of the catalytic Zn(II) and the role of the structural Zn(II) ion in preserving the integrity of the enzyme active site and the overall enzyme-substrate complex synergy with the dynamics of the enzyme and the substrate. Notably, both Zn(II) sites participate in diverse networks of long-range correlated motions that involve the CAT and HPX domains and the THP substrate, thus exercising a complex role in the stability and functionality of the MMP-1•THP complex. Both the Zn(II) ions have a distinct impact on the structural stability and dynamics of the MMP-1•THP complex. The study shifts the paradigm from the "local role" of the Zn(II) ions with knowledge about their essential role in the long-range dynamics and stability of the overall enzyme•substrate (ES) complex.

Comparative QM/MM Molecular Dynamics and Umbrella Sampling Simulations: Interaction of the Zinc-Bound Intermediate Gem-Diolate Trapoxin A Inhibitor and Acetyl-l-lysine Substrate with Histone Deacetylase 8

Targeting the genetic material without destruction is a priority to develop safe anticancer drugs. Histone deacetylase 8 (HDAC8), which is proved to be involved in carcinogenesis, is an enzyme associated with the chromatin for post-translational deacetylation of acetylated lysine. In this study, HDAC8 co-crystallized with the intermediate state tetrapeptide Trapoxin A (TA) inhibitor and the holoenzyme are utilized to find their conformational ensembles. Furthermore, the co-crystallized intermediate gem-diolate TA was used to find optimum interactions with the active site residues by conventional molecular dynamics (MD) simulation and QM/MM umbrella sampling. Finally, the intermediate state of the acetyl-l-lysine substrate was explored by QM/MM steered MD and compared to the binding of the intermediate state of the inhibitor. This research showed that HDAC8 is flexible and exists in conformational ensembles in its holoenzyme state. Binding of the intermediate state TA stabilizes its conformation. The optimum binding to the active site of HDAC8 for structures of gem-diolate TA (intermediate state) and acetyl-l-lysine (intermediate state) was determined according to the corresponding energy profiles. The use of these models will aid in the design of potentially reversible, potent, and selective inhibitors of HDAC8 for cancer treatment.

Activity-dependent interdomain dynamics of matrix metalloprotease-1 on fibrin

The roles of protein conformational dynamics and allostery in function are well-known. However, the roles that interdomain dynamics have in function are not entirely understood. We used matrix metalloprotease-1 (MMP1) as a model system to study the relationship between interdomain dynamics and activity because MMP1 has diverse substrates. Here we focus on fibrin, the primary component of a blood clot. Water-soluble fibrinogen, following cleavage by thrombin, self-polymerize to form water-insoluble fibrin. We studied the interdomain dynamics of MMP1 on fibrin without crosslinks using single-molecule Forster Resonance Energy Transfer (smFRET). We observed that the distance between the catalytic and hemopexin domains of MMP1 increases or decreases as the MMP1 activity increases or decreases, respectively. We modulated the activity using (1) an active site mutant (E219Q) of MMP1, (2) MMP9, another member of the MMP family that increases the activity of MMP1, and (3) tetracycline, an inhibitor of MMP1. We fitted the histograms of smFRET values to a sum of two Gaussians and the autocorrelations to an exponential and power law. We modeled the dynamics as a two-state Poisson process and calculated the kinetic rates from the histograms and autocorrelations. Activity-dependent interdomain dynamics may enable allosteric control of the MMP1 function.

Thermodynamic and Mechanistic Insights into Coupled Binding and Unwinding of Collagen by Matrix Metalloproteinase 1

Local unwinding of the collagen triple helix is a necessary step for initiating the collagen degradation cascade in extracellular matrices. A few matrix metalloproteinases (MMPs) are known to support this key process, but its energetic aspects remain unknown. Here, we captured the thermodynamics of the triple helix unwinding by monitoring interactions between a collagen peptide and MMP-1(E200A) - an active-site mutant of an archetypal vertebrate collagenase - at increasing temperatures, using isothermal titration calorimetry (ITC). Coupled binding and unwinding manifests as a curved relationship between the total enthalpy change and temperature of the reaction, producing increasingly negative heat capacity change (ΔΔC_p ≈ -36.3 kcal/molK²). A specially designed solid-phase binding and cleavage assay (SPBCA) reported strain in the catalytically relevant unwound state, suggesting that this state is distinct from the horizon of sampled conformations of the collagenase-susceptible site. MMP-1 appears to blend selected fit with induced fit mechanisms to catalyse collagen unwinding prior to cleavage of individual collagen chains.

Allosteric Communications between Domains Modulate the Activity of Matrix Metalloprotease-1

An understanding of the structure-dynamics relationship is essential for understanding how a protein works. Prior research has shown that the activity of a protein correlates with intradomain dynamics occurring at picosecond to millisecond timescales. However, the correlation between interdomain dynamics and the function of a protein is poorly understood. Here, we show that communications between the catalytic and hemopexin domains of matrix metalloprotease-1 (MMP1) on type 1 collagen fibrils correlate with its activity. Using single-molecule Förster resonance energy transfer, we identified functionally relevant open conformations in which the two MMP1 domains are well separated, which were significantly absent for catalytically inactive point mutant (E219Q) of MMP1 and could be modulated by an inhibitor or an enhancer of activity. The observed relevance of open conformations resolves the debate about the roles of open and closed MMP1 structures in function. We fitted the histograms of single-molecule Förster resonance energy transfer values to a sum of two Gaussians and the autocorrelations to an exponential and power law. We used a two-state Poisson process to describe the dynamics and calculate the kinetic rates from the fit parameters. All-atom and coarse-grained simulations reproduced some of the experimental features and revealed substrate-dependent MMP1 dynamics. Our results suggest that an interdomain separation facilitates opening up the catalytic pocket so that the collagen chains come closer to the MMP1 active site. Coordination of functional conformations at different parts of MMP1 occurs via allosteric communications that can take place via interactions mediated by collagen even if the linker between the domains is absent. Modeling dynamics as a Poisson process enables connecting the picosecond timescales of molecular dynamics simulations with the millisecond timescales of single-molecule measurements. Water-soluble MMP1 interacting with water-insoluble collagen fibrils poses challenges for biochemical studies that the single-molecule tracking can overcome for other insoluble substrates. Interdomain communications are likely important for multidomain proteins.

Interaction of Selected Terpenoids From Dalbergia sissoo With Catalytic Domain of Matrix Metalloproteinase-1: An In Silico Assessment of Their Anti-wrinkling Potential

Matrix metalloproteinase-1 (MMP-1) is a predominant collagenase enzyme that cleaves collagen fibers, contributing to skin wrinkling. Matrix metalloproteinase-1 inhibitors of herbal origin may provide an earnest probability to offer a novel curative approach against MMP-1-mediated collagenolysis, prompted by ultraviolet (UV)-induced overexpression of MMP-1. In this in silico study, we have explored the MMP-1 inhibitory potential of selected terpenoids from Dalbergia sissoo extracts. Two triterpenoids (lupeol and betulin), 1 diterpenoid (phytol), and 1 ester derivative of lupeol (lupeol acetate) were studied along with a reference inhibitor (doxycycline) using molecular docking approach. Non covalent interaction between the target ligands was found. Lupeol was found interacting with amino acid (AA) residues in the catalytic domain of MMP-1 with 3 hydrogen bonds (H-bond) formation, phytol with 1 and doxycycline with 2 H-bonds, whereas betulin and lupeol acetate were not able to form any H-bond with the AA residues in the catalytic site of the target protein. However, hydrophobic interaction between these ligands and protein was evident with select residues. The binding affinity of lupeol was highest (binding free energy, ΔG = −8.24 kcal/mol), which was greater than reference drug, doxycycline (ΔG = −8.05 kcal/mol). Lupeol acetate and phytol displayed a ΔG value of −7.12 and −7.06 kcal/mol, respectively, whereas betulin holds less binding affinity for the target receptor (ΔG = −4.66 kcal/mol). In silico pharmacokinetic studies demonstrated drug-like properties of the ligand compounds. This study shows that hydroxyl groups present in the ligands play a substantial role in establishing protein ligand interaction via hydrogen bonding.

Conformational flexibility influences structure-function relationships in nucleic acid N-methyl demethylases

N-Methylation of DNA/RNA bases can be regulatory or damaging and is linked to diseases including cancer and genetic disorders. Bacterial AlkB and human FTO are DNA/RNA demethylases belonging to the Fe(ii) and 2-oxoglutarate oxygenase superfamily. Modelling studies reveal conformational dynamics influence structure-function relationships of AlkB and FTO, e.g. why 1-methyladenine is a better substrate for AlkB than 6-methyladenine. Simulations show that the flexibility of the double stranded DNA substrate in AlkB influences correlated motions, including between the core jelly-roll fold and an active site loop involved in substrate binding. The FTO N- and C-terminal domains move in respect to one another in a manner likely important for substrate binding. Substitutions, including clinically observed ones, influencing catalysis contribute to the network of correlated motions in AlkB and FTO. Overall, the calculations highlight the importance of the overall protein environment and its flexibility to the geometry of the reactant complexes.